Why an Integrated Event History Matters
A concentration prediction depends on more than the prescribed regimen. It depends on the administrations that actually occurred and on how the patient’s physiology changed between them. When a dose may have been skipped and renal function has been measured only intermittently, the analyst must make both histories explicit before calculating exposure. Otherwise, a mismatch between observations and predictions can be incorrectly attributed to the pharmacokinetic model.
This case uses hypothetical data and a deliberately simple one-compartment model. Its purpose is to show how to construct and compare internally consistent event histories, not to prescribe a dose or establish an individual patient’s true history from limited observations. The population model parameters are treated as fixed for the calculations; no parameter estimation is performed.
For a dose at time \(t_j\), the concentration contribution remaining at time \(t\geq t_j\) is the dose divided by volume, multiplied by the exponential of the negative cumulative elimination since administration. Contributions from all doses are added. The cumulative elimination is the integral of the time-varying rate constant; it is not generally equal to one constant rate multiplied by elapsed time.
This expression follows by solving \(dC/dt=-k(t)C\) between events: separating variables gives \(dC/C=-k(t)\,dt\), and integration from \(t_j\) to \(t\) gives the exponential survival factor. At a bolus event, concentration increases by \(D_j/V\). Thus the concentration immediately after the dose is the concentration immediately before it plus \(D_j/V\). An observation recorded at exactly the dose time must be identified as pre-dose or post-dose; the two values differ by this increment.
Worked Case: Missed Dose With Stepwise Renal Function
Suppose the hypothetical patient is prescribed 100 mg intravenous boluses every 12 hours, beginning at hour 0. The recorded administrations are 100 mg at hours 0, 12, 36, and 48. The dose scheduled for hour 24 is confirmed to have been missed. eGFR measurements are 60 mL/min/1.73 m² at hour 0, 40 at hour 24, and 30 at hour 48. Use the stepwise convention in the model definition: 60 applies from hour 0 up to 24, 40 from hour 24 up to 48, and 30 from hour 48 onward. The resulting elimination rates are 0.100, 0.06667, and 0.05000 h\(^{-1}\), respectively.
Worked Example: Predicting Concentration Under a Confirmed Missed Dose
Calculate the predicted concentration immediately after the hour-48 dose and at hour 60. Each 100 mg dose contributes \(100/20.0=5.00\) mg/L at the instant it is given.
At hour 48, before the dose: the dose at hour 0 has experienced 24 hours at \(k=0.100\) and 24 hours at \(k=0.06667\). Its cumulative elimination is \(2.40+1.60=4.00\), so its contribution is \(5.00e^{-4.00}=0.09158\) mg/L. The hour-12 dose has experienced 12 hours at 0.100 and 24 hours at 0.06667, giving \(5.00e^{-2.80}=0.3041\) mg/L. The hour-36 dose has experienced 12 hours at 0.06667, giving \(5.00e^{-0.8000}=2.247\) mg/L. Therefore, the pre-dose prediction is approximately \(2.642\) mg/L using unrounded contributions. The hour-48 dose adds 5.00 mg/L, so the post-dose concentration is approximately \(7.642\) mg/L.
At hour 60: the hour-0 contribution has additional elimination \(12(0.05000)=0.6000\), so its total cumulative elimination is 4.60 and its contribution is \(5.00e^{-4.60}=0.05018\) mg/L. For the hour-12 dose the total is \(1.20+1.60+0.6000=3.40\), giving \(5.00e^{-3.40}=0.1669\) mg/L. The hour-36 dose has total \(0.8000+0.6000=1.400\), giving \(5.00e^{-1.400}=1.233\) mg/L. The hour-48 dose contributes \(5.00e^{-0.6000}=2.744\) mg/L. Adding the four contributions gives \(C(60)\approx0.05018+0.1669+1.233+2.744=4.194\) mg/L.
The calculation makes the event history visible: there is no hour-24 dose term. A missing record and a confirmed missed administration are not interchangeable facts. If the record is uncertain rather than confirmed, it is more defensible to calculate a separate history that includes the possible administration than to silently insert a dose into the primary history.
Comparing Alternative Histories
An alternative history should change only the assumption under investigation. Here, one plausible alternative is that the hour-24 dose was actually administered but was absent from the record. Keep the same renal-function history and all other doses. At hour 60, this extra dose has experienced 24 hours at \(k=0.06667\) and 12 hours at \(k=0.05000\), for cumulative elimination \(1.60+0.6000=2.20\). Its predicted contribution is \(5.00e^{-2.20}=0.5540\) mg/L. The alternative-history prediction is therefore \(4.194+0.5540=4.748\) mg/L. Relative to the missed-dose prediction, this is about \(0.5540/4.194=13.21\%\) higher.
That difference is conditional on the stipulated dose, volume, and renal-effect model. It is not a probability that the dose was taken, nor does a higher prediction establish which history occurred. The value of the comparison is that it quantifies how much the uncertain event matters under explicit assumptions.
Worked Example: Checking Histories Against a Pre-Dose Observation
Suppose a hypothetical concentration measured just before the hour-48 administration is 2.70 mg/L. Compare it with each history’s prediction at that same event position. The confirmed-missed history predicts approximately 2.642 mg/L before the hour-48 dose. Its observed-minus-predicted residual is approximately \(0.058\) mg/L using the unrounded prediction. The history that includes an hour-24 dose predicts an additional \(5.00e^{-1.60}=1.009\) mg/L at hour 48, or approximately \(3.652\) mg/L before the hour-48 dose. Its residual is \(2.70-3.652=-0.952\) mg/L.
For this observation, the prediction under the missed-dose history is closer. That is evidence of relative predictive agreement for this measurement under the specified model, not definitive proof of the event history. Measurement error, uncertain parameters, model misspecification, or other unrecorded events could also explain the discrepancy. A formal NLME analysis would account for its specified observation-error and random-effects structure; this arithmetic comparison does not estimate those quantities.
Covariate History Is Also an Assumption
The preceding calculations use a stepwise eGFR history. A different, also explicit, choice is linear interpolation between the recorded measurements. Interpolation assumes the covariate changed smoothly between measurement times; it does not assert that the unobserved renal function is known. Since clearance is proportional to eGFR in this example, interpolation changes the integrated elimination and therefore changes the concentration contributions from earlier doses.
Worked Example: Effect of Interpolating eGFR
Keep the missed-dose history, but linearly interpolate eGFR from 60 at hour 0 to 40 at hour 24, then from 40 at hour 24 to 30 at hour 48. Hold it at 30 after hour 48. Since \(k=eGFR/600\), the rate declines linearly over each interval. From hour 0 to 24, the average rate is \((0.100+0.06667)/2=0.08333\) h\(^{-1}\), and integrated elimination is \(24(0.08333)=2.000\). From hour 24 to 48, the average is \((0.06667+0.05000)/2=0.05833\) h\(^{-1}\), giving \(24(0.05833)=1.400\). From hour 48 to 60, integrated elimination is \(12(0.05000)=0.6000\).
The hour-0 dose therefore has cumulative elimination 4.000 by hour 60 and contributes \(5.00e^{-4.000}=0.09158\) mg/L. For the hour-12 dose, the first-interval cumulative elimination from 12 to 24 is 0.9000: the linearly declining rate has average \((0.08333+0.06667)/2=0.07500\) over those 12 hours. Its cumulative elimination to hour 60 is \(0.9000+1.400+0.6000=2.900\), giving \(5.00e^{-2.900}=0.2751\) mg/L. The hour-36 dose experiences 12 hours in the second interpolated interval, with average rate \((0.05833+0.05000)/2=0.05417\), for 0.6500, then 0.6000 after hour 48. Its contribution is \(5.00e^{-1.250}=1.433\) mg/L. The hour-48 dose contributes \(5.00e^{-0.6000}=2.744\) mg/L. Thus the interpolated-history prediction is approximately \(0.09158+0.2751+1.433+2.744=4.543\) mg/L.
The stepwise prediction was 4.194 mg/L, about 0.349 mg/L lower at hour 60. The direction is consistent with the histories: stepwise carry-forward uses the higher eGFR—and thus higher clearance—until each new measurement, whereas interpolation lowers eGFR progressively. This is a comparison of assumptions, not a test that determines the unobserved renal-function path.
Exposure Over an Interval
Concentration at one time is not the only useful comparison. For an interval with constant elimination rate \(k\), concentration declines exponentially, so the area under the concentration-time curve (AUC) from \(t_0\) to \(t_1\) is \(C(t_0)\{1-e^{-k(t_1-t_0)}\}/k\). This follows by integrating \(C(t_0)e^{-k(t-t_0)}\) over the interval. The units are concentration multiplied by time, here mg·h/L.
Under the missed-dose, stepwise history, the post-dose concentration at hour 48 is approximately 7.642 mg/L. The rate from hour 48 to 60 is constant at 0.05000 h\(^{-1}\), so AUC over that interval is approximately \(7.642273(1-e^{-0.6000})/0.05000=68.96\) mg·h/L. If the hour-24 dose had instead been administered, it would add \(5.00e^{-1.60}=1.009483\) mg/L to the hour-48 pre-dose concentration and therefore to the post-dose concentration. The corresponding AUC increment from hour 48 to 60 is \(1.009483(1-e^{-0.6000})/0.05000=9.109\) mg·h/L. The alternative history’s AUC is thus about 78.07 mg·h/L over this interval. This AUC comparison, like the concentration comparison, is conditional on the same structural assumptions.
Diagnostics, Limitations, and Common Pitfalls
A reliable reconstruction begins with a timeline containing both administrations and observations. Check the event time, dose amount, units, and whether each concentration is pre-dose or post-dose. Confirm that each covariate value is applied over the intended interval and that cumulative elimination includes every portion of the interval since the relevant dose. A simple numerical check is that a dose’s contribution must decline as its cumulative elimination increases, and contributions from doses at the same time add linearly in this model.
- Do not convert an absent record automatically into a missed dose. Distinguish confirmed omission, uncertain administration, and missing documentation. Calculate alternative histories when the record cannot resolve the event.
- Do not apply a new covariate value before its measurement time without stating an interpolation rule. Stepwise carry-forward and interpolation are different assumptions and can produce different predictions.
- Do not compare a pre-dose observation with a post-dose prediction. At a bolus time, the model has a jump of \(D/V\); align observation and prediction on the same side of the event.
- Do not treat a closer prediction as proof of history. Residuals also reflect measurement error and model inadequacy. Repeated residual patterns, event-record review, and sensitivity to plausible histories are more informative than a single comparison.
- Do not overinterpret the renal relationship. The proportional eGFR-clearance relation and parameter values here are hypothetical. Real applications require an appropriate model, relevant units, and evaluation of whether renal measurements adequately represent the physiology affecting clearance.
An alternative to a single assumed history is a sensitivity analysis that reports predictions over plausible event and covariate histories. A richer model could represent gradual renal change, uncertainty in parameters, or random effects, but those additions do not remove the need to establish event ordering and define how covariates are handled between measurements. This workshop’s calculations isolate the impact of those histories while holding the structural parameters fixed.
Check Your Understanding
Answer each question using the model and hypothetical case above; concise solutions follow.
- Under the stepwise history, what elimination rate applies from hour 24 up to hour 48, and why?
- How much does a possible 100 mg dose at hour 24 contribute to concentration at hour 60 under the stepwise history?
- Why does the interpolated history predict a higher hour-60 concentration than the stepwise history in the worked calculation?
- What must be specified when an observation and a bolus dose share the same time?
- Does the closer pre-dose prediction in the worked comparison establish that the dose at hour 24 was missed? Explain.
Solutions. (1) The rate is 0.06667 h\(^{-1}\), because eGFR measured as 40 at hour 24 is carried forward until the hour-48 measurement under the stated stepwise convention. (2) Its remaining cumulative elimination is \(24(0.06667)+12(0.05000)=2.20\), so its contribution is \(5.00e^{-2.20}=0.5540\) mg/L. (3) Interpolation lowers eGFR progressively within each interval, reducing integrated elimination for several earlier doses and leaving more concentration at hour 60. (4) State whether the observation is pre-dose or post-dose, since concentration jumps by \(D/V\) at the bolus. (5) No. It shows better agreement for this single observation under the specified parameter values and model; measurement error, model limitations, and other uncertainties prevent treating it as proof of the event history.