Separating Lag Time From Slow Absorption

Why the Onset Pattern Matters

A concentration profile that rises later than expected can reflect different input processes. With a lag time, drug input is delayed: under the model, no drug enters the systemic circulation before the lag ends. With slow first-order absorption, input begins immediately after dosing but accumulates gradually because the absorption rate constant $k_a$ is small. If measurements begin only after concentrations have risen, these explanations can produce profiles that look similar. Early observations are valuable because the models make different predictions about when input begins and how quickly concentration initially rises.

The aim is not to identify a mechanism from one low measurement. It is to determine whether the data, under a stated model and error structure, support a period of delayed input, gradual input from time zero, or insufficient information to distinguish them.

Model and assumptions: Consider a single oral dose $D$, a one-compartment disposition model with first-order elimination rate constant $k$, and first-order absorption with rate constant $k_a$. Let $T_L$ be the lag time and $A=FD/V$, where $F$ is bioavailability and $V$ is apparent volume of distribution. Assume linear disposition, complete first-order absorption of the bioavailable fraction, and fixed $k$ when comparing the absorption profiles. For a lag-time model, $$ C(t)= \begin{cases} 0, & t\leq T_L,\\[4pt] A\frac{k_a}{k_a-k}\left(e^{-k(t-T_L)}-e^{-k_a(t-T_L)}\right), & t>T_L. \end{cases} $$ For slow absorption without a lag, set $T_L=0$ and use a relatively small $k_a$. The formula assumes $k_a\ne k$; when $k_a=k$, its continuous limit is $C(t)=Ak_a(t-T_L)e^{-k_a(t-T_L)}$ after the lag.

What the Models Predict at Early Times

Before $T_L$, the lag model predicts exactly zero concentration from the dose. Just after the lag, concentration rises from zero. In contrast, the no-lag model predicts positive concentration for every $t>0$ under this idealized input model, although it may be very small and difficult to measure. Thus the distinction is not that slow absorption produces an immediate large concentration; it is that its input starts at time zero rather than after a period of no input.

The early rise also depends on $k_a$. Writing $\tau=t-T_L$ for time since input begins, expansion of the exponentials around $\tau=0$ gives the initial slope and its leading curvature:

e^{-k\tau}-e^{-k_a\tau} =(k_a-k)\tau-\frac{(k_a-k)(k_a+k)}{2}\tau^2+O(\tau^3), \qquad C(t)=A k_a\tau-\frac{A k_a(k_a+k)}{2}\tau^2+O(\tau^3).

Therefore, immediately after the lag, the slope is $A k_a$. A smaller $k_a$ gives a shallower initial rise if $A$ is held fixed. A lag instead shifts the start of that rise to $T_L$. Since $A$ may also be estimated, an early slope alone does not uniquely determine $k_a$: a change in $A$ can partly compensate for a change in $k_a$. Timing and multiple observations are important.

Worked Example: Comparing Predictions Before and After a Lag

Suppose, hypothetically, that $k=0.20\ \mathrm{h}^{-1}$ and $A=10\ \mathrm{mg/L}$ are known. Compare a lag model with $T_L=0.50\ \mathrm{h}$, $k_a=1.0\ \mathrm{h}^{-1}$, with a no-lag slow-absorption model having $k_a=0.25\ \mathrm{h}^{-1}$.

For the lag model, the coefficient is $10(1.0)/(1.0-0.20)=12.5\ \mathrm{mg/L}$. At $t=0.25$ and $0.50\ \mathrm{h}$, the predicted concentration is zero. At $t=0.75\ \mathrm{h}$, time since input begins is $\tau=0.25\ \mathrm{h}$, so $$ C(0.75)=12.5(e^{-0.20(0.25)}-e^{-1.0(0.25)}) =12.5(0.9512-0.7788)=2.155\ \mathrm{mg/L}. $$ At $t=1.00\ \mathrm{h}$, $\tau=0.50\ \mathrm{h}$, giving $$ C(1.00)=12.5(e^{-0.10}-e^{-0.50}) =12.5(0.9048-0.6065)=3.729\ \mathrm{mg/L}. $$

For the no-lag model, the coefficient is $10(0.25)/(0.25-0.20)=50\ \mathrm{mg/L}$. The corresponding predictions at $t=0.25,\ 0.50,\ 0.75,$ and $1.00\ \mathrm{h}$ are $0.591,\ 1.117,\ 1.584,$ and $1.997\ \mathrm{mg/L}$, respectively. For example, $$ C(0.25)=50(e^{-0.20(0.25)}-e^{-0.25(0.25)}) =50(0.9512-0.9394)=0.591\ \mathrm{mg/L}. $$

