Tutorials › Adaptive Clinical Trial Design › Data-Dependent Critical Values

Type I Error Under Data-Dependent Changes · Tutorial 150 of 799

Data-Dependent Critical Values

You will learn to compare an adaptive critical value with the original test’s conditional rejection probability and verify error control point by point in the interim data.

Advanced 12 min read

What You'll Learn

  • Define conditional error as the original test’s rejection probability given interim data.
  • Derive the conditional error for a two-stage normal test with independent increments.
  • Check whether a data-dependent final critical value respects that conditional bound.
  • Translate a threshold on the combined test statistic to one on the final-stage statistic.
  • Diagnose invalid adaptive rules and distinguish conditional checks from an unconditional error calculation.

Why a Critical Value Chosen at Interim Needs a Conditional Check

An interim review can lead investigators to choose a different final critical value depending on the observed data. The difficulty is that the chosen threshold is no longer fixed before the trial begins. To assess its type I error, it is not enough to look only at the nominal significance level attached to each possible threshold. Instead, condition on the interim data and compare the revised test’s rejection probability with the rejection probability the original test had remaining at that same interim result.

That remaining rejection probability is the conditional error. A revised critical value preserves the original test’s error control if, for every relevant interim result, the revised procedure rejects under the null with probability no greater than the original conditional error. This is a pointwise check on the rule for selecting the threshold, not a check performed only after averaging across interim outcomes.

Definition: Let \(X\) denote the interim data, and let \(R_0\) be the rejection event for a prespecified test. The conditional error at interim data \(X=x\) is \(A(x)=P_{H_0}(R_0\mid X=x)\). A revised procedure with rejection event \(R_1\) respects the original conditional-error bound if \(P_{H_0}(R_1\mid X=x)\leq A(x)\) for almost every \(x\) under the null distribution.

A Two-Stage Normal Model

Use a simple one-sided test to make the calculation explicit. Under the null hypothesis, suppose the standardized interim and final-stage increments \(Z_1\) and \(Z_2\) are independent standard normal random variables. The combined statistic is \(Z=(Z_1+Z_2)/\sqrt{2}\), so \(Z\sim N(0,1)\) under the null. The prespecified test rejects when \(Z>1.96\), giving a one-sided significance level of approximately 0.025. This model assumes independent increments, known information weights of one half at each stage, and the stated null distribution. If those assumptions do not describe the trial, the conditional distribution and bound must be recalculated for its actual design.

After observing \(Z_1=x\), the original test rejects exactly when the final increment exceeds a threshold \(q(x)\). Rearranging the original rejection rule gives

$$ \frac{x+Z_2}{\sqrt{2}}>1.96 \quad\Longleftrightarrow\quad Z_2>\sqrt{2}(1.96)-x=q(x). $$

Since \(Z_2\) remains standard normal conditional on \(Z_1=x\), the original conditional error is \(A(x)=1-\Phi(q(x))\), where \(\Phi\) is the standard normal cumulative distribution function. The interim value changes the amount of type I error remaining: a low interim statistic makes the original rejection harder to achieve at the final stage, while a high one makes it easier.

Key result: If the revised procedure rejects when \(Z_2>d(x)\), its conditional rejection probability is \(1-\Phi(d(x))\). Because this tail probability decreases as its threshold increases, a sufficient and exact condition for respecting the original conditional error in this model is \(d(x)\geq q(x)\) for almost every possible interim value \(x\).

Why the Pointwise Bound Preserves Overall Error

The conditional check works because an unconditional rejection probability is the average of the conditional rejection probabilities over the interim data. Under the null, the law of total probability gives

$$ P_{H_0}(R_1) =E_{H_0}\!\left[P_{H_0}(R_1\mid X)\right] \leq E_{H_0}\!\left[A(X)\right] =P_{H_0}(R_0). $$

The last equality is the same law of total probability applied to the original rejection event. Thus, if the original test has type I error at most \(\alpha\), the revised test also has type I error at most \(\alpha\). This short justification depends on the conditional inequality holding across the interim data, not merely on the revised probabilities averaging to a particular value in one calculation.

Worked Example: A Very Unfavorable Interim Result

