Tutorials › Group Sequential Methods and Sample Size Re-Estimation › Sensitivity to Interim Timing

Unblinded Re-Estimation and the Promising Zone · Tutorial 432 of 699

Sensitivity to Interim Timing

Compare alternative interim timings by calculating how information fraction changes promising-zone entry and expected sample size under fixed assumptions and sample limits.

Advanced 12 min read

What You'll Learn

  • Relate interim information fraction to the standard error of the interim treatment-effect estimate.
  • Calculate promising-zone entry probabilities at alternative interim timings.
  • Translate zone-entry probabilities into expected total sample sizes.
  • Recognize why timing effects can be nonmonotonic and depend on the true treatment effect.
  • Diagnose timing comparisons that inadvertently change other design assumptions.
  • Distinguish expansion probability and expected sample size from power and type I error.

Why Interim Timing Changes the Rule’s Behavior

A promising-zone rule uses the interim treatment-effect estimate to decide whether to increase the planned sample size. Moving the interim analysis earlier or later changes how noisy that estimate is. Even when the true treatment effect, zone cutoffs, and sample-size limits stay the same, the chance of landing inside the zone can therefore change. An early look generally gives a noisier estimate, but it does not follow that the probability of entering the zone must always increase or decrease: the result depends on where the true effect lies relative to both cutoffs.

A useful timing comparison isolates this one design choice. Keep the effect assumptions and sample-size limits fixed, vary the interim information fraction, and recalculate the probability of a zone decision and its consequences for expected sample size. These are operating characteristics of the re-estimation rule. They are not, by themselves, the probability of rejecting the null hypothesis or a demonstration of type I error control.

Definition: The interim information fraction \(t\) is the information available at the interim analysis divided by the information planned for the original final analysis, with \(0<t<1\). Under a model in which information is proportional to the number of observed participants, \(t\) is also the fraction of the original planned sample observed at interim. A promising-zone entry probability is the probability, under a specified true treatment effect, that the interim estimate falls between the rule’s lower and upper cutoffs.

Model and Assumptions for a Timing Comparison

Let \(\Delta\) denote the true treatment effect, measured in the outcome’s units, and let \(\hat{\Delta}_t\) be its interim estimate at information fraction \(t\). For the calculations below, use the normal approximation \(\hat{\Delta}_t\sim N(\Delta,SE_t^2)\), where \(SE_t=SE_0/\sqrt{t}\), and \(SE_0\) is the standard error at the original planned final sample size. This relation assumes information accumulates proportionally with sample size and that the variance used in planning is fixed. For example, it follows for a difference in means with known common variance and balanced allocation when the number of participants is proportional to information.

Assume the rule expands the sample when \(L\leq\hat{\Delta}_t\leq U\), where \(L\) and \(U\) are fixed effect-scale cutoffs with \(L<U\). If the estimate is outside that interval, participants continue to the original planned total \(N_0\); inside it, the total is increased to a fixed maximum \(N_{\max}\). Thus the comparison does not include early stopping, changing the zone cutoffs, or changing either sample-size limit. The estimate is assumed to follow the stated normal model at each timing, and the timing comparison holds \(\Delta\), \(SE_0\), \(L\), \(U\), \(N_0\), and \(N_{\max}\) fixed.

This is an operating-characteristic calculation for the zone decision. If the variance is estimated with substantial uncertainty, information does not accrue proportionally to sample size, or the interim estimate has a different distribution, the normal model may be inadequate. Those features should be represented in an appropriate design-level calculation or simulation rather than silently treated as satisfying these assumptions.

Probability of Entering the Zone and Expected Sample Size

Standardizing the zone boundaries gives a direct expression for the entry probability. Because \(\hat{\Delta}_t\) is normal with mean \(\Delta\) and standard deviation \(SE_0/\sqrt{t}\), the probability is the standard normal probability between the standardized cutoffs:

$$ P(\text{zone}\mid\Delta,t) = \Phi\!\left(\frac{U-\Delta}{SE_0/\sqrt{t}}\right) - \Phi\!\left(\frac{L-\Delta}{SE_0/\sqrt{t}}\right), $$

where \(\Phi\) is the cumulative distribution function of a standard normal random variable. Once the entry probability is known, the sample-size calculation follows from the two possible totals. The total is \(N_{\max}\) with probability \(P(\text{zone}\mid\Delta,t)\), and \(N_0\) otherwise. Therefore:

$$ E(N\mid\Delta,t) = N_0+(N_{\max}-N_0)P(\text{zone}\mid\Delta,t). $$

This expression concerns the final total under the stated rule. It does not mean every participant is recruited at the interim: recruitment continues after the interim in both cases, with the larger continuation used only after zone entry.

