Question: How to adjust p-values obtained by pairwise comparison of treatments, but I have multiple treatments?
gravatar for dakterl79
2.6 years ago by
dakterl7920 wrote:

My experiment have several treatments (e.g. A, B, and C).

If I have interest in differentially expressed genes (DEGs) over all treatments, ANOVA will be good enough.

However, my major focus is identification of DEGs between two treatments, but came into trouble with adjusting p-values. What will be the best solution for adjusting p-values?

  1. just perform multiple test correction (e.g. p.adjust function of R) using p-values from each pairwise comaprisons (e.g. A vs B)?

  2. collect all p-values from all pairwise comparisons (A vs B, A vs C, B vs C) and run single multiple testing correction and report q-values?

  3. any best solution for this?

Thank you in advance.

ADD COMMENTlink modified 2.6 years ago by dariober11k • written 2.6 years ago by dakterl7920
gravatar for dariober
2.6 years ago by
WCIP | Glasgow | UK
dariober11k wrote:

Hi- Loosely related, I wouldn't use a straight ANOVA applied to each gene independently but rather use the methods implemented packages like edgeR or DESeq since you get more power and better normalization.

If you perform several comparisons (A vs B, A vs C, B vs C) I think in theory you should correct for all the tests you have performed by collecting all the p values from A vs B, A vs C and B vs C and adjust on that long list (using FDR, probably). However, this may be overly stringent since some tests are not independent.

In practice, I think most of us in this situation correct p-values within each comparison, regardless of how many comparisons have been done. So +1 for taking care of this complication!

ADD COMMENTlink written 2.6 years ago by dariober11k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1136 users visited in the last hour