Should I adjust my q-values in edgeR/DESeq2 in my example?
Entering edit mode
SmallChess ▴ 10
Last seen 11 months ago

Let's say I have 3 groups for differential gene analysis. I want to make two pairwise comparisons, A vs B and B vs C. I have 1000 genes.

I can provide contrasts to edgeR/DESeq2 for each of the pairwise comparison, thus I will have two contrasts.

The first time I test A vs B, I will have a list of q-values for my 1000 genes. Then next time I run B vs C, I will have a new set of q-values for my 1000 genes.

Here's my question: the q-values are adjustment for my 1000 genes. However, I'm actually testing 2000 genes (A vs B and B vs C) simultaneously. Should I further adjust my q-values?



edger deseq2 • 1.2k views
Entering edit mode
Last seen 4 hours ago
WEHI, Melbourne, Australia

You could do an F-like-test for each gene, then you would only be doing 1000 tests. In edgeR, you simply give both pairwise comparisons to glmLRT() or glmQLFTest() simultaneously and it will do a combined test.

Alternatively, if you want to keep the two comparisons separate, see ?decideTests for the various ways in which you can further adjust the p-values across comparisons. Usually people just adjust each comparison separately.


Login before adding your answer.

Traffic: 480 users visited in the last hour
Help About
Access RSS

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

Powered by the version 2.3.6