Question: Best method/package for Gene Set Enrichment Analysis in R?
1
gravatar for peter pfand
3.8 years ago by
peter pfand80
Germany
peter pfand80 wrote:

Dear community,

I recently started working with GSEA of microarray Data in Bioconductor and after a quick search I am quite overwhelmed because the wide supply of different packages to compute GSEA for a given list of differentially expressed probe sets / genes (goseq, topGO, gega, gsea, GOstats...). In many cases, every method claims to be the best, and it's getting hard for me to choose a proper method.

Do you know if there's any benchmark or rather which method is the best one? and for the analysis of KEGG pathways enrichment?

 

go gsea bioconductor • 19k views
ADD COMMENTlink modified 3.7 years ago by bigmawen310 • written 3.8 years ago by peter pfand80

I personally use gage, because it's super easy to make your own custom gene-sets and ranked lists for GSEA. It's also really easy to use it in conjunction with pathview, which is a nice R package for pathway visualization.

ADD REPLYlink written 3.8 years ago by informatics bot570

Hi Lando,

I have a table of DE genes with their respective P- and Q-values calculated with LIMMA package and I would like to start from this data to make the GSEA. I've been reading the gega vignette and it seems it must start from the beginning, with the expression matrix. Do you know whether there's another way to make the GSEA starting from a list of genes with their p-values but using gega?

ADD REPLYlink modified 3.7 years ago • written 3.8 years ago by peter pfand80
3
gravatar for bigmawen
3.7 years ago by
bigmawen310
United States
bigmawen310 wrote:

One thing needs to be pointed our first, the package is called gage (method called GAGE) instead of gega.

Gage/pathview can take outputs from all major RNA-Seq/expression analysis tool, including limma, deseq2/deseq, edger or cufflinks. You should use log2 fold change or their test statistics columns (if available). P-values would also work (although not preferred) when you set rank.test=T when calling gage.

Please check the joint workflow sections in the RNA-Seq pathway analysis tutorial: http://bioconductor.org/packages/release/bioc/vignettes/gage/inst/doc/RNA-seqWorkflow.pdf

ADD COMMENTlink written 3.7 years ago by bigmawen310
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: 1593 users visited in the last hour