Tutorials › Multi-Arm Multi-Stage (MAMS) Designs › Predictive Probability of Final Success

Time-to-Event MAMS With Intermediate Outcomes · Tutorial 323 of 499

Predictive Probability of Final Success

Derive and calculate predictive probabilities of final success, and assess how they depend on the intermediate-to-final model and assumptions about future information.

Advanced 12 min read

What You'll Learn

  • Define predictive probability of final success conditional on interim intermediate outcomes.
  • Derive the normal posterior-predictive calculation for a one-sided final test.
  • Distinguish predictive probability from conditional power at a fixed assumed effect.
  • Integrate over mixture models and uncertain future event information.
  • Check model assumptions, multiplicity boundaries, and common sources of overconfident predictions.

From Intermediate Outcomes to a Probability of Final Success

An intermediate outcome can be informative about a final time-to-event endpoint without being the endpoint itself. At an interim look, the practical question may be: given the intermediate outcomes observed so far, how likely is this arm to satisfy its prespecified final success criterion? The answer depends both on what the intermediate data imply about the final treatment effect and on how much final-endpoint information remains to be collected.

Predictive probability makes both sources of uncertainty explicit. It averages over uncertainty about the effect, as updated by the intermediate outcomes, and over the distribution of future final-endpoint data. It is therefore not simply the probability of success calculated at one chosen effect value. The result is conditional on the interim data and on the stated models and design assumptions; it is not a guarantee that the arm will succeed.

Definition: Let \(D_I\) denote the observed interim intermediate outcomes, and let \(S\) denote the event that an arm meets its prespecified final-endpoint success criterion. The predictive probability of final success is \(P(S\mid D_I)\), evaluated under an explicit model for the final effect given \(D_I\) and for future final-endpoint data given that effect. In a Bayesian analysis, the effect uncertainty is represented by its posterior distribution given \(D_I\). A model-based predictive calculation can also use another explicitly justified distribution for that uncertainty.

For illustration, write \(\theta\) for an arm’s effect on the final endpoint, with larger values defined to be more beneficial. For a time-to-event endpoint, one possible convention is \(\theta=-\log(HR)\), where \(HR\) is the hazard ratio comparing the arm with control. Under that convention, a hazard ratio below 1 corresponds to \(\theta>0\). The final success boundary may be a one-sided critical value for a test of no benefit, or a different prespecified boundary that accounts for the design’s testing procedure.

A Normal Posterior-Predictive Calculation

Consider a simple normal approximation. Suppose the interim intermediate data provide an estimate \(X\) that is linked to the final effect \(\theta\) by the model \(X\mid\theta\sim N(\theta,s_I^2)\). Suppose the analysis assigns \(\theta\) a normal prior distribution \(N(m_0,v_0)\). At the final analysis, let \(Y\) be an estimator of \(\theta\), with \(Y\mid\theta\sim N(\theta,V_F)\). Here \(V_F\) is the variance of the final estimator under the planned final information. Assume, for this calculation, that \(X\) and \(Y\) are conditionally independent given \(\theta\), and that the variances are known.

If final success requires \(Y/\sqrt{V_F}\ge c\), where \(c\) is the prespecified one-sided critical value, then success is equivalent to \(Y\ge c\sqrt{V_F}\). Normal updating gives a posterior distribution for \(\theta\mid X\); the predictive distribution for \(Y\mid X\) then averages the future sampling variation over that posterior uncertainty.

$$ v=\left(\frac{1}{v_0}+\frac{1}{s_I^2}\right)^{-1}, \qquad m=v\left(\frac{m_0}{v_0}+\frac{X}{s_I^2}\right), \qquad \theta\mid X\sim N(m,v). $$

To see the predictive variance, write \(Y=\theta+\varepsilon_F\), where \(\varepsilon_F\sim N(0,V_F)\) is independent of \(\theta\mid X\). Thus \(E(Y\mid X)=m\), and \(\operatorname{Var}(Y\mid X)=\operatorname{Var}(\theta\mid X)+\operatorname{Var}(\varepsilon_F)=v+V_F\). Consequently, the predictive probability is

$$ P(S\mid X) = 1-\Phi\left(\frac{c\sqrt{V_F}-m}{\sqrt{V_F+v}}\right), $$

where \(\Phi\) is the standard normal cumulative distribution function. The denominator contains two distinct uncertainties: \(V_F\), the future sampling variance, and \(v\), the remaining uncertainty about the effect after seeing the intermediate data. This calculation relies on the assumed link between intermediate outcomes and the final effect, as well as on the normal, variance, and independence assumptions stated above.

Worked Example: Normal Updating From an Intermediate Estimate