The key contrast is the prediction before $0.50\ \mathrm{h}$: zero for the lag model, but positive for slow absorption. The later values also differ because the chosen $k_a$ values differ. This comparison is conditional on the stated $A$ and $k$; estimating them from data could change the fitted profiles.

Comparing Profiles With Observations

A fair comparison starts by stating what is held common across the candidate models. For example, one may use the same disposition model and residual-error model, estimate or fix $k$ using suitable information, and compare a lag-time absorption profile with a no-lag profile in which $k_a$ is allowed to be small. In a nonlinear mixed-effects analysis, the comparison concerns the population model and its individual predictions, with the specified random-effects and residual-error structures. A single subject's profile can illustrate the shape contrast, but it does not by itself establish a population-level lag.

An objective profile can make parameter trade-offs visible. To profile $T_L$, fix it at a sequence of values and re-estimate the other parameters at each value; record the minimized objective, such as the negative twice log-likelihood. A similar profile can be formed over $k_a$. A flat region or a long diagonal valley in the joint $T_L$-versus-$k_a$ objective surface indicates that multiple combinations fit nearly as well. This is evidence of weak practical identifiability, not proof that the parameter is exactly unidentifiable.

Worked Example: Early Measurements Favor a Lag Under a Simple Error Model

Consider hypothetical concentrations at $0.25,\ 0.50,\ 0.75,$ and $1.00\ \mathrm{h}$: $0.05,\ 0.12,\ 2.00,$ and $3.40\ \mathrm{mg/L}$. Use the same fixed parameters as the preceding example for each candidate profile. To illustrate the comparison, assume independent additive Gaussian measurement errors with a common known variance. Under these assumptions, comparing the sum of squared residuals (SSE) compares the likelihoods because the variance and number of observations are the same.

The lag-model predictions are $0,\ 0,\ 2.155,$ and $3.729\ \mathrm{mg/L}$. Its residuals, defined as observed minus predicted, are $0.05,\ 0.12,\ -0.155,$ and $-0.329\ \mathrm{mg/L}$. Thus $$ \mathrm{SSE}_{\mathrm{lag}} =(0.05)^2+(0.12)^2+(-0.155)^2+(-0.329)^2 \approx 0.149\ (\mathrm{mg/L})^2. $$

The slow-absorption predictions are $0.591,\ 1.117,\ 1.584,$ and $1.997\ \mathrm{mg/L}$. The residuals are $-0.541,\ -0.997,\ 0.416,$ and $1.403\ \mathrm{mg/L}$, so $$ \mathrm{SSE}_{\mathrm{slow}} =(-0.541)^2+(-0.997)^2+(0.416)^2+(1.403)^2 \approx 3.43\ (\mathrm{mg/L})^2. $$

For these fixed candidate profiles and this assumed error model, the lag profile fits the observations better, chiefly because the first two measurements are close to zero while the slow model predicts appreciable concentrations. This is an illustration, not a general test: with parameters estimated from the observations, the objective values could change. A real comparison must use an appropriate error model and account for the number of estimated parameters.

If candidate models have different parameter counts or are not nested, a raw reduction in residual sum of squares is not, by itself, a calibrated significance test. Likelihood-based criteria or predictive checks may be appropriate, depending on the analysis and assumptions. In particular, do not automatically treat the difference in objective values as having a chi-squared distribution. The usual likelihood-ratio reference distribution requires specific conditions, including an appropriate nesting structure and regularity; those conditions need not hold for competing lag and slow-absorption formulations.

Why Sparse Early Sampling Creates Ambiguity

The initial expansion shows that at very early times a no-lag profile is approximately $A k_a t$. If only such early observations are available, different combinations of $A$ and $k_a$ can give similar concentrations. A lag adds another way to shift the observed rise: a concentration measured after the lag has been active for only a short time reflects $t-T_L$, not the full time since dosing. If sampling misses the period around onset, the data may provide little information about which feature produced the observed rise.

Worked Example: Early Slope Confounding and the Value of Timing

Suppose hypothetically that $k=0.20\ \mathrm{h}^{-1}$. Compare two no-lag profiles: profile A has $A=10\ \mathrm{mg/L}$ and $k_a=0.20\ \mathrm{h}^{-1}$; profile B has $A=20\ \mathrm{mg/L}$ and $k_a=0.10\ \mathrm{h}^{-1}$. The initial slopes are $A k_a=2.0\ \mathrm{mg/L/h}$ for both.

