Tutorials › Group Sequential Methods and Sample Size Re-Estimation › Conditional Error Across Multiple Interim Looks

Type I Error Control in Re-Estimation · Tutorial 450 of 699

Conditional Error Across Multiple Interim Looks

You will derive and apply a backward recursion for the conditional probability of eventual rejection after multiple interim looks.

Advanced 14 min read

What You'll Learn

  • State the canonical null distribution and its assumptions for multiple looks.
  • Derive the transition distribution between successive standardized statistics.
  • Calculate remaining Type I error recursively from future crossing probabilities.
  • Distinguish the conditional error at a particular interim history from the design’s overall Type I error.
  • Diagnose when a recursion is invalid or its conditioning history is incomplete.
  • Compare recursive calculation with direct multivariate-normal probability evaluation.

Why Conditional Error Must Be Updated Across Looks

After an interim analysis, the probability of a later rejection depends on what has already been observed. If the trial continues through several more looks, that probability includes all remaining chances to cross a rejection boundary—not just the chance of crossing at the final analysis. The conditional error is this remaining probability under \(H_0\), conditional on the interim data and on having reached the interim without rejection.

This calculation is useful when a design is being assessed or changed after interim data become available. It makes explicit the Type I error still available at each history. The introductory material on conditional error establishes the principle; here the focus is its recursive calculation across multiple future looks.

Canonical Model and Conditional Transitions

Assume a one-sided test with \(K\) planned looks at strictly increasing information times \(0<t_1<\cdots<t_K=1\). Under the null model, the standardized statistics \(Z_1,\ldots,Z_K\) have the canonical joint distribution: each has a standard normal marginal distribution, and \(\operatorname{Corr}(Z_i,Z_j)=\sqrt{t_i/t_j}\) for \(i<j\). The design rejects at look \(i\) if \(Z_i\ge c_i\), where \(c_i\) is the prespecified boundary. The information times and boundaries are fixed for this calculation.

This model is appropriate when the standardized score process has independent increments under \(H_0\), with information indexing its variance. It is also necessary that the stated information times, statistic, and boundaries describe the design being evaluated. The canonical joint distribution is not automatic for arbitrary statistics or data-dependent information; if its assumptions fail, the recursion below need not give valid conditional error.

Definition: At look \(j\), after observing \(Z_j=z<c_j\) and having had no earlier rejection, the conditional error is the null probability of rejection at one or more of looks \(j+1,\ldots,K\), conditional on the interim history. Under the canonical model and the fixed-boundary rule, the future conditional distribution depends on that history through \(Z_j=z\).

For adjacent looks, the canonical model gives the transition \(Z_{i+1}\mid Z_i=z\sim N\!\left(\sqrt{t_i/t_{i+1}}\,z,\;1-t_i/t_{i+1}\right)\). Thus the transition density is normal with mean \(\sqrt{t_i/t_{i+1}}\,z\) and variance \(1-t_i/t_{i+1}\). This Markov property is what makes a backward recursion possible: once the current statistic is known, earlier statistics do not provide additional information about future statistics under this model.

Backward Recursion for Remaining Error

Write \(E_j(z)\) for the conditional probability of rejection at some later look, given that the trial has reached look \(j\), has not rejected, and \(Z_j=z\). Set \(E_K(z)=0\), because there are no looks after \(K\). For \(j<K\), let \(f_{j+1}(x\mid z)\) be the density of the transition distribution from look \(j\) to \(j+1\). The recursion is

$$ E_j(z) = \int_{c_{j+1}}^\infty f_{j+1}(x\mid z)\,dx + \int_{-\infty}^{c_{j+1}} E_{j+1}(x)f_{j+1}(x\mid z)\,dx . $$

The first integral is the probability of rejecting immediately at the next look. In the second, the next statistic falls below its boundary, so the trial continues; the remaining error from that point is \(E_{j+1}(x)\). These cases are disjoint and cover all possible outcomes at look \(j+1\). Starting at \(j=K-1\) and working backwards therefore accounts for rejection at every future look without counting any path more than once.

The design’s overall Type I error is different from a conditional error at a realized interim value. It is the probability, under \(H_0\), of crossing at any look from the start. Conditional error is calculated at a particular history, after conditioning on survival and the observed statistic. In an error-controlled redesign, the relevant conditional budget is the original design’s conditional error at that history; it is not automatically the overall \(\alpha\).

Worked Example: One Remaining Look and a Correct Tail Calculation

Consider a two-look canonical design with \(t_1=0.60\), \(t_2=1\), and final boundary \(c_2=2\). Suppose the first look has not rejected and \(Z_1=1\). Calculate the conditional error.

There is only one remaining look, so the recursion reduces to the transition tail probability. Under \(H_0\), \(Z_2\mid Z_1=1\sim N(\sqrt{0.60},0.40)\). Therefore,

$$ E_1(1) = P(Z_2\ge 2\mid Z_1=1) = 1-\Phi\!\left(\frac{2-\sqrt{0.60}}{\sqrt{0.40}}\right). $$

The conditional mean is \(\sqrt{0.60}\approx0.77460\), and the conditional standard deviation is \(\sqrt{0.40}\approx0.63246\). Substitution gives \((2-\sqrt{0.60})/\sqrt{0.40}\approx1.93753\), so the upper-tail probability is approximately \(0.02634\), or 2.634%. This is the probability of crossing the final boundary under \(H_0\), conditional on the stated first-look result; it is not the probability that \(H_0\) is true.

Worked Example: Recursing Across Two Future Looks

Suppose there are three looks at \(t_1=0.40\), \(t_2=0.70\), and \(t_3=1\), with boundary 2 at looks 2 and 3. At look 1, the observed statistic is \(Z_1=0.50\), below its boundary. Set up the conditional error for rejection at look 2 or look 3.

