Tutorials › Adaptive Clinical Trial Design › Expected Length Versus Coverage of Adaptive Intervals

Estimation and Confidence Intervals After Adaptation · Tutorial 306 of 799

Expected Length Versus Coverage of Adaptive Intervals

Learn how to compare adaptive intervals by coverage and expected length, and why requiring target coverage on every adaptation path can cost precision.

Advanced 10 min read

What You'll Learn

  • Distinguish overall coverage from coverage conditional on an adaptation route.
  • Calculate expected interval length by weighting route-specific lengths by route probabilities.
  • Compare a 95% overall-coverage procedure with one targeting 95% coverage on every route.
  • Identify the model assumptions needed for route-specific normal intervals to be calibrated.
  • Diagnose when nominal expected-length comparisons may be misleading.
  • Explain why stronger pathwise guarantees can require longer intervals, but do not always do so.

Why Coverage and Length Must Be Compared Together

An adaptive trial can follow different routes, such as using a smaller or larger final sample depending on an earlier decision. A confidence interval may perform differently on each route: its coverage can vary, and its expected length can vary. Reporting only one overall coverage figure can conceal a route on which the interval misses its target often. Requiring the target on every route can protect against that weakness, but may increase the average interval length.

The comparison depends on what is guaranteed. Unconditional coverage averages over the random adaptation route and the final data. Route-specific coverage evaluates coverage among trials that followed a particular route. These are different properties. Neither coverage nor length alone summarizes the trade-off: a very wide interval can have high coverage, while an interval that is short on average can have inadequate coverage on an important route.

Definition: Let \(R\) denote the adaptation route, taking values \(r=1,\ldots,J\), and let \(C_r\) be the confidence interval used on route \(r\). Its unconditional coverage is \(P(\theta\in C_R)\). Its route-specific coverage on route \(r\) is \(P(\theta\in C_r\mid R=r)\). If the interval length is \(|C_r|\), its expected length is \(E(|C_R|)\), averaging over routes and the data.

The expected length is not the length seen in one trial. It is a repeated-trial average under the specified design and model. Likewise, unconditional coverage describes repeated trials before knowing which route will occur. Route-specific coverage describes repeated trials restricted to a route. This tutorial compares those criteria in a model where each route’s precision is known and the final estimator has a normal distribution conditional on that route.

A Model for Comparing the Criteria

Suppose the target parameter is \(\theta\), and the final estimate is \(\widehat\theta\). On route \(r\), assume \(\widehat\theta\mid R=r\sim N(\theta,s_r^2)\), where \(s_r\) is a known standard error for that route. Assume also that every route has positive probability, the route probabilities are specified, and the conditional normal model is correctly calibrated. One way these assumptions can hold is for an independent pilot to determine the route, followed by final data whose estimator has the stated normal distribution on each route.

For a central interval with conditional coverage \(c_r\), define \(z_{(1+c_r)/2}\) as the standard normal quantile with cumulative probability \((1+c_r)/2\). Then the route-specific interval and its length are

$$ C_r(c_r)=\left[\widehat\theta-z_{(1+c_r)/2}s_r,\ \widehat\theta+z_{(1+c_r)/2}s_r\right], \qquad |C_r(c_r)|=2z_{(1+c_r)/2}s_r. $$

Under the stated normal model, the standardized error \((\widehat\theta-\theta)/s_r\), conditional on route \(r\), has a standard normal distribution. The probability that it lies between the two symmetric endpoints is therefore \(c_r\). If \(\pi_r=P(R=r)\), the law of total probability gives

$$ P(\theta\in C_R)=\sum_{r=1}^J \pi_r c_r, \qquad E(|C_R|)=\sum_{r=1}^J \pi_r\,2z_{(1+c_r)/2}s_r. $$

The first expression averages coverage; the second averages length. These identities are a direct justification for comparing procedures route by route before summarizing them overall. They do not justify using the normal intervals if the conditional distribution is not normal or the standard error is not correctly specified.

Worked Example: Overall 95% Coverage With Unequal Route Guarantees

Consider two equally likely routes, \(A\) and \(B\), with known standard errors \(s_A=0.5\) units and \(s_B=2\) units. Suppose a procedure uses a central 99% interval on route \(A\), and a central 91% interval on route \(B\). The route and final estimate satisfy the model above.

The unconditional coverage is \(0.5(0.99)+0.5(0.91)=0.495+0.455=0.950\). Thus the procedure has 95% unconditional coverage, although its route-specific coverage is 99% on \(A\) and 91% on \(B\).

For route \(A\), the central 99% normal interval uses \(z_{0.995}\approx2.575829\), so its length is \(2(2.575829)(0.5)=2.575829\) units. For route \(B\), the central 91% interval uses \(z_{0.955}\approx1.695398\), so its length is \(2(1.695398)(2)=6.781592\) units. The expected length is therefore \(0.5(2.575829)+0.5(6.781592)=4.678711\) units.

This procedure meets the overall 95% target, but it does not give 95% coverage on every route. In particular, trials following route \(B\) have only 91% coverage under the model.

The Precision Cost of a Route-by-Route Requirement

Now require at least 95% coverage on each route. In this model, using the central 95% interval on both routes gives exact route-specific coverage \(0.95\): its critical value is \(z_{0.975}=\Phi^{-1}(0.975)\approx1.959964\). Its unconditional coverage is also 95%, because it averages two 95% values. The lengths are \(2(1.959964)(0.5)=1.959964\) units on route \(A\), and \(2(1.959964)(2)=7.839856\) units on route \(B\). Hence expected length is \(0.5(1.959964)+0.5(7.839856)=4.899910\) units.

