Question: GSVA on enrichment scores about R project
0
gravatar for 1466956114
17 months ago by
14669561140
14669561140 wrote:

Dear all,

I am trying to use GSVA for enrichment scores. I obtain data of gene expression matrix form TCGA, and the the genesets are created by myself as needed. I am a beginner in the R project and have done a lot of hard work, but I am still very confused. Input data format: Gene expression matrix:(for example) genename sample1 sample2 sample3 A 1 2 0 B 2 3 1 C 1 0 2 Genesets data: genesetname description gene geneset1 NA A B C geneset2 NA B D E

Can anyone please give me some suggestions how to use GSVA for my data analysis? Thanks very much. best XX

gene <- read.csv("C:\Users\86188\Desktop\expression.csv",header = TRUE,sep = ",") genesets <- getGmt("C:\Users\86188\Desktop\symbols.gmt") score <- gsva(gene, genesets,method= "ssgsea", + rnaseq=TRUE,abs.ranking=FALSE, min.sz=1, max.sz=Inf, no.bootstraps=0, bootstrap.percent = .632, parallel.sz=0, parallel.type="SOCK", mx.diff=TRUE, tau=switch(method, gsva=1, ssgsea=0.25, NA),kernel=TRUE,ssgsea.norm=TRUE,verbose=TRUE) Error in (function (classes, fdef, mtable) : unable to find an inherited method for function ‘gsva’ for signature ‘"data.frame", "GeneSetCollection"’

rna-seq gsva R • 776 views
ADD COMMENTlink modified 6 weeks ago by DBScan20 • written 17 months ago by 14669561140
0
gravatar for DBScan
6 weeks ago by
DBScan20
DBScan20 wrote:

The first error is because "gene" is a data.frame instead of a matrix (I had the same error same days ago too). The second error is because the GMT file has a wrong format I think. You can check the GMT format here: http://software.broadinstitute.org/cancer/software/gsea/wiki/index.php/Data_formats

ADD COMMENTlink written 6 weeks ago by DBScan20
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 2.3.0
Traffic: 696 users visited in the last hour