Tutorials › Bayesian Clinical Trials › Inverse Probability Weighting for External Controls

External and Real-World Controls · Tutorial 370 of 696

Inverse Probability Weighting for External Controls

Learn to derive and calculate mean-one stabilized weights for external controls, use them to estimate trial-population control outcomes, and assess the consequences of extreme weights.

Advanced 12 min read

What You'll Learn

  • Distinguish a source-membership propensity score from a treatment propensity score.
  • Derive the external-control density ratio from pooled source-membership odds.
  • Calculate stabilized weights and use them to estimate a trial-population control mean.
  • Identify the assumptions required for weighting to support a valid comparison.
  • Evaluate weight truncation as a variance–bias tradeoff rather than an automatic fix.
  • Recognize diagnostics and limitations of inverse probability weighting.

Why Reweight External Controls?

External controls may differ from trial participants in baseline characteristics that predict outcomes. Comparing their unadjusted outcomes can therefore mix differences in patient composition with differences in outcomes under the control strategy. Inverse probability weighting addresses this problem by assigning more influence to external controls whose measured covariate profiles resemble those of trial participants, and less influence to profiles that are relatively common in the external cohort.

The weights in this tutorial concern source membership: they estimate how the distribution of baseline covariates in the external cohort must be reweighted to represent the trial population. They are not treatment-propensity weights. The earlier material on estimating propensity scores provides the basic idea of modeling a probability from covariates; here the outcome of that model is trial versus external-cohort membership.

Definition: Let \(S=1\) identify a trial participant and \(S=0\) identify an external-control participant, and let \(X\) denote measured baseline covariates. The source-membership propensity score is \(e(X)=P(S=1\mid X)\), estimated in the pooled trial and external-control data. Write \(q_T=P(S=1)\) and \(q_E=P(S=0)\) for the trial and external proportions in that pooled data. For an external control, the unstabilized inverse probability weight is the density ratio \(f_T(X)/f_E(X)\), where \(f_T\) and \(f_E\) are the covariate densities in the trial and external cohorts.

This density ratio is large for covariate profiles that are more prevalent in the trial than in the external cohort. Weighting by it changes the external cohort's covariate distribution toward the trial's. Stabilization multiplies every external-control weight by the same factor so that, under the population model, its mean among external controls is one.

Deriving the Weight

The pooled source-membership probability depends both on the covariate distributions and on the proportions of trial and external participants in the pooled fitting sample. Bayes' rule gives the source-membership odds as follows:

$$ \frac{e(x)}{1-e(x)} = \frac{q_T f_T(x)}{q_E f_E(x)}. $$

Rearranging isolates the density ratio needed to reweight external controls:

$$ \frac{f_T(x)}{f_E(x)} = \frac{q_E}{q_T}\frac{e(x)}{1-e(x)}. $$

Thus the mean-one stabilized weight for external control \(i\) is \(w_i=(q_E/q_T)\hat e(X_i)/(1-\hat e(X_i))\), using the fitted source-membership score \(\hat e\). The factor is \(q_E/q_T\), not its inverse. The correction matters because the fitted odds reflect the pooled sample's source proportions as well as the covariate-density ratio.

Justification: For any integrable function \(g(X)\), weighting the external covariate distribution by the density ratio gives
$$ E_E\left[\frac{f_T(X)}{f_E(X)}g(X)\right] = \int f_E(x)\frac{f_T(x)}{f_E(x)}g(x)\,dx = \int f_T(x)g(x)\,dx = E_T[g(X)]. $$

Setting \(g(X)=1\) shows that the expected unstabilized weight among external controls is one. The stabilized expression above equals that same density ratio, so it also has expected value one. In finite samples, its observed average need not be exactly one, particularly when scores are estimated or the model is imperfect.

The derivation requires the pooled source model to represent the relevant sampling process, and the estimated score to capture the covariate differences needed for reweighting. In practice, a logistic regression or another suitable model can estimate \(e(X)\). The model should allow important nonlinearities and interactions when needed; a simple model is not automatically adequate merely because it produces scores between zero and one.

Worked Example: Calculating a Stabilized Weight

Suppose a hypothetical pooled fitting sample is 40% trial participants and 60% external controls, so \(q_T=0.40\) and \(q_E=0.60\). A particular external control has an estimated trial-membership score of \(\hat e(X_i)=0.75\). Calculate the mean-one stabilized weight.