Suppose the observed interim statistic is \(x=-1.00\). Calculate the original conditional error, then assess two proposed thresholds on \(Z_2\): 3.80 and 3.70.

First calculate the original threshold:

$$ q(-1.00)=\sqrt{2}(1.960)-(-1.00) =\frac{1.960}{\sqrt{0.50}}+1.00 =3.7718586\ldots\approx3.772. $$

The original conditional error is \(1-\Phi(3.772)\), approximately 0.0000810. A proposed final threshold of 3.80 is above 3.772, so its conditional rejection probability, \(1-\Phi(3.80)\approx0.0000724\), is below the bound. It respects the conditional error.

A threshold of 3.70 is below the required 3.772. Its conditional rejection probability is \(1-\Phi(3.70)\approx0.0001078\), which exceeds 0.0000810. That choice does not preserve the original conditional-error bound, even though both proposed thresholds may look stringent in isolation.

Checking a Candidate Threshold at a Given Interim Value

A practical assessment has three parts: identify the original rejection rule, condition on the observed interim statistic to find the remaining null rejection probability, and compute the revised rule’s conditional rejection probability using the same conditional distribution. In the normal model, both probabilities are upper-tail probabilities for \(Z_2\), so comparing their thresholds is sufficient. When the revised test uses a different statistic or distribution, compare the probabilities themselves rather than assuming that a numerical threshold alone determines which procedure is more conservative.

Worked Example: A Threshold That Exceeds the Conditional Bound

Suppose \(Z_1=0.950\). Assess a proposed final-stage threshold of 1.75, then compare it with a threshold of 1.90.

The original conditional threshold and error are

$$ q(0.950)=\sqrt{2}(1.960)-0.950 =2.7718586\ldots-0.950 =1.8218586\ldots\approx1.822, $$

so \(A(0.950)=1-\Phi(1.822)\approx0.03423\). For the proposed threshold 1.75, the conditional rejection probability is \(1-\Phi(1.75)\approx0.04006\). Since \(0.04006>0.03423\), this threshold exceeds the available conditional error and is invalid under the stated bound. With a threshold of 1.90, the probability is \(1-\Phi(1.90)\approx0.02872\), which is below 0.03423, so it respects the bound.

The comparison is conditional on the observed value \(Z_1=0.950\). It does not claim that the observed interim data make the null hypothesis more or less likely; it describes rejection probabilities under the null model.

When the Adaptive Cutoff Is Written on the Combined Statistic

A protocol may specify a data-dependent critical value \(k(x)\) for the combined statistic \(Z\), rather than a threshold \(d(x)\) for \(Z_2\). These forms are equivalent in the two-stage model: conditional on \(Z_1=x\), the rule \(Z>k(x)\) requires \(Z_2>\sqrt{2}k(x)-x\). The original critical value is 1.96. Therefore, a combined-statistic cutoff no smaller than 1.96 gives a final-stage threshold no smaller than the original conditional threshold. A cutoff below 1.96 gives a lower final-stage threshold and exceeds the original conditional error in this model.

Worked Example: Translating an Adaptive Combined-Statistic Cutoff

At \(x=0.950\), suppose the revised procedure rejects when the combined statistic exceeds \(k(x)=1.980\). Translate this cutoff to the final-stage scale and check the conditional bound.

The required final increment is

$$ d(0.950)=\sqrt{2}(1.980)-0.950 =2.8001428\ldots-0.950 =1.8501428\ldots\approx1.850. $$

This exceeds the original conditional threshold \(q(0.950)\approx1.822\). Consequently, the revised conditional rejection probability is \(1-\Phi(1.850)\approx0.0322\), which is below the original bound of approximately 0.03423. The cutoff of 1.980 therefore respects the bound at this interim value. This calculation checks one value; a data-dependent rule must satisfy the condition throughout the range of interim values where it can be applied.

Assumptions, Diagnostics, and Alternatives