Suppose a hypothetical arm has final effect \(\theta=-\log(HR)\), and its interim intermediate data yield \(X=0.180\) with standard error \(s_I=0.120\). Assume \(\theta\sim N(0,0.200^2)\), and that the planned final estimator has standard error \(\sqrt{V_F}=0.100\). For this example, assume the intermediate estimate is unbiased for \(\theta\) under the stated normal model, the final estimator is unbiased conditional on \(\theta\), and the two estimates are conditionally independent given \(\theta\). Let final success require \(Y/0.100\ge1.96\).

The prior variance is \(v_0=0.0400\), and the intermediate estimate variance is \(s_I^2=0.0144\). The posterior variance and mean are

$$ v=\left(\frac{1}{0.0400}+\frac{1}{0.0144}\right)^{-1} =\left(25.000+69.4444\right)^{-1} =0.0105882, $$
$$ m=0.0105882\left(\frac{0}{0.0400}+\frac{0.180}{0.0144}\right) =0.0105882(12.500) =0.132353. $$

Success requires \(Y\ge1.96(0.100)=0.196\). The predictive standard deviation of \(Y\) is

$$ \sqrt{V_F+v} = \sqrt{0.0100+0.0105882} = \sqrt{0.0205882} \approx 0.1435. $$

Therefore,

$$ P(S\mid X) = 1-\Phi\left(\frac{0.196-0.132353}{0.1435}\right) = 1-\Phi(0.4436) \approx 0.329. $$

Under these assumptions, the arm has about a 32.9% predictive probability of meeting the stated final criterion. This is a model-conditional probability, not a claim that there is a 32.9% chance the effect is beneficial or that the final result is already determined. The conversion of the intermediate estimate into information about \(\theta\) is a substantive assumption, not a consequence of the arithmetic.

Why Predictive Probability Differs From Conditional Power

Conditional power calculates the chance of final success under a specified value of the true effect, conditional on the interim information included in its calculation. Predictive probability instead averages over uncertainty about that effect. In the normal example, using the posterior mean \(m=0.132353\) as though it were the true effect gives a conditional power of

$$ 1-\Phi\left(\frac{0.196-0.132353}{0.100}\right) = 1-\Phi(0.6365) \approx 0.262. $$

That is lower than the predictive probability of about 0.329. The predictive calculation includes the possibility that the true effect is above the posterior mean, as well as the possibility that it is below. Averaging success probabilities over this uncertainty need not equal the success probability at the mean because the probability function is nonlinear. The two quantities answer different questions; neither is universally preferable. Their usefulness depends on the decision being considered and on whether the assumptions about the effect are credible.

In multi-arm designs, this calculation is initially an arm-specific probability for a specified final boundary. A marginal probability for one arm does not by itself describe the probability that any arm succeeds, the probability of selecting the correct arm, or a familywise error rate. If final success uses an adjusted or otherwise design-specific boundary, use that boundary in the predictive calculation. Dependence between arms, for example from a shared control group, also matters when calculating joint probabilities.

Worked Example: Averaging Over a Mixture of Plausible Effects

Consider a second hypothetical analysis in which the intermediate data support two possible regions for the final effect. Suppose the model gives posterior probability \(0.80\) to \(\theta\sim N(0.100,0.040^2)\) and probability \(0.20\) to \(\theta\sim N(0.300,0.050^2)\). Assume the final estimator has standard error \(0.120\), conditionally independently of the intermediate data given \(\theta\), and final success requires \(Y\ge0.196\). These mixture weights and component distributions are assumptions of this example, not facts inferred from the numbers alone.

Within the first component, the predictive standard deviation is \(\sqrt{0.120^2+0.040^2}=\sqrt{0.0160}=0.1265\). Its success probability is \(1-\Phi((0.196-0.100)/0.1265)=1-\Phi(0.759)\approx0.224\).

Within the second component, the predictive standard deviation is \(\sqrt{0.120^2+0.050^2}=\sqrt{0.0169}=0.1300\). Its success probability is \(1-\Phi((0.196-0.300)/0.1300)=1-\Phi(-0.800)=\Phi(0.800)\approx0.788\).

Weight the component probabilities by their posterior probabilities:

$$ P(S\mid D_I) = 0.80(0.224)+0.20(0.788) \approx 0.337. $$

The predictive probability is approximately 33.7%. Replacing the mixture by one normal distribution with the same mean, or evaluating success only at the mixture’s mean effect, is not generally equivalent: the success probability must be averaged over the specified effect distribution.

Uncertain Future Information

The final information may itself be uncertain at the interim look. For a time-to-event endpoint, for example, the eventual number of observed events can vary with accrual, follow-up, and event occurrence. If future information \(I\) has a specified distribution given the interim data, predictive success should average over that distribution as well. Under the normal approximation \(Y\mid\theta,I\sim N(\theta,1/I)\), the final statistic is \(Z=Y\sqrt{I}\), and the success event is \(Z\ge c\). Conditional on \(I\) and the posterior \(\theta\mid D_I\sim N(m,v)\), \(Z\) has mean \(m\sqrt I\) and variance \(1+vI\). Thus the probability conditional on \(I\) is

