Tutorials › Group Sequential Methods and Sample Size Re-Estimation › Stopping-Induced Bias in Treatment Effect Estimates

Bias-Adjusted Estimation and Confidence Intervals · Tutorial 286 of 699

Stopping-Induced Bias in Treatment Effect Estimates

Learn why estimates reported after data-dependent stopping can be systematically extreme even when each prespecified-look estimator is unbiased.

Advanced 10 min read

What You'll Learn

  • Distinguish fixed-sample unbiasedness from bias induced by selecting a look using interim results.
  • Express the stopped estimate’s bias as a sum over the possible stopping looks.
  • Work through a two-look normal example showing upward bias at the null.
  • Interpret an efficacy-stopped estimate without treating one observed trial as proof of bias.
  • Recognize when the direction of stopping-induced bias depends on the stopping rule.
  • Identify why fixed-sample inference may not describe estimates selected by a sequential rule.

Why the Reported Estimate Can Be Extreme

A treatment-effect estimator can be unbiased when the sample size is fixed in advance and still be biased when the study reports the estimate from a look chosen in response to interim results. The distinction is the selection: an efficacy rule tends to stop when the observed estimate is unusually favorable, so estimates from early stops are selected values, not a random sample of all possible estimates at that look.

This does not mean every estimate from a sequential trial is necessarily biased in the same direction, or that a particular early estimate is known to be too large. The direction and size of the bias depend on the stopping rule, the true effect, and which outcomes are included in the expectation. The central point is that the rule for choosing when to report can change the repeated-sampling distribution of the reported estimator.

Definition: Let \(\hat{\theta}_k\) be the treatment-effect estimator at look \(k\), and let \(\tau\) be the random look at which the trial stops. The reported estimator is \(\hat{\theta}_\tau\). Its stopping-induced bias at parameter value \(\theta\) is \(E_\theta(\hat{\theta}_\tau)-\theta\). This is distinct from the fixed-look bias \(E_\theta(\hat{\theta}_k)-\theta\), which refers to an estimator evaluated at a prespecified look \(k\).

Fixed-Look Unbiasedness and Selection

For a simple reference model, suppose that at each prespecified look \(k\), the estimator satisfies \(\hat{\theta}_k\sim N(\theta,1/I_k)\), where \(I_k\) is the information available at that look. Assume the treatment-effect model is correctly specified, the variance or information scale is known, and the estimator is unbiased at each fixed look. The normal model is a useful approximation in many large-sample settings; its assumptions should not be treated as automatic for every trial.

A fixed-look estimator has expectation \(\theta\) before the data are observed. But a stopping rule makes the event \(\{\tau=k\}\) depend on interim data. Consequently, even if \(E_\theta(\hat{\theta}_k-\theta)=0\), the selected error \((\hat{\theta}_k-\theta)\) need not have mean zero among trials that stop at \(k\). In general,

$$ E_\theta(\hat{\theta}_\tau-\theta) = \sum_k E_\theta\left[(\hat{\theta}_k-\theta)\mathbf{1}\{\tau=k\}\right]. $$

Here \(\mathbf{1}\{\tau=k\}\) equals 1 when stopping occurs at look \(k\), and 0 otherwise. Fixed-look unbiasedness says that the expectation of the error without this indicator is zero. It does not say that the error multiplied by an outcome-dependent selection indicator has expectation zero. The indicator can give more weight to favorable estimation errors.

Worked Example: A Fixed-Sample Estimate Is Unbiased Before Selection

Suppose a hypothetical fixed-sample trial estimates a mean treatment effect \(\theta\), measured in milligrams per deciliter, with \(\hat{\theta}=\bar{X}_T-\bar{X}_C\). Assume independent randomized groups, finite outcome variances, and a prespecified sample size in each group. Then \(E_\theta(\bar{X}_T)=\mu_T\) and \(E_\theta(\bar{X}_C)=\mu_C\), where \(\theta=\mu_T-\mu_C\).

Linearity of expectation gives \[ E_\theta(\hat{\theta})=E_\theta(\bar{X}_T)-E_\theta(\bar{X}_C) =\mu_T-\mu_C=\theta. \] Thus the estimator is unbiased under these assumptions when the sample size and analysis are fixed in advance.

For a hypothetical observed difference of \(1.2\) milligrams per deciliter with standard error \(0.4\) milligrams per deciliter, the standardized statistic is \(1.2/0.4=3.00\). That observed value does not, by itself, show that the estimator is biased. Bias is a property of the estimator’s expectation across repeated trials under a specified procedure, not a label assigned to one unusually large result.

