Šidák correction for t-test

One of the application of Student's t-test is to test the location of one sequence of independent and identically distributed random variables. If we want to test the locations of multiple sequences of such variables, Šidák correction should be applied in order to calibrate the level of the Student's t-test. Moreover, if we want to test the locations of nearly infinitely many sequences of variables, then Šidák correction should be used, but with caution. More specifically, the validity of Šidák correction depends on how fast the number of sequences goes to infinity.

Introduction

Suppose we are interested in m different hypotheses, ${\displaystyle H_{1},...,H_{m}}$, and would like to check if all of them are true. Now the hypothesis test scheme becomes

${\displaystyle H_{null}}$: all of ${\displaystyle H_{i}}$ are true;

${\displaystyle H_{alternative}}$: at least one of ${\displaystyle H_{i}}$ is false.

Let ${\displaystyle \alpha }$ be the level of this test (the type-I error), that is, the probability that we falsely reject ${\displaystyle H_{null}}$ when it is true.

We aim to design a test with certain level ${\displaystyle \alpha }$.

Suppose when testing each hypothesis ${\displaystyle H_{i}}$, the test statistic we use is ${\displaystyle t_{i}}$.

If these ${\displaystyle t_{i}}$'s are independent, then a test for ${\displaystyle H_{null}}$ can be developed by the following procedure, known as Šidák correction.

Step 1, we test each of m null hypotheses at level ${\displaystyle 1-(1-\alpha )^{\frac {1}{m}}}$.

Step 2, if any of these m null hypotheses is rejected, we reject ${\displaystyle H_{null}}$.

Finite case

For finitely many t-tests, suppose ${\displaystyle Y_{ij}=\mu _{i}+\epsilon _{ij},i=1,...,N,j=1,...,n,}$ where for each i, ${\displaystyle \epsilon _{i1},...,\epsilon _{in}}$ are independently and identically distributed, for each j ${\displaystyle \epsilon _{1j},...,\epsilon _{Nj}}$ are independent but not necessarily identically distributed, and ${\displaystyle \epsilon _{ij}}$ has finite fourth moment.

Our goal is to design a test for ${\displaystyle H_{null}:\mu _{i}=0,\forall i=1,...,N}$ with level α. This test can be based on the t-statistic of each sequences, that is,

${\displaystyle t_{i}={\frac {{\bar {Y}}_{i}}{S_{i}/{\sqrt {n}}}},}$

where:

${\displaystyle {\bar {Y}}_{i}={\frac {1}{n}}\sum _{j=1}^{n}Y_{ij},\qquad S_{i}^{2}={\frac {1}{n}}\sum _{j=1}^{n}(Y_{ij}-{\bar {Y}}_{i})^{2}.}$

Using Šidák correction, we reject ${\displaystyle H_{null}}$ if any of the t-tests based on the t-statistics above reject at level ${\displaystyle 1-(1-\alpha )^{\frac {1}{N}}.}$ More specifically, we reject ${\displaystyle H_{null}}$ when

${\displaystyle \exists i\in \{1,\ldots ,N\}:|t_{i}|>\zeta _{\alpha ,N},}$

where

${\displaystyle P(|Z|>\zeta _{\alpha ,N})=1-(1-\alpha )^{\frac {1}{N}},\qquad Z\sim N(0,1)}$

The test defined above has asymptotic level α, because

${\displaystyle {\begin{aligned}{\text{level}}&=P_{null}\left({\text{reject }}H_{null}\right)\\&=P_{null}\left(\exists i\in \{1,\ldots ,N\}:|t_{i}|>\zeta _{\alpha ,N}\right)\\&=1-P_{null}\left(\forall i\in \{1,\ldots ,N\}:|t_{i}|\leq \zeta _{\alpha ,N}\right)\\&=1-\prod _{i=1}^{N}P_{null}\left(|t_{i}|\leq \zeta _{\alpha ,N}\right)\\&\to 1-\prod _{i=1}^{N}P\left(|Z_{i}|\leq \zeta _{\alpha ,N}\right)&&Z_{i}\sim N(0,1)\\&=\alpha \end{aligned}}}$

