Why Integration Errors Matter
An integrated evidence chain can be numerically tidy and still be misleading. A model may pass an exposure summary to a response model, which passes an estimate to a trial-outcome model. At each handoff, a mismatch in what is being estimated, an unrecorded dependence, or discarded uncertainty can make the final prediction look more certain—or more coherent—than the underlying evidence warrants. The concern is not that different models must agree. It is that any comparison or combination must preserve the meaning and uncertainty of the quantities being connected.
This tutorial focuses on five recurring errors: combining incompatible exposure metrics, treating correlated estimates as independent, reporting unjustified parameter precision, counting the same data more than once, and losing uncertainty when passing estimates between models. These failures can overlap. For example, reusing data can create dependence, while passing only a point estimate can conceal both its uncertainty and its dependence on other inputs.
Check Whether the Quantities Match
An exposure metric is part of the definition of a model input, not merely a label for a concentration summary. AUC, measured in concentration-by-time units, and \(C_{\max}\), measured in concentration units, describe different aspects of a concentration-time profile. A response model using AUC and one using \(C_{\max}\) therefore do not have interchangeable slopes. Their coefficients have different units and generally answer different questions.
Worked Example: Incompatible Exposure Metrics
Suppose, hypothetically, two response models describe change in the same biomarker, measured in units. Model A uses AUC and estimates a slope of \(0.20\) biomarker units per \(\mathrm{mg\cdot h/L}\). Model B uses \(C_{\max}\) and estimates a slope of \(1.5\) biomarker units per \(\mathrm{mg/L}\). For a hypothetical regimen with AUC \(30\ \mathrm{mg\cdot h/L}\) and \(C_{\max}=4\ \mathrm{mg/L}\), the model-based changes are:
The predictions happen to be numerically equal for this regimen. The slopes are not comparable: averaging \(0.20\) and \(1.5\) would combine quantities with different units. Even equal predictions do not show that the models represent the same exposure-response relationship. A meaningful comparison requires a shared target, such as both models’ predicted biomarker change for the same regimen, and appropriate uncertainty for each prediction. If one model must be translated to the other exposure metric, that translation needs a justified relationship between the metrics rather than an assumption that they are interchangeable.
A useful diagnostic is to write the target, population, regimen or scenario, endpoint, and units beside every quantity at an interface. Then check whether the upstream model’s output is exactly the input the downstream model requires. If the quantities differ, the integration needs an explicit transformation or separate predictions on a common scale. Sensitivity to the choice of exposure definition can reveal fragility, but it does not itself establish which metric is correct.
Account for Dependence When Combining Estimates
Two model estimates may share participants, measurements, fitted parameters, or other evidence. Their errors can therefore be correlated. For a weighted sum of two estimates \(T_1\) and \(T_2\), the variance is:
The covariance term is omitted only when independence is justified, or when the term is known to be zero for another reason. Positive covariance often means that the estimates contribute less independent information than a calculation assuming independence would suggest. Negative covariance can have the opposite effect. Dependence should be determined from how the estimates were generated, not chosen to make a combined result more convenient.
For two estimates, the result follows by minimizing the variance of \(wT_1+(1-w)T_2\). The weights sum to one, so the combination remains unbiased for their shared target. Differentiating that variance with respect to \(w\), setting the derivative to zero, and solving gives the covariance-aware weights above. If the covariance matrix is misspecified, the resulting precision is also misspecified.
Worked Example: Ignored Covariance Changes the Combination
Suppose two hypothetical unbiased model estimates target the same quantity, in outcome units. They are \(T_1=10\), with \(SE=2\), and \(T_2=12\), with \(SE=3\). Their correlation is \(0.60\). Assume the estimates are approximately jointly normal and that these standard errors and their correlation are known for this calculation. The covariance is \(0.60(2)(3)=3.6\), so:
Normalizing the two entries gives weights \(5.4/5.8=0.9310\) and \(0.4/5.8=0.0690\). The combined estimate and its standard error are:
If the correlation were incorrectly set to zero, inverse-variance weights would be \(9/13\) and \(4/13\). That calculation gives an estimate of \(10.62\) and \(SE=\sqrt{36/13}=1.664\). Here, ignoring positive covariance makes the combined result appear more precise. The difference is not evidence that one calculation is universally preferable: the correct calculation depends on the actual shared evidence and a valid covariance estimate.
Do Not Manufacture Precision or Reuse Evidence
A parameter estimate passed from one model to another is uncertain unless its value is genuinely fixed by design or known without relevant uncertainty. Substituting an estimate as though it were a constant usually removes variation from the downstream result. Precision can also be overstated when a likelihood or summary is used twice as if each use came from new data. These are distinct errors: the first drops uncertainty at a handoff; the second inflates the apparent amount of information.
Worked Example: Uncertainty Lost at a Model Handoff
Suppose, hypothetically, an upstream model estimates an exposure scaling parameter as \(\hat\theta=0.80\), with \(SE=0.12\). A downstream model estimates a response slope \(\hat\beta=2.0\), with \(SE=0.30\), in outcome units per unit of scaled exposure. For illustration, assume the two estimates are independent, approximately unbiased, and sufficiently well behaved for a first-order delta-method approximation. The predicted outcome change is \(g(\theta,\beta)=\theta\beta\), with estimate \(1.60\) outcome units. The first-order variance is:
Thus the propagated standard error is approximately \(\sqrt{0.1152}=0.3394\) outcome units. Treating the upstream estimate as exact would keep only the second variance contribution, giving \(SE\approx\sqrt{(0.80)^2(0.30)^2}=0.2400\). That understates uncertainty because it excludes the estimated exposure scaling. If the two estimates are dependent, the first-order variance also includes \(2(\partial g/\partial\theta)(\partial g/\partial\beta)\operatorname{Cov}(\hat\theta,\hat\beta)\); independence must not be assumed without justification.
A point estimate is not automatically a fixed input. A sound handoff carries forward the needed uncertainty and dependence information, for example through joint parameter draws or a suitable joint covariance matrix. For nonlinear models, propagating joint draws can retain distributional features that a local approximation may miss. The choice depends on the model and available information; merely reporting more decimal places does not improve precision.
Worked Example: Duplicated Data Create False Precision
Imagine a hypothetical dataset yields an estimate of \(10.0\) outcome units with \(SE=2.0\). Two downstream summaries are produced using the same estimator on the same participants and data; they therefore produce identical estimates of \(10.0\), and in this stipulated example each reports \(SE=2.0\). If someone incorrectly treats these as independent estimates and combines them by inverse-variance weighting, the result is \(10.0\) with:
The apparent gain in precision is fictitious: no new observations were added. With perfect dependence from exact reuse of the same estimator and data, the combined value is still \(10.0\) and its \(SE\) remains \(2.0\), not \(1.414\). Using the same data does not always force different analyses or estimators to give identical estimates. The identity here follows specifically from the stated same-estimator, same-data setup. In a real integration, trace each likelihood, dataset, and summary to its source before treating evidence contributions as independent.
Diagnostics, Alternatives, and Limits
An integration audit begins with the information at each handoff, not with the final standard error. Check the target and units; identify overlapping data and shared model components; record whether each input is an estimate, a fixed quantity, or a distribution; and inspect how uncertainty is propagated. Compare an integrated prediction with relevant model-specific predictions and examine whether conclusions change under plausible interface or dependence assumptions. Residual patterns and predictive checks can help locate a failing interface, but a lack of obvious discrepancy does not prove that data were independent or uncertainty was fully propagated.
Write what the integrated quantity predicts, for which population, regimen or scenario, endpoint, and units.
Map shared participants, data, priors, fitted parameters, and summaries. Mark any repeated contribution before combining results.
Verify that the upstream output is the downstream input required, including its definition, scale, units, and uncertainty.
Where dependence or interface choices are uncertain, evaluate plausible alternatives and report whether the substantive conclusion changes.
The main alternatives address different problems. Aligning models on a common predicted outcome can allow comparison when their internal exposure metrics differ, but does not make their coefficients interchangeable. A covariance-aware combination is preferable to an independence calculation when shared information is quantified; if dependence is unknown, sensitivity analysis is more honest than silently setting covariance to zero. Joint simulation or another propagation method can preserve input uncertainty across a handoff, whereas passing only point estimates is simpler but generally less complete. These approaches cannot repair a mismatched target, duplicated likelihood, or incorrect model specification.
Some limitations remain even after these checks. Covariance may be difficult to estimate when models were fit separately or when only summary outputs are available. Approximate propagation can be unreliable for highly nonlinear transformations, bounded parameters, or strongly skewed distributions. Model agreement is not independent confirmation when models share data or assumptions. A careful analysis therefore reports what was combined, what dependence and interface assumptions were used, and which uncertainties could not be represented.
Common pitfall: “The two models give similar numbers, so they can be averaged and the combined standard error will be smaller.” Similar numbers do not establish a common estimand, and averaging does not create independent information. First verify the target and units, then account for covariance and any reused data. If those checks cannot be made, present the outputs separately or show how the conclusion varies under explicit assumptions.
Check Your Understanding
For each question, identify the relevant integration risk before checking the solution.
- A model estimates a response per unit AUC and another estimates a response per unit \(C_{\max}\). Is averaging the slopes a defensible way to combine them? What could be compared instead?
- Two unbiased estimates of the same target have variances \(1\) and \(4\), and covariance \(0\). Calculate their inverse-variance combined estimate if their values are \(8\) and \(10\), and give its standard error.
- An upstream parameter estimate has standard error \(0.20\), but a downstream analysis treats it as fixed. What uncertainty is at risk of being omitted?
- The same dataset is summarized twice with the same estimator, producing the same estimate and standard error. Why does an independence-based combination overstate precision?
Solutions. (1) No. The slopes have different exposure units and correspond to different inputs. Compare predictions for a shared regimen and endpoint, with uncertainty, or justify a transformation to a common metric. (2) The weights are \(1\) and \(0.25\), normalized to \(0.8\) and \(0.2\). The estimate is \(0.8(8)+0.2(10)=8.4\), and its variance is \(1/(1+0.25)=0.8\), so \(SE=\sqrt{0.8}=0.8944\) target units. (3) The contribution of uncertainty in that estimated parameter, and possibly its covariance with downstream estimates, is omitted; the resulting prediction can be too precise. (4) The summaries contain no new independent observations. Their errors are dependent because they come from the same data and estimator, so treating them as independent falsely adds information and reduces the calculated standard error.