Question: GSEA using clusterProfiler
0
gravatar for dxc008263
3 months ago by
dxc0082630
dxc0082630 wrote:

Hi, I am trying to do GSEA with clusterProfiler using gseKEGG function. However, I keep getting this error:

preparing geneSet collections... --> Expected input gene ID: 18639,11997,74551,18642,104112,18640 Error in check_gene_id(geneList, geneSets) : --> No gene can be mapped....

This is what I did. I sorted my gene list first:

GSEAList <- sort(WTvsKoBasal,decreasing = TRUE) head(GSEAList) Gene.Id112 Gene.Id264 Gene.Id156 Gene.Id121 Gene.Id177 Gene.Id403 "99899" "99543" "98365" "97114" "97086" "94284"

Then I ran this command:

WTvsKOBasalGSEA<-gseKEGG(geneList = GSEAList, organism = "mmu", keyType = "kegg", exponent = 1).

Then I will get an error as I mentioned in the beginning.

Any suggestion will be appreciated!

Thanks,

Ben

ADD COMMENTlink modified 3 months ago • written 3 months ago by dxc0082630

Maybe you have your IDs as numbers - clusterProfiler expects them to be characters.

ADD REPLYlink written 3 months ago by mkulecka300
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: 1598 users visited in the last hour