Why Covariate Imbalance Can Look Like Conflict
Historical and current trials may enroll patients with different covariate distributions. If treatment benefit varies across covariate values, each trial’s average treatment effect depends partly on which patients it enrolled. The resulting marginal effects can differ even when the treatment effects within comparable covariate groups are the same. Conversely, if the treatment effect is constant across covariates, a prognostic covariate imbalance alone does not change a within-trial treatment contrast, although it can change the outcome levels in both arms.
The distinction matters when historical data inform a current treatment effect. A commensurate model that compares unadjusted marginal effects may interpret a difference in patient composition as evidence that the treatment effects themselves conflict. Covariate adjustment can clarify the comparison, but it does not automatically make the studies comparable: it depends on measuring relevant covariates, having adequate overlap, and specifying an appropriate effect model.
A Model for Composition and Treatment-Effect Modification
For a continuous outcome, consider the following mean model in each trial:
Here \(\alpha_s\) allows the baseline outcome level to differ by trial, \(\beta\) describes the covariate’s prognostic association, \(\theta_s\) is the treatment effect at \(X=0\), and \(\gamma\) describes how the treatment effect changes with \(X\). Thus the conditional treatment effect is \(\theta_s+\gamma X\). For binary \(X\), let \(p_s=P(X=1\mid s)\). Averaging over the covariate distribution in trial \(s\) gives:
The trial-specific intercept and the prognostic covariate term cancel in this within-trial contrast. If the conditional effect parameters are shared across trials, so that \(\theta_C=\theta_H\) and \(\gamma\) is common, then:
This identity isolates the source of apparent conflict: when \(\gamma\ne0\), different covariate prevalences produce different marginal effects even without a difference in the conditional effect at either covariate value. If \(\gamma=0\), the marginal treatment effects agree under this model, regardless of \(p_C\) and \(p_H\).
The model makes substantive assumptions. Treatment is randomized within each trial, outcomes and covariates use compatible definitions and scales, and the specified conditional mean form describes the relevant treatment-effect modification. The common-\(\gamma\) assumption says that the way treatment effect changes with \(X\) is the same in both trials. If that assumption is implausible, a shared treatment effect at \(X=0\) does not by itself imply agreement at other values of \(X\). For more general covariates, standardization averages the fitted conditional effect over the chosen target distribution rather than over \(p_s\).
Worked Example: Different Marginal Effects, Matching Conditional Effects
Suppose a hypothetical outcome is change in systolic blood pressure in mmHg, with negative values favouring treatment. Let \(X=1\) indicate a high-risk baseline subgroup. In both trials, the conditional treatment effect is \(-2\) mmHg for \(X=0\) and \(-8\) mmHg for \(X=1\). Thus \(\theta_H=\theta_C=-2\) and \(\gamma=-6\) mmHg. Suppose \(p_H=0.20\) and \(p_C=0.80\).
The historical marginal effect is:
The current marginal effect is:
The marginal difference is \(-6.8-(-3.2)=-3.6\) mmHg, but the conditional treatment effects agree exactly in both subgroups. For a target population with \(P(X=1)=0.50\), the standardized effect in either trial is \(-2+(-6)(0.50)=-5.0\) mmHg. In this hypothetical example, the apparent marginal discrepancy is entirely attributable to the different subgroup proportions under the stated model.
Adjusting and Standardizing Before Borrowing
A direct way to compare studies is to estimate treatment effects conditional on the same covariates, then average those effects over one common target distribution. For a binary \(X\), define the target prevalence as \(p_*\). Under the model above, the standardized effect for trial \(s\) is:
Choosing the target distribution is part of the estimand, not a cosmetic calculation. A current-trial target answers what the treatment effect would be for the current trial’s covariate mix. A prespecified clinical target may answer a broader question. Whichever target is used, both trials must be standardized to the same distribution for their effects to be directly comparable.
A covariate-adjusted commensurate model places the borrowing relationship on corresponding adjusted treatment-effect parameters rather than on raw marginal contrasts. For example, with a fixed commensurate precision \(\kappa>0\), a simple specification is:
This prior expresses a preference for the current and historical effects at the same reference covariate value to be close, with the strength of that preference controlled by \(\kappa\). Its interpretation depends on the chosen reference value, covariate coding, and outcome units. If treatment-effect modification is present, borrowing on \(\theta_s\) alone does not determine agreement in \(\gamma\); that parameter must be modeled and assessed as well. A commensurate prior on adjusted effects is not a substitute for checking whether the adjustment model is credible.
Worked Example: Standardizing to a Common Patient Mix
In a hypothetical pair of trials, the estimated conditional mean differences are \(-3\) mmHg for \(X=0\) and \(-7\) mmHg for \(X=1\) in each trial. The historical trial has \(p_H=0.25\), and the current trial has \(p_C=0.75\). Calculate both marginal effects and standardize to a target population with \(p_*=0.40\).
The implied conditional model has \(\theta=-3\) mmHg and \(\gamma=-7-(-3)=-4\) mmHg. The historical marginal effect is:
The current marginal effect is:
These marginal effects differ by \(-2.0\) mmHg, despite identical conditional effects. Standardizing both trials to \(p_*=0.40\) gives:
The calculation supports attributing the marginal difference to composition only if the conditional effects and covariate definitions are comparable and the fitted model is adequate. It does not establish that the studies are alike in every clinically relevant respect.
What Covariate Adjustment Can and Cannot Resolve
Adjustment separates composition from conditional effect differences only to the extent that the model captures relevant effect modification. A useful analysis therefore examines both the adjusted effect and the assumptions behind it. Randomization supports an unbiased within-trial comparison, but it does not guarantee that the historical and current trials have the same conditional treatment effect, nor does it repair differences in outcome measurement, adherence, follow-up, or co-interventions.
A concise justification for the standardization step follows directly from averaging over the target distribution. For binary \(X\), the conditional effects are \(\theta_s\) when \(X=0\) and \(\theta_s+\gamma\) when \(X=1\). The target population assigns probabilities \(1-p_*\) and \(p_*\) to those strata, so:
Because \(p_*\) is identical for both trials, a difference in their standardized effects under this model is:
Thus standardization removes the contribution of differing covariate prevalence, while retaining a difference in the reference conditional effect. This result relies on the common-\(\gamma\) model; if treatment-effect modification differs by trial, the standardized contrast also depends on that difference.
Worked Example: A Conditional Commensurate Update
Suppose a hypothetical covariate-adjusted historical treatment-effect estimate is \(-5.0\) mmHg. For illustration, treat this estimate as fixed when centring the prior. Let the adjusted current-trial likelihood be approximately normal, with estimate \(-4.0\) mmHg and standard error \(2.0\) mmHg. Use the conditional commensurate prior \(\theta_C\sim N(-5.0,4)\), corresponding to \(\kappa=1/4=0.25\) mmHg\(^{-2}\).
The current likelihood has variance \(SE^2=2^2=4\) mmHg\(^2\), so its precision is \(1/4=0.25\) mmHg\(^{-2}\). The prior precision is also \(1/4=0.25\) mmHg\(^{-2}\). Combining two normal distributions with known variances gives posterior precision equal to the sum of precisions and posterior mean equal to the precision-weighted mean:
The posterior standard deviation is \(\sqrt{2}=1.414\) mmHg, and a normal 95% credible interval is approximately \(-4.5\pm1.96(1.414)\), or \((-7.27,-1.73)\) mmHg. Conditional on this model and the plug-in historical estimate, the posterior probability that \(\theta_C\) lies in this interval is 0.95. This illustrative update borrows on an adjusted effect; it does not propagate uncertainty in the historical estimate, which a full joint model would need to do.
Diagnostics, Alternatives, and Limitations
Before interpreting adjusted commensurability, inspect the covariate distributions and the model’s fit. In particular, check whether the covariates are measured in the same way, whether treatment-effect modification is clinically plausible, and whether the trials contain patients in the covariate regions used for standardization. If one trial has virtually no patients in a target stratum, its standardized effect there is driven mainly by extrapolation rather than direct information. Large changes in the adjusted effect under reasonable model specifications are a warning that the apparent resolution of conflict is fragile.
The main alternatives answer related but different questions. Comparing raw trial-level marginal effects is simple, but confounds effect differences with covariate composition when effects vary across covariates. Separate analyses avoid borrowing, but do not use historical information to estimate the current effect. Stratified comparisons can make effect modification visible when strata are sufficiently populated; regression adjustment can use covariates more efficiently but relies on the chosen model. Standardization makes the target population explicit. None of these approaches resolves unmeasured differences or guarantees transportability.
A Bayesian commensurate prior also adds assumptions beyond covariate adjustment: it specifies how strongly adjusted historical and current parameters should be related. The posterior may appear more precise because information is borrowed, but that precision is conditional on the prior and the likelihood model. A substantial residual difference between adjusted effects, poor overlap, or sensitivity to the adjustment specification should not be dismissed merely because the model produces a pooled estimate.
- Wrong reasoning: “The historical and current marginal effects differ, so their treatment effects must conflict.”
- Correction: First check whether treatment-effect modification and differing covariate distributions can explain the marginal difference. Compare conditional effects or effects standardized to a common target before judging commensurability.
- Another common error: “The baseline covariate distribution differs, so the within-trial randomized treatment contrast must be biased.”
- Correction: Under a constant additive treatment effect, a prognostic covariate imbalance changes outcome levels but cancels in the within-trial contrast. It changes the marginal treatment effect when treatment effects vary with the covariate, or can mislead when investigators compare outcome levels rather than within-trial contrasts.
Check Your Understanding
Work through each question, then compare with the solutions.
- In the binary-covariate model, let \(\theta=-1\) mmHg, \(\gamma=-5\) mmHg, \(p_H=0.20\), and \(p_C=0.60\). Calculate both marginal effects and their difference.
- If \(\gamma=0\), can different covariate prevalences alone produce different within-trial marginal treatment effects under the stated model? Explain.
- Two trials have conditional effects of \(-2\) mmHg at \(X=0\) and \(-6\) mmHg at \(X=1\). What is the standardized effect for a target population with \(p_*=0.25\)?
- Why can standardizing to a target distribution with little or no representation in one trial be unreliable?
- What does a commensurate prior on \(\theta_C\) and \(\theta_H\) fail to establish if treatment-effect modification differs by trial?
Solutions. (1) \(\Delta_H=-1+(-5)(0.20)=-2.0\) mmHg and \(\Delta_C=-1+(-5)(0.60)=-4.0\) mmHg, so \(\Delta_C-\Delta_H=-2.0\) mmHg. (2) No. If \(\gamma=0\), both effects equal \(\theta\), independent of \(p_s\). (3) Here \(\theta=-2\) and \(\gamma=-6-(-2)=-4\) mmHg, so \(\Delta(p_*)=-2+(-4)(0.25)=-3.0\) mmHg. (4) The estimate relies on extrapolating the fitted conditional effect into covariate regions with little or no data, making it sensitive to model assumptions. (5) It does not establish agreement in the treatment-effect modification parameter or in standardized effects away from \(X=0\); those require their own modeling and assessment.