Tutorials › Group Sequential Methods and Sample Size Re-Estimation › Likelihood-Based Intervals After Early Stopping

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

Likelihood-Based Intervals After Early Stopping

Learn how to construct a likelihood-ratio interval from a stopped Gaussian trial and why its nominal cutoff need not provide the target sequential coverage.

Advanced 12 min read

What You'll Learn

  • Write the canonical Gaussian likelihood for the observed increments up to the stopping look.
  • Derive the likelihood-ratio statistic and its usual interval under known variance.
  • Distinguish a likelihood interval from one obtained by inverting stopping-aware tests.
  • Calculate and compare both intervals for a hypothetical early efficacy stop.
  • Identify when conditioning on the stopping stage changes the likelihood being used.
  • Recognize why nominal fixed-sample cutoffs do not guarantee sequential coverage.

What a likelihood interval uses after a stop

An early-stopped trial supplies fewer observations than planned, and the observed estimate may be unusually far from the null because the stopping rule selected the time of analysis. A likelihood-based interval starts with the probability model for the data actually observed. Its construction is direct; its coverage, however, must be assessed under the sequential design rather than assumed from a fixed-sample formula.

We use a canonical Gaussian model with a known variance and deterministic information at each planned look. Write \(I_k\) for the cumulative information at look \(k\), and \(S_k\) for the cumulative score. The independent increments satisfy \(\Delta S_j=S_j-S_{j-1}\sim N(\theta\Delta I_j,\Delta I_j)\), where \(\Delta I_j=I_j-I_{j-1}\), \(I_0=S_0=0\), and \(\theta\) is the treatment effect in the outcome's units. Thus \(S_k\sim N(\theta I_k,I_k)\). The stopping rule is fixed in advance and depends on the observed data, not on unknown \(\theta\).

Definition: For a realized path observed through stopping look \(K\), the Gaussian likelihood, apart from factors not involving \(\theta\), is \(L(\theta)\propto\exp\{\theta S_K-\tfrac12\theta^2 I_K\}\). Its maximum-likelihood estimate is \(\hat{\theta}=S_K/I_K\). A likelihood-ratio confidence set with cutoff \(c\) contains the values of \(\theta\) satisfying \(2\{\ell(\hat{\theta})-\ell(\theta)\}\le c\), where \(\ell\) is the log likelihood.

This likelihood uses all observed increments through the stop. The event that the realized path met a stopping boundary is not an extra parameter-dependent factor in the likelihood of that full path: the rule determines whether that path is observed, but does not alter its sampling density. This statement does not mean that the stopping rule is irrelevant to coverage. It shapes the distribution of the estimate and information at the time of stopping.

Deriving the likelihood-ratio interval

Completing the square in the log likelihood gives a particularly simple result:

$$ \ell(\theta)=\theta S_K-\frac{1}{2}\theta^2 I_K+C =\ell(\hat{\theta})-\frac{1}{2}I_K(\theta-\hat{\theta})^2, \qquad 2\{\ell(\hat{\theta})-\ell(\theta)\}=I_K(\theta-\hat{\theta})^2. $$

For a nominal 95% likelihood-ratio interval, the usual one-degree-of-freedom cutoff is \(c=\chi^2_{1,0.95}\approx3.841\). Since \(\sqrt{3.841}\approx1.96\), solving the likelihood-ratio inequality gives:

$$ I_K(\theta-\hat{\theta})^2\le3.841 \quad\Longleftrightarrow\quad \hat{\theta}-\frac{1.96}{\sqrt{I_K}} \le\theta\le \hat{\theta}+\frac{1.96}{\sqrt{I_K}}. $$

This is exactly the familiar fixed-information Gaussian interval evaluated at the observed information. The algebra is exact for the stated Gaussian model. What is not automatic is 95% coverage across repetitions of a sequential trial: \(K\) and the estimate at \(K\) are jointly affected by the stopping rule. In particular, a cutoff justified by a fixed-sample reference distribution is not, by itself, a guarantee of sequential coverage.