$$ P(S\mid D_I,I) = 1-\Phi\left(\frac{c-m\sqrt I}{\sqrt{1+vI}}\right). $$

Worked Example: Averaging Over Future Event Information

Use the hypothetical posterior \(m=0.132353\) and \(v=0.0105882\) from the normal-updating example, but suppose final information is uncertain: \(I=100\) with probability \(0.70\), or \(I=64\) with probability \(0.30\), conditional on the interim data. Assume the normal information model above applies in either case and use \(c=1.96\).

For \(I=100\), the standardized cutoff is \((1.96-0.132353\sqrt{100})/\sqrt{1+0.0105882(100)} =(1.96-1.32353)/\sqrt{2.05882} =0.4436\). The conditional success probability is \(1-\Phi(0.4436)\approx0.329\).

For \(I=64\), the standardized cutoff is \((1.96-0.132353\sqrt{64})/\sqrt{1+0.0105882(64)} =(1.96-1.05882)/\sqrt{1.67765} \approx0.696\). The conditional success probability is \(1-\Phi(0.696)\approx0.243\).

Averaging over future information gives \(0.70(0.329)+0.30(0.243)\approx0.303\). The predictive probability is about 30.3%, rather than the 32.9% obtained under fixed information \(I=100\). This difference follows from the specified information distribution; a different accrual or event model could produce a different predictive probability.

Assumptions, Diagnostics, and Alternatives

The calculation is only as credible as its ingredients. Check that the final success criterion is stated on the correct scale and direction; that the intermediate-to-final relationship is supported by an appropriate model; and that the model reflects the trial’s accrual, follow-up, and endpoint maturity. For the normal calculations, also check the normal approximation, variance assumptions, and any claimed conditional independence. In survival settings, censoring patterns, competing events, or changing risk over follow-up can make a simple fixed-information approximation inadequate.

Useful diagnostics include comparing predicted and observed final outcomes in relevant historical or simulated settings, checking calibration across ranges of intermediate results, and examining sensitivity to plausible alternative effect distributions and future-information assumptions. A prediction that changes sharply under small reasonable changes to the intermediate-to-final link or prior distribution is fragile. Predictions made when few outcomes are mature may also be strongly driven by assumptions rather than by observed evidence.

The main alternatives answer related but distinct questions. Conditional power evaluates a fixed assumed effect; predictive probability averages over specified uncertainty about effects and, if needed, future information. A plug-in forecast based only on a point estimate is simpler but omits uncertainty in that estimate. Simulation can handle more complex time-to-event mechanisms, but it still requires a defensible model and must implement the actual final criterion. None of these methods removes the need to prespecify the final testing rule or to account for the design’s multiplicity.

Predictive probability is not automatically a frequentist operating characteristic. It depends on the distribution used for effect uncertainty and on the future-data model. It should not be interpreted as the probability that the null hypothesis is true, nor as a replacement for error-rate evaluation of the complete design. It is a conditional forecast useful for reasoning about an arm given interim information, subject to explicit assumptions.

Common Mistake

  • Wrong reasoning: “The posterior mean effect is \(0.132\), so the predictive probability is the conditional power computed at \(0.132\).” This treats the uncertain effect as known and omits its posterior variance.
  • Correction: Integrate the final success probability over the effect distribution. In the normal model, this adds \(v\) to \(V_F\) in the predictive variance; with mixtures or uncertain information, average over the relevant components or information values as well.
  • Another warning: A high predictive probability is not proof that the arm is effective. It is conditional on the stated model and does not replace the prespecified final test.
Key takeaway: Predictive probability of final success is obtained by combining what interim intermediate outcomes imply about the final effect with the distribution of future final data. State the success boundary, effect model, future-information assumptions, and remaining uncertainty before interpreting the resulting probability.

Check Your Understanding

Work through each question, then compare with the solutions below.

  1. In a normal posterior-predictive model, what two sources of variance appear in the predictive variance of the final estimator?
  2. Suppose a one-sided final criterion is \(Y/\sqrt{V_F}\ge c\). Write the equivalent criterion for \(Y\).
  3. Why can predictive probability differ from conditional power calculated at the posterior mean effect?
  4. In the uncertain-information example, why must the two information-specific success probabilities be weighted and averaged?

Solutions. (1) The predictive variance is \(V_F+v\): \(V_F\) is future sampling variance, and \(v\) is posterior variance in the effect. (2) Multiplying by the positive standard error gives \(Y\ge c\sqrt{V_F}\). (3) Conditional power fixes the effect at one value; predictive probability averages over its uncertainty, and the success-probability function is nonlinear. (4) The final information is random under the stated model, so the overall probability is the sum of conditional probabilities weighted by the probabilities of the corresponding information values.

Clinical Biostats