Tutorials › Nonlinear Mixed-Effects Models (NLME) › Diagnosing Absorption Misspecification From Early-Time Patterns

Goodness-of-Fit Diagnostics and Their Pitfalls · Tutorial 565 of 797

Diagnosing Absorption Misspecification From Early-Time Patterns

Learn how early post-dose residual patterns and individual overlays can suggest absorption misspecification, and what checks help separate a model problem from timing or assay artifacts.

Intermediate 12 min read

What You'll Learn

  • Interpret early residual signs relative to a model’s predicted concentration-time curve.
  • Recognize patterns consistent with delayed, rapid, or prolonged absorption.
  • Use individual overlays alongside residuals without treating either as conclusive.
  • Check whether recorded dose and sampling times could explain an apparent pattern.
  • Consider assay and residual-error effects before changing the structural model.
  • Understand why early-time plots suggest, but do not prove, absorption misspecification.

Why Early-Time Patterns Matter

After an oral dose, absorption shapes the rise in concentration and the timing and sharpness of the peak. A model that describes absorption poorly can therefore show systematic disagreement soon after dosing even when its later predictions look reasonable. Early-time residuals and individual concentration overlays help locate that disagreement in time and suggest whether the model’s absorption is too delayed, too rapid, or insufficiently prolonged.

These plots are diagnostic evidence, not a stand-alone test of a particular absorption mechanism. A similar pattern can arise if dose or sample times were recorded inaccurately, if early samples were handled differently by an assay, or if the residual-error model gives the wrong scale to those observations. The useful question is not simply whether early residuals are large, but whether their direction and timing are coherent across individuals and remain after plausible data and measurement explanations have been checked.

Definition: Let \(y_{ij}\) be the observed concentration for individual \(i\) at sample time \(t_{ij}\), and let \(\hat{y}_{ij}\) be the corresponding model prediction. The raw residual is \(e_{ij}=y_{ij}-\hat{y}_{ij}\). A positive residual means the observation is above the prediction; a negative residual means it is below. A standardized residual divides this difference by an appropriate model-based residual scale. In this tutorial, “early-time” means the post-dose portion of the profile in which absorption materially shapes the concentration curve, not one universal fixed time interval.

Reading the Sign and Timing of Residuals

Suppose the fitted model predicts a concentration curve \(C(t)\), but the actual curve is shifted later by a small time \(\delta>0\). The delayed curve can be represented approximately as \(C(t-\delta)\). A first-order Taylor expansion gives:

$$ C(t-\delta)-C(t)\approx -\delta C'(t). $$

On the rising part of the predicted curve, \(C'(t)>0\), so a delayed observed curve tends to lie below the prediction: residuals are negative. After the peak, while the curve is falling, \(C'(t)<0\), so the delayed curve tends to lie above the prediction: residuals are positive. A model that absorbs too rapidly relative to the data can show the reverse sequence: positive residuals during the rise and negative residuals after the model has peaked too soon.

This argument describes a time shift, not every possible difference in absorption. A model that absorbs too briefly may produce a profile that rises and peaks too narrowly, while data with more prolonged input can stay below that prediction early and above it later. The exact pattern depends on the dose, sampling times, structural model, and other parameters. It is therefore more reliable to describe what the residuals do than to label a mechanism from a single point or plot.

Diagnostic principle: An absorption interpretation is more credible when the residual direction changes in a way that aligns with the observed and predicted curve shapes, appears in multiple individuals, and is supported by individual overlays. A lone large early residual, or a pattern confined to one assay batch or one sampling time, is weak evidence for structural absorption misfit.

Worked Example: A Pattern Consistent With Delayed Absorption

Consider a hypothetical individual with samples after one oral dose. The fitted model predicts a rising concentration through 2 hours and a decline by 4 hours. The table gives concentrations in mg/L. For illustration, suppose the residual scale is 1.5 mg/L at each time, so the standardized residual here is \(e_{ij}/1.5\).

Time (h)Observed (mg/L)Predicted (mg/L)Raw residual (mg/L)Standardized residual
0.524-2-1.33
168-2-1.33
2111010.667
49721.33

