Tutorials › Statistical Theory III - Theories of Inference, Multivariate Methods, and Contingency Tables › Power of Bonferroni, Scheffé, and Max-t Procedures

Simultaneous Inference and Power · Tutorial 293 of 972

Power of Bonferroni, Scheffé, and Max-t Procedures

Learn how critical values, contrast correlations, and the scope of simultaneous protection determine power under Bonferroni, Scheffé, and max-t procedures.

Advanced 13 min read

What You'll Learn

  • Express simultaneous contrast statistics using a joint normal vector and a shared variance estimate
  • Calculate exact finite-family power from the joint noncentral t distribution
  • Compare Bonferroni, max-t, and Scheffé critical values and their protection guarantees
  • Distinguish power for at least one rejection from power for every contrast
  • Identify why Scheffé protection can cost power for a small prespecified family
  • Evaluate how correlation among contrasts affects max-t power

Power Depends on Both the Alternative and the Protection Goal

Simultaneous procedures are designed to control the probability of making one or more false rejections across a family of contrasts. That protection affects power: a larger critical value makes rejection harder. But comparing procedures requires more than comparing critical values. One must also specify the family of contrasts, their dependence, and what counts as success under the alternative.

For a specified finite family, a max-t procedure uses the joint distribution of its contrast statistics, whereas Bonferroni uses a bound that does not require their dependence structure. Scheffé protection is broader: it covers every contrast in a contrast space, not just a prespecified finite list. That wider guarantee can be costly when the scientific target is a small finite family.

Definition: Let \(C_1,\ldots,C_m\) be a specified family of contrasts, and let \(R_j\) denote rejection of the null hypothesis for \(C_j\). The any-rejection power at an alternative is \(P(\bigcup_{j=1}^m R_j)\). The all-rejections power is \(P(\bigcap_{j=1}^m R_j)\). These are different operating characteristics; which one matters depends on the scientific objective.

Use the usual exact normal-theory setting: the contrast estimates are jointly normal, and the common error variance is estimated independently by \(s^2\), with \(\nu\) error degrees of freedom. Let \(q\) be the dimension of the contrast space, and write \(v_j\) for the variance factor of contrast \(C_j\), so that \(\operatorname{Var}(\hat C_j)=\sigma^2v_j\). Under an alternative, define the standardized effect \(\delta_j=(C_j-\theta_{j,0})/(\sigma\sqrt{v_j})\), where \(\theta_{j,0}\) is the null value. Let \(R\) be the correlation matrix of the standardized contrast errors.

Theorem: In this setting, the vector of studentized statistics has the representation \(T_j=(Z_j+\delta_j)/\sqrt{W/\nu}\), \(j=1,\ldots,m\), where \(Z\sim N_m(0,R)\), \(W\sim\chi^2_\nu\), and \(Z\) and \(W\) are independent. Thus all the statistics share the same random denominator.

This representation gives an exact route from an alternative to family power. For a common two-sided cutoff \(c\), conditional on \(W=w\), rejection occurs if at least one coordinate satisfies \(|Z_j+\delta_j|>c\sqrt{w/\nu}\). Therefore

$$ P(\text{any rejection}) =\int_0^\infty P\!\left(\max_{1\le j\le m}|Z_j+\delta_j|> c\sqrt{w/\nu}\right) f_{\chi^2_\nu}(w)\,dw. $$

The formula retains the dependence among contrasts and the shared variance estimate. Treating the studentized statistics as independent t variables generally gives the wrong family power.

Critical Values and What They Protect

For two-sided tests at familywise level \(\alpha\), the Bonferroni cutoff is \(t_{\nu,1-\alpha/(2m)}\). The max-t cutoff is the \(1-\alpha\) quantile of \(\max_j |Z_j|/\sqrt{W/\nu}\) under the joint null, using the family’s correlation matrix. Scheffé’s cutoff for the full \(q\)-dimensional contrast space is \(\sqrt{qF_{q,\nu;1-\alpha}}\), where \(F_{q,\nu;1-\alpha}\) is the \(1-\alpha\) quantile of the central \(F_{q,\nu}\) distribution.

These choices have different guarantees. Bonferroni controls the familywise error rate for the specified \(m\) tests, regardless of their dependence. Max-t uses the joint distribution for that finite family and controls its familywise error rate. Scheffé controls the probability of any false rejection over the entire \(q\)-dimensional contrast space. Consequently, a finite-family max-t cutoff can be smaller than the Scheffé cutoff: the maximum over a finite set cannot exceed the maximum over all directions.

