Tutorials › Adaptive Clinical Trial Design › Stagewise Error Rates and Final Test Calibration

The Combination Test Principle · Tutorial 213 of 799

Stagewise Error Rates and Final Test Calibration

Learn to distinguish stage-specific thresholds from the overall type I error level and calculate a two-stage weighted inverse-normal test’s final rejection rule.

Advanced 10 min read

What You'll Learn

  • State the assumptions that give an inverse-normal combination statistic an exact standard normal null distribution.
  • Derive the final rejection cutoff for a two-stage weighted inverse-normal test.
  • Calculate a combined statistic and p-value from two stagewise p-values.
  • Find the stage-two evidence needed to cross the final combination boundary.
  • Distinguish a final combination test from testing each stage separately.
  • Identify when discrete or dependent stagewise p-values require different calibration.

Why the Final Cutoff Is Not a Stagewise Cutoff

A two-stage adaptive trial produces evidence at two points in time. Each stage may have its own nominal significance threshold, but those thresholds do not automatically define the type I error rate of the final test. The final decision depends on the rule applied to the evidence across stages: for example, whether success means crossing a cutoff at either stage, or whether the two stagewise p-values are combined and compared with one final cutoff.

This distinction matters because type I error belongs to the whole decision rule. A test that rejects whenever either of two independent stagewise p-values is below 0.05 does not have an overall type I error of 0.05. By contrast, a properly calibrated combination test uses a single final threshold chosen to control the probability of rejecting the global null across the two stages.

Definition: A stagewise significance threshold is a cutoff applied to evidence from one stage, such as rejecting at that stage when \(p_i \leq \alpha_i\). The final combination-test threshold is the cutoff applied to a statistic formed from the stages’ evidence, chosen so that the probability of final rejection under the null is the prespecified overall level \(\alpha\), under the assumptions used for calibration.

Two-Stage Inverse-Normal Calibration

Consider a one-sided test with null hypothesis \(H_0\) and alternative hypothesis \(H_a\) in the prespecified direction. Let \(p_1\) and \(p_2\) be the stagewise p-values. For the exact normal calibration below, assume that under \(H_0\), each \(p_i\) has the continuous Uniform(0,1) distribution and that \(p_1\) and \(p_2\) are independent. The test procedures and direction must be specified so that small p-values consistently indicate evidence for \(H_a\). These assumptions commonly correspond to valid continuous stagewise tests based on independent data from the two stages.

Transform each stagewise p-value to a normal score using \(Z_i=\Phi^{-1}(1-p_i)\), where \(\Phi\) is the standard normal cumulative distribution function. Let \(w_1\) and \(w_2\) be fixed, prespecified nonnegative weights satisfying \(w_1^2+w_2^2=1\). The two-stage inverse-normal combination statistic is \(Z_{\mathrm{comb}}=w_1Z_1+w_2Z_2\).

Theorem: Under the stated assumptions, \(Z_1\) and \(Z_2\) are independent standard normal random variables. Therefore \(Z_{\mathrm{comb}}\) has a standard normal distribution under \(H_0\). For an upper-tailed test at overall level \(\alpha\), reject \(H_0\) if \(Z_{\mathrm{comb}} \geq z_{1-\alpha}\), where \(z_{1-\alpha}=\Phi^{-1}(1-\alpha)\). The combined p-value is \(1-\Phi(Z_{\mathrm{comb}})\).

The normalization of the weights is essential to this reference distribution. The weights may reflect a chosen combination rule, but once selected they must be fixed independently of the observed stage results for this calibration. The final cutoff is determined by the overall \(\alpha\), not by applying each \(w_i\) or each stagewise threshold as a separate significance level.

Why the Calibration Works

The key facts follow directly from the assumed null distributions. If \(p_i\) is Uniform(0,1), then for any real \(z\),

$$ P\{\Phi^{-1}(1-p_i)\leq z\} =P\{p_i\geq 1-\Phi(z)\} =\Phi(z). $$

