1 3 Multiple Hypothesis Testing
1 3 Multiple Hypothesis Testing
Statistics Refresher
Lecture 3: Multiple Hypothesis Testing
Caroline Uhler (MIT) MITx: Statistics, Computation & Applications Lecture 3 1/9
Some quotes and research findings
Caroline Uhler (MIT) MITx: Statistics, Computation & Applications Lecture 3 2/9
Some quotes and research findings
”Orange cars are less likely to have serious damages that are discovered
only after the purchase.”
Caroline Uhler (MIT) MITx: Statistics, Computation & Applications Lecture 3 2/9
Jelly Beans and Acne
Caroline Uhler (MIT) MITx: Statistics, Computation & Applications Lecture 3 3/9
Problematic of selective inference
http://imgs.xkcd.com/comics/significant.png
Caroline Uhler (MIT) MITx: Statistics, Computation & Applications Lecture 3 4/9
Wonder-syrup
measure 100 variables before and after taking the syrup: weight,
blood pressure, etc.
Caroline Uhler (MIT) MITx: Statistics, Computation & Applications Lecture 3 5/9
Wonder-syrup
measure 100 variables before and after taking the syrup: weight,
blood pressure, etc.
Caroline Uhler (MIT) MITx: Statistics, Computation & Applications Lecture 3 5/9
Wonder-syrup
measure 100 variables before and after taking the syrup: weight,
blood pressure, etc.
Caroline Uhler (MIT) MITx: Statistics, Computation & Applications Lecture 3 5/9
Different protection levels
Caroline Uhler (MIT) MITx: Statistics, Computation & Applications Lecture 3 6/9
Corrections for multiple testing
Bonferroni correction:
Reject H0 when: m · p-value ≤ α
where m is the total number of hypothesis tests performed
Bonferroni correction implies FWER ≤ α
Caroline Uhler (MIT) MITx: Statistics, Computation & Applications Lecture 3 7/9
Corrections for multiple testing
Bonferroni correction:
Reject H0 when: m · p-value ≤ α
where m is the total number of hypothesis tests performed
Bonferroni correction implies FWER ≤ α
Holm-Bonferroni correction:
Sort p-values in increasing order: p(1) ≤ · · · ≤ p(m)
Reject H0 when: (m − i + 1)p(i) ≤ α (more power than Bonferroni)
Holm-Bonferroni correction implies FWER ≤ α
Caroline Uhler (MIT) MITx: Statistics, Computation & Applications Lecture 3 7/9
Corrections for multiple testing
Bonferroni correction:
Reject H0 when: m · p-value ≤ α
where m is the total number of hypothesis tests performed
Bonferroni correction implies FWER ≤ α
Holm-Bonferroni correction:
Sort p-values in increasing order: p(1) ≤ · · · ≤ p(m)
Reject H0 when: (m − i + 1)p(i) ≤ α (more power than Bonferroni)
Holm-Bonferroni correction implies FWER ≤ α
Benjamini-Hochberg correction:
Sort p-values in increasing order: p(1) ≤ · · · ≤ p(m)
Reject H0 when: mp(i) /i ≤ α
Benjamini-Hochberg correction implies FDR ≤ α
Caroline Uhler (MIT) MITx: Statistics, Computation & Applications Lecture 3 7/9
Commonly accepted practice
Caroline Uhler (MIT) MITx: Statistics, Computation & Applications Lecture 3 8/9
References
http://simons.berkeley.edu/talks/yoav-benjamini-2013-12-11a
Caroline Uhler (MIT) MITx: Statistics, Computation & Applications Lecture 3 9/9