Question: GSVA for single sample
0
gravatar for jaison75
19 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 • 456 views
ADD COMMENTlink modified 19 months ago • written 19 months ago by jaison7520
Answer: GSVA for single sample
1
gravatar for Robert Castelo
19 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 19 months ago by Robert Castelo2.3k
Answer: GSVA for single sample
0
gravatar for jaison75
19 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 19 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: 134 users visited in the last hour