Tutorials › Group Sequential Methods and Sample Size Re-Estimation › Classical Efficacy Boundaries And Effect Size

Classical Boundaries · Tutorial 117 of 699

Classical Efficacy Boundaries And Effect Size

Learn how null-calibrated efficacy boundaries combine with effect size and information to determine the probability of stopping for efficacy.

Intermediate 10 min read

What You'll Learn

  • Distinguish null calibration of efficacy boundaries from power under an alternative.
  • Express crossing and first-crossing probabilities for a sequence of looks.
  • Relate treatment effect and information to the mean of a group sequential statistic.
  • Interpret how effect size changes the chance of crossing a fixed boundary.
  • Identify why marginal look-by-look probabilities are not overall group sequential power.
  • Diagnose assumptions that can undermine a calibrated design.

Why the Same Boundary Has Different Roles

A classical efficacy boundary is chosen to limit the chance of a false efficacy conclusion when the null hypothesis is true. Its calibration therefore uses the distribution of the test statistics under \(H_0\). Once the boundary is fixed, however, the chance of crossing it in a trial with a real treatment effect is a different quantity: it depends on how far the alternative distribution is shifted toward the boundary. That probability is part of the design’s power.

This distinction matters because a critical value is not a promise that a trial will have a particular probability of success. It is a decision threshold. The null distribution determines how often a sequence of statistics crosses the threshold when there is no effect; the alternative distribution determines how often it crosses when an effect of a specified size is present.

Definition: In a one-sided group sequential efficacy design, let \(Z_k\) be the test statistic at planned look \(k\), for \(k=1,\ldots,K\), and let \(c_k\) be the upper efficacy boundary at that look. The design rejects \(H_0\) for efficacy if \(Z_k\ge c_k\) at any look. Its overall type I error is \(P_0(\text{cross at least one boundary})\); its power at a specified alternative is \(P_\Delta(\text{cross at least one boundary})\).

Model, Assumptions, and the Null Calibration

A standard setting has a prespecified sequence of looks, with information \(I_1<\cdots<I_K\). For a treatment effect \(\Delta\), the statistic at look \(k\) is approximately normal with unit variance and mean \(\Delta\sqrt{I_k}\). Thus, under the null effect \(\Delta=0\), \(Z_k\) has mean zero; under a positive alternative, its mean is shifted upward. Information is measured on the scale that makes \(\Delta\sqrt{I_k}\) dimensionless.

For cumulative statistics based on nested data, the standard canonical model also has \(\operatorname{Corr}(Z_j,Z_k)=\sqrt{I_j/I_k}\) for \(j<k\). The statistics are therefore generally dependent: observations contributing to an earlier statistic also contribute to later ones. In applications, a normal or asymptotic normal model is used when appropriate, and the information sequence and dependence structure used for calibration must represent the planned design.

For upper boundaries, define \(A_k\) as the event that the trial first crosses at look \(k\): it has not crossed at any earlier look, and \(Z_k\ge c_k\). The events \(A_1,\ldots,A_K\) are disjoint. Therefore the null error rate is the sum of their null probabilities:

$$ \alpha_{\text{design}} = P_0\left(\bigcup_{k=1}^{K}\{Z_k\ge c_k\}\right) = \sum_{k=1}^{K} P_0(A_k). $$

Classical boundaries are calibrated by choosing the \(c_k\)'s so that this joint probability equals the target significance level \(\alpha\), or is no greater than it. Because the statistics are dependent, adding the marginal null tail probabilities \(P_0(Z_k\ge c_k)\) generally double-counts paths that cross at more than one look. Calibration instead accounts for the joint path and, in particular, for the first look at which a crossing occurs.

Key distinction: The boundary values are selected using the null distribution to control type I error. Power is calculated after those values are fixed, using the alternative distribution. Changing the assumed effect size changes power, not the null-calibrated boundary.

How Effect Size Changes Crossing Probability

Under the canonical normal model, an alternative effect shifts the mean of each statistic by \(\Delta\sqrt{I_k}\). For a fixed look, the marginal probability of crossing its upper boundary is \(P_\Delta(Z_k\ge c_k)=1-\Phi(c_k-\Delta\sqrt{I_k})\), where \(\Phi\) is the standard normal cumulative distribution function. Increasing a positive \(\Delta\), or increasing information while holding \(\Delta\) fixed, increases this marginal probability.

