Tutorials › Bayesian Clinical Trials › Worked Continuous Conflict Analysis

Robustness to Prior-Data Conflict · Tutorial 448 of 696

Worked Continuous Conflict Analysis

Work through how historical-current conflict changes posterior estimates and uncertainty under robust mixture and commensurate models, including a careful comparison of borrowing.

Advanced 12 min read

What You'll Learn

  • Set up a normal likelihood for a current continuous-outcome estimate and state the assumptions needed to borrow historical information.
  • Calculate robust-mixture component probabilities, posterior means, and total posterior variance.
  • Derive the normal-normal posterior used to evaluate commensurate precision choices.
  • Compare borrowing with the current-only estimate using an explicitly defined descriptive measure.
  • Recognize diagnostics and limitations that can make a borrowing analysis misleading.
  • Distinguish sensitivity analysis from evidence that one model is automatically correct.

Why Continuous-Outcome Conflict Needs Careful Analysis

A historical trial may estimate a treatment effect close to zero, while a current trial estimates a clinically important benefit. For a continuous outcome, the size of that discrepancy must be judged relative to the standard errors, the historical information's relevance, and the borrowing model. A difference of one outcome unit can be substantial when estimates are precise and modest when they are noisy.

This tutorial works with hypothetical trial summaries to show how conflict affects posterior effects, uncertainty, and effective borrowing. The calculations assume that both trials target the same treatment-effect estimand, measured in the same units. If the populations, endpoint definitions, follow-up, or analysis methods differ in ways that change the estimand, a model that borrows across the estimates may not repair that mismatch.

Definition: Let \(\hat\theta_C\) be the current-trial estimate of a continuous treatment effect \(\theta_C\), with known standard error \(s_C\). The normal summary-data likelihood is \(\hat\theta_C\mid\theta_C\sim N(\theta_C,s_C^2)\). Let \(D_H\) denote the historical data, and let \(\theta_H\) denote the corresponding historical effect. Conflict means that the historical and current information are discrepant relative to their uncertainty and the assumed relationship between \(\theta_H\) and \(\theta_C\); it is not, by itself, proof that either estimate is wrong.

The normal likelihood is a suitable summary model when the estimate is approximately normal and its standard error is reliable. We treat the standard errors as known for the calculations below, as is common when working from sufficiently large-sample trial summaries. We also assume that the historical and current data are conditionally independent given their respective effects. These assumptions should be reconsidered when standard errors are estimated imprecisely, outcomes are highly skewed, or the estimates share participants or controls.

A Robust Mixture Prior for the Current Effect

A robust mixture prior represents two possible degrees of agreement with historical information. Its informative component transfers the historical estimate and uncertainty to the current effect. Its protective component is deliberately broader, allowing the current estimate more room to move away from history. The mixture does not guarantee that borrowing will be small: the prior component weights and the predictive fit of each component both matter.

Definition: A two-component mixture prior for \(\theta_C\) is \(\pi(\theta_C\mid D_H)=\omega\,\pi_I(\theta_C\mid D_H)+(1-\omega)\,\pi_P(\theta_C)\), where \(0\leq\omega\leq1\), \(\pi_I\) is the informative component, and \(\pi_P\) is the protective component. After observing the current estimate, the posterior probability of component \(j\) is proportional to its prior probability times its predictive density for \(\hat\theta_C\). The posterior for \(\theta_C\) is the corresponding probability-weighted mixture of the component posteriors.

For a normal component \(j\), suppose its prior is \(N(a_j,v_j)\). Combining it with the current normal likelihood gives a normal posterior. The result follows by multiplying the two densities and collecting the terms involving \(\theta_C\):

$$ -\frac{1}{2}\left\{\frac{(\hat\theta_C-\theta_C)^2}{s_C^2} +\frac{(\theta_C-a_j)^2}{v_j}\right\} = -\frac{1}{2V_j}(\theta_C-m_j)^2+\text{constant}, $$

where \(V_j=(s_C^{-2}+v_j^{-1})^{-1}\) and \(m_j=V_j(s_C^{-2}\hat\theta_C+v_j^{-1}a_j)\). The predictive distribution of \(\hat\theta_C\) under that component is \(N(a_j,s_C^2+v_j)\). Thus the posterior component weight \(w_j\) is found by normalizing the component's prior weight times this predictive density.

Worked Example: A Robust Mixture Under Continuous-Outcome Conflict

Suppose the historical summary yields a historical-effect distribution \(N(0.20,0.04)\): its mean is 0.20 outcome units and its variance is 0.04 outcome-units squared. The current estimate is \(\hat\theta_C=1.20\) units with \(s_C=0.50\) units. Use an informative component \(N(0.20,0.04)\) with prior weight 0.80 and a protective component \(N(0,4)\) with weight 0.20. For this illustration, zero is a reference value and the protective variance is intentionally broad; neither choice is universal.