For profile A, $k_a=k$, so use the limiting formula $C(t)=A k_a t e^{-k_a t}$. At $t=0.10\ \mathrm{h}$, $$ C_A(0.10)=10(0.20)(0.10)e^{-0.020} =0.1960\ \mathrm{mg/L}. $$ For profile B, the coefficient $A k_a/(k_a-k)$ is $20(0.10)/(0.10-0.20)=-20\ \mathrm{mg/L}$. Therefore, $$ C_B(0.10)=-20(e^{-0.020}-e^{-0.010}) =0.1970\ \mathrm{mg/L}. $$ At $0.20\ \mathrm{h}$, the corresponding values are $C_A=10(0.20)(0.20)e^{-0.040}=0.3843\ \mathrm{mg/L}$ and $C_B=-20(e^{-0.040}-e^{-0.020})=0.3882\ \mathrm{mg/L}$. These profiles are very similar at these early times despite their different $A$ and $k_a$.

Now compare profile B with a lag version having the same $A$, $k_a$, and $k$, but $T_L=0.10\ \mathrm{h}$. At $t=0.10\ \mathrm{h}$, that model predicts zero; at $t=0.20\ \mathrm{h}$, it predicts $C_B(0.10)=0.1970\ \mathrm{mg/L}$, because only $0.10\ \mathrm{h}$ has elapsed since the lag ended. A well-measured sample at $0.10\ \mathrm{h}$ could distinguish these particular predictions, whereas observations only after the onset may leave more room for parameter trade-offs. The example demonstrates why sampling time matters, not that these parameter values are generally interchangeable.

Diagnostics, Limitations, and Common Mistakes

Check the fitted predictions against the actual early observations, not just a smooth line through later concentrations. Useful checks include plotting observed and predicted concentrations on a time scale that makes the onset visible, examining residuals around the first measurable concentrations, and inspecting the objective profile over $T_L$ and $k_a$. In a mixed-effects fit, compare individual predictions as well as the population prediction: averaging subjects with different onset patterns can obscure the shape seen in any one subject.

Interpret early measurements in light of the assay and residual model. A model prediction of zero before a lag does not imply that every recorded measurement must equal zero: measurement error, baseline signal, or other sources of concentration may produce nonzero observations. Conversely, a result below the limit of quantification does not establish an exact zero. The likelihood or observation model must handle such measurements appropriately; taking logarithms of a zero prediction is not valid. If early concentrations are too small to measure reliably, an apparent lag may simply reflect limited sensitivity.

The comparison also depends on the assumptions. If the disposition rate $k$ is poorly informed, uncertainty in disposition can affect the absorption profile. If $A$ is uncertain, it can trade off with $k_a$, especially early in the profile. A fitted lag should therefore be reported with its uncertainty and interpreted in the context of parameter correlations, sampling coverage, assay limits, and the assumed structural and residual models. A concentration pattern alone does not prove a particular physiological cause.

A common mistake is to say, “The first observed concentration is low, so there must be a lag.” A low early concentration can also result from slow first-order absorption, a small $A$, measurement error, or an observation taken before concentrations are reliably quantifiable. The correction is to compare the explicit predictions of both models at the actual sampling times and examine whether the data contain observations capable of distinguishing their onset patterns.

The main practical alternative to relying on one fitted profile is to improve information about the onset: obtain adequately timed early samples when feasible, use appropriate information about disposition, and assess the fitted profile across a range of plausible $T_L$ and $k_a$ values. If the profiles remain similar over that range, report the distinction as uncertain rather than selecting a mechanism solely because one optimizer returned a particular estimate.

Key takeaway: A lag model predicts no input before $T_L$; a no-lag first-order model predicts input from time zero, with an early slope proportional to $A k_a$. Early observations can distinguish these patterns only when their timing, measurement quality, and model assumptions provide enough information.

Check Your Understanding

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

Under the stated model, what does the lag-time profile predict for concentration before $T_L$, and how does that differ from slow absorption without a lag?
For $A=8\ \mathrm{mg/L}$ and $k_a=0.10\ \mathrm{h}^{-1}$, what is the approximate initial slope after input begins?
Why can two early no-lag profiles with different values of $A$ and $k_a$ be similar?
A fitted lag model predicts zero before $T_L$, but a measured concentration there is small and positive. Does this alone rule out the lag model? Explain.
What does a flat profile of the objective over several lag-time values suggest?

Solutions. (1) It predicts zero before $T_L$; a no-lag first-order model predicts positive, potentially very small, concentrations for $t>0$. (2) The initial slope is $A k_a=(8\ \mathrm{mg/L})(0.10\ \mathrm{h}^{-1})=0.80\ \mathrm{mg/L/h}$. (3) Their early concentrations are approximately $A k_a t$, so different parameter combinations with similar products $A k_a$ can give similar predictions. (4) No. Measurement error, baseline signal, or other sources can produce a positive observation; assess it using the specified observation model and the full profile. (5) The data may not determine $T_L$ precisely under the fitted model, indicating weak practical identifiability over that range.