The difference from the procedure in the worked example is \(4.899910-4.678711=0.221199\) units, or about \(4.51\%\) of the route-by-route procedure’s expected length. In this particular comparison, the procedure requiring 95% on both routes is longer on average. Its benefit is a clear guarantee for route \(B\); the cost is not a universal constant, but follows here from how coverage was allocated across routes.

The calculation also shows why stronger coverage requirements do not logically imply that every possible procedure must be longer. If a procedure already has at least 95% coverage on every route, it meets both criteria without any increase. If route-specific intervals have very different lengths, an unconditional procedure might allocate coverage in a way that is inefficient, and another procedure could satisfy a stronger requirement at similar or even lower expected length. The comparison must specify both procedures, the route probabilities, and the precision on each route.

Worked Example: Interpreting Intervals on the Observed Routes

Use the two-route setup above. Suppose, for illustration, that a hypothetical trial follows route \(A\) and produces \(\widehat\theta=1.20\) units. Under the 99% route-\(A\) interval, the margin of error is \(2.575829(0.5)=1.287915\) units. The interval is \([1.20-1.287915,\ 1.20+1.287915]=[-0.087915,\ 2.487915]\) units.

If instead a hypothetical trial follows route \(B\) and produces \(\widehat\theta=0.40\) units, its 91% interval has margin of error \(1.695398(2)=3.390796\) units. The interval is \([0.40-3.390796,\ 0.40+3.390796]=[-2.990796,\ 3.790796]\) units.

These realized intervals illustrate that the less precise route produces a wider interval, even though its assigned coverage is lower. Their endpoints do not determine whether either interval contains the unknown \(\theta\). The coverage statements concern the long-run behavior of the procedure under the assumed model, not a probability attached to the fixed parameter after observing one interval.

Coverage Allocation When Routes Have Different Probabilities

The cost also depends on how likely each route is. Consider two routes with probabilities \(0.8\) and \(0.2\), and the same standard errors \(0.5\) and \(2\) units, respectively. Assign 96% coverage to the more frequent, more precise route and 91% to the other. The unconditional coverage is \(0.8(0.96)+0.2(0.91)=0.768+0.182=0.950\).

Worked Example: Weighting Lengths by Route Probability

For the 96% interval, \(z_{0.98}\approx2.053749\), so route 1 has length \(2(2.053749)(0.5)=2.053749\) units. For the 91% interval, \(z_{0.955}\approx1.695398\), so route 2 has length \(2(1.695398)(2)=6.781592\) units. The expected length is \(0.8(2.053749)+0.2(6.781592)=1.642999+1.356318=2.999317\) units.

With 95% coverage on both routes, the expected length is \(0.8[2(1.959964)(0.5)]+0.2[2(1.959964)(2)]\) \(=0.8(1.959964)+0.2(7.839856)=1.567971+1.567971=3.135942\) units. Here the route-by-route target increases expected length by \(3.135942-2.999317=0.136625\) units. The overall target is the same, but route probabilities and route-specific precision determine the length comparison.

Assumptions, Diagnostics, and Limitations

The calculations rely on more than the stated overall coverage target. Before using them, verify that the parameter and estimator are defined consistently across routes; the route-specific standard errors are appropriate; and the conditional distribution of the estimator is adequately represented by the normal model. The route probabilities must also correspond to the adaptive design being evaluated. If the route is chosen using data that are also used in the final estimator, conditional normality and the stated standard error do not follow automatically. Selection can change the estimator’s conditional distribution, so a simple normal interval may not have its claimed route-specific coverage.

Useful diagnostics include simulating the full adaptation rule under parameter values relevant to the trial, then tabulating coverage and length separately by route as well as overall. Compare empirical route frequencies with the probabilities used in the expected-length calculation. Check whether nominal route-specific coverage is close to its target and whether interval length changes in the direction predicted by the route-specific standard errors. Large discrepancies can indicate a miscalibrated conditional distribution, incorrect route probabilities, or an implementation error. A satisfactory overall coverage estimate can hide a route-specific shortfall.

Expected length is only one precision summary. Two procedures with the same expected length can have different route-specific lengths, tail behavior, or probabilities of producing very wide intervals. Moreover, the numerical comparison depends on the chosen coverage allocations; the examples do not establish a universally best allocation. Alternatives include imposing a minimum coverage on every route, reporting both overall and route-specific operating characteristics, or using a common critical value across routes. Each choice answers a different design or reporting question.

Common pitfall: “The interval has 95% overall coverage, so it has 95% coverage whichever route the trial follows.” This does not follow. Overall coverage is the probability-weighted average of route-specific coverages; one route may be below 95% while another is above it. Calculate or simulate route-specific coverage when that guarantee matters.

Check Your Understanding

Work through the questions, then compare your reasoning with the solutions.

  1. Two routes are equally likely. Their route-specific coverages are 98% and 92%. What is unconditional coverage? Does the procedure achieve at least 95% on each route?
  2. Two equally likely routes have interval lengths 3 and 7 units. What is expected length?
  3. A procedure has route-specific coverages 95% and 95%, with route probabilities 0.7 and 0.3. Find its unconditional coverage.
  4. Why can a normal interval based on a route-specific standard error fail when the same interim data both select the route and contribute to the final estimate?

Solutions. 1. Unconditional coverage is \(0.5(0.98)+0.5(0.92)=0.95\), or 95%; it does not achieve 95% on each route because the second route has 92%. 2. Expected length is \(0.5(3)+0.5(7)=5\) units. 3. It is \(0.7(0.95)+0.3(0.95)=0.95\), or 95%. 4. Selection using the same data can alter the estimator’s distribution conditional on the selected route, so the assumed normal distribution and standard error may no longer give the claimed conditional coverage.

Clinical Biostats