Correction for multiple comparison by DeSeq2 - treatment1 vs. control, treatment2 vs. control, ... treatment n vs. control
Entering edit mode
ncorak • 0
Last seen 7 weeks ago

Hi everynone!

I am using DeSeq2 to generate a list of differentially expressed genes when different treatments are applied (A vs. con, B vs. con and C vs. con). I am interested what kind of p-value correction method should I use. I know that padj calculated by DeSeq adjusts the value based on the number of tested genes, but should the correction method consider that I also have 3 different comparisons (in fact, for each gene I have 3 tests, one for each comparison)? I would like to emphasize that treatments weren't compared at any point, just treatment vs. control.

Also, I would appreciate it if you could provide examples of scientific papers where this issue was addressed and the correct correction method was applied to a real scientific problem and data.

Thank you all very much!

DESeq2 MultipleComparison r • 179 views
Entering edit mode
Last seen 4 hours ago
United States

I don't typically correct across a few pairwise comparisons, but instead just consider that each results table controls its own FDR.

You can see how often, under the null, would BH give you any false rejections across multiple results tables, it's not a huge concern:

> replicate(10, sum(p.adjust(runif(1000)) < .1))
 [1] 0 0 0 0 0 0 1 0 0 0

Login before adding your answer.

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