The fitted source odds are \(0.75/(1-0.75)=0.75/0.25=3\). The stabilization factor is \(q_E/q_T=0.60/0.40=1.5\). Therefore,

$$ w_i = \frac{0.60}{0.40}\frac{0.75}{0.25} = 1.5(3) = 4.500. $$

This control receives 4.5 times the influence of a control with weight one in a weighted sum. The weight is not a 4.5-fold probability of being in the trial, nor is it a posterior probability. It reflects the estimated trial-to-external covariate-density ratio for this patient's measured profile. Here, the profile is estimated to be much more prevalent in the trial than in the external cohort.

Using Weights to Estimate a Control Outcome

For a measured outcome \(Y\), a normalized weighted mean among external controls is

$$ \bar Y_{E,w} = \frac{\sum_{i:S_i=0}w_iY_i}{\sum_{i:S_i=0}w_i}. $$

Under the density-ratio identity, this estimates the mean control outcome standardized to the trial's measured covariate distribution, provided the outcome is comparable across sources conditional on the adjusted covariates. The normalized form also makes the estimate invariant to multiplying every weight by the same positive constant. That invariance does not make the correct stabilization factor irrelevant: the mean-one convention remains useful for interpreting and reporting weights, and other analyses may depend on their scale.

Worked Example: Reweighting Two Covariate Strata

Consider hypothetical data with 100 trial participants and 150 external controls. Let \(X\) indicate a baseline risk factor. Among trial participants, 40 have \(X=1\) and 60 have \(X=0\). Among external controls, 30 have \(X=1\) and 120 have \(X=0\). Suppose the mean control outcomes in the two external strata are 10 and 6 points, respectively. The goal is to estimate the external-control mean standardized to the trial's distribution of \(X\).

The pooled proportions are \(q_T=100/250=0.40\) and \(q_E=150/250=0.60\). For \(X=1\), the fitted source membership proportion is \(e(1)=40/(40+30)=4/7\), so its odds are \(4/3\). Its stabilized weight is

$$ w(1)=\frac{0.60}{0.40}\frac{4/7}{3/7} =1.5\left(\frac{4}{3}\right)=2. $$

For \(X=0\), \(e(0)=60/(60+120)=1/3\), giving odds \(1/2\) and stabilized weight

$$ w(0)=\frac{0.60}{0.40}\frac{1/3}{2/3} =1.5\left(\frac{1}{2}\right)=0.75. $$

The weighted external counts are \(30(2)=60\) for \(X=1\) and \(120(0.75)=90\) for \(X=0\). Their proportions, \(60/150=0.40\) and \(90/150=0.60\), match the trial proportions. The weighted outcome mean is

$$ \bar Y_{E,w} = \frac{60(10)+90(6)}{60+90} = \frac{1140}{150} = 7.600\text{ points}. $$

The unweighted external mean would be \(\{30(10)+120(6)\}/150=6.800\) points. The difference arises because the external cohort overrepresents the lower-outcome \(X=0\) stratum relative to the trial. This calculation verifies balance for the single binary covariate in this example; with multiple covariates, matching all relevant features is not guaranteed by fitting a source model alone.

Assumptions and Interpretation

Weighting can align measured covariate distributions, but it does not by itself establish that trial and external outcomes are exchangeable. A causal comparison using a trial treatment group and weighted external controls needs several substantive conditions:

  • Consistency and comparable strategies: the control strategy and outcome have the same relevant meaning in both sources. Differences in endpoint definitions, follow-up, or treatment versions can undermine this condition.
  • Conditional source exchangeability for control outcomes: given the measured covariates used for adjustment, the expected control outcome for eligible trial participants is represented by the external controls. Important unmeasured prognostic differences violate this condition.
  • Positivity: covariate profiles in the target trial population must have representation in the external cohort. Formally, the external density must be positive wherever the trial density is positive, \(f_E(x)>0\) whenever \(f_T(x)>0\).
  • Adequate source-score estimation: the fitted score must capture the covariate distribution differences relevant to the weighting. A misspecified model can leave important imbalance.

Under these assumptions, the weighted external mean can serve as an estimate of the control outcome mean in the trial population. A contrast with the trial's treatment-group outcome mean may then estimate a treatment effect for that population, subject to the trial's treatment-assignment design and outcome analysis. The estimate is not automatically causal merely because weights were calculated.

Extreme Weights: A Practical Tradeoff

