Question: Best method/package for Gene Set Enrichment Analysis in R?
1
gravatar for peter pfand
4.6 years ago by
peter pfand90
Germany
peter pfand90 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 • 24k views
ADD COMMENTlink modified 4.6 years ago by bigmawen340 • written 4.6 years ago by peter pfand90

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 4.6 years ago by informatics bot610

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 4.6 years ago • written 4.6 years ago by peter pfand90
3
gravatar for bigmawen
4.6 years ago by
bigmawen340
United States
bigmawen340 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 4.6 years ago by bigmawen340
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: 1251 users visited in the last hour