The current likelihood variance is \(s_C^2=0.25\). The predictive variances for the informative and protective components are \(0.04+0.25=0.29\) and \(4+0.25=4.25\), respectively. Their predictive densities at 1.20 have ratio

$$ \frac{f_I(1.20)}{f_P(1.20)} = \sqrt{\frac{4.25}{0.29}}\, \exp\left\{-\frac12\left[ \frac{(1.20-0.20)^2}{0.29} -\frac{(1.20-0)^2}{4.25}\right]\right\} \approx 0.809. $$

The observed estimate is less predictive under the informative component than under the broad protective component. But the informative component began with four times the prior weight. Updating the prior odds gives

$$ w_I=\frac{0.80 f_I(1.20)} {0.80 f_I(1.20)+0.20 f_P(1.20)} \approx 0.7639,\qquad w_P=1-w_I\approx0.2361. $$

The informative-component posterior has variance \(V_I=(1/0.25+1/0.04)^{-1}=0.03448\) units squared and mean \(m_I=V_I(1.20/0.25+0.20/0.04)=0.3379\) units. The protective-component posterior has variance \(V_P=(1/0.25+1/4)^{-1}=0.2353\) units squared and mean \(m_P=V_P(1.20/0.25+0/4)=1.129\) units.

The mixture posterior mean uses the unrounded component probabilities and means:

$$ E(\theta_C\mid D_H,\hat\theta_C) =w_I m_I+w_Pm_P \approx 0.5248\text{ units}. $$

For a mixture, total variance includes within-component variance and variation between component means. With the same unrounded values,

$$ \operatorname{Var}(\theta_C\mid D_H,\hat\theta_C) =\sum_{j\in\{I,P\}}w_j\{V_j+(m_j-0.524831)^2\} \approx0.1949\text{ units}^2. $$

The posterior standard deviation is therefore approximately \(0.4415\) units. The between-component term matters: reporting only an average of \(V_I\) and \(V_P\) would omit uncertainty about which component best represents the current effect.

Commensurate Models and Precision

A commensurate model relates the current effect to a historical effect through a precision parameter. Larger precision means stronger concentration of \(\theta_C\) around \(\theta_H\), hence a stronger potential for borrowing; smaller precision permits greater separation. Unlike a fixed informative prior alone, the model makes the degree of agreement an explicit part of the relationship. Here we evaluate fixed precision choices to show the consequences clearly; a full analysis may instead assign a prior to the precision and integrate over it.

Model: Suppose the historical posterior is \(\theta_H\mid D_H\sim N(m_H,v_H)\), and, conditional on \(\theta_H\) and commensurate precision \(\tau\), \(\theta_C\mid\theta_H,\tau\sim N(\theta_H,1/\tau)\). The current estimate satisfies \(\hat\theta_C\mid\theta_C\sim N(\theta_C,s_C^2)\). Assume the current and historical data are conditionally independent. Integrating over the historical posterior gives \(\theta_C\mid D_H,\tau\sim N(m_H,v_H+1/\tau)\).

Worked Example: Comparing Two Commensurate Precision Choices

Use the same hypothetical summaries: \(m_H=0.20\) units, \(v_H=0.04\) units squared, \(\hat\theta_C=1.20\) units, and \(s_C^2=0.25\) units squared. Compare \(\tau=4\) with \(\tau=0.25\), where precision is measured in inverse outcome-units squared. The resulting prior variances for \(\theta_C\) are \(0.04+1/4=0.29\) and \(0.04+1/0.25=4.04\).

For a normal prior \(N(m_H,v)\), the posterior variance and mean are \(V=(1/s_C^2+1/v)^{-1}\) and \(m=V(\hat\theta_C/s_C^2+m_H/v)\). With \(\tau=4\), \(v=0.29\), so

$$ V=\left(\frac1{0.25}+\frac1{0.29}\right)^{-1} =\frac{0.25(0.29)}{0.54}=0.1343\text{ units}^2, \qquad m=V\left(\frac{1.20}{0.25}+\frac{0.20}{0.29}\right) =0.7370\text{ units}. $$

The posterior standard deviation is \(\sqrt{0.1343}=0.3664\) units. With \(\tau=0.25\), \(v=4.04\), giving

$$ V=\left(\frac1{0.25}+\frac1{4.04}\right)^{-1} =\frac{0.25(4.04)}{4.29}=0.2354\text{ units}^2, \qquad m=V\left(\frac{1.20}{0.25}+\frac{0.20}{4.04}\right) =1.141\text{ units}. $$

The posterior standard deviation is \(0.4852\) units. The higher commensurate precision pulls the estimate further toward historical information and reduces uncertainty in this example. These are conditional results for two specified precision values, not a data-driven conclusion that either value is correct. If \(\tau\) is uncertain, its posterior distribution and the resulting mixture over \(\theta_C\) must also be included in the reported uncertainty.

Comparing Posterior Effects and Effective Borrowing