Worked Example: Likelihood Interval at an Early Stop

Suppose a hypothetical trial has known-variance Gaussian information \(I_1=64\) at its first look and stops there. Its observed score is \(S_1=11.52\), so the estimated treatment effect is \(11.52/64=0.1800\) outcome units. Assume the Gaussian model and known variance specified above.

The log-likelihood difference from its maximum is \(2\{\ell(\hat{\theta})-\ell(\theta)\}=64(\theta-0.1800)^2\). Using the nominal 95% cutoff, the interval is:

$$ 0.1800\pm\frac{1.96}{\sqrt{64}} =0.1800\pm\frac{1.96}{8} =0.1800\pm0.2450 =[-0.0650,\ 0.4250]. $$

On the stated outcome scale, the likelihood-ratio interval ranges from \(-0.0650\) to \(0.4250\) outcome units. It describes the values retained by this likelihood-ratio rule. Calling it a 95% sequential confidence interval requires a separate coverage justification under the stopping design.

Why inversion of stopping-aware tests can differ

An interval obtained by inverting tests is built from a different question: for each candidate value \(\theta_0\), would the prespecified test reject \(H_0:\theta=\theta_0\) given the observed stopped path? A stopping-aware test calibrates its rejection rule to the sequential sampling plan. It accounts for the opportunities to cross boundaries at multiple looks, rather than applying a fixed-sample reference distribution at whichever look happened to stop the trial.

For example, suppose a two-sided test rejects a candidate \(\theta_0\) if a standardized score crosses a specified boundary at any look. At look \(j\), the statistic is \(Z_j(\theta_0)=(S_j-I_j\theta_0)/\sqrt{I_j}\). If the design's early boundary is \(b_1\), then at an early stop the candidate is retained when \(|Z_1(\theta_0)|<b_1\). Inverting this condition produces an interval centered at the observed estimate with half-width \(b_1/\sqrt{I_1}\), provided the test uses symmetric boundaries and no other observed-look condition further restricts the set.

The boundary is a feature of the test design, not a universal likelihood-ratio constant. A sequence of boundaries must be calibrated to the design's desired error rate; the resulting confidence set inherits that test calibration. By contrast, the usual likelihood interval uses the fixed cutoff 3.841. Consequently, likelihood-based and test-inversion intervals can have different widths even when they use the same observed estimate and information.

Worked Example: Inverting an Early-Look Test

Consider a hypothetical Gaussian trial with known variance, early information \(I_1=36\), and observed score \(S_1=16.2\). The trial stops for efficacy at this look under a prespecified boundary. The estimate is \(\hat{\theta}=16.2/36=0.4500\) outcome units, and its standardized score under the null value zero is \(16.2/6=2.70\). Suppose the stopping-aware two-sided test has early boundary \(b_1=2.50\); its later boundary is calibrated as part of the design, but is not needed for this early-stop inversion.

A candidate \(\theta_0\) is retained by the early-look test when \(\left|(16.2-36\theta_0)/6\right|<2.50\). Rearranging gives \(\left|0.4500-\theta_0\right|<2.50/6\), so the inverted interval is:

$$ 0.4500\pm\frac{2.50}{6} =0.4500\pm0.4167 =[0.0333,\ 0.8667]\ \text{outcome units}. $$

The observed null statistic, 2.70, exceeds the specified boundary 2.50, so zero is excluded by this test inversion. The interval is wider than the nominal likelihood interval, whose half-width is \(1.96/6=0.3267\) outcome units. This comparison is valid for the stipulated boundary test; it does not establish that 2.50 alone guarantees a particular overall error rate. That guarantee depends on calibration of the complete sequential test.

What stopping does—and does not—change in the likelihood