At 0.5 hours, \(e=2-4=-2\) mg/L and the standardized residual is \(-2/1.5=-1.33\). At 4 hours, \(e=9-7=2\) mg/L and the standardized residual is \(2/1.5=1.33\). The residuals move from negative during the rise to positive after the predicted peak. This is consistent with delayed absorption, but the example is only one individual: a wrong sample time or an early assay bias could also contribute. The next check is whether the same phase-related pattern recurs across individuals and whether their observed overlays show a later rise or peak than the model predicts.

Individual Overlays and the Population Pattern

An individual overlay places that person’s observations and predicted concentration-time curve on the same time axis. For this diagnosis, the plot should make dose time, actual sample times, and the early rising and peak regions easy to see. Where useful, display the population prediction and individual prediction distinctly: they answer different questions. A population prediction shows the typical model profile, whereas an individual prediction incorporates the person-specific parameter estimates used by the fitted model.

Individual predictions can help reveal whether the model can describe the observed profile after accounting for estimated between-individual differences. But they are not independent confirmation: the observations helped estimate those individual parameters. A good-looking individual curve does not by itself rule out structural misspecification, particularly with sparse early samples or substantial parameter shrinkage. Compare overlays with residuals and inspect whether discrepancies recur in subjects with enough early observations to show the curve shape.

Worked Example: A Model That Absorbs Too Slowly

In a second hypothetical profile, the fitted model predicts a slow rise. Assume the model-based residual scale is 2 mg/L at each listed time. The observations are unusually high early and then fall below the prediction after the model’s later peak.

Time (h)Observed (mg/L)Predicted (mg/L)Raw residual (mg/L)Standardized residual
0.55231.50
19542.00
210910.500
468-2-1.00

For example, at 1 hour, \(e=9-5=4\) mg/L and \(e/2=2.00\). At 4 hours, \(e=6-8=-2\) mg/L and \(e/2=-1.00\). The early positive and later negative residuals are compatible with the model absorbing too slowly: the observed profile rises sooner than the prediction. If this sequence recurs across individuals, it supports checking the absorption component. It does not establish a specific absorption rate or mechanism; sparse sampling may not locate the peak precisely enough for that conclusion.

Distinguishing Absorption Misfit From Timing and Assay Effects

Before changing the structural model, verify the information that determines where each observation belongs on the time axis. Check that dose times, sample times, time zones or clock conventions, and dose-to-sample intervals are correct. If a sample was recorded later than it was actually collected, plotting it at the recorded time can create an apparent phase shift. A dose-time error can affect every post-dose time for that occasion. If the pattern disappears after correcting a documented timing error, it was not evidence that the absorption model needed to change.

Next, look for clustering by assay, analytical run, site, or study. An additive assay offset can move early observations in one direction without changing the underlying profile shape. A scale difference can make discrepancies increase with concentration. If early samples were disproportionately measured in one batch, an assay-related pattern may mimic time-dependent structural misfit. Earlier material on residuals by study, site, and assay and residual diagnostics under proportional and combined error provides relevant checks: examine whether the trend remains after accounting for measurement scale and grouping.

Worked Example: An Early Assay Shift That Mimics a Profile Pattern

Suppose a hypothetical assay run measures early samples about 1.2 mg/L higher than the reference method, while later samples use the reference method. The model predicts 4 mg/L at 0.5 hours and 7 mg/L at 1 hour; the corresponding observations from the shifted run are 5.2 and 8.2 mg/L.

At 0.5 hours, \(e=5.2-4=1.2\) mg/L. At 1 hour, \(e=8.2-7=1.2\) mg/L. The residuals are positive at both early times, but their shared 1.2 mg/L offset matches the stated assay shift. That pattern does not, by itself, indicate faster absorption. Check assay and run labels and compare observations measured by different runs at similar times, if such data exist. If the assay shift is not supported, this explanation should not simply be assumed; the example illustrates why assay information matters before a structural change is made.

When Absorption Is Too Brief in the Model

A model may also fail because its input ends too quickly relative to the data. In that case, observed concentrations can remain lower than the fitted curve during the early rise and higher than it later, as the observed profile continues to receive input after the model’s predicted rise has largely passed. The transition may be gradual rather than a clean reversal at a single time. Individual overlays can show whether the observed rise or peak is broader than the prediction, while early and near-peak residuals indicate which parts of the curve are missed.

Worked Example: A Profile Suggesting More Prolonged Absorption

