Tutorials › Group Sequential Methods and Sample Size Re-Estimation › Futility Spending With Unequal Information Fractions

Futility Boundaries and Beta Spending · Tutorial 181 of 699

Futility Spending With Unequal Information Fractions

Calculate lower futility thresholds at uneven information fractions, accounting for paths that have already crossed a boundary and for the distinction between information time and calendar time.

Advanced 12 min read

What You'll Learn

  • Express interim timing as information fractions rather than calendar-time fractions.
  • Translate cumulative beta spending into the probability target for a new futility crossing.
  • Calculate a two-look lower boundary while accounting for the dependence between Z statistics.
  • Distinguish a marginal tail probability from the probability of a new crossing.
  • Diagnose thresholds based on calendar time or on unadjusted marginal quantiles.
  • Recognize assumptions and limitations of information-based calibration.

Why Unequal Information Fractions Matter

Interim analyses are often scheduled at unequal intervals in calendar time. That fact alone does not tell us how much information has accumulated at each look. Enrollment may be faster or slower than expected, event accrual may vary, or measurements may arrive on a different schedule. Futility boundaries are calibrated using the statistical information available at each analysis, not the number of days since the trial began.

This distinction matters because the standardized test statistics at different looks are dependent, with dependence determined by their information fractions. When a later boundary is calibrated, its target is the probability of a new futility crossing among paths that did not cross earlier. Unequal information fractions affect both that dependence and the conditional distribution used in the calculation.

Definition: Let \(I_k\) denote the information at interim look \(k\), and let \(I_K\) be the planned final information. The information fraction at look \(k\) is \(t_k=I_k/I_K\), with \(0<t_1<\cdots<t_K=1\). A beta-spending function \(B(t)\) specifies the cumulative target probability of futility crossing by information fraction \(t\), under the chosen design alternative. The incremental target at look \(k\) is \(\Delta B_k=B(t_k)-B(t_{k-1})\), where \(B(t_0)=0\).

For the calculations below, assume the usual canonical group sequential model: information is correctly represented by \(I_k\), the standardized statistics have a joint normal distribution, and their increments are independent when expressed on the information scale. Under the design alternative, write \(\theta=\delta\sqrt{I_K}\), where \(\delta\) is the alternative effect on the scale used to standardize the statistic. Then \(Z_k\) has mean \(\theta\sqrt{t_k}\) and variance 1, while for \(j<k\),

$$ \operatorname{Corr}(Z_j,Z_k)=\sqrt{\frac{t_j}{t_k}}. $$

These assumptions are the basis for the standard joint-normal calibration. In practice, the statistic and information measure must be appropriate for the trial; if the canonical model is a poor description, the calculated boundaries may not achieve their nominal operating characteristics.

Calibrating the New Crossing at Each Look

Let \(l_k\) be the lower futility boundary at look \(k\). A path crosses for futility at look \(k\) if \(Z_k\leq l_k\) and it has not crossed at any earlier look. Thus the beta-spending calibration is

$$ P_A\left(Z_k\leq l_k,\ Z_j>l_j\text{ for every }j<k\right)=\Delta B_k, $$

where \(P_A\) denotes probability under the design alternative. Summing these mutually exclusive first-crossing probabilities gives \(B(t_k)\), the cumulative target through look \(k\). At the first look, the equation is simply \(P_A(Z_1\leq l_1)=B(t_1)\). At a later look, the equation must exclude paths that already crossed.

For two looks, put \(\rho=\sqrt{t_1/t_2}\), \(\mu_k=\theta\sqrt{t_k}\), and \(s=\sqrt{1-\rho^2}\). The joint normal model implies

$$ P_A(Z_2\leq l_2,\ Z_1>l_1) = \int_{-\infty}^{l_2-\mu_2} \phi(y)\, \Phi\left(\frac{\mu_1+\rho y-l_1}{s}\right)\,dy, $$