The calculations above are valid only under the stated probability model. In particular, the final increment must have the specified conditional distribution given the interim data. Independent standard normal increments are a convenient model, not a universal feature of adaptive trials. Unequal information at the two stages changes the weights in the combined statistic. Estimated variance, discrete outcomes, dependence, or nuisance parameters can also change the conditional distribution. In those settings, derive or bound the conditional rejection probability using the actual design assumptions; do not reuse \(q(x)=\sqrt{2}(1.96)-x\) automatically.

A useful diagnostic is to evaluate the proposed rule over the full range of interim values it may encounter, paying particular attention to values where the adaptive cutoff is lowest. Check that the rule is defined at boundary values and that any numerical approximations to tail probabilities are accurate enough to determine which side of the bound they fall on. If nuisance parameters are present, the bound must be valid for the relevant null parameter values; a calculation at only one convenient parameter value may not establish control for a composite null.

The main alternative is to keep the original, prespecified critical value. That avoids recalculating conditional bounds, but it may not permit the desired design changes. Another approach is to specify a valid adaptive procedure in advance and evaluate its unconditional type I error under the full null model. Conditional-error checking instead provides a direct way to judge whether a post-interim choice spends more rejection probability than the original test had remaining at each interim result. These approaches answer related but distinct questions: a nominal cutoff by itself does not show that a data-dependent selection rule is valid.

A common mistake is to argue, “The revised cutoff is stringent, and the average rejection rate appears to be at most \(\alpha\), so the change is safe.” The correction is to compute the original conditional error and compare it with the revised conditional rejection probability for each relevant interim result. A threshold can be numerically large yet still be too permissive at a particular interim value, and an average calculation alone does not establish the pointwise conditional bound required here.

The method also has limits. It does not make an incorrect conditional model correct, nor does it by itself validate every aspect of an adaptive design. Its conclusion is specifically about preserving the stated null rejection bound under the assumed conditional distribution. If a design changes the statistic, endpoint, or hypotheses, the revised rejection event must be expressed precisely before its conditional probability can be assessed.

1
Write the original rejection rule.
State the test statistic, null distribution, and original critical value before using interim data.
2
Condition on the interim result.
Use the model to find the original test’s conditional rejection probability at that result.
3
Evaluate the revised rule under the same condition.
Compute its conditional rejection probability or compare thresholds when the same monotone distribution applies.
4
Check the whole adaptive rule.
Verify the inequality for all relevant interim values and account for any nuisance parameters or changed distributions.
Key takeaway: A critical value chosen using interim data is assessed against the original test’s remaining conditional error. Under the two-stage normal model, a final-increment threshold must be at least \(q(x)=\sqrt{2}(1.96)-x\) at interim value \(x\) to preserve that bound.

Check Your Understanding

Use the two-stage normal model above, with original combined-statistic cutoff 1.96. Solutions are included after the questions.

  1. At \(Z_1=0\), calculate the original final-increment threshold \(q(0)\). Is a proposed threshold of 2.80 valid under the conditional-error bound?
  2. At \(Z_1=1.20\), calculate \(q(1.20)\). Would a final-increment cutoff of 1.50 respect the bound?
  3. Explain why checking only the average rejection probability of a data-dependent rule does not establish the pointwise conditional-error inequality.
  4. If a revised procedure uses combined-statistic cutoff \(k(x)=2.00\), what is its conditional final-increment threshold at \(x=0.50\), and does it exceed the original conditional threshold?

Solutions. (1) \(q(0)=\sqrt{2}(1.96)=2.7718586\ldots\), approximately 2.772. Since 2.80 is greater than 2.772, its conditional rejection probability is smaller than the original bound, so it is valid at \(x=0\). (2) \(q(1.20)=2.7718586\ldots-1.20=1.5718586\ldots\), approximately 1.572. A cutoff of 1.50 is smaller, so its upper-tail rejection probability exceeds the conditional bound. (3) An average can conceal interim values at which the revised conditional rejection probability exceeds the original conditional error; the pointwise inequality must be checked over the relevant interim values. (4) The revised threshold is \(\sqrt{2}(2.00)-0.50=2.3284271\ldots\), approximately 2.328. The original threshold is \(q(0.50)=2.7718586\ldots-0.50=2.2718586\ldots\), approximately 2.272. The revised threshold is larger, so it respects the conditional bound at \(x=0.50\).

Clinical Biostats