Overall power is not just the marginal probability at the final look. It is the probability of crossing at any look, with earlier crossing paths counted once. In full form,

$$ \operatorname{Power}(\Delta) = P_\Delta\left(\bigcup_{k=1}^{K}\{Z_k\ge c_k\}\right) = \sum_{k=1}^{K}P_\Delta(A_k). $$

The joint alternative distribution is needed to calculate this probability. For example, an early crossing and a later crossing are not independent events. The first-crossing probabilities must use the joint distribution of the sequence, just as they do under the null.

Worked Calculations

Worked Example: A Single-Look Boundary Is Calibrated Under the Null

Suppose a one-sided, single-look design tests \(H_0:\Delta=0\) against a positive treatment effect, using a normally distributed test statistic with variance 1 under \(H_0\). The target significance level is \(\alpha=0.025\). Find the boundary and state what it guarantees.

With one look, the overall type I error is simply the probability of exceeding the boundary under the null. The standard normal upper-tail cutoff is

$$ c_1=\Phi^{-1}(1-\alpha) =\Phi^{-1}(0.975) \approx 1.960. $$

Under \(H_0\), \(Z_1\sim N(0,1)\), so \(P_0(Z_1\ge 1.960)\approx 0.0250\). The boundary therefore controls the chance of a false efficacy declaration at 2.5% for this one-sided test. This calculation says nothing by itself about the probability of crossing when \(\Delta>0\); that requires specifying the effect and information.

Worked Example: Effect Size Changes an Interim Crossing Probability

Consider a hypothetical two-look design with a first-look boundary \(c_1=2.50\). Assume this boundary is part of a complete set of boundaries calibrated under the null, and that the first-look information is \(I_1=0.0625\) per \((\text{mmHg})^2\). At the first look, compare two positive treatment effects: \(\Delta=2\) mmHg and \(\Delta=4\) mmHg. Use the canonical normal model.

The first-look statistic has variance 1 and mean \(\Delta\sqrt{I_1}\). Since \(\sqrt{I_1}=0.25\) per mmHg, the mean is \(0.50\) for the 2 mmHg effect and \(1.00\) for the 4 mmHg effect. For the smaller effect,

$$ P_{\Delta=2}(Z_1\ge 2.50) =1-\Phi(2.50-2(0.25)) =1-\Phi(2.00) \approx 0.0228. $$

For the larger effect,

$$ P_{\Delta=4}(Z_1\ge 2.50) =1-\Phi(2.50-4(0.25)) =1-\Phi(1.50) \approx 0.0668. $$

The interim boundary is identical in both calculations because it belongs to the fixed design. The probability of crossing it is higher for the larger effect. These are first-look crossing probabilities, not the overall power of the two-look design: overall power also includes paths that do not cross at look 1 but cross at look 2.

Worked Example: Separating Effect Size From Information

Suppose a hypothetical one-look design has an upper boundary \(c_1=1.960\), calibrated for a one-sided null error rate of 0.025. A treatment effect is measured in mmHg. Compare the power for \(\Delta=0.5\) mmHg at information \(I_1=4\) per \((\text{mmHg})^2\) with the power for the same effect at \(I_1=16\) per \((\text{mmHg})^2\).

At \(I_1=4\), the alternative mean of the statistic is \(\Delta\sqrt{I_1}=0.5\sqrt{4}=1.00\). The one-look power is

$$ P_{\Delta=0.5}(Z_1\ge 1.960) =1-\Phi(1.960-1.00) =1-\Phi(0.960) \approx 0.1685. $$

At \(I_1=16\), the alternative mean is \(0.5\sqrt{16}=2.00\). The boundary remains 1.960, since the null-calibrated design has not changed. Power is now

$$ P_{\Delta=0.5}(Z_1\ge 1.960) =1-\Phi(1.960-2.00) =1-\Phi(-0.040) \approx 0.5160. $$

Thus, for the same effect in mmHg, greater information increases the noncentrality and the chance of crossing. This is a one-look calculation, so the probability shown is also the overall power. With multiple looks, the same mean-shift principle applies, but overall power must account for their joint distribution and the possibility of crossing earlier.

Why Power Increases With a Positive Effect

For upper efficacy boundaries, the dependence of crossing probability on a positive effect can be justified directly. Couple the statistics under different effects by using the same null-centered random path \(Z_k(0)\), and write \(Z_k(\Delta)=Z_k(0)+\Delta\sqrt{I_k}\). If \(\Delta_2\ge\Delta_1\ge0\), then for every look \(k\),

