Tutorials › Multi-Arm Multi-Stage (MAMS) Designs › Interpreting Estimates Augmented With Nonconcurrent Controls

Platform and Umbrella Trials · Tutorial 386 of 499

Interpreting Estimates Augmented With Nonconcurrent Controls

Learn to distinguish the population an augmented estimate targets from the historical data used to improve its precision, and to interpret its assumptions, uncertainty, and sensitivity.

Intermediate 10 min read

What You'll Learn

  • Define the target-population treatment contrast that an augmented estimate is intended to estimate.
  • Distinguish control-evidence weights from weights that define a target population.
  • Calculate an augmented estimate and its standard error under stated independence assumptions.
  • Interpret confidence intervals while recognizing uncertainty not captured by the standard error.
  • Assess how the estimate changes when the contribution from nonconcurrent controls changes.
  • Compare augmented estimates with concurrent-only and pooled-control estimates.

What an Augmented Estimate Says

A treatment-effect estimate that uses nonconcurrent controls is not interpreted just by reading off its point estimate. One must first ask which population’s treatment effect is being estimated, then determine what the nonconcurrent data contribute and which assumptions make that contribution relevant. A narrow confidence interval does not answer those questions: it describes sampling uncertainty under the analysis model, not whether historical and current patients are sufficiently comparable.

For a clear starting point, consider a continuous outcome measured at a fixed follow-up time, with lower values preferred. Let \(Y(1)\) and \(Y(0)\) denote the outcomes a patient would have under treatment and control, respectively. Let \(T\) denote the population of patients eligible for the current treatment comparison. The target effect is the mean treatment effect in that population. The same interpretation applies to other outcome scales, provided the estimand and model are specified for that scale.

Definition: The target-population effect is \(\Delta_T=E_T[Y(1)-Y(0)]\), the mean difference between the potential outcome under treatment and the potential outcome under control for patients in the target population \(T\). For the examples below, a negative value favors treatment because lower outcomes are better. A nonconcurrent control is a control observation collected in a different calendar period from the current treatment arm.

An augmented analysis combines information from concurrent controls, enrolled during the same period as the treatment arm, and nonconcurrent controls, enrolled earlier or otherwise in a different period. To interpret the combination as an estimate of \(\Delta_T\), the control means must represent the target population’s expected outcome under control, either directly or after suitable adjustment. In the examples, assume that each adjusted control estimate is intended to represent \(E_T[Y(0)]\).

The weight assigned to each source of control information is not automatically the proportion of patients in the target population who came from that period. It may instead be an analysis weight describing how much each estimate contributes to the control reference. If raw data are pooled, or the analysis explicitly defines a mixture of populations, the resulting target may differ from \(T\). The estimand must be read from the actual standardization and weighting scheme, not inferred from the word “augmented.”

Assumptions Behind the Interpretation

In a randomized concurrent comparison, randomization supports a causal comparison for the patients randomized in the current period, subject to the trial’s treatment, outcome, and follow-up definitions. Adding nonconcurrent controls requires further justification: their untreated outcomes must inform the untreated outcomes that would have been observed in the current target population.

One way to state the key requirement is conditional exchangeability across periods. Given the measured covariates used in the adjustment, calendar period must not leave a systematic difference in the control outcome that is relevant to the target contrast. The data must also have adequate covariate overlap: the target population’s covariate patterns should be represented among the nonconcurrent controls used to estimate its control outcome. The treatment, outcome measurement, eligibility criteria, follow-up, and standard care must be sufficiently aligned or explicitly accounted for.

These assumptions are substantive, not consequences of a small standard error. If unmeasured changes over time affect outcomes, an adjusted historical mean can remain biased even with a very large historical sample. Earlier material on sensitivity analysis for unmeasured temporal confounding and on trimming and restriction when historical controls lack overlap addresses ways to probe or limit these concerns; the resulting estimate still needs a target-specific interpretation.

Interpretation rule: An augmented estimate can be described as an estimate of the target-population treatment effect only if each control contribution is relevant to that target under the stated adjustment and transport assumptions. Its conventional standard error and confidence interval quantify sampling uncertainty conditional on the analysis model; they do not, by themselves, include uncertainty from unmeasured temporal differences or a misspecified adjustment.

Estimate, Weight, and Uncertainty

Suppose \(\hat\mu_{C}\) is the concurrent-control estimate of the target control mean and \(\hat\mu_{N}\) is the nonconcurrent-control estimate after adjustment to that same target. Let \(\bar Y_T\) estimate the target mean under treatment. For a fixed weight \(0\leq\lambda\leq1\), one simple augmented estimator is

$$ \hat\Delta_{\lambda} = \bar Y_T-\{\lambda\hat\mu_C+(1-\lambda)\hat\mu_N\}. $$

Here \(\lambda\) is the contribution from the concurrent estimate; \(1-\lambda\) is the contribution from the nonconcurrent estimate. The formula makes the sign convention explicit: a negative estimate means the estimated outcome under treatment is lower than the weighted control reference. The weight can affect both the point estimate and its standard error.

