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.
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:
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.
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.5 | 2 | 4 | -2 | -1.33 |
| 1 | 6 | 8 | -2 | -1.33 |
| 2 | 11 | 10 | 1 | 0.667 |
| 4 | 9 | 7 | 2 | 1.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.5 | 5 | 2 | 3 | 1.50 |
| 1 | 9 | 5 | 4 | 2.00 |
| 2 | 10 | 9 | 1 | 0.500 |
| 4 | 6 | 8 | -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.
Check Your Understanding
Answer each question, then compare your reasoning with the solutions below.
- On the predicted rising limb, residuals are mostly negative and become positive after the predicted peak. What absorption difference is one plausible explanation?
- 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?
- Positive residuals occur only for early samples from one assay run. Name one check that should come before changing the absorption model.
- Why should an individual overlay not be treated as independent proof that the absorption model is correct?
- 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.