Tutorials › Group Sequential Methods and Sample Size Re-Estimation › Rounding Lower Boundaries In Practice

Futility Boundaries and Beta Spending · Tutorial 188 of 699

Rounding Lower Boundaries In Practice

See how small changes to a lower futility boundary affect crossing probabilities and how to verify the rounded design before implementation.

Advanced 12 min read

What You'll Learn

  • Explain why rounding a lower futility boundary upward can increase the chance of stopping for futility.
  • Calculate the change in a conditional crossing probability when a threshold is rounded.
  • Distinguish a conditional probability at one observed history from an overall operating characteristic.
  • Verify operating characteristics under the null and a prespecified alternative.
  • Choose a rounding direction that respects a stated futility-probability target.
  • Identify diagnostics for errors caused by stored or implemented boundaries.

Why Boundary Rounding Needs a Design Check

A computed lower futility boundary is often stored or reported to fewer decimal places than the software used to calculate it. The numerical change may look negligible, but a lower boundary determines which interim results trigger stopping: the statistic stops for futility when it is at or below the boundary. Moving the boundary upward therefore includes additional outcomes in the stopping region. Whether this change matters depends on the probability mass near the boundary and on the design’s joint monitoring rule.

The relevant check is not merely whether the rounded number is close to the computed number. It is whether using the rounded boundary changes the operating characteristics enough to affect the design’s intended error rates, power, or futility-spending targets. This tutorial focuses on that check. Earlier material on binding and nonbinding futility, attained beta, and simulation provides the decision-rule and operating-characteristic context.

Definition: At look \(i\), let \(l_i\) be a lower futility boundary and \(Z_i\) the monitoring statistic. The rule stops for futility at that look when \(Z_i\leq l_i\), provided the trial has not already stopped. If the computed boundary is \(\tilde l_i\), rounding replaces it with a stored value \(l_i^{(r)}\). The rounded rule is evaluated using \(l_i^{(r)}\), not the unrounded value.

For clarity, suppose \(l_i^{(r)}> \tilde l_i\), as happens when a negative boundary is rounded upward toward zero. Then the rounded rule stops for futility on the additional interval \(\tilde l_i<Z_i\leq l_i^{(r)}\). Rounding downward has the reverse effect: it removes outcomes from the futility-stopping region. These statements concern the boundary rule itself; overall effects can also depend on whether earlier efficacy or futility decisions prevent the trial from reaching look \(i\).

Model and Conditional Probability Change

For the calculations below, use the standard canonical group sequential model at increasing information fractions \(0<t_1<\cdots<t_K=1\). Under a fixed standardized effect \(\theta\), assume the vector of standardized statistics is jointly normal, each statistic has variance one, and

$$ E_\theta(Z_i)=\theta\sqrt{t_i}, \qquad \operatorname{Corr}(Z_i,Z_j)=\sqrt{\frac{t_i}{t_j}}\quad (i<j). $$

This model is appropriate when the information-based normal approximation and the planned analysis schedule are adequate. It is not a substitute for checking the actual analysis model when, for example, the statistic has a materially different distribution or the information fractions do not represent the planned information. At each look, assume the stopping rule is applied to the statistic and boundary actually specified in the monitoring plan.

In this model, the conditional distribution at a later look, given the immediately preceding statistic, is also normal. If \(t_{i-1}<t_i\), then

$$ Z_i\mid Z_{i-1}=z \sim N\left( \theta\sqrt{t_i} + \sqrt{\frac{t_{i-1}}{t_i}} \left[z-\theta\sqrt{t_{i-1}}\right], \; 1-\frac{t_{i-1}}{t_i} \right). $$

Write the conditional mean as \(m\) and the conditional standard deviation as \(s_c\), the square root of the variance shown above. For an upward-rounded boundary \(l^{(r)}>l\), the conditional increase in the chance of crossing the futility boundary, given that history, is exactly

$$ P(l<Z_i\leq l^{(r)}\mid\text{history}) = \Phi\left(\frac{l^{(r)}-m}{s_c}\right) - \Phi\left(\frac{l-m}{s_c}\right). $$

This expression gives a direct local sensitivity calculation. It is not, by itself, the change in the overall probability of stopping for futility: the overall calculation must also account for the probability of reaching this look and the distribution of earlier statistics among trials that continue.

Worked Boundary Changes

Worked Example: A One-Look Rounding Change

Hypothetical data. Suppose a one-look design has a computed lower futility boundary of \(-0.554\), but the monitoring document records it to two decimal places as \(-0.55\). Under a specified alternative, assume \(Z\sim N(1,1)\). Compare the futility-crossing probabilities.

