Tutorials › Bayesian Clinical Trials › Posterior Futility Rules for Arm Removal

Adaptive and Platform Designs · Tutorial 557 of 696

Posterior Futility Rules for Arm Removal

Learn to calculate and interpret a posterior futility probability for a binary-outcome experimental arm and apply a prespecified arm-removal threshold.

Advanced 10 min read

What You'll Learn

  • Define futility using the posterior probability that an experimental arm’s risk difference reaches a target effect.
  • Update independent beta priors with interim binary-outcome data.
  • Calculate the posterior target-effect probability by integrating over both response-rate distributions.
  • Apply a prespecified probability threshold consistently across experimental arms.
  • Diagnose how prior choice, interim sample size, and outcome-model assumptions can affect an arm-removal decision.
  • Distinguish posterior futility from superiority rules and predictions about future trial outcomes.

Why Futility Is About the Target, Not Just Superiority

At an interim analysis, an experimental arm may look unlikely to be better than control. But “better” can mean any positive difference, including one too small to matter for the trial’s prespecified objective. A futility rule instead asks whether the arm still has a sufficiently plausible chance of reaching a specified target effect. If that posterior probability is very small, the protocol may stop assigning new participants to the arm.

This tutorial uses a binary endpoint, such as response by a defined assessment time. Let \(p_E\) denote the response probability in an experimental arm and \(p_C\) the response probability in the control arm. Define the risk difference as \(p_E-p_C\), measured in probability units. A positive target \(\delta\) represents the effect level the rule asks the arm to achieve. The futility threshold \(\gamma\) is a small probability chosen before applying the rule.

Definition: For an experimental arm, its posterior futility probability is \(q=P(p_E-p_C\geq\delta\mid D)\), where \(D\) is the interim data and \(\delta\) is the prespecified target risk difference. A posterior futility rule removes the arm from further enrollment when \(q<\gamma\), for a prespecified threshold \(\gamma\). The rule concerns the probability of reaching the target effect, not the probability that the null hypothesis is true.

The target \(\delta\) and threshold \(\gamma\) answer different questions. The target says what effect counts as reaching the goal; the threshold says how little posterior support for that goal is enough to trigger the rule. Neither should be selected after looking at the interim result. A value of \(\gamma=0.10\), for example, means the decision criterion is a posterior probability below 10%; it does not assert a 10% frequentist type I error rate.

Binary-Outcome Model and Calculation

Suppose an interim analysis contains \(x_E\) responses among \(n_E\) participants assigned to the experimental arm and \(x_C\) responses among \(n_C\) participants assigned to control. Conditional on the response probabilities, model the counts as \(X_E\sim\mathrm{Binomial}(n_E,p_E)\) and \(X_C\sim\mathrm{Binomial}(n_C,p_C)\). A convenient specification assigns independent beta priors, \(p_E\sim\mathrm{Beta}(a_E,b_E)\) and \(p_C\sim\mathrm{Beta}(a_C,b_C)\). With these assumptions, the interim posteriors are:

$$ p_E\mid D\sim\mathrm{Beta}(a_E+x_E,\ b_E+n_E-x_E),\qquad p_C\mid D\sim\mathrm{Beta}(a_C+x_C,\ b_C+n_C-x_C). $$

The beta prior parameters are positive. They describe prior distributions on probabilities, not extra observed participants unless a separate interpretation has been explicitly justified. The binomial likelihood assumes a common response probability within each arm for the endpoint and assessment definition being analyzed. The model also treats participants’ outcomes as conditionally independent given their arm-specific probabilities. The independent priors encode prior independence of the two arm probabilities; they are not an automatic requirement of all Bayesian models.

Write \(f_E(u)\) for the experimental-arm posterior density and \(F_C(v)\) for the control-arm posterior cumulative distribution function. Conditional posterior independence lets us integrate over the control probability for each possible experimental probability:

$$ q=P(p_E-p_C\geq\delta\mid D) =\int_{\delta}^{1} f_E(u)F_C(u-\delta)\,du. $$

This expression follows because, for a fixed \(p_E=u\), the event \(p_E-p_C\geq\delta\) is equivalent to \(p_C\leq u-\delta\). If \(u-\delta\leq 0\), that conditional probability is zero; if \(u-\delta\geq 1\), it is one. Numerical integration or posterior simulation can evaluate \(q\). The resulting posterior probability is compared with \(\gamma\); it is not itself a test statistic with a frequentist reference distribution.