For independent treatment, concurrent-control, and nonconcurrent-control samples, and fixed weights, the variance follows from the variance of a linear combination. If the three estimates have variances \(V_T\), \(V_C\), and \(V_N\), then

$$ \operatorname{Var}(\hat\Delta_{\lambda}) = V_T+\lambda^2V_C+(1-\lambda)^2V_N. $$

The squared weights appear because \(\operatorname{Var}(aX)=a^2\operatorname{Var}(X)\). Independence makes the covariance terms zero; if estimates share observations or are otherwise dependent, the relevant covariance terms must be included. When adjustment parameters or weights are estimated, uncertainty from estimating them must also be represented, for example through an appropriate variance estimator or resampling procedure. The simple formula is not automatically valid for every augmented analysis.

Worked Example: A Continuous Outcome With an Adjusted Historical Mean

In a hypothetical trial, the outcome is a symptom score in points, with lower scores better. The treatment group has mean 7.0 points, \(n=90\), and sample standard deviation \(s=4.0\). The concurrent-control group has mean 10.0 points, \(n=100\), and \(s=4.0\). After adjustment to the current target population, the nonconcurrent-control estimate is 12.0 points, with standard error 0.559 points. Use \(\lambda=0.60\). Assume the groups are independent, the adjusted historical standard error includes uncertainty from its adjustment, and the target and transport assumptions stated above are considered plausible for this illustration.

The treatment mean standard error is \(4/\sqrt{90}\), so its estimated variance is \(16/90=0.1778\) square points. The concurrent-control mean variance is \(16/100=0.1600\) square points. The historical estimate variance is \(0.559^2=0.3125\) square points.

First calculate the weighted control reference:

$$ 0.60(10.0)+0.40(12.0)=10.8\text{ points}. $$

Thus the treatment-effect estimate is \(7.0-10.8=-3.8\) points. Its estimated variance and standard error are

$$ \widehat{\operatorname{Var}}(\hat\Delta_{\lambda}) = 0.1778+(0.60)^2(0.1600)+(0.40)^2(0.3125) = 0.2854, \qquad SE=\sqrt{0.2854}=0.5342\text{ points}. $$

Using a normal 95% interval, \(-3.8\pm1.96(0.5342)\) gives approximately \([-4.85,-2.75]\) points. Conditional on the model and its assumptions, the interval estimates that the mean symptom score under treatment in the target population is about 2.75 to 4.85 points lower than under control. It is not a 95% probability statement about the fixed parameter, and the interval does not protect against bias from an unmeasured difference between periods.

How Dependence on Historical Data Changes the Result

An estimate can be precise and still depend materially on historical data. A useful interpretation therefore reports the weight, the adjusted historical estimate, and how the result compares with the concurrent-only estimate. The concurrent-only analysis uses \(\lambda=1\); it may be less precise, but it does not rely on the historical-control transport assumption. Full reliance on the adjusted historical estimate corresponds to \(\lambda=0\), which is generally more exposed to temporal differences.

Worked Example: Changing the Evidence Weight

Suppose, hypothetically, the treatment mean is 8.0 points, the concurrent-control mean is 11.0 points, and the adjusted nonconcurrent-control mean is 13.0 points. These estimates target the same current population. Let their sample sizes and standard deviations be \(n_T=90,s_T=4\), \(n_C=100,s_C=4\), and \(n_N=80,s_N=5\), respectively. Assume independent groups and fixed weights. With \(\lambda=0.75\), the weighted control mean is

$$ 0.75(11.0)+0.25(13.0)=11.5\text{ points}, \qquad \hat\Delta_{0.75}=8.0-11.5=-3.5\text{ points}. $$

The estimated variance is

$$ \frac{4^2}{90} +(0.75)^2\frac{4^2}{100} +(0.25)^2\frac{5^2}{80} = 0.2873, \qquad SE=\sqrt{0.2873}=0.5360\text{ points}. $$

The corresponding normal 95% interval is \(-3.5\pm1.96(0.5360)\), or approximately \([-4.55,-2.45]\) points. With concurrent controls alone, the estimate would be \(8.0-11.0=-3.0\) points. With \(\lambda=0.50\), the weighted control mean would be 12.0 points and the estimate \(-4.0\) points. These changes show how the historical data pull the estimate; they do not show which weight is scientifically correct. Because both control estimates were assumed to target the same population, changing \(\lambda\) changes the evidence combination, not the intended target population.

If instead the weights were used to define a mixture of different populations, changing them could change the estimand itself. That is a different interpretation and should be stated explicitly rather than described as merely borrowing more historical information.

Binary Outcomes and Scale-Specific Meaning

For a binary outcome, an augmented estimate may be a risk difference, a risk ratio, or another contrast. The scale determines the units and interpretation. The following example uses a risk difference, where a negative value indicates a lower event risk under treatment. The approximate standard error below uses independent binomial groups and treats the adjusted historical risk estimate as a binomial proportion; more complex adjustment requires a variance calculation that accounts for that adjustment.

Worked Example: Augmented Risk Difference