Two-Look Derivation: How the Stopping Rule Changes the Mean

The following canonical normal example isolates the mechanism. Let the true standardized treatment effect be \(\theta=0\). At the first look, let the information be \(I_1=1\), so the score \(Z_1=\sqrt{I_1}\hat{\theta}_1\) has distribution \(N(0,1)\). The second look has total information \(I_2=2\). Write its cumulative score as \(Z_2=Z_1+\varepsilon\), where \(\varepsilon\sim N(0,1)\) is independent of \(Z_1\). Then \(\hat{\theta}_1=Z_1\) and \(\hat{\theta}_2=Z_2/2\); both estimates are unbiased at their respective fixed looks.

Suppose the rule stops early if \(Z_1\geq c\), with \(c=1.96\), and otherwise continues to the second look. This is a hypothetical one-sided efficacy rule. It assumes a prespecified boundary, the stated independent normal increments, and that no other stopping rule changes which estimate is reported.

1
Write the estimate selected by the rule.
If \(Z_1\geq c\), report \(Z_1\). If \(Z_1<c\), report \((Z_1+\varepsilon)/2\). Therefore the stopped estimate is \(Z_1\mathbf{1}\{Z_1\geq c\} + (Z_1+\varepsilon)\mathbf{1}\{Z_1<c\}/2\).
2
Take its expectation.
Independence and \(E(\varepsilon)=0\) imply that \(E[\varepsilon\mathbf{1}\{Z_1<c\}]=0\). Also \(E(Z_1)=0\), so \(E[Z_1\mathbf{1}\{Z_1<c\}]=-E[Z_1\mathbf{1}\{Z_1\geq c\}]\).
3
Simplify and evaluate.
The first-look contribution has weight 1 on early stops and weight \(1/2\) on continuations. The resulting expectation is one-half of the upper-tail first-look score contribution, which is positive.
$$ \begin{aligned} E_0(\hat{\theta}_\tau) &= E[Z_1\mathbf{1}\{Z_1\geq c\}] + \frac{1}{2}E[(Z_1+\varepsilon)\mathbf{1}\{Z_1<c\}]\\ &= \frac{1}{2}E[Z_1\mathbf{1}\{Z_1\geq c\}]\\ &= \frac{1}{2}\phi(c) = \frac{1}{2}\phi(1.96) \approx \frac{1}{2}(0.05844) = 0.02922. \end{aligned} $$

