Trustworthiness Is a Decision-Specific Judgment
An individual forecast can be numerically precise yet unhelpful, or uncertain yet adequate for a low-stakes decision. Its trustworthiness therefore cannot be judged by the width of an interval alone. The relevant question is whether the forecast, given its data and assumptions, is dependable enough for the intended decision: for example, deciding whether a concentration is likely to fall outside a specified range or whether another measurement is needed.
Let \(y_i\) denote the observed concentrations for individual \(i\), \(D_i\) the recorded dosing and sampling history, and \(M\) the population PK model, including its structural, variability, and residual-error assumptions. Let \(\theta_i\) denote the individual parameters. A Bayesian forecast uses the posterior distribution \(p(\theta_i\mid y_i,D_i,M)\) together with the model to predict a future concentration or other quantity of interest. The assessment in this tutorial concerns whether that forecast is sufficiently supported for its stated use; it does not certify that the model is true.
A useful assessment keeps several questions distinct. Are the inputs accurate? Is there enough informative data to update the prior? How uncertain is the quantity that drives the decision? Do the observations agree reasonably with the prior and the model’s predictive implications? Are there clinical or design changes that make the historical model inappropriate for the forecast period? A reassuring answer to one question cannot compensate automatically for a serious problem in another.
What the Forecast Uncertainty Contains
A future concentration is uncertain for more than one reason. The individual parameters are uncertain after observing the available data; even if those parameters were known, residual variability can make an individual concentration differ from its model expectation. Uncertainty in dose amounts, dose times, adherence, or future dosing can contribute as well. The posterior predictive distribution integrates over these uncertainties that are represented in the model.
Here \(D_{\mathrm{future}}\) is the assumed future dose history. If that history is uncertain, a fuller prediction also averages over its plausible values. For a forecast conditional on a fixed future dose plan, the displayed distribution does not account for uncertainty that the plan will actually be followed. That distinction matters: a precise conditional forecast is not a guarantee about an uncertain real-world dosing course.
The law of total variance makes the two main sources of predictive spread explicit. Conditional on the individual parameters, the first term represents residual or process variability; the second represents variation in the model’s expected future value across plausible individual parameters.
This decomposition explains why a narrow posterior for a parameter does not necessarily imply a narrow prediction interval: residual variability may still be substantial. Conversely, a posterior mean forecast alone conceals how much probability lies on either side of a clinical threshold. The uncertainty measure must match the question. A prediction interval describes a range for a future observation under the model; a credible interval for a parameter describes uncertainty about that parameter conditional on the model and data. They are not interchangeable.
Worked Example: Is the Forecast Precise Enough for the Decision?
Suppose a hypothetical patient’s future concentration under a specified dose plan is modelled by a posterior predictive distribution approximated as normal, with mean 10.0 mg/L and standard deviation 2.00 mg/L. For illustration, suppose the decision is whether the concentration will lie within a desired range of 8.00 to 12.0 mg/L. The normal approximation is assumed adequate for this calculation; in practice, its fit and the model’s predictive assumptions also need assessment.
A central 90% prediction interval is approximately the mean plus or minus 1.645 standard deviations:
The probability of falling in the desired range is, under this normal predictive distribution,
Thus the model assigns about a 68.3% probability to the concentration being within the specified range and about a 31.7% probability to it being outside. Whether that is adequate depends on the consequences and available alternatives. The interval alone does not decide that question. Nor does this calculation establish that the forecast is trustworthy: the dose record, data support, model fit, and applicability to the patient still require review.
Assess the Information in the Individual Data
Data adequacy is not just a count of concentrations. Information depends on when samples were taken relative to dosing, the parameters that affect the observed profile, assay and residual error, and whether the concentrations span regions that distinguish competing parameter values. Several samples clustered in a region where the model predictions barely change with clearance may provide little information about clearance. A smaller number of well-timed samples can sometimes be more informative, though sparse data demand caution.
The practical question is whether the likelihood from the individual’s observations materially informs the parameters relevant to the forecast, rather than merely producing a posterior estimate. This can be assessed by comparing prior and posterior distributions, examining uncertainty in the forecast-relevant parameter combinations, and considering whether the data provide leverage at the relevant times. Eta shrinkage, influential observations, and sparse-data limitations can help identify weak individual information, but none should be interpreted in isolation.
Worked Example: Prior-Data Agreement and the Size of the Update
Consider a hypothetical scalar parameter \(\eta\), expressed as a dimensionless log-clearance deviation from a reference value. Suppose its population prior is \(\eta\sim N(0,0.30^2)\). A patient’s data are summarized, under a normal likelihood approximation, by an estimate \(z=0.30\) with standard error \(0.10\), so \(z\mid\eta\sim N(\eta,0.10^2)\). Assume the prior and likelihood are appropriate and the summary is unbiased with the stated standard error.
For normal prior and likelihood, posterior precision is the sum of prior and data precisions. The posterior variance and mean are therefore:
The posterior standard deviation is \(\sqrt{0.00900}=0.0949\). A normal 95% posterior interval is approximately \(0.270\pm1.96(0.0949)\), or \([0.0840,0.456]\). The data have shifted the estimate away from the prior mean and reduced its uncertainty; the posterior is not simply a report of the population typical value.
A simple prior-predictive discrepancy for the summary is its standardized difference from the prior mean. Under the prior and likelihood assumptions, \(z\) has prior-predictive variance \(0.30^2+0.10^2=0.100\), giving:
This value is not an alarming discrepancy by itself under a standard normal reference. It supports the limited conclusion that this summary is not strongly at odds with the specified prior predictive distribution. It does not show that the prior is correct, that the model is adequate, or that other observations are well predicted. The posterior interval is also not a prediction interval for a future concentration.
Check Predictive Agreement, Not Just Parameter Agreement
Prior-data agreement asks whether the observed patient information is surprising under the population prior and observation model. Predictive diagnostics ask whether the model can reproduce relevant features of the observed concentrations. These address related but different concerns. A patient can have a plausible parameter update while the model still fails to reproduce the concentration-time pattern, for example because of a structural mismatch or a recording error.
Posterior predictive checks compare observed data with data replicated from the model using the patient’s posterior distribution and relevant sampling and dosing conditions. Useful summaries may include residual patterns over time, the maximum deviation, or whether observations fall in predictive intervals. A check is more informative when its summary targets a plausible failure mode. A single observation inside a broad predictive interval is weak evidence; a systematic pattern or extreme discrepancy across relevant observations deserves investigation. Because the data also informed the posterior, posterior predictive checks can be less sensitive to misspecification than checks based on genuinely held-out observations.
Worked Example: A Forecast Fit With a Concentration-Time Pattern
Suppose four hypothetical concentrations are compared with posterior predictive expectations at their recorded sampling times. For a simple diagnostic, assume predictive standard deviations of 1.00, 1.00, 0.80, and 1.00 mg/L, respectively. The model expectations are 8.00, 6.00, 4.00, and 3.00 mg/L; the observed values are 8.40, 6.30, 4.20, and 7.00 mg/L. This diagnostic assumes the dose and sample times are correct and that the stated predictive distributions include the relevant variability.
The standardized differences \((\text{observed}-\text{expected})/\text{predictive standard deviation}\) are:
The first three observations are close to their expectations, while the fourth is 4.00 predictive standard deviations above its expectation. This is a notable discrepancy to investigate, not an automatic proof that the model is wrong. The sampling and dose times, dose amount, assay result, and clinical changes should be checked first. If the record is sound, the pattern may indicate that the model does not capture the patient’s concentration trajectory. A forecast based on the same model should then be treated cautiously, even if its posterior parameter estimates appear precise.
Dose-History Quality and Model Applicability
A forecast conditions on a dose history. Incorrect dose amounts, administration times, missed doses, infusion durations, or sample times can distort the inferred parameters and the predicted concentration. Record quality should therefore be assessed as part of the forecast, not as clerical detail. Distinguish what is documented from what is assumed. When adherence or timing is uncertain, consider whether plausible alternatives change the forecast enough to change the decision.
A forecast also assumes that the model remains applicable over the forecast period. Changes in renal or hepatic function, interacting medicines, body size, disease state, or other modelled covariates can matter if they alter relevant parameters. If the model assumes stable clearance but a patient’s clearance is changing, a narrow interval conditional on stable clearance may give false reassurance. A changing covariate or nonlinear process should be represented only if the model and available information support it; otherwise, the limitation should remain explicit.
A practical sensitivity assessment asks whether reasonable alternative assumptions about uncertain inputs or model components alter the clinically relevant prediction. For example, calculate the forecast under two plausible dose-time histories. If one implies a low probability of crossing a threshold and the other implies a high probability, the forecast is not robust to the record uncertainty. This is more informative than reporting one precise result while hiding an uncertain input.
A Structured Assessment and Its Limits
A compact assessment can be organized into linked checks. It is a reasoning aid, not a validated numerical score. The conclusion should say what is supported, what remains uncertain, and which decision the forecast can or cannot support.
Check dose amount and timing, sample times, concentration units, and relevant covariates. Identify assumptions about adherence or missing records.
Review sample timing and informativeness, and determine whether the data update parameters relevant to the forecast rather than simply returning a population-based estimate.
Use the posterior predictive distribution and summarize probabilities or intervals in relation to the clinical threshold or action under consideration.
Check prior-data agreement and predictive diagnostics; investigate influential observations and test whether plausible alternative inputs materially change the result.
Specify the model, dose plan, time horizon, and known limitations. If a key check fails, qualify the forecast or obtain additional information before relying on it.
The main alternative to this structured synthesis is to rely on a single indicator: a narrow interval, low shrinkage, a good fit statistic, or agreement with the population prior. Each can be useful, but each answers only part of the question. A narrow interval can be conditional on a wrong dose record; low shrinkage does not rule out structural misspecification; a good in-sample fit does not guarantee future calibration; and close prior agreement may mean that sparse data contributed little. Conversely, a visible prior-data difference is not inherently a failure: it may represent real individual variation when the data and model support it.
Calibration across patients can reveal whether forecasts are systematically too high, too low, or overconfident in a population. It cannot establish that a particular patient’s dose record is correct or that this patient’s forecast is adequate for a specific decision. Individual diagnostics and population-level calibration are complementary. Neither removes uncertainty about unrepresented model errors or clinical changes.
Common Mistake: Treating Precision as Trustworthiness
- Wrong reasoning: “The forecast interval is narrow, so the patient’s forecast is reliable.”
- Correction: The interval may be narrow only under the assumed model and recorded history. Verify those assumptions and check predictive agreement before interpreting its width.
- Wrong reasoning: “The observed concentrations agree with the prior, so the model is validated.”
- Correction: Agreement is only absence of an obvious conflict for the data examined. It does not establish predictive adequacy, especially when the observations are sparse or weakly informative.
- Wrong reasoning: “One extreme residual proves that the model is misspecified.”
- Correction: First check data, dose, and sample records and consider expected variability. A substantial unexplained pattern raises concern, but a diagnostic discrepancy is evidence to investigate rather than a standalone diagnosis.
Check Your Understanding
Answer each question, then compare your reasoning with the solutions.
- A forecast has a narrow posterior interval for clearance but a wide prediction interval for the next concentration. Give one model-based reason this can occur.
- In the normal-prior example, what does the standardized prior-predictive discrepancy of 0.949 establish, and what does it not establish?
- A concentration is far outside its predictive expectation. Name two input checks to make before concluding that the model is inadequate.
- Why is a forecast conditional on a fixed future dose plan not necessarily a reliable prediction of what will occur in practice?
Solutions. (1) Residual or process variability can remain substantial even when parameter uncertainty is small; the predictive distribution includes both sources. (2) It indicates that the observed summary is not strongly discrepant from the specified prior predictive distribution under the stated assumptions. It does not validate the prior or model, establish good prediction of other observations, or show that the future forecast is accurate. (3) Check, for example, the dose amount and administration time, and the sample time and concentration units; also review assay and adherence records. (4) The patient may not follow the assumed plan, or the actual dose history may differ. The conditional prediction does not integrate over that uncertainty unless it is explicitly modelled.