Question: Error in GSEA function of clusterProfiler on R
0
gravatar for jimmy0958073736
7 weeks ago by
Taiwan
jimmy09580737360 wrote:

Hello everyone I just met some problem with GSEA when I run the clusterProfiler. Can anyone help me how to fix the problem, although the data still could be processed (I'm not so sure the data could be used under those error message.

Follow my code

geneList_organism = "mmu"                                                                                                                                            
           kk2 <- gseKEGG(organism     = geneList_organism,
           geneList     = geneList,
           nPerm        = 10000,
           minGSSize    = 3,
           maxGSSize    = 800,
           pvalueCutoff = 0.05,
           pAdjustMethod = "none")******

Output error message:

> preparing geneSet collections...
GSEA analysis...
leading edge analysis...
done...
Warning messages:

 1.  In .GSEA(geneList = geneList, exponent = exponent, minGSSize = minGSSize,  :
  We do not recommend using nPerm parameter incurrent and future releases
 2.  In fgsea(pathways = geneSets, stats = geneList, nperm = nPerm, minSize = minGSSize,  :
  You are trying to run fgseaSimple. It is recommended to use fgseaMultilevel. To run fgseaMultilevel, you need to remove the nperm argument in the fgsea function call.
 3. In preparePathwaysAndStats(pathways, stats, minSize, maxSize, gseaParam,  :
  There are ties in the preranked stats (16.25% of the list).
 4. In serialize(data, node$con) :  'package:stats' may not be available when loading
 5. In serialize(data, node$con) :  'package:stats' may not be available when loading
clusterprofiler gsea rna-seq R • 251 views
ADD COMMENTlink modified 6 weeks ago • written 7 weeks ago by jimmy09580737360

Check this post

ADD REPLYlink written 7 weeks ago by brunobsouzaa490

I already try that method than before. But it’s just show the other error message said It’s not suitable for hsa GSEA( but it’s really so wired that I do not type any associating with human species). So I wonder to know have any method to fix that.

ADD REPLYlink written 7 weeks ago by jimmy09580737360
kk2 <- gseKEGG(geneList=geneList, 
             keyType = "ENTREZID", 
             nPerm = 10000, 
             minGSSize = 3, 
             maxGSSize = 800, 
             pvalueCutoff = 0.05, 
             verbose = TRUE, 
             OrgDb = org.Mm.eg.db, 
             pAdjustMethod = "none")

It will process the following error message.

Reading KEGG annotation online:

Reading KEGG annotation online:

Reading KEGG annotation online:

fail to download KEGG data... Error in KEGG_convert("kegg", keyType, species) : ENTREZID is not supported for has ...

Then I try to exchange the code of following

gsepathway <- gseKEGG(geneList=geneList, 
             nPerm = 10000, 
             minGSSize = 3, 
             maxGSSize = 800, 
             pvalueCutoff = 0.05, 
             verbose = TRUE, 
             OrgDb = org.Mm.eg.db, 
             pAdjustMethod = "none")

It will process this error message

Error in .GSEA(geneList = geneList, exponent = exponent, minGSSize = minGSSize, : unused argument (OrgDb = new("OrgDb", .xData = <environment>))

How should I adjust my code?

ADD REPLYlink written 7 weeks ago by jimmy09580737360

Hi Jimmy, first of all you are not getting an error - this is a warning message. Error would be a break in the code. Here the code performs but tells you that something is fishy. Here is what I can say to your warnings:

Warning 1 - 2: I am still trying for myself how to change the fgseaSimple to fgseaMultilevel within clusterProfiler's gseKEGG function. However unfortunately with little success so far ... i think the problem is that fgsea was further developed recently and the fgseaMultilevel was implemented and is now considered to perform better than fgseaSimple while clusterProfiler did not integrate this novel development. Again, idk if that is true or not but I think it would explain the warning comming from the fgsea function. Maybe someone else in this community can help out here.

Warning 3: Exactly what the warning message says. Check here for more info.

Warning 4-5: As you nPerm is relatively large this could be the problem. Try some smaller nPerm like 1000 or 2000 and see if the warning still appears.

ADD REPLYlink written 7 weeks ago by hannes.reinwald20
0
gravatar for jimmy0958073736
6 weeks ago by
Taiwan
jimmy09580737360 wrote:

Thanks, Hannes

Warning 1/2 I think warning one and two seems like a long term problem with those involved in this field. Exactly I'm not so sure that I could put this figure on my review. In fact, I also ran out the same data to GSEA tools, but the data looks so different in both. I'm so worried about that.

Warning 3 Thanks for your hint. I think that won't be my problem again.

Warning 4/5 I also checked some user code and found that seems not so necessary to type this parameter. But I will take your advance to adjust THX.

ADD COMMENTlink written 6 weeks ago by jimmy09580737360
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: 2215 users visited in the last hour
_