Assumptions to check: The endpoint and follow-up window are consistently defined; the analyzed participants and responses are correctly classified; the arm-specific binomial likelihood is appropriate; the priors are stated and applied as specified; and the interim analysis follows the protocol’s rules for available outcomes and arm comparisons. If these conditions are doubtful, the computed probability may not answer the intended question.

Worked Interim Decisions

Worked Example: A Small Probability of Reaching the Target

Suppose a hypothetical trial uses independent \(\mathrm{Beta}(1,1)\) priors for both response probabilities. At an interim analysis, 8 of 20 experimental-arm participants and 10 of 20 control participants have responded. The target is \(\delta=0.10\), and the futility threshold is \(\gamma=0.15\). Assume the response endpoint is measured at the same time in both arms and the binomial model is appropriate.

The posterior distributions are \(p_E\mid D\sim\mathrm{Beta}(1+8,1+20-8)=\mathrm{Beta}(9,13)\) and \(p_C\mid D\sim\mathrm{Beta}(1+10,1+20-10)=\mathrm{Beta}(11,11)\). The target-event probability is therefore

$$ q=\int_{0.10}^{1} f_{\mathrm{Beta}(9,13)}(u) F_{\mathrm{Beta}(11,11)}(u-0.10)\,du\ \approx\ 0.102. $$

Numerical integration gives a posterior probability of about 10.2% that the experimental response rate exceeds the control rate by at least 10 percentage points. Since \(0.102<0.15\), this arm meets the stated futility rule. The conclusion is not that the treatment cannot work; it is that, under this model and these priors, the posterior probability of reaching this particular target is below the prespecified decision threshold.

Worked Example: A Promising Arm Does Not Meet the Futility Criterion

Consider a second hypothetical interim analysis with 15 responses among 20 experimental participants and 10 among 20 control participants. Use the same independent \(\mathrm{Beta}(1,1)\) priors, target \(\delta=0.10\), threshold \(\gamma=0.15\), and outcome-model assumptions as in the first example.

The posteriors are \(p_E\mid D\sim\mathrm{Beta}(16,6)\) and \(p_C\mid D\sim\mathrm{Beta}(11,11)\). Substitution into the target-probability expression gives

$$ q=\int_{0.10}^{1} f_{\mathrm{Beta}(16,6)}(u) F_{\mathrm{Beta}(11,11)}(u-0.10)\,du\ \approx\ 0.81. $$

The posterior probability of a risk difference of at least 0.10 is about 81%, which is not below 15%. The arm therefore does not meet this futility criterion. This does not establish that the target will be reached or that the arm is superior; it says only that the rule does not remove the arm at this analysis.

Worked Example: Sparse Interim Data

In a third hypothetical trial, suppose 1 of 8 experimental participants and 2 of 8 control participants have responded. With independent \(\mathrm{Beta}(1,1)\) priors, let the target be \(\delta=0.05\) and the threshold be \(\gamma=0.25\). Assume the response counts and endpoint definitions are reliable, while recognizing that the small sample leaves substantial uncertainty.

The posteriors are \(p_E\mid D\sim\mathrm{Beta}(2,8)\) and \(p_C\mid D\sim\mathrm{Beta}(3,7)\). Evaluating the integral gives

$$ q=\int_{0.05}^{1} f_{\mathrm{Beta}(2,8)}(u) F_{\mathrm{Beta}(3,7)}(u-0.05)\,du\ \approx\ 0.22. $$

The estimated probability is about 22%, just below the 25% threshold, so the stated rule would remove the arm. The result should be read in light of the sparse interim data and the specified prior: it is a rule-based decision under the model, not a claim of certainty about the arm’s eventual effect. In practice, the protocol should specify whether this analysis is permitted at such an information level.

Derivation, Application, and Diagnostics

The posterior update and the futility decision have two distinct steps. First, the observed counts update the distribution of each response probability. Second, the joint posterior distribution is used to calculate the probability of the target event. Keeping these steps separate helps prevent a common error: comparing posterior means alone with the target.

1
Update each arm’s posterior.
For each arm, add its observed responses to the beta prior’s first parameter and its nonresponses to the second parameter.
2
Calculate the target probability.
Evaluate \(P(p_E-p_C\geq\delta\mid D)\) from the joint posterior, using the same target and endpoint as the protocol.
3
Apply the prespecified threshold.
Remove the arm under this rule if and only if \(q<\gamma\). Record the analysis data, prior, target, threshold, and decision.