Thus \(Z_i\) is standard normal. Independence of the p-values implies independence of their transformed scores. A weighted sum of independent normal variables is normal; its mean is zero and its variance is \(w_1^2+w_2^2=1\). Hence \(Z_{\mathrm{comb}}\sim N(0,1)\) under \(H_0\). The upper-tail probability at \(z_{1-\alpha}\) is \(\alpha\), which gives the stated overall type I error.

This argument also shows why exact uniformity is not a cosmetic assumption. A discrete or conservative p-value is generally not exactly Uniform(0,1); its transformed score need not be standard normal. If independent stagewise p-values are only superuniform—that is, \(P(p_i\leq t)\leq t\) for every \(t\) between zero and one—the normal-score combination can be conservative under suitable independent-stage conditions, but the exact standard-normal claim no longer follows from the argument above. Its actual null calibration should be established for the p-value distributions and decision rule being used.

Worked Example: Combining Two Nonsignificant Stagewise Results

Suppose a hypothetical two-stage trial uses one-sided stagewise p-values \(p_1=0.10\) and \(p_2=0.12\). Assume each is exactly uniform under \(H_0\), the stagewise p-values are independent, and the prespecified weights are \(w_1=0.6\) and \(w_2=0.8\). These weights satisfy \(0.6^2+0.8^2=1\). Test at overall \(\alpha=0.05\).

The normal scores are \(Z_1=\Phi^{-1}(0.90)\approx1.2816\) and \(Z_2=\Phi^{-1}(0.88)\approx1.1750\). The combined statistic is

$$ Z_{\mathrm{comb}} =(0.6)(1.2816)+(0.8)(1.1750) =0.76896+0.94000 =1.70896\approx1.709. $$

The one-sided final critical value is \(z_{0.95}\approx1.6449\). Since \(1.709>1.6449\), reject \(H_0\) at the overall 0.05 level. The combined p-value is \(1-\Phi(1.70896)\approx0.0437\). Neither stagewise p-value is at most 0.05, yet the prespecified weighted combination crosses the final boundary. The result is a rejection by the combination rule, not by either stage’s separate threshold.

Stagewise Thresholds Define Different Rules

To see why threshold labels alone do not determine overall error, consider a different rule: reject if either of two independent, exactly uniform stagewise p-values is at most 0.05. Under \(H_0\), the probability of no rejection is the probability that both exceed 0.05. Independence gives

$$ P(\text{no rejection}\mid H_0) =(1-0.05)^2=0.9025, \qquad P(\text{at least one rejection}\mid H_0) =1-0.9025=0.0975. $$

So applying a 0.05 cutoff at each stage and rejecting at either one yields a 0.0975 overall type I error in this example, not 0.05. One alternative is to use 0.025 at each stage with the same “either stage” rule. Its overall type I error is \(1-(1-0.025)^2=0.049375\), which is slightly below 0.05. This is a separate multiple-testing rule; its thresholds should not be confused with the final cutoff of a combination statistic.

A trial could also have a formal interim efficacy boundary that permits early rejection and a later combination test for participants who continue. In that case, the global decision rule includes both the early-rejection event and the later rejection event. The final combination cutoff alone does not automatically calibrate the probability of rejection across both opportunities. The early and final parts must be calibrated as one rule.

Worked Example: Finding the Stage-Two Boundary

Use the same weights \(w_1=0.6\), \(w_2=0.8\), and overall one-sided level \(\alpha=0.05\). Suppose the first-stage p-value is \(p_1=0.10\). Find how small the second-stage p-value must be for the final combination test to reject.

The first-stage normal score is \(Z_1=\Phi^{-1}(0.90)\approx1.2816\), and the final cutoff is \(1.6449\). The rejection condition is

$$ 0.6(1.2816)+0.8Z_2\geq1.6449. $$

Solving for \(Z_2\) gives

$$ Z_2\geq \frac{1.6449-0.6(1.2816)}{0.8} =\frac{1.6449-0.76896}{0.8} \approx1.0949. $$