Infinite case

In some cases, the number of sequences, ${\displaystyle N}$, increase as the data size of each sequences, ${\displaystyle n}$, increase. In particular, suppose ${\displaystyle N(n)\rightarrow \infty {\text{ as }}n\rightarrow \infty }$. If this is true, then we will need to test a null including infinitely many hypotheses, that is

${\displaystyle H_{null}:{\text{ all of }}H_{i}{\text{ are true, }}i=1,2,....}$

To design a test, Šidák correction may be applied, as in the case of finitely many t-test. However, when ${\displaystyle N(n)\rightarrow \infty {\text{ as }}n\rightarrow \infty }$, the Šidák correction for t-test may not achieve the level we want, that is, the true level of the test may not converges to the nominal level ${\displaystyle \alpha }$ as n goes to infinity. This result is related to high-dimensional statistics and is proven by Fan, Hall and Yao (2007).[1] Specifically, if we want the true level of the test converges to the nominal level ${\displaystyle \alpha }$, then we need a restraint on how fast ${\displaystyle N(n)\rightarrow \infty }$. Indeed,

• When all of ${\displaystyle \epsilon _{ij}}$ have distribution symmetric about zero, then it is sufficient to require ${\displaystyle \log N=o(n^{1/3})}$ to guarantee the true level converges to ${\displaystyle \alpha }$.
• When the distributions of ${\displaystyle \epsilon _{ij}}$ are asymmetric, then it is necessary to impose ${\displaystyle \log N=o(n^{1/2})}$ to ensure the true level converges to ${\displaystyle \alpha }$.
• Actually, if we apply bootstrapping method to the calibration of level, then we will only need ${\displaystyle \log N=o(n^{1/3})}$ even if ${\displaystyle \epsilon _{ij}}$ has asymmetric distribution.

The results above are based on Central Limit Theorem. According to Central Limit Theorem, each of our t-statistics ${\displaystyle t_{i}}$ possesses asymptotic standard normal distribution, and so the difference between the distribution of each ${\displaystyle t_{i}}$ and the standard normal distribution is asymptotically negligible. The question is, if we aggregate all the differences between the distribution of each ${\displaystyle t_{i}}$ and the standard normal distribution, is this aggregation of differences still asymptotically ignorable?

When we have finitely many ${\displaystyle t_{i}}$, the answer is yes. But when we have infinitely many ${\displaystyle t_{i}}$, the answer some time becomes no. This is because in the latter case we are summing up infinitely many infinitesimal terms. If the number of the terms goes to infinity too fast, that is, ${\displaystyle N(n)\rightarrow \infty }$ too fast, then the sum may not be zero, the distribution of the t-statistics can not be approximated by the standard normal distribution, the true level does not converges to the nominal level ${\displaystyle \alpha }$, and then the Šidák correction fails.

See also

• Šidák correction
• Multiple comparisons
• Bonferroni correction
• Familywise error rate
• Closed testing procedure

Notes

1. Fan, Jianqing; Hall, Peter; Yao, Qiwei (2007). "To How Many Simultaneous Hypothesis Tests Can Normal, Student's t or Bootstrap Calibration Be Applied". Journal of the American Statistical Association. 102 (480): 1282–1288. arXiv:math/0701003. doi:10.1198/016214507000000969.CS1 maint: ref=harv (link)

References

• Fan, Jianqing; Hall, Peter; Yao, Qiwei (2007). "To How Many Simultaneous Hypothesis Tests Can Normal, Student's t or Bootstrap Calibration Be Applied". Journal of the American Statistical Association. 102 (480): 1282–1288. arXiv:math/0701003. doi:10.1198/016214507000000969.CS1 maint: ref=harv (link)