Question: GSA for comparison of geneset enricment analysis between multiple class
0
20 months ago by
Seymoo0
Oslo
Seymoo0 wrote:

I have a gene expression data AllExp and a data.frame info with sample name of AllExp as row names and first column in it specifies groups that each sample belongs to.

I am trying to perform a a gene set enrichment analysis using GSA package. but I get following error

Error in cut.default(sd, br, labels = F) : invalid number of intervals

GS <- Sys.glob("E:/gmts/h.all.v6.1.symbols.gmt")
GSA(AllExp, y=factor(info[,1]), genesets=GS\$genesets, genenames = row.names(AllExp),method="maxmean", resp.type="Multiclass", minsize=10)

what is wrong with my code?

gsa • 244 views