Result: For the same finite family and level \(\alpha\), the exact max-t critical value is no greater than the Bonferroni critical value. The finite-family max-t critical value is also no greater than the Scheffé critical value when the finite contrasts lie in the \(q\)-dimensional space protected by Scheffé.

For the first comparison, under the null each statistic has a t distribution with \(\nu\) degrees of freedom. The union bound gives \(P(\max_j |T_j|>c_B)\leq\sum_jP(|T_j|>c_B)=m(\alpha/m)=\alpha\), where \(c_B=t_{\nu,1-\alpha/(2m)}\). The max-t cutoff is the smallest cutoff achieving the specified null probability \(\alpha\), so it cannot be larger than \(c_B\).

For the second comparison, consider a contrast-space estimate represented in \(q\) orthonormal coordinates. The squared maximum absolute studentized contrast over all unit directions is \(qF\), with \(F\sim F_{q,\nu}\) under the null. Restricting the maximum to a finite set of directions cannot increase it. Hence the finite-family maximum has a critical value no larger than the full-space Scheffé critical value. These are ordering results for critical values, not a universal ordering between Bonferroni and Scheffé; their relative sizes depend on \(m\), \(q\), \(\nu\), and \(\alpha\).

Worked Power Comparisons

Worked Example: Comparing Cutoffs for Two Orthogonal Contrasts

Consider two orthogonal contrasts in a two-dimensional contrast space, with known error variance, familywise level \(\alpha=0.05\), and standardized effects \((\delta_1,\delta_2)=(3,0)\). This known-variance normal example is a transparent benchmark for how the protection scope changes power. Here the two standardized statistics are independent unit-variance normal variables.

For Bonferroni, each two-sided test has size \(0.05/2=0.025\), so its cutoff is \(z_{0.9875}\approx2.241\). The probability that the first statistic does not reject is \(\Phi(2.241-3)-\Phi(-2.241-3)=\Phi(-0.759)-\Phi(-5.241)\approx0.2239\). The null second statistic is inside its interval with probability \(0.975\). Independence therefore gives \(P(\text{any rejection})=1-(0.2239)(0.975)\approx0.782\).

For Scheffé, \(q=2\) and the known-variance cutoff is \(\sqrt{\chi^2_{2;0.95}}=\sqrt{5.991}\approx2.448\). The first coordinate’s nonrejection probability is \(\Phi(\sqrt{5.991}-3)-\Phi(-\sqrt{5.991}-3)\approx0.29036\). The second coordinate is inside its interval with probability \(2\Phi(\sqrt{5.991})-1\approx0.9856\). Thus \(P(\text{any rejection})=1-\left[\Phi(\sqrt{5.991}-3)-\Phi(-\sqrt{5.991}-3)\right]\left[2\Phi(\sqrt{5.991})-1\right]\approx0.7138\). In particular, the power for the first contrast alone is \(\Phi(0.552)+\Phi(-5.448)\approx0.710\), not \(0.709\).

For max-t, independence gives the cutoff from \((2\Phi(c)-1)^2=0.95\), because each coordinate must fall inside \([-c,c]\) for no rejection. Thus \(c=\Phi^{-1}((1+\sqrt{0.95})/2)\approx2.236\). The first-coordinate nonrejection probability is approximately \(\Phi(-0.764)-\Phi(-5.236)\approx0.2224\); the second-coordinate nonrejection probability is \(\sqrt{0.95}\). Therefore the any-rejection power is approximately \(1-(0.2224)\sqrt{0.95}\approx0.783\). The finite-family max-t and Bonferroni powers are close here, while Scheffé is less powerful because it protects more directions than these two.

Worked Example: Exact Finite-Sample Power with a Shared t Denominator

Suppose two standardized contrast errors are independent, \(\nu=10\), and the effects are \((1.5,0)\). Consider the any-rejection power for the Bonferroni procedure at \(\alpha=0.05\). Its cutoff is \(c=t_{10,0.9875}\). Let \(W\sim\chi^2_{10}\), and define \(a(w)=c\sqrt{w/10}\). Conditional on \(W=w\), the probability that the first contrast does not reject is \(\Phi(a(w)-1.5)-\Phi(-a(w)-1.5)\). For the second, whose effect is zero, the nonrejection probability is \(2\Phi(a(w))-1\). Conditional independence of the normal numerators gives

