Question: GSEA using clusterprofileR
2
gravatar for sup230
2.1 years ago by
sup23030
sup23030 wrote:

Hi, 

I am wondering how the result of GSEA using cluster profiler is different from that of GSEA GUI java application. I initially used GSEA GUI desktop application and tried clusterprofileR package in R using gseKEGG function. I assigned latest kegg database available online and pvalue cutoff of 0.05 for cluster profileR. I also assigned the same permutation number and minimum geneset size to be using the same condition as what I used for GSEA GUI software. However, the results are different; some gene sets I got as downregulated in one were upregulated in the other and vice versa. Anyone has an explanation for this and which is preferred for general practice?

Thank you for your input in advance!

clusterprofiler gsea • 638 views
ADD COMMENTlink modified 2.1 years ago by assaron150 • written 2.1 years ago by sup23030
Answer: GSEA using clusterprofileR
2
gravatar for assaron
2.1 years ago by
assaron150
assaron150 wrote:

clusterProfiler/fgsea differs from Broad GSEA in the multiple-hypotheses correction procedures. The latter uses an ad-hoc procedure while the former uses any standard method with Benjamini-Hochberg being the default one. I would recommend using at least 10000 permuations for  clusterProfiler/fgsea (which is still pretty fast).

In details this is covered in the fgsea manuscript (https://www.biorxiv.org/content/early/2016/06/20/060012), comments on biorXiv, and also here: https://bioinformatics.stackexchange.com/questions/149/are-fgsea-and-broad-institute-gsea-equivalent

 

 

ADD COMMENTlink written 2.1 years ago by assaron150
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: 360 users visited in the last hour