where \(\phi\) and \(\Phi\) are the standard normal density and distribution functions. To see why, condition on \(Z_2=\mu_2+y\). The conditional distribution of \(Z_1\) has mean \(\mu_1+\rho y\) and variance \(1-\rho^2\). The factor involving \(\Phi\) is therefore the conditional probability that the first-look statistic is above its futility boundary. Integrating over values of \(Z_2\) below \(l_2\) yields the probability of a new crossing.

Worked Example: A Two-Look Design With Unequal Information Fractions

Suppose the information fractions are \(t_1=0.50\) and \(t_2=1.00\), and under the design alternative \(\theta=1.549193\). The planned cumulative beta spending is \(B(t_1)=0.030\) and \(B(t_2)=0.080\). Find the lower boundaries.

The means are \(\mu_1=1.549193\sqrt{0.50}=1.095445\) and \(\mu_2=1.549193\). At the first look the marginal target is 0.030, so

$$ l_1=\mu_1+\Phi^{-1}(0.030) =1.095445-1.88079 \approx -0.7853. $$

The new-crossing target at look 2 is not 0.080; it is the increment \(0.080-0.030=0.050\). Here \(\rho=\sqrt{0.50}=0.707107\) and \(s=0.707107\). Substituting into the conditional probability gives

$$ \frac{\mu_1+\rho y-l_1}{s} \approx y+2.660, \qquad P_A(Z_2\leq l_2,\ Z_1>l_1) = \int_{-\infty}^{l_2-1.549193} \phi(y)\Phi(y+2.660)\,dy. $$

Numerically solving this equation with the right-hand side set to 0.050 gives \(l_2\approx0.047\). At this value the integration limit is \(0.047-1.549193=-1.5022\), and the integral is approximately 0.050. The cumulative futility-crossing probability is therefore approximately \(0.030+0.050=0.080\), as targeted. The calculation accounts for the fact that some paths with \(Z_2\leq l_2\) would already have crossed at look 1.

Unequal Calendar Spacing Is Not Unequal Beta Spending

The spending function is evaluated at the observed or planned information fractions. It is not evaluated at the fraction of calendar time elapsed. Nor does an uneven calendar schedule automatically imply a front-loaded or back-loaded beta-spending schedule. Those are distinct design choices: the calendar determines when analyses occur, while \(B(t)\) determines the cumulative crossing target as information accrues.

Worked Example: Using Information Rather Than Calendar Time

Suppose a planned interim analysis occurs halfway through the calendar duration of a trial, but only 30 of the planned 100 information units have accrued. Thus \(t=30/100=0.30\), not 0.50. Assume \(\theta=2\), there has been no earlier futility look, and the spending function is \(B(t)=0.20t\). Find the lower boundary.

At the actual information fraction, the target crossing probability is \(B(0.30)=0.20(0.30)=0.060\). The mean under the design alternative is \(\mu=2\sqrt{0.30}=1.0954\). Since this is the first futility look, the marginal equation is appropriate:

$$ l=\mu+\Phi^{-1}(0.060) =1.0954-1.5548 \approx -0.4594. $$

Using the calendar fraction instead would assign \(B(0.50)=0.100\), which is not the specified cumulative spending at the information actually available. It would also misstate the mean if the statistic were evaluated at \(t=0.50\): \(2\sqrt{0.50}=1.4142\), rather than 1.0954. Calendar time is relevant to scheduling and monitoring logistics, but it does not replace the information fraction in either the spending target or the canonical mean.

Marginal Probability Is Not a New-Crossing Probability

A common shortcut is to choose \(l_k\) so that \(P_A(Z_k\leq l_k)=\Delta B_k\). That is generally incorrect after the first look. The marginal probability includes paths that crossed at earlier analyses. The target \(\Delta B_k\), by contrast, applies only to surviving paths. A marginal quantile can be useful for a first look or for a look with no earlier possibility of crossing, but later boundaries require the joint crossing calculation.

