Question: GSVA for single sample
0
gravatar for jaison75
14 months ago by
jaison7510
United States
jaison7510 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 • 332 views
ADD COMMENTlink modified 14 months ago • written 14 months ago by jaison7510
Answer: GSVA for single sample
1
gravatar for Robert Castelo
14 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 14 months ago by Robert Castelo2.3k
Answer: GSVA for single sample
0
gravatar for jaison75
14 months ago by
jaison7510
United States
jaison7510 wrote:

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

This solved the issue.

ADD COMMENTlink written 14 months ago by jaison7510
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: 263 users visited in the last hour