limma symbols2indices from gmt files versus rdata files
1
0
Entering edit mode
biominer ▴ 10
@biominer-7701
Last seen 7.0 years ago
European Union

Hello,

I'd like to do an enrichment using romer/roast/camera and my input are HGNC official human gene symbols.

For the symbols2indices() function I need a list of vectors for the reference gene sets, either from the rdata files provided by you or generated from the gmt-files from Broad institute.

As discussed previously in this post:

Romer and symbols2indices query

you had to curate the gene symbols in the gmt-files heavily to make them consistent, which is why you nicely provide these downloads. Does this all apply to the entrez-id gmts only or also to the gene-symbol gmts downloadable from Broad?

I wonder what is the better option, to use the gene symbol gmts or to map my ids to entrez first and use the rdata files. To my experience enrichments can be pretty vulnerable to how good the mappings work.

Do you have any experience with the gene-symbol version of the gmt-files? Or what would you recommend to do starting from HGNC-symbols (or ensembl gene ids)?

Best,

Maik

 

 

limma romer roast • 1.5k views
ADD COMMENT
2
Entering edit mode
@gordon-smyth
Last seen 3 hours ago
WEHI, Melbourne, Australia

Either way is fine.

Mapping official human symbols to Entrez Gene IDs is easy and pretty much one-to-one. That's what I would do.

Or working from the Broad Institute gene symbol gmts is also fine. It was different 6 years ago, but these days the symbol gmts from MSigDb contain official human symbols in a consistent format.

ADD COMMENT

Login before adding your answer.

Traffic: 532 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6