Question: GSVA for single sample
0
gravatar for jaison75
21 months ago by
jaison7520
United States
jaison7520 wrote:

Hello,

I am trying to use GSVA to calculate an enrichment score for a single sample.  I thought the ssgsea method implemented in GSVA can perm such a calculation.  Do I have this wrong?  My searches have returned empty handed.  I am pasting an example code below I used to test the concept (modifed from the help page). 

 

## consider three disjoint gene sets
geneSets <- list(set1=paste("g", 1:3, sep=""),
                 set2=paste("g", 4:6, sep=""),
                 set3=paste("g", 7:10, sep=""))

y1 <- as.matrix(runif( n = 10, min = 5.3, max = 15),ncol =1)
colnames(y1) <- c("s1")
row.names(y1) <- paste("g", 1:10, sep="")
gsva(y1, geneSets, method == "ssgsea", ssgsea.norm = FALSE, mx.diff = FALSE, kernel = F)

Any help would be appreciated.

 

Thanks

-Jaison

gsva gsea ssgsea • 479 views
ADD COMMENTlink modified 21 months ago • written 21 months ago by jaison7520
Answer: GSVA for single sample
1
gravatar for Robert Castelo
21 months ago by
Robert Castelo2.3k
Spain/Barcelona/Universitat Pompeu Fabra
Robert Castelo2.3k wrote:

hi,

you have a typo in your code, the '==' operator is for comparison, not for assignment. please replace

method == "ssgsea"

by

method = "ssgsea"

and it will work.

cheers,

robert.

ADD COMMENTlink written 21 months ago by Robert Castelo2.3k
Answer: GSVA for single sample
0
gravatar for jaison75
21 months ago by
jaison7520
United States
jaison7520 wrote:

Thank you very much Robert.  I cannot believe, I did not catch that. 

This solved the issue.

ADD COMMENTlink written 21 months ago by jaison7520
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 205 users visited in the last hour