$$ P(\text{any rejection})= 1-\int_0^\infty \left[\Phi(a(w)-1.5)-\Phi(-a(w)-1.5)\right] \left[2\Phi(a(w))-1\right] f_{\chi^2_{10}}(w)\,dw. $$

This is an exact finite-sample power expression: the cutoff, degrees of freedom, effects, and distribution in the integral are all specified. The product is valid conditional on \(W\), but it cannot be replaced by a product of unconditional marginal t nonrejection probabilities. Both statistics use the same \(W\), so they are dependent even though their normal numerators are independent.

Worked Example: A Finite Family Versus the Entire Contrast Space

Suppose \(q=3\) and the scientific questions concern only two contrasts, with an estimated common variance and \(\nu\) error degrees of freedom. Bonferroni uses \(t_{\nu,1-\alpha/4}\), since \(m=2\). Max-t uses the \(1-\alpha\) quantile of the maximum of those two correlated absolute t statistics. Scheffé instead uses \(\sqrt{3F_{3,\nu;1-\alpha}}\), protecting all directions in the three-dimensional space.

The max-t test’s null event of no rejection is the event that both specified absolute statistics are at most its cutoff. The Scheffé no-rejection event requires the maximum over every unit direction to be at most its cutoff. Because the two specified directions are contained in the full space, the finite maximum is no larger than the full-space maximum for every realized data set. Its null quantile is therefore no larger than Scheffé’s. At a particular alternative, however, power depends on where the effect vector points and on the correlation of the two selected contrasts; the cutoff ordering alone does not determine power for every possible success criterion.

How to Compare Power Correctly

A fair comparison holds the alternative, significance level, and target family fixed, then computes the relevant rejection probability under each rule. For the any-rejection criterion, the target is \(P(\bigcup_j R_j)\). For all-rejections power, the target is \(P(\bigcap_j R_j)\), which can be much smaller, especially when some contrasts have effects near zero. A procedure can have respectable any-rejection power while rarely rejecting every member of a family.

The contrast correlation matters because max-t calibrates to the distribution of the maximum. Strong positive dependence can make simultaneous extreme values less likely than under independence, often lowering the max-t critical value. Negative or more complex dependence can change the joint tail probability differently. Bonferroni does not adapt to this dependence, while Scheffé pays for protecting all directions rather than just the chosen contrasts.

1
Specify the family and goal.
State whether the family is a finite list or the full contrast space, and whether success means any rejection or all rejections.
2
Set the joint model.
Record the standardized effects, contrast correlation matrix, and error degrees of freedom. Retain the common variance denominator when using t statistics.
3
Calibrate each procedure.
Use the Bonferroni cutoff, the finite-family max-t quantile, or the Scheffé full-space cutoff appropriate to the stated protection.
4
Evaluate power under the same alternative.
Compute the probability of the desired rejection event under the joint distribution; do not compare procedures using different families or success criteria.

Common Pitfall

A common mistake is to say, “Scheffé is more conservative, so it always has lower power than Bonferroni.” This does not follow. Scheffé’s cutoff may be larger or smaller than the Bonferroni cutoff, depending on the number of tests and the contrast-space dimension. Moreover, power depends on the alternative’s direction and on the criterion being measured. The reliable comparison is to calculate the cutoffs for the same setting and evaluate the same joint rejection event under the same alternative.

A second error is to infer family power from one contrast’s marginal power. In the first example, the Scheffé marginal power for the shifted contrast is about \(0.710\), but any-rejection power is about \(0.714\) because the null second contrast can also reject. For all-rejections power, the same null contrast instead makes success unlikely. The event in the probability statement must match the scientific question.

Key takeaway: Bonferroni and max-t protect a specified finite family, while Scheffé protects every contrast in a space. Their power comparison is meaningful only after the family, joint dependence, alternative, and rejection objective have all been made explicit.

Check Your Understanding

Use the definitions and joint-distribution representation above to answer these questions.

  1. Why can two t statistics with independent normal numerators still be dependent?
  2. For a fixed finite family, why is the max-t critical value no larger than the Bonferroni critical value?
  3. What additional set of contrasts does Scheffé protect beyond a specified finite family?
  4. Why does any-rejection power differ from all-rejections power under an alternative that shifts only one contrast?
  5. Which information, in addition to the critical value, is required to compare max-t power under two different alternatives?
Clinical Biostats