With the computed boundary,

$$ P(Z\leq -0.554) = \Phi(-0.554-1) = \Phi(-1.554) \approx 0.0601. $$

With the rounded boundary,

$$ P(Z\leq -0.55) = \Phi(-0.55-1) = \Phi(-1.55) \approx 0.0606. $$

Thus, rounding upward increases the futility-crossing probability by approximately \(\Phi(-1.55)-\Phi(-1.554)\approx 0.00048\), or \(0.048\) percentage points. The probability of continuing past this look falls by the same amount, from approximately \(0.93991\) to \(0.93943\). This calculation checks the effect at this one look under the stated normal model; in a multi-look design, the overall operating characteristics require the joint stopping rule.

Worked Example: Conditional Crossing at a Later Look

Hypothetical design. Consider two looks at \(t_1=0.50\) and \(t_2=1\), under \(\theta=2\). Suppose the first-look statistic is \(Z_1=-0.20\), and compare a computed second-look lower boundary of \(-0.60\) with an upward-rounded value of \(-0.55\). This is a conditional calculation for the specified first-look result, not an overall trial probability.

The means are \(E(Z_1)=2\sqrt{0.50}=1.4142\) and \(E(Z_2)=2\); the correlation is \(\sqrt{0.50}=0.7071\). The conditional mean and standard deviation of \(Z_2\) are therefore

$$ m = 2+\sqrt{0.50}\left[-0.20-2\sqrt{0.50}\right] = 0.8586, \qquad s_c=\sqrt{1-0.50}=\sqrt{0.50}=0.7071. $$

At the computed boundary, the conditional crossing probability is

$$ P(Z_2\leq -0.60\mid Z_1=-0.20) = \Phi\left(\frac{-0.60-0.858579}{\sqrt{0.50}}\right) = \Phi(-2.063) \approx 0.01957. $$

At the rounded boundary, it is

$$ P(Z_2\leq -0.55\mid Z_1=-0.20) = \Phi\left(\frac{-0.55-0.858579}{\sqrt{0.50}}\right) = \Phi(-1.992) \approx 0.02318. $$

The conditional increase is approximately \(0.02318-0.01957=0.00361\), or \(0.361\) percentage points, for this specified earlier statistic. The increase is larger than in the one-look example because the conditional distribution’s standard deviation is \(0.7071\), and the relevant probability mass between these two thresholds is correspondingly different. It would be incorrect to report \(0.003615\) as the design’s unconditional increase: that requires averaging over the distribution of first-look results that reach look 2.

From a Local Change to Operating Characteristics

An operating characteristic is a probability defined by the complete stopping rule under a specified parameter value. Examples include the probability of rejecting \(H_0\) under the null, the probability of rejecting under an alternative (power), and the probability of stopping for futility. A rounded boundary can affect these probabilities differently. In particular, futility crossing probability is not automatically the same as overall type II error: the decision rule and whether futility is binding determine how the events relate.

For a binding futility rule, an exact calculation partitions the joint normal distribution into the regions that correspond to stopping for efficacy, stopping for futility, and continuing. For example, with lower boundary \(l_1\) and upper efficacy boundary \(u_1\) at the first look, continuation requires \(l_1<Z_1<u_1\). If the lower boundary is rounded upward, the first-look futility region expands, and the continuation region contracts. At later looks, only histories that have not already stopped contribute to that look’s unconditional stopping probability.

A direct verification compares the original and rounded designs under the same model, schedule, parameter values, and efficacy boundaries. Evaluate the full stopping rule once with the computed \(l_i\) values and once with the stored \(l_i^{(r)}\) values. Check the planned quantities under \(H_0\) and under the design alternative, as well as any futility-spending targets used to construct the boundaries. A check only at the alternative can miss a change in type I error; a check only under the null may not reveal a material change in power or futility behavior.

1
Record both boundary sets.
Keep the full-precision computed values and the exact values to be implemented. State the rounding rule and confirm that each rounded lower boundary remains below its corresponding upper efficacy boundary.
2
Recalculate the complete rule.
Use the planned joint distribution, information fractions, and stopping decisions. Recompute the operating characteristics for both boundary sets rather than adjusting a reported probability by a local conditional difference.
3
Compare against design goals.
Report the absolute changes in the relevant probabilities, in percentage points where helpful. Check that the rounded design still meets the prespecified tolerances for error rates, power, and futility behavior.
4
Confirm implementation.
Verify that the boundary used by the analysis code matches the value in the monitoring plan, including its sign, comparison convention (\(\leq\) versus \(<\)), and precision.