Because \(Z_2=\Phi^{-1}(1-p_2)\), this is equivalent to \(p_2\leq1-\Phi(1.0949)\approx0.1368\). For instance, \(p_2=0.13\) gives \(Z_2\approx1.1264\) and \(Z_{\mathrm{comb}}\approx0.76896+0.90112=1.6701\), so the test rejects. With \(p_2=0.14\), \(Z_2\approx1.0803\) and \(Z_{\mathrm{comb}}\approx0.76896+0.86424=1.6332\), so it does not reject. The stage-two boundary depends on the first-stage result: it is not a single fixed stage-two significance threshold for every possible \(p_1\).

Applying the Rule Reliably

A calibrated decision requires more than plugging two p-values into a formula. Before analysis, specify the one-sided direction, the stagewise tests, the combination weights, the overall \(\alpha\), and whether any separate interim rejection rule is allowed. Check that the two stagewise p-values have the null properties assumed for the chosen calibration and that their dependence structure supports the reference distribution.

1
Verify the inputs.
Confirm that both p-values test the same prespecified direction and that the null assumptions give the intended joint calibration.
2
Apply the fixed transformation and weights.
Calculate \(Z_i=\Phi^{-1}(1-p_i)\) for each stage and form \(Z_{\mathrm{comb}}=w_1Z_1+w_2Z_2\), with \(w_1^2+w_2^2=1\).
3
Compare with the global boundary.
Reject for \(Z_{\mathrm{comb}}\geq z_{1-\alpha}\), or report the corresponding combined p-value \(1-\Phi(Z_{\mathrm{comb}})\).
4
Account for any additional rejection opportunity.
If the design also permits early rejection, evaluate the full decision rule rather than assuming the final combination cutoff controls the error across all stages.

Diagnostics include checking that the weights square to one, that the final cutoff matches the intended one-sided overall level, and that the combination statistic and combined p-value lead to the same decision. If a stagewise p-value is discrete, heavily rounded, or computed from data that overlap with another stage, the exact independent-uniform derivation may not apply. In those cases, use a calibration justified for the actual null distributions and dependence rather than treating the normal reference as automatic.

Common Mistakes and Limits

  • Calling each stage’s cutoff the overall level. A 0.05 cutoff at each of two opportunities does not generally produce a 0.05 overall error rate. State the complete rejection rule and evaluate its null rejection probability.
  • Comparing the combined statistic with a stagewise cutoff. The statistic \(Z_{\mathrm{comb}}\) is compared with \(z_{1-\alpha}\) for the chosen overall level, not with a cutoff selected independently for either stage.
  • Assuming every valid p-value is exactly uniform. Validity may mean only that a p-value is conservative. Exact normal-score calibration requires exact uniformity and independence under the stated setup; other cases require appropriate null calibration.
  • Changing weights after seeing the results. The standard calibration above assumes prespecified weights. Choosing them in response to observed evidence changes the decision procedure and is not covered by this derivation.

The inverse-normal rule is one way to combine stagewise evidence, not a universal description of every adaptive trial. Its exact cutoff follows from its specified assumptions and statistic. Alternative combination rules can have different statistics and reference distributions, so their final thresholds must be calibrated for those rules rather than borrowed from the inverse-normal test.

Check Your Understanding

Answer each question using the stated rule and assumptions.

  1. For weights \(w_1=0.6\), \(w_2=0.8\), overall one-sided \(\alpha=0.05\), and \(p_1=0.10\), what is the approximate maximum \(p_2\) that still leads to rejection?
  2. Two independent, exactly uniform null p-values are each tested at 0.05, with rejection if either passes its threshold. What is the overall type I error?
  3. Why does exact standard-normal calibration of the inverse-normal statistic require more than merely valid stagewise p-values?
  4. With two independent, exactly uniform null p-values and a rule that rejects if either is at most 0.025, calculate the overall type I error.

Solutions. (1) The boundary is \(Z_2\geq(1.6449-0.6(1.2816))/0.8\approx1.0949\), or \(p_2\leq1-\Phi(1.0949)\approx0.1368\). (2) The error is \(1-(0.95)^2=0.0975\). (3) A valid p-value may be discrete or conservative rather than exactly uniform, so its inverse-normal transform need not be standard normal; dependence can also invalidate the weighted-sum reference distribution. (4) The error is \(1-(0.975)^2=0.049375\), approximately 0.0494.

Clinical Biostats