Why Use a Bayesian Rule for Enrichment?
At an interim analysis, a trial team may need to decide whether to keep enrolling a biomarker-defined population, restrict enrollment to another population, or stop a population’s enrollment. A Bayesian decision rule expresses this choice in terms of the evidence about the population-specific treatment effect after combining interim data with a prior distribution. It provides a direct probability statement for a prespecified decision, but its meaning depends on the model, prior, and decision threshold.
Let \(\theta_j\) denote the treatment effect for population \(j\), with positive values defined as favoring treatment. The effect must use the population-specific estimand and scale specified for the trial. For example, if the outcome is a change score where a larger value is better, positive effects may represent improvement; for an outcome where smaller values are better, the effect should be coded or interpreted accordingly.
The threshold \(q\) is not the significance level \(\alpha\). It is an operational decision criterion: a larger \(q\) requires stronger posterior evidence of benefit before retaining a population. Its choice should reflect the consequences of continuing enrollment in a population that is not beneficial, and of excluding one that is beneficial. The interim rule also needs to be specified before examining the data on which it will act.
A Normal Model for an Interim Effect Estimate
A useful working model assumes that the interim estimate \(\hat{\theta}_j\) has an approximately normal sampling distribution with known standard error \(SE_j\), conditional on the true effect:
Here \(m_{0j}\) is the prior mean and \(\tau_j\) is the prior standard deviation. The likelihood approximation is most appropriate when the estimator is reasonably well approximated by a normal distribution and its standard error is reliable. In practice, uncertainty in nuisance parameters, small samples, or complex endpoints may require a more suitable likelihood rather than treating \(SE_j\) as fixed.
For the calculation below, assume the prior is proper, the estimate and standard error correspond to the prespecified population-specific estimand, and the interim data are used once in the likelihood. If multiple population estimates share participants or controls, their estimates may be correlated; separate calculations that assume independence would then be inappropriate. A joint model should represent that dependence.
Posterior Calculation and the Decision Boundary
Under the normal likelihood and prior, the posterior is also normal. Its variance is determined by adding prior and likelihood precisions, and its mean is the precision-weighted average of the prior mean and the interim estimate:
The posterior probability of a beneficial effect is therefore \(P(\theta_j>0\mid D)=\Phi(m_j/\sqrt{v_j})\), where \(\Phi\) is the standard normal cumulative distribution function. This follows by standardizing the posterior: if \(\theta_j\mid D\sim N(m_j,v_j)\), then \((\theta_j-m_j)/\sqrt{v_j}\) is standard normal, so \(P(\theta_j>0\mid D)=P((\theta_j-m_j)/\sqrt{v_j}>-m_j/\sqrt{v_j})=\Phi(m_j/\sqrt{v_j})\).
The same rule can be written as a cutoff for the interim estimate. Define \(w=\tau_j^2/(\tau_j^2+SE_j^2)\). Then \(m_j=w\hat{\theta}_j+(1-w)m_{0j}\), and \(\sqrt{v_j}=\sqrt{\tau_j^2SE_j^2/(\tau_j^2+SE_j^2)}\). Since \(w>0\), the rule is equivalent to:
This gives a concrete way to see how changing the prior or \(q\) changes the data threshold for retention. It does not mean that an interim estimate alone determines the decision: the cutoff also depends on the prior and standard error.
Worked Decisions
Worked Example: Retaining a Population With Strong Posterior Support
Suppose a hypothetical biomarker-positive population has an interim treatment-effect estimate of \(0.30\) units with \(SE=0.15\) units. Use a prior \(\theta\sim N(0,0.25^2)\) and retain the population if \(P(\theta>0\mid D)\geq0.90\). Assume the normal likelihood is adequate for this estimate and standard error, and that the population-specific effect is defined so positive values favor treatment.
The prior precision is \(1/0.25^2=16\), and the likelihood precision is \(1/0.15^2=44.444\). Therefore:
The posterior standard deviation is \(\sqrt{0.01654}=0.1286\) units. Hence \(P(\theta>0\mid D)=\Phi(0.2206/0.1286)=\Phi(1.715)\approx0.9568\). Because \(0.9568\geq0.90\), the prespecified rule retains this population. This is a posterior probability under the stated model and prior, not a probability that the null hypothesis is true.
Worked Example: A Prior Choice Changes the Selection Decision
In a second hypothetical population, suppose the interim estimate is \(0.12\) units with \(SE=0.20\) units, and the decision cutoff is \(q=0.70\). First use a prior \(N(0,0.15^2)\). Assume the normal sampling model is adequate and the estimate targets the same defined effect scale.
The prior precision is \(1/0.15^2=44.444\), and the likelihood precision is \(1/0.20^2=25\). The posterior variance is \(1/(44.444+25)=0.01440\), its standard deviation is \(0.1200\), and its mean is \(25(0.12)/69.444=0.04320\). Thus \(P(\theta>0\mid D)=\Phi(0.04320/0.1200)=\Phi(0.360)=0.6406\), below \(0.70\); the rule excludes the population.
Now use a less concentrated prior \(N(0,0.50^2)\), leaving the data and cutoff unchanged. The prior precision is \(4\), so the posterior variance is \(1/(4+25)=0.03448\), the posterior standard deviation is \(0.1857\), and the posterior mean is \(25(0.12)/29=0.1034\). The probability is \(\Phi(0.1034/0.1857)=\Phi(0.557)=0.7112\). This just exceeds \(0.70\), so the decision changes to retention. The data did not change; the prior assumption did.
Worked Example: Selecting Between Two Biomarker Populations
Consider a hypothetical interim analysis with two distinct biomarker populations. For the positive population, use the estimate \(0.30\) units and standard error \(0.15\) units from the first example. For the negative population, suppose the estimate is \(0.08\) units with \(SE=0.18\) units. Use independent \(N(0,0.25^2)\) priors and retain each population if its posterior probability of benefit is at least \(0.90\). For simplicity, assume these estimates are independent; if the trial’s sampling design induces dependence, this calculation would need a joint model.
The positive population has posterior probability \(0.9568\), so it meets the cutoff. For the negative population, prior precision is \(16\), likelihood precision is \(1/0.18^2=30.864\), and:
Its posterior standard deviation is \(\sqrt{0.02134}=0.1461\) units, giving \(P(\theta>0\mid D)=\Phi(0.05268/0.1461)=\Phi(0.361)\approx0.641\). Under this rule, retain the positive population and exclude the negative population. The conclusion is limited to the stated model, prior, data, and cutoff; it does not establish that the treatment effect differs between the two populations.
Choosing and Checking the Rule
A probability cutoff should reflect the decision, not simply be selected because it is a familiar number. One way to make the trade-off explicit is to assign utilities to continuing or stopping. Suppose continuing has utility \(U_B\) if the effect is beneficial and \(U_N\) if it is not, while stopping has utility \(U_S\). If \(U_B>U_N\), continuing has greater expected utility than stopping when \(pU_B+(1-p)U_N>U_S\), where \(p=P(\theta>0\mid D)\). When \(U_N<U_S<U_B\), this is equivalent to \(p>(U_S-U_N)/(U_B-U_N)\). This simple calculation illustrates how costs and benefits can motivate a threshold; real trial decisions may involve more than two possible states and actions.
Prior choice deserves direct scrutiny because it can alter the posterior probability and the selected population, as the second example shows. Examine whether the prior is scientifically defensible for the specific estimand, whether its spread is appropriate, and how the decision changes under plausible alternative priors. A prior centered at zero is not automatically neutral: its concentration determines how strongly it pulls estimates toward zero. Borrowing information from other populations or earlier evidence also requires justification that the information is relevant to the population and outcome at issue.
A posterior threshold does not, by itself, guarantee a chosen frequentist false-positive rate or other long-run operating characteristic. The rule should be evaluated over plausible trial scenarios, including no-effect settings, beneficial effects, differing subgroup effects, and realistic accrual and variance patterns. Check for unstable estimates, unexpected standard errors, missing biomarker data, and discrepancies between the assumed likelihood and observed interim data. If a model or data-quality problem is material, the numerical posterior probability may not be a sound basis for an adaptation.
Common Pitfalls and Alternatives
- Equating \(q\) with \(\alpha\). A posterior cutoff is not a frequentist significance level. It is a decision rule conditional on the model and prior; long-run error behavior must be assessed separately.
- Ignoring prior sensitivity. A posterior probability can look precise while depending strongly on a prior. Report or examine the decision under plausible prior alternatives.
- Using separate subgroup calculations when estimates are dependent. Shared controls or nested populations can induce covariance. Use a joint likelihood or model that represents it.
- Updating twice with the same interim data. When incorporating later data, use a coherent likelihood that includes each observation once; do not treat interim results as an independent new data set.
- Interpreting selection as proof of a subgroup difference. Retaining one population and excluding another under separate rules does not establish that their treatment effects differ. The decision is an adaptation, not a direct comparison of effects.
A conditional-power rule, as used in frequentist design calculations, evaluates the chance of meeting a specified future statistical criterion under assumptions about the effect and remaining data. A posterior-probability rule instead evaluates the current posterior probability of a state such as \(\theta_j>0\). They answer different questions, and neither removes the need to prespecify the adaptation and assess the trial’s operating characteristics.
Check Your Understanding
Answer each question, then check the brief solutions below.
- If \(P(\theta_j>0\mid D)=0.82\) and the prespecified cutoff is \(q=0.85\), what does the rule prescribe?
- Why does changing a prior while keeping the interim estimate fixed sometimes change the enrichment decision?
- For a normal posterior \(N(m,v)\), write the posterior probability of a beneficial effect and the equivalent mean-based decision condition.
- Does excluding one biomarker population and retaining another prove that their treatment effects differ? Why or why not?
Solutions. (1) Exclude population \(j\), because \(0.82<0.85\). (2) The prior contributes precision and a prior mean to the posterior, changing both its center and uncertainty, and therefore changing \(P(\theta_j>0\mid D)\). (3) The probability is \(\Phi(m/\sqrt{v})\); retention requires \(m\geq\Phi^{-1}(q)\sqrt{v}\). (4) No. Separate decisions compare each population’s posterior probability with a threshold; they are not a test or direct posterior comparison of the two treatment effects.