Why the Cause of Disagreement Matters
When a QSP model is calibrated to multiple data sources, disagreement between observations and predictions can have at least two different explanations. The sources may measure the same biological quantity differently, or the model may fail to represent a biological pattern that appears across sources. These explanations imply different next steps: a measurement problem may call for calibration or assay review, whereas shared structured residuals may motivate reconsidering the mechanistic representation.
Residuals alone do not prove which explanation is correct. A source-specific pattern is evidence about measurement only when the study design and controls support that interpretation. Conversely, a pattern repeated across assays is consistent with shared model misfit, but can also arise from a measurement bias common to all assays. The objective is therefore to compare explicit predictions of the two explanations against controls, replicates, and matched experimental conditions.
Assumptions and a Diagnostic Decomposition
This distinction is most informative when sources have measurements at overlapping biological conditions, controls that identify how instruments respond, and replicates that characterize within-condition variability. The model prediction must correspond to the quantity actually measured, or a stated measurement mapping must connect the prediction to the observed scale. For example, a model may predict free concentration while an assay reports total concentration; treating those as interchangeable would manufacture apparent residual structure.
A simple diagnostic representation separates a source pattern from a condition pattern:
Here \(a_s\) represents a source-associated residual component, \(h_c\) a condition-associated component common across sources, and \(\varepsilon_{sc}\) remaining variation. This is a diagnostic decomposition, not automatically a validated biological model. A constraint such as \(\sum_s a_s=0\) is needed to define the reference level of the source components. If every source measures a different condition, \(a_s\) and \(h_c\) cannot be separated from the observations alone: source and condition are confounded.
A useful justification comes from comparing sources at the same condition. Suppose the model prediction is shared, and observations on the comparable scale satisfy \(y_{sc}=f_c(\theta)+a_s+\varepsilon_{sc}\). Subtracting observations from two sources cancels the common model prediction:
Thus a consistent contrast across matched conditions can point to a source difference without requiring the model prediction to be correct. But this cancellation is valid only if the sources measure the same underlying quantity under comparable conditions. With a source-specific scale transformation, the raw difference need not be constant; controls are then needed to estimate or assess the transformation.
Residual size should be judged relative to uncertainty, not just by visual distance. If an observation has standard error \(SE\), a simple standardized residual is \(r/SE\). For correlated observations, the residual vector’s covariance matters: treating repeated measurements from one sample as independent can make a pattern appear more precise than it is. Replicates estimate measurement variability under their experimental conditions; they do not by themselves establish that the mechanistic model is correct.
Use Controls to Test Measurement Explanations
A reference control has a known or independently characterized value, so it can reveal whether sources differ even when the biological model is not at issue. For a control value \(x\), a simple affine measurement model is \(y_s=a_s+b_s x+\epsilon\), where \(a_s\) is an offset, \(b_s\) a scale factor, and \(\epsilon\) measurement error. One control level can inform an offset only if the scale is already known; multiple levels are needed to assess a linear scale difference. A visibly curved control response would call the affine assumption into question.
Controls should span the range relevant to the biological samples and, where appropriate, use a matrix comparable to those samples. A control measured only at a low concentration may not reveal a source difference at high concentrations. Likewise, a reference material in buffer may not expose a matrix effect in plasma. Control-based adjustment also carries uncertainty: it does not turn noisy measurements into exact values.
Worked Example: A Control Reveals an Offset
Two hypothetical assay sources measure a reference material whose value is 10 concentration units. Source A reports 10.2, 10.0, and 9.8 units; source B reports 12.1, 11.9, and 12.0 units. Assess whether this control supports an offset difference.
The sample means are \(\bar{x}_A=(10.2+10.0+9.8)/3=10.0\) units and \(\bar{x}_B=(12.1+11.9+12.0)/3=12.0\) units. The estimated difference is \(12.0-10.0=2.0\) units.
For A, the squared deviations sum to \(0.2^2+0^2+(-0.2)^2=0.08\), so \(s_A=\sqrt{0.08/(3-1)}=0.200\) units. For B, the squared deviations sum to \(0.1^2+(-0.1)^2+0^2=0.02\), so \(s_B=\sqrt{0.02/2}=0.100\) units. If the replicate measurements are independent within and across sources, the standard error of the difference in means is
The observed two-unit difference is large relative to replicate variability, and the reference value is known to be 10 units. This supports a source-associated measurement offset for this control. It does not establish that the same offset applies to every biological matrix or concentration; that requires relevant controls and an adequate measurement model.
Look for Replicated Residual Structure
After considering the measurement process, inspect residuals by condition, time, dose, and source. A source pattern that is approximately stable across matched conditions, and is also visible in controls, is more consistent with measurement conflict. A residual pattern that changes with biological condition in a similar way across sources is more consistent with shared model misfit. These are evidential patterns, not rules that uniquely identify cause.
Replicate patterns help distinguish a stable shift from ordinary measurement scatter. If replicate values within a condition vary widely, a difference between source means may be poorly determined. If replicates cluster tightly but their mean repeatedly departs from the model prediction, random measurement noise alone is less persuasive as an explanation. However, shared sample handling, common reagents, or common preprocessing can induce correlated errors that replicate counts do not remove.
Worked Example: A Condition Pattern Shared Across Sources
Suppose a hypothetical model predicts biomarker values of 10, 8, and 5 units at three times. Two sources measure each time twice. Source A reports 10.0 and 10.0; 9.0 and 9.2; 7.0 and 7.2 units. Source B reports 10.1 and 9.9; 9.0 and 9.2; 7.1 and 6.9 units.
The source means at the three times are A: 10.0, 9.1, and 7.1 units; B: 10.0, 9.1, and 7.0 units. Averaging the two source means gives 10.0, 9.1, and 7.05 units. The corresponding residuals relative to predictions are \(10.0-10=0.00\), \(9.1-8=1.10\), and \(7.05-5=2.05\) units.
The source contrast in the means is \(0.0\), \(0.0\), and \(-0.1\) units, while the residual grows from zero to about two units over time. At the second time, each source has replicate deviations of \(0.1\) units from its mean; at the third time, each source’s deviations are also \(0.1\) units. The repeated increasing residual pattern is therefore not accompanied by a comparable source separation or large within-source replicate scatter. It is consistent with shared model misfit over time, such as a process missing from the model’s temporal response. It does not identify which mechanism is missing, and shared assay bias remains possible without suitable controls.
When a Measurement Scale Difference Is Not an Offset
Source disagreement need not be constant across the measurement range. An offset-only correction can leave residual structure if sources differ in scale. Controls at multiple known levels can reveal this distinction. The following example uses an exact hypothetical calibration pattern to make the arithmetic transparent; real control data would have uncertainty that should be carried into subsequent interpretation.
Worked Example: Control Data Reveal a Scale Difference
Assay A reports the reference concentrations 0, 10, and 20 units as 0, 10, and 20 units. Assay B reports those same reference values as 1, 12, and 23 units. Assume, provisionally, that each assay has a linear measurement relationship over this range.
For assay B, the change in reported value from the first to second control is \(12-1=11\) units for a 10-unit change in reference concentration. The estimated slope is \(11/10=1.1\). The intercept is 1 unit, since a reference value of zero is reported as 1. Thus the control relationship is \(y_B=1+1.1x\).
Suppose a biological sample is reported as 6.0 units by A and 7.6 units by B. Transform B’s result back to the reference scale: \(x=(y_B-1)/1.1=(7.6-1)/1.1=6.0\) units. The two results agree after this control-based mapping. An offset-only correction would subtract 1 from B’s result and give 6.6 units, leaving a 0.6-unit disagreement. The example supports a scale difference over the calibrated range, subject to the linearity assumption; it does not establish that the mapping is valid outside that range.
Practical Diagnostic Sequence
Check units, sample matrix, endpoint definition, and any mapping from the model-predicted quantity to the assay result. A mismatch here can mimic either source conflict or model failure.
Compare controls across sources and across the relevant measurement range. Summarize replicate variability and check whether replicates share samples, batches, or other sources of dependence.
For conditions measured by multiple sources, inspect source contrasts and model residuals separately. A consistent source contrast and a common condition trend are different patterns and can coexist.
If a source appears only at one condition, the data cannot distinguish a source effect from a condition effect without additional assumptions or external controls.
Report whether the evidence is more consistent with a measurement explanation, shared model misfit, both, or neither. Do not treat a residual pattern as proof of a specific cause.
Alternatives, Failure Modes, and Limits
A residual plot is a useful first check because it displays direction and structure, but it does not identify measurement behavior on its own. Control-based calibration directly investigates assay response, but depends on controls being relevant to the biological samples and on the chosen calibration form being adequate. Replicate summaries quantify repeatability, but cannot detect a bias shared by every replicate. Comparing source contrasts at matched conditions helps cancel a common model prediction, but fails if the sources measure different quantities or if their measurement mappings vary with condition.
A common mistake is to say, “Source B has larger residuals, so the biological model is wrong for source B.” This reasoning confuses the observed scale with the underlying mechanism. If controls show a source-specific offset or gain, the source pattern may arise before the model’s biological prediction is compared with the data. The correction is to inspect measurement behavior and matched conditions first, then assess whether residual structure remains across sources on a comparable scale.
The reverse mistake is also possible: applying a source correction until residuals look small and concluding that the model is adequate. A correction can conceal a shared biological discrepancy, especially when source and condition are confounded or when the correction is estimated from too few controls. Flexible corrections can also absorb real differences between experimental conditions. Preserve the distinction between what controls identify and what the mechanistic data support.
The strongest conclusions require crossed evidence: controls measured by multiple sources, biological conditions shared across sources, and replicates sufficient to assess variability. If these are absent, report the ambiguity and identify the additional measurement that would be most informative. A control spanning the disputed range can test a scale explanation; a second source measured at the discrepant biological condition can help separate source from condition effects.
Check Your Understanding
For each problem, state the relevant diagnostic reasoning. Solutions are included after the questions.
- Two sources measure the same control and the same biological condition. The control differs by 1.5 units between sources, and the biological sample also differs by about 1.5 units. What explanation is supported, and what limitation remains?
- A model underpredicts observations increasingly over time in two assays. Their matched-condition means are nearly equal, and controls show no obvious source difference. What pattern is present, and what can it not establish?
- Source A measures only low-dose samples and source B only high-dose samples. Residuals are higher for B. Can this design distinguish source bias from dose-related model misfit? Explain.
- Assay C has replicate standard deviation 0.1 units and assay D has replicate standard deviation 0.8 units. Why is comparing only their raw residual magnitudes potentially misleading?
Solutions. (1) The control and sample contrasts are consistent with a source-associated offset of about 1.5 units. The control supports that explanation for the control conditions, but does not prove the offset is identical in the biological matrix or across the measurement range. (2) The increasing residual is a condition-associated pattern shared across sources, consistent with shared model misfit. It does not identify a missing mechanism, and common measurement bias is not excluded solely by similar assays. (3) No. Source and dose are confounded, so a source effect and dose-related residual pattern cannot be separated from these data alone; overlapping doses or relevant controls are needed. (4) The same raw residual has different scale relative to measurement variability. Residuals should be interpreted with their standard errors and dependence structure, not by magnitude alone.