Choosing a Rounding Direction

Worked Example: Protecting a Crossing-Probability Target

Hypothetical one-look calculation. Suppose \(Z\sim N(1,1)\) under the alternative and the target futility-crossing probability at this look is \(0.05\). The boundary that attains this target solves \(\Phi(l-1)=0.05\), so

$$ l=1+\Phi^{-1}(0.05) \approx 1-1.6449 = -0.6449. $$

If a two-decimal boundary must be used and the goal is not to exceed the target, round downward to \(-0.65\). Then

$$ P(Z\leq -0.65) = \Phi(-1.65) \approx 0.04947, $$

which is below \(0.05\). Rounding upward to \(-0.64\) instead gives

$$ P(Z\leq -0.64) = \Phi(-1.64) \approx 0.05050, $$

which exceeds the target. This illustrates the direction of the effect and a conservative choice for this particular probability constraint. In a multi-look design, rounding every boundary downward is not a general guarantee that all operating-characteristic constraints are met; the complete design still needs to be recalculated.

Diagnostics, Alternatives, and Limitations

A useful diagnostic is to check the probability mass between the full-precision and implemented boundaries. If that mass is tiny under the relevant conditional distributions, the local effect is likely to be small; if substantial density lies near a boundary, even a small numerical change may matter. The overall impact also depends on how often the design reaches the look and on the prior history of trials that continue. Consequently, a tiny absolute change in a boundary is not, on its own, evidence that the operating characteristics are unchanged.

Exact joint-normal integration is the direct verification approach when the canonical model is appropriate and the calculation is available. It avoids Monte Carlo sampling error and permits close comparisons of the original and rounded designs. Simulation is a useful independent check and can accommodate more complex implementation details, but finite simulation produces random error. When comparing two designs by simulation, use the same generated paths for both where possible; this paired comparison reduces noise in their difference. Report the number of simulated trials and the uncertainty of estimated probabilities rather than treating simulation estimates as exact.

Common warning signs include a rounded lower boundary that is numerically below the computed value when the documented rule says “round upward,” an unexpected decrease in futility stopping after upward rounding, or an unexplained discrepancy between code and the monitoring plan. Also investigate a rounded boundary that meets or exceeds the efficacy boundary: the intended continuation region may then be empty or incoherent. These checks can catch sign, precision, and inequality errors before they affect decisions.

Rounding is not the only source of discrepancy. The normal model may be an approximation; information fractions or effect assumptions may be misspecified; and numerical integration or simulation may be inaccurate. Rounding checks establish the consequences of using particular thresholds under the model being evaluated. They do not establish that the model itself is appropriate, nor do they replace a prespecified tolerance for acceptable changes in operating characteristics.

Common Mistake

Common mistake: “The boundary changed by only 0.01, so the probability of stopping must have changed by only 0.01.” A boundary is on the scale of \(Z_i\), not the probability scale. The probability change depends on the distribution near that boundary, conditional on reaching the look. Calculate the probability mass between the two thresholds, then verify the complete operating characteristics under the rounded rule.

Check Your Understanding

Answer each question using the stated rule. Solutions follow the questions.

  1. A lower futility boundary changes from \(-0.72\) to \(-0.70\). Which rule stops for futility on more possible statistic values?
  2. Why is a conditional crossing-probability change at one observed interim result not generally the unconditional change in futility probability?
  3. Under \(Z\sim N(0.8,1)\), calculate the futility probability for boundary \(-0.50\). What is it if the boundary is rounded upward to \(-0.48\)?
  4. Name two operating characteristics that should be checked when verifying rounded boundaries, and explain why checking only under the alternative is insufficient.
  5. What implementation details should be checked if the calculated and simulated results disagree?

Solutions. (1) The \(-0.70\) boundary stops on more values because the rule is \(Z\leq l\), and \(-0.70> -0.72\). (2) The unconditional probability also weights by how often the trial reaches the look and by the distribution of histories among trials that continue. (3) At \(-0.50\), \(P(Z\leq-0.50)=\Phi(-1.30)\approx0.0968\). At \(-0.48\), it is \(\Phi(-1.28)\approx0.1003\); upward rounding increases the probability by about \(0.0035\), or \(0.35\) percentage points. (4) Check type I error under \(H_0\) and power under the alternative, along with any planned futility target. An alternative-only check does not establish null error control. (5) Check the boundary values and precision, sign, comparison convention, look schedule, and whether the simulated rule matches the rule used in the calculation.

Clinical Biostats