1
Fix the design assumptions.
Specify the true effect used for evaluation, the planned-final standard error, the zone cutoffs, the original total, and the maximum total.
2
Change the timing only.
For each information fraction \(t\), calculate \(SE_t=SE_0/\sqrt{t}\). A smaller fraction means less interim information and a larger standard error.
3
Calculate and compare operating characteristics.
Standardize both cutoffs, find the zone probability, then use it to calculate the expected total sample size. Keep intermediate precision until the final rounding.

Worked Example: Comparing an Early Look With a Later Look

Suppose a hypothetical trial has an original planned total of \(N_0=200\) participants and a maximum of \(N_{\max}=300\). The planned-final standard error of the estimated treatment effect is \(SE_0=0.20\) outcome units. The fixed promising zone is \(0.20\leq\hat{\Delta}_t\leq0.45\) outcome units, and the true effect for this calculation is \(\Delta=0.30\) outcome units. Compare interim fractions \(t=0.25\) and \(t=0.80\). The normal model and proportional information assumption are taken to hold for this hypothetical calculation.

At \(t=0.25\), \(SE_t=0.20/\sqrt{0.25}=0.40\) outcome units. The standardized cutoffs are \((0.20-0.30)/0.40=-0.25\) and \((0.45-0.30)/0.40=0.375\). Thus \(P(\text{zone}\mid\Delta=0.30,t=0.25)=\Phi(0.375)-\Phi(-0.25)\approx0.64617-0.40129=0.24488\).

The expected total is \(200+(300-200)(0.24488)=224.488\), or about 224.5 participants. The rule expands to 300 participants about 24.5% of the time under the specified true effect and model.

At \(t=0.80\), \(SE_t=0.20/\sqrt{0.80}\approx0.2236068\) outcome units. The standardized cutoffs are \((0.20-0.30)/0.2236068\approx-0.4472136\) and \((0.45-0.30)/0.2236068\approx0.6708204\). Using the normal probabilities before rounding, \(P(\text{zone}\mid\Delta=0.30,t=0.80)\approx0.74883-0.32736=0.42147\).

The corresponding expected total is \(200+100(0.42147)\approx242.147\), or about 242.1 participants. In this example, moving the interim later raises the zone-entry probability by approximately \(0.42147-0.24488=0.17659\), or 17.7 percentage points, and raises expected enrollment by about 17.7 participants. The cutoffs and sample-size limits did not change; the interim estimate became more precise.

Why the Effect’s Location Matters

The standardized interval changes at both ends as timing changes. A later look reduces the standard error, so the standardized cutoffs move farther from zero in the direction determined by their distance from \(\Delta\). The resulting probability depends on both distances, \(L-\Delta\) and \(U-\Delta\). It is not generally a monotone function of \(t\). In particular, a true effect near the middle of the zone can have a high chance of remaining inside it with a precise estimate, while a noisy estimate can either spread probability into the zone or spread it out of the zone. For a true effect outside the zone, precision may instead make estimates more consistently fall on that same outside side.

The next example holds the effect and all other design assumptions fixed across the timing comparison, and shows why an intermediate timing need not behave like a simple step between early and late looks. The sample-size expectation uses the same \(N_0=200\) and \(N_{\max}=300\) as above.

Worked Example: A Nonmonotonic Pattern Across Timings

Use the same hypothetical assumptions: \(\Delta=0.30\) outcome units, \(SE_0=0.20\), and zone \(0.20\) to \(0.45\) outcome units. Compare \(t=0.25\), \(t=0.50\), and \(t=0.80\). At \(t=0.25\), the preceding calculation gives a zone probability of approximately 0.24488 and expected total of 224.5 participants.

At \(t=0.50\), the standard error is \(0.20/\sqrt{0.50}\approx0.2828427\). The standardized cutoffs are approximately \(-0.3535534\) and \(0.5303301\). Therefore, \(P(\text{zone}\mid\Delta=0.30,t=0.50)\approx\Phi(0.5303301)-\Phi(-0.3535534)\approx0.70206-0.36184=0.34022\). The expected total is \(200+100(0.34022)=234.022\), or about 234.0 participants.

At \(t=0.80\), the probability is approximately 0.42147 and the expected total is about 242.1 participants. Across these three timings, the zone-entry probability rises from about 0.2449 to 0.3402 to 0.4215. That pattern is specific to this true effect and these cutoffs; it is not a general rule that later interim analyses always increase zone entry.

The arithmetic also illustrates why probabilities should be calculated from unrounded normal CDF values, then rounded at the end. Rounding each CDF value early and subtracting can change the final reported probability, and that discrepancy carries into the expected sample size.

Worked Example: A Different True Effect Gives a Different Timing Pattern

Keep the same cutoffs, planned-final standard error, original total, and maximum total, but evaluate the rule at a different hypothetical true effect, \(\Delta=0.10\) outcome units. Compare \(t=0.25\) with \(t=0.80\). The purpose is to show that the effect assumption matters; this is not a change made within either timing comparison.