In a hypothetical study of a binary adverse event, the treatment risk is \(0.18\) among 200 patients, the concurrent-control risk is \(0.25\) among 160 patients, and the adjusted nonconcurrent-control risk is \(0.30\) among 150 patients. Assume the adjusted historical risk targets the current population, estimates are independent, and \(\lambda=0.50\).

The weighted control risk is \(0.50(0.25)+0.50(0.30)=0.275\). The estimated risk difference is \(0.18-0.275=-0.095\), or \(-9.5\) percentage points.

Using the binomial variance formula \(p(1-p)/n\) for each group, the estimated variance is

$$ \frac{0.18(0.82)}{200} +(0.50)^2\frac{0.25(0.75)}{160} +(0.50)^2\frac{0.30(0.70)}{150} = 0.001381. $$

Thus \(SE=\sqrt{0.001381}=0.03716\), or 3.716 percentage points. A normal 95% interval is \(-0.095\pm1.96(0.03716)\), approximately \([-0.168,-0.022]\), equivalent to a reduction of about 2.2 to 16.8 percentage points under the model. This interval concerns sampling uncertainty on the risk-difference scale; it does not account for possible bias in transporting the historical risk.

Diagnostics, Alternatives, and Limits

Before interpreting an augmented estimate, check whether its construction and assumptions are credible. Useful questions include:

  • Do the treatment, control, eligibility, outcome definition, follow-up, and standard-care conditions align across periods, or are differences explicitly handled?
  • Do measured baseline covariates overlap between the target patients and the historical controls used? Large weights or a lack of comparable historical patients can make an adjusted estimate unstable or dependent on extrapolation.
  • Do adjusted covariate distributions and outcome patterns look plausible, and are there indications of calendar-time or site-related changes in control outcomes?
  • How much do the estimate and interval change under reasonable alternative weights, adjustment specifications, or a concurrent-only analysis?

These checks can reveal fragility; they cannot prove that all important temporal differences have been removed. A reported standard error usually reflects the chosen sampling model. It does not automatically include uncertainty about the choice of historical period, the correctness of the transport assumptions, or unmeasured confounding. Sensitivity analyses make this dependence more visible, but do not convert untestable assumptions into verified facts.

ApproachWhat it estimates or usesMain interpretive trade-off
Concurrent controls onlyCurrent treatment effect using randomized concurrent comparisonsLess dependent on historical transport, but may be less precise
Augmented controlsA weighted control reference using concurrent and adjusted nonconcurrent informationMay improve precision, but its interpretation depends on transport assumptions and the specified weight
Unadjusted pooled controlsA reference formed by pooling control observations across periodsSimple pooling can mix different populations or standards of care and may not target the current population

The comparison is not simply “more data versus fewer data.” Concurrent-only estimation can be a useful reference for judging how much an augmented result relies on historical information. Pooling all controls without suitable adjustment may alter the target or introduce bias. An augmented analysis is most informative when the target, the role of each control source, and the sensitivity of the result are made explicit.

Common Mistake

Wrong reasoning: “The historical controls make the sample larger, so the augmented estimate is more reliable and its narrow interval includes all relevant uncertainty.” A larger historical sample can reduce the modeled sampling variance, but it cannot ensure that historical and current patients are comparable. If control outcomes changed over time for reasons not captured in the adjustment, the estimate may be biased while its conventional interval remains narrow.

Correction: State the target population and assumptions first, report the concurrent and nonconcurrent contributions, and describe the interval as conditional on the model. Then compare with a concurrent-only estimate and examine sensitivity to the historical contribution. Precision and validity are separate properties.

Check Your Understanding

Answer each question, then compare your reasoning with the solutions.

  1. An estimate uses an adjusted historical-control mean intended to represent the current eligible population. What population does the estimate target, and what must be true for the historical mean to be relevant?
  2. The treatment mean is 6 points, the concurrent-control mean is 9 points, the adjusted nonconcurrent-control mean is 11 points, and \(\lambda=0.50\). Calculate the weighted control mean and treatment-effect estimate. Lower outcomes are better.
  3. Why can an augmented estimate have a smaller standard error than a concurrent-only estimate but still be less trustworthy?
  4. If \(\lambda\) is an evidence-combination weight and both control estimates target the same current population, does changing \(\lambda\) necessarily change the estimand? Explain.

Solutions. (1) It targets the mean treatment effect in the current eligible population. The adjusted historical controls must represent that population’s expected control outcome, requiring adequate measured-covariate overlap and no important residual period differences after adjustment. (2) The weighted control mean is \(0.50(9)+0.50(11)=10\) points, so the estimate is \(6-10=-4\) points, favoring treatment on this scale. (3) Its standard error can decrease because more control information is included, while unmeasured temporal differences can still bias the estimate and are not necessarily represented in that standard error. (4) No. Under the stated condition the target remains the current population; changing \(\lambda\) changes the evidence combination and potentially the estimate and precision, not the estimand. If the weights instead define a population mixture, changing them can change the estimand.

Clinical Biostats