The likelihood calculation above uses the full observed path. If only a summary is available—for example, the stopping stage and endpoint statistic, but not the earlier increments—the likelihood must reflect the probability law of that recorded summary. Likewise, a likelihood conditional on stopping at look \(k\) is a different object from the likelihood of the full path. In general, conditioning on \(K=k\) introduces a parameter-dependent normalizing probability, because \(P_\theta(K=k)\) can vary with \(\theta\). One should not silently replace the full-path likelihood with a conditional likelihood, or the reverse.

The distinction matters because a likelihood formula and a confidence procedure are not the same thing. The likelihood describes relative support for parameter values under the chosen data model. A confidence interval is a procedure with a coverage property under repeated sampling. A nominal likelihood-ratio cutoff may be convenient and informative, but its sequential coverage must be checked for the stopping rule and the data actually recorded. Inverting a stopping-aware test targets coverage through its test calibration; it may produce less symmetric limits or limits with different widths at different stopping stages.

Practical checks and limitations

Before reporting a likelihood-based interval, verify that the assumptions behind the calculation fit the design: the Gaussian increment model, known variance, deterministic information increments, and correctly recorded path through the stop. If the variance is estimated, the displayed likelihood and cutoff no longer give the stated interval; uncertainty in the variance needs its own treatment. If information increments are random or depend on the observed outcomes, the deterministic-information argument also needs revision.

Useful diagnostics include checking that the reported \(S_K\), \(I_K\), and \(\hat{\theta}=S_K/I_K\) agree; confirming which looks and stopping boundaries were actually used; and comparing the interval's claimed coverage with the sequential design's error calibration. Large differences between the nominal likelihood interval and a stopping-aware interval are not an arithmetic error by themselves. They can signal that the selection induced by stopping is consequential.

The principal alternatives answer related but distinct questions. A standard likelihood-ratio interval summarizes likelihood support with a chosen cutoff. An interval from inverting stopping-aware tests is tied directly to a sequential rejection rule and its calibration. Conditional intervals address inference given a specified stopping stage, while unconditional intervals average coverage across stopping stages; neither is interchangeable with an unadjusted fixed-sample interval. The appropriate choice depends on the inferential target, the recorded data, and the error-control properties required.

Common pitfall: “The likelihood is ordinary Gaussian, so the usual 95% likelihood interval must have 95% coverage after stopping.” The Gaussian likelihood calculation can be correct while its fixed-sample calibration is not valid for the sequential procedure. State which likelihood and cutoff were used, and distinguish likelihood support from sequential coverage.

Check Your Understanding

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

  1. In the known-variance canonical model, what is the maximum-likelihood estimate at stopping look \(K\), and which two quantities determine the likelihood-ratio interval's center and width?
  2. A hypothetical stopped trial has \(I_K=100\) and \(\hat{\theta}=0.12\) outcome units. Calculate its nominal 95% likelihood-ratio interval.
  3. Why does the ordinary Gaussian likelihood not, by itself, establish nominal sequential coverage?
  4. At an early look with \(I_1=25\), an inverted symmetric stopping-aware test uses boundary \(b_1=2.40\). If the estimate is \(0.30\) outcome units, calculate the early-look inverted interval.
  5. Explain why conditioning on the stopping stage can change the likelihood compared with using the full observed path.

Solutions. (1) \(\hat{\theta}=S_K/I_K\); the center is this estimate and the usual interval's half-width is \(1.96/\sqrt{I_K}\). (2) \(0.12\pm1.96/10=0.12\pm0.196=[-0.076,\ 0.316]\) outcome units. (3) Stopping changes the joint sampling distribution of the observed estimate and stopping look, whereas the usual cutoff comes from a fixed-information reference calculation; coverage must be justified under the sequential design. (4) The half-width is \(2.40/5=0.48\), giving \(0.30\pm0.48=[-0.18,\ 0.78]\) outcome units. (5) The probability of stopping at that stage can depend on \(\theta\), so a conditional likelihood generally has a parameter-dependent normalization absent from the full-path likelihood.

Clinical Biostats