The calculation has a useful basic check. If \(\delta\) increases while the posterior distributions stay fixed, the event \(p_E-p_C\geq\delta\) becomes narrower, so its posterior probability cannot increase. Similarly, if \(\gamma\) increases, more arms can satisfy \(q<\gamma\). A software result that violates either property may indicate a coding, sign, or threshold error.

Other diagnostics focus on whether the model answers the intended question. Check arm-specific denominators, response counts, outcome timing, and any missing or pending outcomes included in the interim data. Review whether the beta priors are defensible for the setting and whether a common response probability per arm is plausible. If outcomes vary substantially across sites or patient groups, or if dependence is important, an unadjusted binomial model may misrepresent uncertainty. A sensitivity analysis using prespecified alternative priors or a more suitable model can reveal whether the decision is fragile.

At multiple interim looks, repeatedly applying a threshold does not make the decision automatically calibrated to a desired frequentist operating characteristic. The protocol should specify when the rule is evaluated and the decision threshold at each permitted look. This tutorial’s posterior calculation uses observed interim data to describe uncertainty about the arm probabilities; it does not calculate the chance that future enrollment will eventually cross a stopping boundary.

How Futility Differs From Other Arm-Dropping Rules

A posterior superiority rule, as in the earlier material on Posterior Superiority Rules for Dropping Arms, evaluates posterior support for a specified ordering such as \(p_E>p_C\). The futility rule here evaluates support for a target-sized difference, \(p_E-p_C\geq\delta\). These events are not generally the same: evidence for any positive difference need not imply much evidence for reaching a larger target. Conversely, a target can be defined as zero, in which case the probability event reduces to superiority, though the decision threshold still determines the action.

The posterior probability in this rule is also distinct from a probability about future trial outcomes. It conditions on the current data and concerns the unknown arm probabilities. A calculation about future enrollment would require modeling future outcomes and the planned continuation of the trial, which is a different decision quantity.

The rule’s limitations follow from its inputs. A posterior probability depends on the likelihood, prior, target, and interim data; changing any of them can change the decision. A small \(q\) is not proof of no treatment effect, and failing to meet the futility criterion is not proof of success. The rule also does not, by itself, determine whether an arm is removed permanently, how participants already assigned to it are followed, or how randomization changes afterward. Those operational details belong in the trial protocol.

Common Mistake

Common mistake: “The posterior mean risk difference is below the target, so the arm is futile.” This reasoning substitutes a point summary for the probability of the target event. For example, a posterior mean below \(\delta\) can coexist with substantial posterior probability above \(\delta\), depending on uncertainty and the shape of the posterior. Calculate \(P(p_E-p_C\geq\delta\mid D)\) directly, then compare that probability with \(\gamma\). Do not interpret the resulting posterior probability as the probability that the null hypothesis is true.

Check Your Understanding

Answer each question, then compare your reasoning with the solutions.

  1. In a futility rule, what does \(\delta\) specify, and what does \(\gamma\) specify?
  2. With independent \(\mathrm{Beta}(1,1)\) priors, an experimental arm has 6 responses among 12 participants. What is its posterior distribution?
  3. Why does comparing the posterior mean risk difference with \(\delta\) not calculate the futility probability?
  4. If \(q=0.08\) and \(\gamma=0.10\), what decision follows from the rule? State what that decision does not establish.
  5. Holding the posterior distributions fixed, what happens to \(q\) if the target \(\delta\) is increased?

Solutions. (1) \(\delta\) is the target risk difference; \(\gamma\) is the posterior-probability threshold below which the arm is removed. (2) The posterior is \(\mathrm{Beta}(1+6,1+12-6)=\mathrm{Beta}(7,7)\). (3) The mean is one summary of the posterior; the rule requires the probability of the entire event \(p_E-p_C\geq\delta\), which depends on posterior uncertainty and shape. (4) Since \(0.08<0.10\), remove the arm under this rule. This does not prove that the treatment has no effect or that it cannot eventually reach the target. (5) \(q\) cannot increase: raising \(\delta\) makes the event \(p_E-p_C\geq\delta\) smaller or leaves it unchanged.

Clinical Biostats