For a third hypothetical individual, suppose the model predicts a relatively narrow rise. The residual scale is 1 mg/L. The observations are 2, 6, 10, and 8 mg/L at 0.5, 1, 2, and 4 hours; predictions are 3, 7, 10, and 6 mg/L, respectively.

The residuals are \(2-3=-1\), \(6-7=-1\), \(10-10=0\), and \(8-6=2\) mg/L. Dividing by the 1 mg/L scale gives standardized residuals \(-1,-1,0,\) and \(2\). The observed profile is below prediction during the rise and remains above it at 4 hours. This is compatible with a broader or more prolonged input than the model predicts. However, the 4-hour point alone cannot establish prolonged absorption: other model components, a timing error, or measurement variability could produce it. Repeated profiles and adequate sampling around the rise and peak are needed to assess the shape.

Assumptions, Failure Modes, and Alternatives

These interpretations assume that the dose and sample times are reliable, observations are assigned to the correct dose and occasion, concentrations are measured on a consistent basis, and residuals are interpreted using a scale appropriate to the fitted residual-error model. They also assume that the model predictions being compared refer to the same individual, dose, and covariate conditions as the observations. When these conditions are doubtful, a sign pattern may describe data recording or measurement rather than absorption.

Useful diagnostics complement rather than replace one another. A residual-versus-time plot summarizes direction and spread across observations; individual overlays show whether those residuals correspond to a coherent mismatch in profile shape. A plot of residuals against dose or regimen can reveal that a pattern is confined to a particular dosing condition. Review of assay and sampling records targets nonstructural explanations. A more complex absorption model is an alternative to the current model, but it should be considered only after these checks and judged with the same diagnostic scales and relevant data, not selected simply because it can follow a few points more closely.

  • Too few early samples: the data may not distinguish a delayed rise from a different peak time. Do not infer a precise absorption process from an unobserved curve segment.
  • Unequal sampling: a pooled residual plot can be dominated by individuals with many samples. Check individual patterns and how observations are distributed over time.
  • Large individual differences: a population-level pattern may reflect a subgroup or a few influential profiles. Check whether it is broadly repeated, rather than relying on the pooled trend alone.
  • Residual scale effects: a concentration-dependent spread can make raw residual magnitudes misleading. Interpret residual direction and use the appropriate standardized scale, as discussed in residual diagnostics under proportional and combined error.
  • Confounding with other model features: early residuals do not uniquely identify absorption. The same observed concentration can be affected by other model components and by measurement error. Avoid using this early-time diagnosis to settle questions that require examining other portions of the profile.
Common pitfall: “The early residuals are positive, so absorption must be too fast.” This reverses the comparison without checking the prediction. Positive residuals mean observations exceed predictions. If they occur on the early rising limb and later turn negative, they are consistent with the model absorbing too slowly; if the sign sequence is reversed, overly rapid model absorption is one possibility. Always state which curve is above the other, at what phase, and whether timing or assay patterns could explain it.
Key takeaway: Early residual signs and individual overlays can indicate the direction and timing of absorption misfit, but they do not identify its cause uniquely. Interpret a repeated, phase-coherent pattern only after checking dose and sample times, assay groupings, residual scale, and the amount of information available around the rise and peak.

Check Your Understanding

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

  1. On the predicted rising limb, residuals are mostly negative and become positive after the predicted peak. What absorption difference is one plausible explanation?
  2. A model predicts 3 mg/L at 1 hour, and the observed value is 5 mg/L. What is the raw residual, and what does its sign mean?
  3. Positive residuals occur only for early samples from one assay run. Name one check that should come before changing the absorption model.
  4. Why should an individual overlay not be treated as independent proof that the absorption model is correct?
  5. There is only one sample before the peak for most individuals. What limitation does this create for diagnosing absorption?

Solutions. (1) Delayed absorption relative to the model is plausible: the observed curve is below prediction during the rise and above prediction later. (2) \(e=5-3=2\) mg/L; the observation is 2 mg/L above the prediction. (3) Check assay-run assignments and whether a run-specific offset or scale difference could explain the pattern. (4) The individual parameter estimates used in the overlay were estimated using that person’s observations, so the displayed fit is not independent confirmation. (5) The early curve shape and peak timing are weakly observed; a delayed rise, different peak shape, or random variation may be hard to distinguish.