Question: GSEA using clusterprofileR
gravatar for sup230
5 months ago by
sup23010 wrote:


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!

ADD COMMENTlink modified 5 months ago by assaron140 • written 5 months ago by sup23010
gravatar for assaron
5 months ago by
assaron140 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 (, comments on biorXiv, and also here:



ADD COMMENTlink written 5 months ago by assaron140
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.2.0
Traffic: 337 users visited in the last hour