At \(t=0.25\), \(SE_t=0.40\). The standardized cutoffs are \((0.20-0.10)/0.40=0.25\) and \((0.45-0.10)/0.40=0.875\). Thus the zone probability is \(\Phi(0.875)-\Phi(0.25)\approx0.80921-0.59871=0.21051\), giving an expected total of \(200+100(0.21051)\approx\) 221.1 participants.

At \(t=0.80\), \(SE_t\approx0.2236068\). The standardized cutoffs are approximately \(0.4472136\) and \(1.5652476\), so the zone probability is about \(\Phi(1.5652476)-\Phi(0.4472136)\approx0.9412-0.6726=0.2686\). The expected total is approximately \(200+100(0.2686)=\) 226.9 participants.

In both timings the true effect is below the lower zone cutoff, yet the later look has a higher probability of entering the zone under this particular model and interval. That result is possible because the zone is not centered on the true effect and its two boundaries are at different distances from it. It reinforces the need to calculate, rather than infer, the direction of a timing effect.

Diagnostics, Alternatives, and Limitations

A clean sensitivity analysis changes the information fraction while holding the other inputs fixed. Check that the interim standard error follows the assumed information relationship; that the effect, cutoffs, and sample-size limits are identical across rows; and that the expected total lies between \(N_0\) and \(N_{\max}\). If a calculation yields an expected total outside those limits, or a zone probability outside \([0,1]\), there is an arithmetic or implementation error. Also inspect the standardized endpoints: if they have been calculated using the wrong sign or the wrong standard error, the resulting probability can look plausible while being wrong.

In a real design, compare analytic calculations with simulation when the normal approximation, known variance, or proportional information assumption is questionable. Simulation can represent the planned data-generating model, interim timing, actual re-estimation rule, and sample-size caps together. Monte Carlo uncertainty should be reported when simulated operating characteristics are used. A conditional-power rule is another alternative: unlike fixed effect cutoffs, it bases the decision on a conditional probability of the final result under specified assumptions. That rule can also be timing-sensitive, but it is a different decision rule and should not be substituted into a fixed-zone calculation without changing the design being evaluated.

Zone-entry probability and expected sample size do not determine power. To evaluate the probability of a final rejection, one must specify the final testing procedure and calculate its behavior under the adaptive design. In particular, one should not assume that a conventional fixed-sample final test automatically retains its nominal type I error after the sample size depends on an unblinded interim estimate. That inferential question requires separate evaluation; a timing comparison of expansion probabilities does not resolve it.

Finally, the normal calculations here are conditional on a specified true treatment effect. They do not account for uncertainty about which effect is the right planning value, variance misspecification, deviations from the assumed information growth, or operational delays and recruitment constraints. Sensitivity to timing is useful because it reveals how the rule behaves under its stated assumptions, not because it removes those limitations. Report the assumptions and the rule’s full operating characteristics together.

Key takeaway: With fixed treatment-effect assumptions, zone cutoffs, and sample-size limits, interim timing affects the precision of the interim estimate and therefore the probability of sample-size expansion. Recalculate the zone probability at each timing; do not assume the effect is monotonic, and do not equate expansion probability with power or error-rate control.

Check Your Understanding

Use the normal model and expected-sample-size formula in this tutorial. Solutions are included after the questions.

  1. For \(\Delta=0.30\), \(SE_0=0.20\), zone \(0.20\) to \(0.45\), and \(t=0.50\), what are the interim standard error and approximate zone-entry probability?
  2. If \(N_0=200\), \(N_{\max}=300\), and the zone-entry probability is 0.42147, calculate the expected total sample size.
  3. With the same cutoffs, \(SE_0=0.20\), and \(t=0.80\), calculate the zone-entry probability when \(\Delta=0.20\).
  4. Why can a later interim analysis either raise or lower the probability of zone entry?
  5. Does a higher zone-entry probability by itself establish higher power or controlled type I error? Explain briefly.

Solutions. (1) \(SE_t=0.20/\sqrt{0.50}\approx0.28284\). The standardized cutoffs are approximately \(-0.35355\) and \(0.53033\), giving \(\Phi(0.53033)-\Phi(-0.35355)\approx0.3402\). (2) \(E(N)=200+100(0.42147)=242.147\), or about 242.1 participants. (3) The standard error is \(0.2236068\); the standardized cutoffs are \(0\) and \(0.25/0.2236068\approx1.11803\). Thus the probability is \(\Phi(1.11803)-\Phi(0)\approx0.8682-0.5000=0.3682\). (4) Timing changes the estimate’s standard error, while the true effect can be at different distances from the two zone boundaries. Consequently, probability mass can move into or out of the zone; the direction is not universal. (5) No. Zone entry describes a sample-size decision. Power and type I error depend on the final test and its behavior under the adaptive design.

Clinical Biostats