FUNCTIONAL ENRICHMENT ANALYSIS
1
0
Entering edit mode
nia ▴ 30
@nia-12707
Last seen 4.1 years ago

Dear fellows,

I need your suggestion i have around 200 degs which i retrieved after doing its linear modeling , i submit the genes in (genecodis3) tool for the biological and molecular enrichment analysis and i got around 250 GO annotations is it fine i just take top 10 for table (fdr<0.05) (test=hyper geometrical test)

genecodis3 fdr • 1.2k views
ADD COMMENT
0
Entering edit mode
Robert Castelo ★ 3.3k
@rcastelo
Last seen 6 days ago
Barcelona/Universitat Pompeu Fabra

hi,

It's unclear to me to what Bioconductor package are your referring to when you talk about the "genecodis3 tools" as a text-free search among the software packages does not give any hit.

About the general question on how to select enriched GO terms, there are many options, just think for instance that you're talking about the top-10 GO terms without defining what is the metric you use to rank them.

If you're going to select GO terms on the basis of multiple test correction on a one-tailed Fisher's exact test p-value, which is a tricky business given the overlapping nature of GO terms and the consequently lack of independence of the raw p-values assumed by correction methods such as FDR, I would avoid testing GO terms annotated to fewer than 5 and more than 300 genes in your data set. This should alleviate the burden of the multiple test correction, while likely discarding GO terms for which either the representation in your data is too small (< 5 genes) to do a reliable inference, or too large (> 300) to draw useful conclusions. You may of course use different thresholds than the ones I'm proposing.

Regarding the way in which you may rank GO terms, once you've selected them, i'd use the odds-ratio of enrichment.

Cheers,

robert.

ADD COMMENT

Login before adding your answer.

Traffic: 503 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