# Šidák correction

In statistics, the Šidák correction, or Dunn–Šidák correction, is a method used to counteract the problem of multiple comparisons. It is a simple method to control the familywise error rate. When all null hypotheses are true, the method provides familywise error control that is exact for tests that are stochastically independent, is conservative for tests that are positively dependent, and is liberal for tests that are negatively dependent. It is credited to a 1967 paper [1] by the statistician and probabilist Zbyněk Šidák.[2]

## UsageEdit

• Given m different null hypotheses and a familywise alpha level of ${\displaystyle \alpha }$ , each null hypotheses is rejected that has a p-value lower than ${\displaystyle \alpha _{SID}=1-(1-\alpha )^{\frac {1}{m}}}$ .
• This test produces a familywise Type I error rate of exactly ${\displaystyle \alpha }$  when the tests are independent from each other and all null hypotheses are true. It is less stringent than the Bonferroni correction, but only slightly. For example, for ${\displaystyle \alpha }$  = 0.05 and m = 10, the Bonferroni-adjusted level is 0.005 and the Šidák-adjusted level is approximately 0.005116.
• One can also compute confidence intervals matching the test decision using the Šidák correction by using 100(1 − α)1/m% confidence intervals.

## ProofEdit

The Šidák correction is derived by assuming that the individual tests are independent. Let the significance threshold for each test be ${\displaystyle \alpha _{1}}$ ; then the probability that at least one of the tests is significant under this threshold is (1 - the probability that none of them are significant). Since it is assumed that they are independent, the probability that all of them are not significant is the product of the probabilities that each of them are not significant, or ${\displaystyle 1-(1-\alpha _{1})^{m}}$ . Our intention is for this probability to equal ${\displaystyle \alpha }$ , the significance level for the entire series of tests. By solving for ${\displaystyle \alpha _{1}}$ , we obtain ${\displaystyle \alpha _{1}=1-(1-\alpha )^{1/m}.}$