First, \(Z_2\mid Z_1=0.50\) has mean \(\sqrt{0.40/0.70}(0.50)\approx0.37796\) and variance \(1-0.40/0.70=3/7\), with standard deviation approximately \(0.65465\). The chance of rejecting at look 2 is \(1-\Phi((2-0.37796)/0.65465)\), where the standardized boundary is approximately \(2.478\).

If look 2 does not reject and its observed value is \(x<2\), then \(Z_3\mid Z_2=x\sim N(\sqrt{0.70}\,x,0.30)\). The probability of a final rejection from that state is \(1-\Phi((2-\sqrt{0.70}\,x)/\sqrt{0.30})\). Substituting both transition distributions into the recursion gives the exact conditional error:

$$ \begin{aligned} E_1(0.50) ={}&1-\Phi\!\left(\frac{2-\sqrt{0.40/0.70}(0.50)} {\sqrt{1-0.40/0.70}}\right)\\ &+\int_{-\infty}^{2} \left[ 1-\Phi\!\left(\frac{2-\sqrt{0.70}\,x}{\sqrt{0.30}}\right) \right] \frac{1}{\sqrt{3/7}}\, \phi\!\left(\frac{x-\sqrt{0.40/0.70}(0.50)}{\sqrt{3/7}}\right)\,dx . \end{aligned} $$

The first term covers rejection at look 2. The integral covers all look-2 values below 2, weighting each value’s chance of final rejection by its conditional density. This one-dimensional integral can be evaluated by numerical integration; retaining it as an integral also specifies the conditional error exactly, without rounding a numerical approximation. Omitting the integral would calculate only the chance of crossing at look 2 and would miss rejection at the final look.

Worked Example: Updating at the Last Interim Look

In the same three-look setting, suppose the trial reaches look 2 without rejection and observes \(Z_2=1\). The only remaining look is at \(t_3=1\), and its boundary is 2. Calculate the conditional error from this updated history.

The transition distribution is \(Z_3\mid Z_2=1\sim N(\sqrt{0.70},0.30)\). Thus the standardized boundary is \((2-\sqrt{0.70})/\sqrt{0.30}\approx(2-0.83666)/0.54772\approx2.124\).

$$ E_2(1) = 1-\Phi\!\left(\frac{2-\sqrt{0.70}}{\sqrt{0.30}}\right) \approx 0.01684 . $$

The remaining Type I error is about 1.684% for this look-2 history. It is smaller than the conditional error at a different history only if that history and its future crossing chances imply a smaller probability; conditional error is history-specific, not a single amount carried unchanged from look to look.

Calculation Checks, Alternatives, and Limitations

A practical calculation proceeds backwards from the last look. At each step, compute the probability of crossing the next boundary, then integrate the conditional error beyond that look over values that do not cross. Numerical integration can be used when the recursion has no closed form. The resulting conditional error is between 0 and 1, and it cannot be smaller than the probability of crossing at the immediate next look, because that probability is included in the total.

A direct alternative is to calculate a multivariate-normal probability for the union of all future crossing events, conditional on the interim data. It expresses the same probability but can be computationally demanding as the number of remaining looks grows. The recursion instead uses one-dimensional integrals and the Markov transition. For a single remaining look, both approaches reduce to a normal tail probability.

  • Check the conditioning history. Use the observed statistic at the current look and condition on survival through that look. A rejection at an earlier look means the future recursion is irrelevant for that path.
  • Check the transition variance. It must be \(1-t_i/t_{i+1}\), not \(t_i/t_{i+1}\). The transition variance is positive when information strictly increases.
  • Check that every future look is included. A final-look tail alone is not the conditional error if there is an earlier future boundary.
  • Check the model and rule. The displayed recursion assumes the canonical transitions and fixed upper boundaries. A different stopping rule, statistic, or conditional distribution requires a corresponding calculation.

A common mistake is to treat the conditional error after an interim look as the original \(\alpha\), or to add separate marginal crossing probabilities as if the future looks were independent. The remaining error is conditional on the observed history, and future statistics are correlated through the canonical process. The recursion respects both facts by using conditional transition distributions and by partitioning paths into crossing now versus continuing.

The method has limits. It does not establish that a chosen canonical model is appropriate, nor does it fix invalid information-time estimates or an incorrectly specified boundary. If the transition law does not depend on the observed history only through the current statistic, the scalar recursion is insufficient: the conditioning state must contain all relevant information. Numerical integration should also be checked for stability, especially when boundaries are extreme or many looks remain.

Key takeaway: Under the canonical model, remaining Type I error is calculated by a backward recursion: count rejection at the next look, then add the conditional error from paths that survive that look. Each update is specific to the observed interim history.

Check Your Understanding

Try these questions, then compare your reasoning with the solutions.

  1. For a canonical design with \(t_i=0.50\) and \(t_{i+1}=0.80\), state the conditional variance of \(Z_{i+1}\mid Z_i=z\).
  2. Why does the recursion integrate the later conditional error only over values below the next boundary?
  3. A design has just one future look, with conditional mean 0.6, conditional standard deviation 0.8, and boundary 2. What is the conditional error in terms of \(\Phi\)?
  4. Explain why adding marginal crossing probabilities at two future looks generally does not give the probability of crossing at either look.

Solutions. (1) The variance is \(1-0.50/0.80=0.375\). (2) Values at or above the boundary are already counted as immediate rejection in the first integral; values below it continue and contribute the later conditional error. (3) It is \(1-\Phi((2-0.6)/0.8)=1-\Phi(1.75)\), approximately \(0.04006\). (4) The crossing events are dependent under the canonical joint distribution, so their marginal probabilities cannot generally be added without accounting for paths that cross at both looks.

Clinical Biostats