Why Renal Function Needs a Time Index
A renal-function value recorded at study entry is not necessarily a suitable description of renal function throughout follow-up. Acute illness, recovery, fluid shifts, and changes in treatment can accompany changes in kidney function. If clearance depends on renal function, a model that uses only the baseline value may therefore give systematically inaccurate predictions for later concentrations. The key modelling decision is not merely which renal metric to use, but when each value is assumed to apply.
A time-varying covariate is a covariate whose value for an individual can change during observation. It is useful to distinguish the underlying renal function at time \(t\), the times at which it is measured, and the rule used to represent it between measurements. These are related but not identical: a laboratory result is an observation of renal function at a particular time, while the model needs a covariate value over the dosing and sampling timeline.
The renal metric and its units should be stated explicitly. For example, an estimated filtration measure might be recorded in mL/min or normalized to body surface area. A value in one unit convention cannot be substituted into a model parameterized for another without the appropriate conversion. The covariate model also needs a reference value \(R_{\mathrm{ref}}\) in the same units.
Putting a Changing Covariate Into Clearance
Earlier material on renal and nonrenal contributions describes clearance models in which the renal component varies with renal function while the nonrenal component remains unchanged. One useful representation is
Here \(CL_{\mathrm{ref}}\) is clearance at the reference renal function for an individual with \(\eta_i=0\), \(\eta_i\) is the individual random effect on clearance, \(f_R\) is the renal fraction of reference clearance, and \(\theta\) describes the renal-function relationship. At \(R_i(t)=R_{\mathrm{ref}}\), the bracket equals one. For \(0\leq f_R<1\) and \(\theta>0\), the bracket tends to \(1-f_R>0\) as renal function approaches zero, representing a nonrenal component. At \(f_R=1\), no nonrenal component remains. Other supported covariate forms may be used; the essential time-varying feature is that the model evaluates the renal relationship using \(R_i(t)\), not automatically \(R_{i,0}\).
The covariate value enters the concentration-time model at the time it applies. In a one-compartment model with constant volume \(V_i\), no input over the interval, and time-varying clearance, the amount \(A_i(t)\) satisfies \(dA_i(t)/dt=-CL_i(t)A_i(t)/V_i\). Solving this differential equation gives
Thus, concentration depends on the accumulated clearance over the interval, not just clearance at the endpoint. With dosing or another input, the input is included in the differential equation and the model is solved across each dose and covariate-update time. In a multicompartment model, the time-varying clearance enters the corresponding system of amount equations.
Representing Renal Function Between Measurements
A population PK model typically receives renal-function values at discrete times, whereas its concentration predictions evolve continuously. The analysis must therefore define how observed values determine \(R_i(t)\) between measurements. A common choice is a piecewise-constant covariate: each measured value applies from its specimen-collection time until the next measurement. This is a modelling assumption, not a claim that biological renal function changes instantaneously at each blood draw.
Linear interpolation is another possible representation when gradual change between measurements is scientifically plausible and the data support it. The choice can matter when renal function changes quickly, when measurements are far apart, or when the covariate-clearance relationship is nonlinear. A result becoming available in the clinical record is not necessarily the time at which the biological measurement applies: the specimen-collection time is generally the relevant time for assigning the observed value.
Worked Example: Clearance Changes During a 24-Hour Interval
Consider a hypothetical individual with \(V_i=30.0\) L and \(C_i(0)=8.00\) mg/L. The renal-function covariate is measured as 80 mL/min at time 0 and 50 mL/min at hour 12. For this individual, suppose the fitted covariate relationship gives clearance of 5.00 L/h at 80 mL/min and 3.378414 L/h at 50 mL/min. Assume the measured value applies from its collection time until the next measurement, so clearance is piecewise constant: 5.00 L/h during hours 0–12 and 3.378414 L/h during hours 12–24. There is no dose or other input during the interval.
The accumulated clearance is the sum of clearance multiplied by the duration of each segment:
Applying the one-compartment concentration equation,
The calculation checks the timing assumption directly: the earlier renal value governs the first 12 hours, and the later value governs the next 12 hours. If the later clearance of 3.378414 L/h were incorrectly applied to the full 24 hours, the prediction would instead be
That alternative is higher because it assumes lower clearance for the entire interval, including the first 12 hours when the model assigns higher clearance. It is not a harmless change of bookkeeping; it changes the predicted concentration.
Measurement Timing and Update Rules
The preceding calculation follows a general result. If clearance is constant within each of \(K\) successive intervals, with clearance \(CL_{i,k}\) applied for duration \(\Delta t_k\), then the total clearance accumulated by the end of the intervals is \(\sum_{k=1}^K CL_{i,k}\Delta t_k\). For a one-compartment model without input,
This follows by splitting the integral of clearance into the intervals on which it is constant: \(\int_0^{t_K}CL_i(u)\,du=\sum_{k=1}^K\int_{t_{k-1}}^{t_k}CL_{i,k}\,du=\sum_{k=1}^KCL_{i,k}\Delta t_k\). The formula makes the timing consequence explicit: changing the time assigned to a renal measurement changes the duration \(\Delta t_k\) over which its associated clearance applies.
Worked Example: Updating Clearance at Specimen Times
In a hypothetical one-compartment example, let \(V_i=25.0\) L and \(C_i(0)=12.0\) mg/L. Suppose the covariate model is \(CL_i(t)=1.50+2.50[R_i(t)/60]\) L/h, where \(R_i(t)\) is in mL/min. Specimens collected at hours 0, 8, and 20 give renal-function values of 60, 90, and 45 mL/min, respectively. Use a piecewise-constant rule: each value applies from its collection time through the next collection time, and the hour-20 value applies to hour 24.
The corresponding clearances are \(1.50+2.50(60/60)=4.00\) L/h for hours 0–8, \(1.50+2.50(90/60)=5.25\) L/h for hours 8–20, and \(1.50+2.50(45/60)=3.375\) L/h for hours 20–24. Their durations are 8, 12, and 4 hours, so
Therefore,
For comparison, using only the baseline renal value of 60 mL/min gives \(CL_i=4.00\) L/h for all 24 hours. The resulting concentration is \(12.0\exp[-(4.00)(24)/25.0]=12.0\exp(-3.84)\approx0.2579\) mg/L. The difference reflects the two later covariate values and their assigned durations. Neither prediction should be interpreted as definitive without considering whether the update rule and model assumptions are reasonable.
Nonlinear Covariate Effects and Averaging
When the clearance relationship is nonlinear in renal function, predicting clearance from an average renal-function value need not equal averaging the clearances predicted from the changing values. For a power relationship, in general, the average of \(R_i(t)^\theta\) is not the \(\theta\)th power of the average of \(R_i(t)\). Consequently, substituting a single mean value for a changing covariate can alter the accumulated clearance even if that mean is calculated correctly.
Worked Example: Why the Mean Renal Value May Mislead
Suppose a hypothetical individual has \(V_i=40.0\) L, \(C_i(0)=10.0\) mg/L, and clearance model \(CL_i(t)=2.00+3.00[R_i(t)/60]^{1.2}\) L/h. Renal function is 30 mL/min for 12 hours and 90 mL/min for the next 12 hours. The power terms are \((30/60)^{1.2}\approx0.4353\) and \((90/60)^{1.2}\approx1.6267\). Thus the two clearances are approximately \(2.00+3.00(0.4353)=3.306\) L/h and \(2.00+3.00(1.6267)=6.880\) L/h.
The accumulated clearance and resulting concentration are
The arithmetic mean renal function is 60 mL/min. Substituting that mean for the full 24 hours gives clearance \(2.00+3.00(60/60)^{1.2}=5.00\) L/h and concentration \(10.0\exp[-(5.00)(24)/40.0]\approx0.4979\) mg/L. The two predictions differ because the covariate relationship is nonlinear: the average of the two clearance predictions is about 5.093 L/h, not 5.00 L/h. This example does not establish that one update rule is universally preferable; it shows why the chosen time representation should be evaluated rather than assumed inconsequential.
Diagnostics, Alternatives, and Limitations
A time-varying covariate is useful only if its timing and relationship to clearance are adequately supported by the data. Review the individual renal-function records against dose times, concentration samples, and relevant clinical events. Look for implausible jumps, inconsistent units, duplicate or misordered specimen times, and long intervals with no measurement. Compare predictions under defensible alternatives, such as baseline-only versus time-varying values or piecewise-constant versus interpolated values. Large changes in predicted concentrations indicate sensitivity to assumptions that should be reported and investigated.
These approaches answer somewhat different questions. A baseline-only covariate is simple and may be appropriate when renal function is stable or repeated measurements are sparse, but it cannot represent within-person change. A categorical occasion effect can capture differences between occasions, but it does not describe a continuous renal-function relationship unless that relationship is explicitly modelled. A piecewise-constant covariate is transparent and easy to align with specimen times, but it treats the assigned value as unchanged until the next measurement. Interpolation represents gradual change between observations, but may imply unsupported precision or a trajectory that was never observed.
Important limitations remain even when the time index is handled correctly. Renal-function measurements are imperfect, may be delayed relative to the underlying physiological change, and may be collected more often when a patient is unwell. A covariate-clearance association in an observational data set does not by itself establish that a change in renal function caused the entire change in clearance. Other time-varying factors, model misspecification, sparse concentration sampling, and uncertainty in the renal metric can also affect estimates. The model should not extrapolate a fitted relationship far beyond the renal-function range that informed it without explicit justification.
Check Your Understanding
Answer each question, then compare your reasoning with the solutions below.
- In a model that uses \(R_{i,0}\) for all times, what distinguishes the renal-function covariate from one that uses \(R_i(t)\)?
- Clearance is 4.0 L/h for 6 hours and 3.0 L/h for the following 6 hours. What is the accumulated clearance over 12 hours?
- Why can clearance predicted at mean renal function differ from the time-average predicted clearance when the covariate effect is a power function?
- A renal-function result is reported at hour 10, but the specimen was collected at hour 8. Under a specimen-time update rule, when does the measured value begin to apply?
Solutions. (1) The baseline model holds the designated initial value \(R_{i,0}\) fixed throughout; the time-varying model allows the assigned value \(R_i(t)\) to change. (2) The accumulated clearance is \((4.0)(6)+(3.0)(6)=42\) L. (3) A nonlinear transformation generally does not commute with averaging, so the average of transformed values need not equal the transformation of the average. (4) The value begins to apply at hour 8 under that rule, because the specimen-collection time, not the later reporting time, defines when the observed value pertains.