In the last line, \(\phi\) is the standard normal density. The identity \(E[Z_1\mathbf{1}\{Z_1\geq c\}]=\phi(c)\) follows by integrating \(z\phi(z)\) from \(c\) to infinity and using \(\phi'(z)=-z\phi(z)\). Thus, although both fixed-look estimators have expectation zero, the reported estimator has positive expectation, about \(0.02922\) standardized treatment-effect units under this particular rule and model.

Worked Example: The Early-Stopped Subgroup

In the same hypothetical model, consider only repetitions that stop at the first look. The early-stop event is \(Z_1\geq1.96\). Its probability under \(\theta=0\) is approximately \(0.02500\). The average estimate among those selected repetitions is \[ E(Z_1\mid Z_1\geq1.96) =\frac{E[Z_1\mathbf{1}\{Z_1\geq1.96\}]}{P(Z_1\geq1.96)} =\frac{0.05844}{0.02500} \approx 2.338. \]

So early-stopped estimates average about \(2.338\) first-look standard-error units in this null model, despite a true effect of zero. This is a conditional expectation among selected early stops, not the unconditional bias of every possible reported estimate. The unconditional bias in the preceding calculation is much smaller because most repetitions continue and contribute a later estimate.

Reading a Stopped Estimate in Context

Consider a separate hypothetical trial with first-look information \(I_1=25\) and an efficacy boundary \(Z_1\geq2.40\). Suppose the observed first-look estimate is \(0.48\) treatment-effect units and its standard error is \(0.20\) units. The statistic is \[ Z_1=\frac{0.48}{0.20}=2.40. \] The boundary is reached, so the protocol reports the estimate \(0.48\) units and stops. The calculation verifies the rule’s application; it does not calculate how much this particular estimate is biased. That would require a repeated-sampling model and an expectation over the stopping procedure.

The appropriate interpretation is that \(0.48\) is the observed estimate at the selected stopping look. It should not be described as though it came from a sample size fixed independently of the observed treatment effect. Nor does the selection argument prove that the true treatment effect is smaller than \(0.48\). It warns that, across repetitions of the same design, early efficacy stops preferentially occur when interim estimates are favorable.

Worked Example: What One Early Stop Does and Does Not Show

Suppose the protocol in the hypothetical trial above specifies the first look at \(I_1=25\), an efficacy boundary of \(2.40\), and a later planned look. At the first look, the estimated treatment effect is \(0.48\) units with standard error \(0.20\) units. The statistic \(0.48/0.20=2.40\) meets the boundary exactly, so the trial stops and reports \(0.48\) units.

The reported estimate is the result of applying the prespecified stopping rule to these observed data. The selection mechanism implies that estimates from early efficacy stops tend to be high relative to their fixed-look sampling distribution. However, this single observation cannot establish the size of the repeated-sampling bias or determine the unobserved estimate that would have resulted at the later look. Those are different questions from calculating the observed estimate.

Assumptions, Diagnostics, and Alternatives

The normal derivation depends on a correctly stated information scale, unbiased fixed-look estimators, independent increments, and a stopping rule based on the interim statistic. Actual analyses may depart from these assumptions because of model misspecification, estimated nuisance parameters, delayed outcomes, changes to the look schedule, or additional futility and safety rules. The direction of bias also need not be upward in every design: an upper efficacy boundary favors positive estimates among early stops, while a lower boundary or a different reporting rule can select in another direction.

Useful checks focus on whether the analysis accounts for how the estimate was selected. Confirm that the actual stopping rule, look timing, and boundary used in the trial match the rule assumed for inference. In repeated simulations under plausible parameter values, compare the mean of the fixed-look estimator with the mean of the reported stopped estimator, and examine early-stop estimates separately from all reported estimates. A large difference is evidence that selection matters under the simulated model; it is not a universal bias correction for the real trial.

The main alternative is not to pretend that the sample size was fixed. One may report the ordinary estimate with clear disclosure that it was selected by a stopping rule, and use inferential procedures calibrated to the sequential design when estimation uncertainty is needed. A fixed-sample analysis can remain a useful descriptive reference, but its usual sampling properties do not automatically carry over after data-dependent selection. Procedures that adjust the estimate or construct intervals accounting for stopping require additional methods; the key prerequisite is to specify the stopping rule and the target of inference.

Common mistake: “The estimator is unbiased at every fixed look, so the estimate reported at the stopping look is unbiased too.” The conclusion does not follow: the stopping look is selected using the data. The correction is to distinguish \(E_\theta(\hat{\theta}_k)=\theta\) at a prespecified \(k\) from \(E_\theta(\hat{\theta}_\tau)\), where \(\tau\) depends on interim results.

Stopping-induced bias matters because an estimate often communicates the magnitude of benefit, not merely whether a boundary was crossed. If the estimate is selected to be favorable, presenting it without the stopping context can overstate how representative it is of the design’s repeated-sampling behavior. The bias is a property of the complete procedure—data collection, interim decisions, and reporting—not simply a defect in the arithmetic of the estimator.

Check Your Understanding

Answer each question, then compare with the solutions below.

  1. Why does fixed-look unbiasedness not imply unbiasedness of the reported estimate \(\hat{\theta}_\tau\)?
  2. In the two-look example, what is the unconditional bias at \(\theta=0\), and what does its sign mean?
  3. Why is the early-stop conditional mean of about \(2.338\) not the unconditional bias of the stopped estimator?
  4. A hypothetical estimate is \(0.30\) units with standard error \(0.15\) units at an upper boundary of \(2.00\). Does it cross the boundary, and what can be concluded about its bias from this one result?
  5. Name two assumptions used in the two-look derivation and one way the actual trial could violate the simplified model.

Solutions. (1) The stopping event depends on interim data, so the selection indicator can be associated with the estimation error; a zero mean without selection need not remain zero after selection. (2) The bias is \(\phi(1.96)/2\approx0.02922\) standardized units, which is upward under the stated null model and rule. (3) The value \(2.338\) is conditional on stopping early; the unconditional mean also includes continuations, whose estimates offset much of the early-stop contribution. (4) The statistic is \(0.30/0.15=2.00\), so it reaches the boundary. This one result gives the observed estimate and stopping decision, not the repeated-sampling bias. (5) Assumptions include unbiased fixed-look estimates and independent normal increments; for example, an unplanned change in look timing or a different safety stopping rule could violate the specified design.

Clinical Biostats