Question: Performing GSEA using MSigDB gene sets in R
1
gravatar for wilsonav
11 months ago by
wilsonav10
wilsonav10 wrote:

I am trying to perform a gene set enrichment analysis in r using the gene sets available from msigdb and a list of gene names from my own data set.

I am able to to use the msigdbr library to import the gene collections from msigdb into r, but I am unsure of how to specifically use a function to compute the overlaps between the genes in my gene set and the gene sets in msigdb and obtain the FDR p-values. Are there any tutorials online for this method or example codes?

Thank you

gsea R msigdb • 1.6k views
ADD COMMENTlink modified 11 months ago by igor8.3k • written 11 months ago by wilsonav10
3
gravatar for igor
11 months ago by
igor8.3k
United States
igor8.3k wrote:

You can try the fgsea package, which is probably most similar to the original GSEA. It can be run in a single command:

fgseaRes <- fgsea(pathways = examplePathways,  stats = exampleRanks)

Check the vignette for more details: https://www.bioconductor.org/packages/devel/bioc/vignettes/fgsea/inst/doc/fgsea-tutorial.html

ADD COMMENTlink written 11 months ago by igor8.3k
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: 1503 users visited in the last hour