Why Combine External and Concurrent Controls?
A randomized trial’s concurrent control group supports a direct comparison with the treated group because treatment assignment is randomized within the trial. External controls may add information about the control outcome, but they come from a different data source and were not randomized alongside the treated patients. Combining the groups can improve precision only if the analysis addresses differences in patient mix, outcome measurement, and care over time. The external data do not acquire the protection of randomization merely by being included in the same model.
The central question is therefore not simply how many control observations are available. It is how much the external controls should inform the control outcome for the trial population, given the evidence that the populations and data-generating processes are comparable. The target population should be specified first. In this tutorial, the target is the eligible population represented by the randomized trial, and the treatment effect is the average difference between assignment to the trial treatment and assignment to concurrent control in that population.
Define the Comparison Before Borrowing
Let \(A_i=1\) indicate assignment to the trial treatment and \(A_i=0\) indicate assignment to concurrent control. Let \(S_i=0\) identify trial participants and \(S_i=1\) identify external controls. Let \(X_i\) denote measured baseline covariates, and let \(Y_i\) be the outcome at the prespecified assessment time. The treatment effect of interest may be an average mean difference, a risk difference, or another clinically defined contrast. The estimand must specify the outcome, time point, treatment strategy, and target population.
A useful outcome model for a continuous endpoint is
Here \(\beta\) is the treatment contrast conditional on the included covariates under the model, while \(\delta\) represents a source difference between external and trial controls after covariate adjustment. The model assumes a common treatment effect across covariate values unless treatment-by-covariate interactions are added. It also assumes that the source difference can be represented adequately by the stated source term. A more elaborate model may be needed if outcomes are clustered, binary, or time-to-event.
For a population-average effect, the analysis should align comparisons to the trial target population. One approach is to model outcomes conditional on \(X_i\) and standardize predicted outcomes to the trial covariate distribution. Another is to construct transport weights for eligible external controls so their measured covariate distribution resembles that of the target trial population. Either approach requires adequate overlap: covariate patterns common in the trial must be represented in the external data. A conditional regression coefficient is not automatically the same as a marginal population-average effect, particularly for nonlinear outcome models.
Adjustment and the Role of the Source Difference
Covariate adjustment is intended to make the external controls informative about the control outcome in the trial population, not to erase every distinction between data sources. Baseline severity, age, prior treatment, calendar time, and eligibility criteria may differ. Adjusting for measured covariates can address observed differences, but residual confounding and source-specific measurement or care remain possible. The earlier material on transport weights and evaluating external-control selection bias provides relevant context for assessing these issues.
A Bayesian analysis can encode uncertainty about the residual source difference. For example, assign \(\delta \sim N(0,\tau^2)\). A small \(\tau\) says that large source differences are considered unlikely and allows stronger borrowing; a large \(\tau\) allows substantial source differences and reduces borrowing. This is a modeling choice that should be clinically justified and examined in sensitivity analyses. It is not evidence that the sources are exchangeable.
The following simplified normal model illustrates the calculation after external-control outcomes have been adjusted to the trial covariate distribution. Let \(T\) be the estimated treatment-group mean, \(C\) the concurrent-control mean, and \(E\) the adjusted external-control mean. Suppose their sampling variances are \(v_T\), \(v_C\), and \(v_E\). Write
Here \(\mu\) is the target-population control mean, and \(\beta\) is the treatment-minus-control mean difference. This model assumes independent sampling errors across the three disjoint groups, known or suitably estimated variances, and a common target-population mean after adjustment. With a flat prior for \(\mu\), integrating over \(\delta\) gives \(E\mid\mu \sim N(\mu,v_E+\tau^2)\). The external estimate thus contributes less precision when the assumed source difference is more uncertain.
Deriving the Combined Posterior
The concurrent and adjusted external control estimates are conditionally independent normal observations of \(\mu\), with variances \(v_C\) and \(v_E+\tau^2\). Multiplying their normal likelihoods and completing the square yields the posterior for \(\mu\). With a flat prior for \(\beta\), the treated-group mean contributes information about \(\beta\), not about \(\mu\): for any \(\mu\), the unrestricted \(\beta\) can account for the observed treated mean.
Conditional on \(\mu\), the treatment contrast is \(T-\mu\) with sampling variance \(v_T\). Averaging over uncertainty in \(\mu\) gives
The posterior variance contains both uncertainty in the treatment mean and uncertainty in the target control mean. The formulas show precisely how borrowing operates: the adjusted external estimate receives precision \(1/(v_E+\tau^2)\), rather than its sampling precision alone. When \(\tau\) is large, its contribution approaches zero.
Worked Example: Posterior Estimate With Partial Borrowing
Suppose a hypothetical randomized trial measures a symptom score at six months, where lower scores are better. The treatment group has \(n=100\), mean \(T=10.2\), and standard deviation \(4\). The concurrent-control group has \(n=50\), mean \(C=12.0\), and standard deviation \(4\). After covariate adjustment and standardization to the trial population, hypothetical external controls have \(E=12.5\), \(n=100\), and standard deviation \(4\). For illustration, treat these standard deviations as known and set \(\tau=1.5\) score points. Assume the adjustment is adequate, the outcome measurement is harmonized, and the independent normal sampling model is reasonable.
The sampling variances are \(v_T=4^2/100=0.16\), \(v_C=4^2/50=0.32\), and \(v_E=4^2/100=0.16\). The external-control variance after allowing for source difference is \(v_E+\tau^2=0.16+2.25=2.41\). Therefore,
The combined target-control mean is
Thus the posterior mean treatment contrast is \(10.2-12.06=-1.859\) score points, and its posterior standard deviation is \(\sqrt{0.16+0.2825}=0.6652\). A normal 95% credible interval is \(-1.859\pm1.96(0.6652)\), or approximately \((-3.163,-0.555)\) points. Under this model and prior, the posterior probability that treatment lowers the mean score is \(P(\beta<0\mid\text{data})\), approximately 0.997. This is a probability statement about the parameter under the specified Bayesian model, not a frequentist coverage claim. The result concerns a mean difference in symptom-score points; its clinical importance requires a clinically meaningful threshold.
Worked Example: Why Naïve Pooling Can Mislead
Consider a different hypothetical setting with a continuous outcome for which lower is better. The randomized trial has 40 concurrent controls with mean 12.0 and 50 treated participants with mean 10.5. Sixty external controls have mean 15.0 because, before adjustment, they are older and more severely ill. Simply pooling all 100 controls gives
The naïve treatment-minus-control contrast is \(10.5-13.8=-3.3\) outcome units. The concurrent randomized comparison is \(10.5-12.0=-1.5\) units. The difference between these estimates is not evidence that external controls reveal a larger treatment effect; it can arise because the unadjusted external group has a different prognosis. Suppose, for illustration, covariate adjustment standardizes the external mean to 12.4 in the trial population. That adjusted value is closer to the concurrent mean, but it does not alone establish exchangeability or justify giving the external group full weight. A borrowing model should still account for residual uncertainty about source differences.