Worked Example: Diagnosing an Unadjusted Second-Look Boundary

Use the two-look model and first boundary from the preceding two-look example: \(\mu_2=1.549193\), \(t_1/t_2=0.50\), and \(l_1\approx-0.7853\). Suppose someone proposes \(l_2=0.09\) because they are comparing a marginal tail probability with the desired new-crossing probability.

First calculate the marginal probability:

$$ P_A(Z_2\leq0.09) =\Phi(0.09-1.549193) =\Phi(-1.459193) \approx0.0723. $$

The new-crossing probability is smaller. Its integral has upper limit \(-1.459193\), so

$$ P_A(Z_2\leq0.09,\ Z_1>l_1) = \int_{-\infty}^{-1.459193}\phi(y)\Phi(y+2.660)\,dy \approx0.055. $$

Consequently, approximately \(0.0723-0.055=0.0173\) of the marginal probability lies on paths already counted at look 1. The proposed boundary does not achieve a new-crossing target of 0.050: its new-crossing probability is about 0.055. The numerical root near \(l_2=0.047\), not the marginal-quantile calculation, is the boundary that meets that target.

Calculation Checks and Design Limitations

A numerical solution is credible only if its inputs and event definition match the design. For a planned schedule, verify that information fractions are strictly increasing and that the cumulative spending values are nondecreasing. Each increment \(\Delta B_k\) must be nonnegative. Also check that the requested new-crossing probability is attainable: it cannot exceed the probability of reaching look \(k\) without a previous futility crossing. A solver failure or a boundary that moves unexpectedly can signal an incoherent spending schedule, incorrect information fractions, or a mistake in the survival event.

For more than two looks, the event at look \(k\) requires survival of every earlier lower boundary. The joint probability can be evaluated using the canonical multivariate normal distribution or by integrating sequential conditional distributions. Using only the immediately preceding statistic is justified when the canonical independent-increment model gives the required Markov conditional structure; the full event must still exclude crossings at all earlier looks.

Compared with direct joint-probability calibration, a marginal-quantile rule is simpler but does not generally spend the specified incremental beta after the first look. Assigning equal beta increments at equally spaced calendar dates is another tempting substitute; it is valid only if those increments are actually the chosen spending schedule at the information fractions, not merely because the dates are evenly spaced. The standard beta-spending approach is flexible about timing, but it depends on a correctly specified information scale and design alternative. If the effect used for calibration is poorly chosen, or if information fractions differ materially from the planned values, actual futility probabilities can differ from the design targets.

Key takeaway: At an interim look, evaluate the spending function at the information fraction, then calibrate the lower boundary to the incremental probability of a first futility crossing. Calendar-time spacing and marginal tail probabilities do not perform those two jobs.

Check Your Understanding

Answer each question, then compare with the solutions below.

  1. Information fractions are \(0.25\) and \(0.75\), with cumulative spending \(B(0.25)=0.02\) and \(B(0.75)=0.07\). What new-crossing probability should the second boundary target?
  2. Why does a later-look boundary usually not satisfy \(P_A(Z_k\leq l_k)=\Delta B_k\)?
  3. An analysis occurs at 40% of the planned calendar duration, but the information fraction is 0.28. At which fraction is the spending function evaluated?
  4. Under the canonical model, what is the correlation between statistics at information fractions \(t_j=0.36\) and \(t_k=0.81\), for \(j<k\)?

Solutions. (1) The increment is \(0.07-0.02=0.05\), so the new-crossing probability target is 0.050. (2) The marginal event \(Z_k\leq l_k\) includes paths that crossed at an earlier look, whereas \(\Delta B_k\) applies only to paths surviving all earlier futility boundaries. (3) Evaluate spending at \(t=0.28\), the information fraction; calendar time does not determine the information-based target. (4) The correlation is \(\sqrt{0.36/0.81}=\sqrt{4/9}=2/3\approx0.6667\).

Clinical Biostats