A score near one for an external control makes \(\hat e/(1-\hat e)\) large. Such a control has a covariate profile estimated to be much more common in the trial than in the external cohort. It may be essential for representing part of the trial population, but a very large weight can also make the result depend heavily on one or a few observations. Stabilization changes the overall scale; it does not remove large relative weights.

Analysts should inspect the estimated score and weight distributions, identify which records carry the largest weights, and examine whether those records have plausible covariate values and reliable outcome measurements. They should also check whether weighted covariate summaries move toward the trial's summaries. These checks can reveal poor overlap, model problems, or data errors; they cannot prove that unmeasured confounding is absent.

Worked Example: Assessing a Weight Cap

Suppose five hypothetical external controls have estimated stabilized weights \(1,1,1,1,24\), with outcomes of \(5,6,7,8,20\) points. For illustration, suppose the first four weights correspond to a score of \(0.50\) when \(q_T=q_E\), while the last corresponds to a score of \(0.96\), since \(0.96/0.04=24\). Calculate the weighted mean before and after capping weights at 10.

Without a cap, the total weight is \(1+1+1+1+24=28\), and the weighted outcome total is \(5+6+7+8+24(20)=506\). Thus the weighted mean is \(506/28=18.07\) points.

After capping the final weight at 10, the total weight is \(1+1+1+1+10=14\), and the weighted outcome total is \(5+6+7+8+10(20)=226\). The capped weighted mean is \(226/14=16.14\) points. Capping reduces the influence of the extreme record, but changes the weighted estimate by about 1.93 points. It trades reduced sensitivity to that observation for a departure from the original density-ratio weighting.

Weight truncation or capping should therefore be treated as a sensitivity analysis, not a universal repair. If used, the rule should be justified and specified transparently; report the uncapped estimate and how the result changes under the chosen rule. When extreme weights indicate that important trial profiles are scarcely represented in the external cohort, no numerical cap can create the missing information. The later material on overlap weighting addresses a different target under limited comparability; it should not be silently substituted for trial-population weighting.

Alternatives and Limitations

Matching, introduced in the earlier material on matching trial participants to external controls, selects or pairs controls to improve comparability. Weighting instead retains external observations and changes their contributions to a summary. Matching can discard observations and may not represent the full trial covariate distribution; weighting can retain them but may be unstable when some profiles receive very large weights.

Outcome regression is another adjustment strategy: it models the outcome conditional on covariates rather than reweighting the control cohort. It relies on its own model assumptions and is treated separately in the later material on outcome regression. These approaches are not interchangeable by default. In particular, weighting does not repair incompatible endpoints, unmeasured confounding, source-specific measurement error, or severe lack of covariate support.

Finally, weights are estimated quantities. Treating them as fixed when calculating uncertainty can understate uncertainty, and treating weighted records as literal replicated patients can misrepresent the data. In Bayesian analyses, weighting must be incorporated in a way appropriate to the model and inferential target; a weight is not itself a likelihood contribution or a probability. The present calculations establish the adjustment construction, not a complete inferential procedure.

Key takeaway: The external-control density-ratio weight is \(f_T(X)/f_E(X)=(q_E/q_T)e(X)/(1-e(X))\). Its mean-one stabilization uses the external-to-trial pooled proportion ratio \(q_E/q_T\). Valid interpretation still depends on measured-covariate exchangeability, positivity, comparable outcomes, and adequate score estimation.

Check Your Understanding

Try each question, then check the solutions below.

  1. If \(q_T=0.30\), \(q_E=0.70\), and an external control has \(\hat e(X)=0.60\), what is the mean-one stabilized weight?
  2. Why does a mean-one stabilized weight use \(q_E/q_T\) rather than \(q_T/q_E\)?
  3. In a finite sample, what does a very large external-control weight suggest, and name one concern it raises.
  4. Does matching measured covariates with weights establish that external controls have no unmeasured confounding? Explain briefly.

Solutions. (1) The odds are \(0.60/0.40=1.5\), so \(w=(0.70/0.30)(1.5)=3.500\). (2) Bayes' rule gives \(e(x)/(1-e(x))=(q_T/q_E)f_T(x)/f_E(x)\); rearrangement therefore requires \(q_E/q_T\) to recover the density ratio. (3) The profile may be much more common in the trial than in the external cohort; the estimate may be dominated by that record, and limited support should be investigated. (4) No. Weighting adjusts measured covariates represented in the model; unmeasured prognostic differences can remain.

Clinical Biostats