A useful reference for describing borrowing is the current-only analysis. Under a flat prior for \(\theta_C\), the current likelihood gives posterior mean 1.20 units and variance 0.25 units squared. Relative to this reference, the robust-mixture mean is 0.5248 units, and the fixed-\(\tau\) commensurate means are 0.7370 and 1.141 units. Each analysis gives a different balance of historical influence and current evidence.

One descriptive summary is the fraction of the distance from the current estimate to the historical mean that the posterior mean has moved:

$$ B=\frac{\hat\theta_C-E(\theta_C\mid\text{data})} {\hat\theta_C-m_H}. $$

For these data, the robust-mixture value is approximately \((1.20-0.524831)/(1.20-0.20)=0.6752\). The commensurate values are approximately \(0.4630\) for \(\tau=4\) and \(0.05897\) for \(\tau=0.25\). This is a descriptive mean-shift measure, not a universal effective-sample-size calculation. It can be misleading if the posterior lies beyond the historical mean, if the historical mean is not the relevant borrowing target, or if the model has several sources of historical information. Always pair it with the posterior uncertainty and model specification.

A concise comparison should report the current-only reference, the posterior mean and uncertainty under each model, and the assumptions that create the borrowing. In the mixture example, the posterior component probability is also informative: it indicates how the data updated the specified component probabilities, not the probability that a scientific hypothesis or model is true.

Diagnostics, Limitations, and a Common Pitfall

Before interpreting a calculation, verify that both estimates refer to the same effect, outcome scale, direction of benefit, and analysis population. Check the reported standard errors, the independence assumption, and whether the historical posterior already includes uncertainty about the historical effect. A commensurate model that treats \(m_H\) as fixed when \(v_H\) is material can overstate historical precision. A protective component that is broad but poorly centered may also behave differently from a genuinely data-driven no-borrowing analysis.

Useful diagnostics include plotting the component prior and predictive distributions against the current estimate, inspecting posterior component probabilities or the posterior for \(\tau\), and repeating the analysis across plausible prior weights, protective variances, and precision specifications. If the posterior effect changes sharply across reasonable choices, that sensitivity is a substantive result to report. Numerical integration or MCMC also requires computational checks appropriate to the method, such as verifying that numerical weights normalize and that chains explore the relevant regions when sampling is used.

The main alternative is no borrowing: analyze the current trial without historical information. It is transparent and avoids cross-trial exchangeability assumptions, but may be less precise. A fixed informative prior can be simpler and more precise when historical relevance is well justified, but it offers less protection against conflict. Robust mixtures and commensurate models allow disagreement to influence borrowing, at the cost of additional choices and sensitivity. None removes the need to assess clinical and design comparability, and none turns a conflict diagnostic into a definitive test of whether the trials are exchangeable.

  • Pitfall: “The informative component had the largest prior weight, so the current posterior must be close to the historical estimate.” This ignores the current data's predictive fit under each component and the posterior component-specific estimates.
  • Correction: Update component probabilities using their predictive densities, then calculate the posterior mean and total variance from the component posteriors. In the worked mixture, the informative component retained posterior weight about 0.7639, but the protective component still contributed both a distinct mean and appreciable between-component uncertainty.
  • Limitation: A numerical borrowing summary depends on the reference and model. It cannot establish that historical and current populations are clinically exchangeable or that a borrowed estimate has causal validity.
Key takeaway: Continuous-outcome conflict analysis is a comparison of explicit assumptions and posterior consequences. Report the estimated effect with uncertainty, show how borrowing was constructed, and describe sensitivity to plausible alternatives rather than presenting one borrowing measure as definitive.

Check Your Understanding

Try the questions first, then compare your reasoning with the solutions.

  1. In a two-component robust mixture, why is the posterior component probability not necessarily equal to its prior probability?
  2. Why does the total variance of a mixture include a term for the distance between component means?
  3. In the commensurate example, what happens to the prior variance of \(\theta_C\) as \(\tau\) decreases?
  4. Using the worked mixture's posterior mean and \(m_H=0.20\), calculate the descriptive borrowing fraction \(B\).
  5. Name two assumptions or diagnostics that should be checked before treating a continuous historical estimate as borrowable information.

Solutions. (1) The current estimate updates the component probabilities according to each component's predictive density; a component that predicts the observed estimate relatively well gains posterior weight. (2) The component identity is uncertain, so posterior draws can differ because they come from distributions with different means, even after accounting for each component's own variance. (3) Since the commensurate variance includes \(1/\tau\), decreasing \(\tau\) increases that variance and weakens concentration around the historical effect. (4) \(B=(1.20-0.524831)/(1.20-0.20)\approx0.6752\). It is a descriptive fraction, not a universal effective sample size. (5) Suitable checks include common estimand and outcome scale, reliable standard errors, conditional independence, historical-effect uncertainty, predictive fit, and sensitivity to mixture or precision choices.

Clinical Biostats