Question: ClusterProfiler GSEA with geneList
0
gravatar for lucap
4 months ago by
lucap0
lucap0 wrote:

Hi everyone, Just a quick question. I want to perform GSEA on my geneList using clusterProfiler. Since I have a lot of genes and for some of them the logFC is small, I would like to perform GSEA only on the genes that have a logFC > 2. How do I do that in R?

I am following the recommended steps to prepare the geneList:

geneList = d[,2] names(geneList) = as.character(d[,1]) geneList = sort(geneList, decreasing = TRUE)

followed by:

gene <- names(geneList)[abs(geneList) > 2]

but then I do not understand how to "filter" the gene list that will be submitted to GSEA? Maybe I am not supposed to do it, and I should only run GSEA with the whole geneList set? The code for GSEA I am using is:

kkGSEA <- gseKEGG(geneList=geneList, organism = "hsa", nPerm = 1000, pvalueCutoff = pvalue_gsea, minGSSize = 5)

Any help would be much much appreciated. Thanks :D

clusterprofiler gsea genelist • 229 views
ADD COMMENTlink modified 4 months ago by Ludwig Geistlinger50 • written 4 months ago by lucap0
Answer: ClusterProfiler GSEA with geneList
0
gravatar for Ludwig Geistlinger
4 months ago by
USA/NYC/CUNY
Ludwig Geistlinger50 wrote:

GSEA works on the full gene vector, testing whether genes of a gene set (here: a KEGG pathway) accumulate at the top or bottom of the full gene vector ordered by direction and magnitude of expression change.

If you are interested only in genes of significant expression change (here: abs(log2FC) > 2), and want to know whether certain gene sets (here: KEGG pathways) contain a disproportional number of these significant genes, you would rather carry out a over-representation analysis (ORA).

In clusterProfiler: use enrichKEGG for that purpose.

ADD COMMENTlink modified 4 months ago • written 4 months ago by Ludwig Geistlinger50

Thanks Ludwig for the clear explanation! I see what you mean. I’ll carry the over-representation test then

Thanks Luca

ADD REPLYlink written 4 months ago by lucap0
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: 468 users visited in the last hour