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

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
ADD COMMENT
3
Entering edit mode
@gordon-smyth
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.

ADD COMMENT

Login before adding your answer.

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