$$ Z_k(\Delta_2) =Z_k(\Delta_1)+(\Delta_2-\Delta_1)\sqrt{I_k} \ge Z_k(\Delta_1). $$

Consequently, any path that crosses an upper boundary under \(\Delta_1\) also crosses an upper boundary under \(\Delta_2\), using this coupling. The crossing event for the smaller effect is contained in the crossing event for the larger effect, so its probability cannot be greater. This monotonicity concerns positive effects and upper boundaries; it does not mean that every individual trial becomes more likely to stop at a particular look.

Interpreting the Design and Diagnosing Problems

A null-calibrated boundary and an alternative power calculation answer different design questions. Calibration asks how often the planned rule rejects \(H_0\) when the null model holds. Power asks how often that same rule rejects \(H_0\) when a specified effect and information sequence hold. It is therefore useful to report the assumed effect size and information alongside a power estimate; “the design has power” is incomplete without them.

A common alternative to group sequential monitoring is a fixed-sample design with only a final analysis. In a one-look design, a marginal normal tail probability is the full type I error or power. In a group sequential design, repeated looks create multiple opportunities to cross, so the nominal cutoff for a single test cannot simply be used at every look while claiming the same overall error rate. Calibrated classical boundaries account for those repeated opportunities under the null. Under an alternative, one must then calculate power for the entire sequence rather than treating the looks as independent tests.

Several checks help identify a mistaken calculation or an unsuitable model:

  • Check which distribution is being used. The boundary calibration uses the null distribution. Power uses an alternative with a stated effect size.
  • Check the statistic’s mean shift. Confirm that \(\Delta\sqrt{I_k}\) is dimensionless and uses the effect and information on compatible scales.
  • Check the whole path. Summing marginal crossing probabilities is generally wrong because crossings at different looks can occur on the same path.
  • Check the planned design. The information times, look schedule, and dependence structure used for calibration should match the analysis plan. Material departures may invalidate the stated error rate.
  • Check the direction and rule. The monotonicity argument above is for positive effects and upper boundaries. A two-sided design uses a different crossing event.

The normal model and canonical dependence are approximations in many applications. If the test statistic is poorly approximated by its assumed distribution, the information estimates are inaccurate, or the actual analysis departs from the planned rule, the calculated type I error and power may not describe the trial well. Calibration does not remove these modeling requirements.

Common pitfall: It is incorrect to say, “The boundary is 2.50, so the probability of crossing it is fixed.” The boundary is fixed, but the crossing probability is not: under a normal model it depends on the statistic’s mean, which changes with effect size and information. The correction is to name the probability’s reference distribution and calculate it under either \(H_0\) or a specified alternative.

Check Your Understanding

Answer each question, then compare your reasoning with the brief solutions below.

  1. A design has three planned looks and calibrated upper boundaries. Which distribution is used to control overall type I error, and which distribution is used to calculate power?
  2. At one look, \(c_k=2.00\), \(\Delta=1.5\) units, and \(I_k=0.64\) per unit squared. Under the canonical normal model, calculate the marginal crossing probability at that look.
  3. Why is the sum of the marginal probabilities \(P_\Delta(Z_k\ge c_k)\) generally not the overall power of a multi-look design?
  4. For positive effects and upper boundaries, what happens to the crossing event when \(\Delta\) increases while the information and boundaries stay fixed?

Solutions. (1) Overall type I error is calculated under the joint null distribution; power is calculated under the joint distribution for a specified alternative. (2) The mean is \(1.5\sqrt{0.64}=1.20\). Thus \(P(Z_k\ge2.00)=1-\Phi(2.00-1.20)=1-\Phi(0.80)\approx0.2119\). This is a marginal probability at one look, not necessarily overall power. (3) The statistics are dependent, and a path may cross at multiple looks; summing marginal probabilities can count that path more than once. Overall power is the probability of the union of crossing events, or equivalently the sum of the disjoint first-crossing probabilities. (4) Under the same null-centered path, every statistic shifts upward by \((\Delta_2-\Delta_1)\sqrt{I_k}\) when the effect increases from \(\Delta_1\) to \(\Delta_2\). The crossing event for the smaller positive effect is therefore contained in the crossing event for the larger effect, so power cannot decrease.

Clinical Biostats