Tutorials › Population Pharmacokinetics › Observed Versus Population-Predicted Concentration Plots

Diagnostics: Goodness of Fit and Shrinkage · Tutorial 342 of 599

Observed Versus Population-Predicted Concentration Plots

You will learn to build observed-versus-population-predicted plots and assess calibration, scatter, and groupwise bias without mistaking visual agreement for evidence of predictive performance.

Intermediate 10 min read

What You'll Learn

  • Define population-predicted concentrations for a subject’s dosing history, covariates, and observation time.
  • Construct a plot with observed concentrations on the vertical axis and population predictions on the horizontal axis.
  • Use the identity line to assess calibration and identify systematic over- or underprediction.
  • Distinguish scatter around the identity line from groupwise bias.
  • Recognize why in-sample visual agreement does not establish predictive performance.

What the Plot Compares

An observed-versus-population-predicted concentration plot compares each measured concentration with the concentration predicted for that observation by the population model. It is a compact way to ask whether the model’s typical predictions are calibrated to the observations across the concentration range and across relevant subject groups.

Let \(C_{ij}\) be the observed concentration for subject \(i\) at observation \(j\), and let \(\widehat C^{\mathrm{pop}}_{ij}\) be the corresponding population prediction. The latter uses the fitted population parameters, the subject’s recorded covariates and dosing history, and the observation time, but sets the subject-specific random effects to their population-typical value of zero. Thus, it is not the same as an overall mean concentration: two observations can have different predictions because their times, regimens, or covariates differ.

Definition: An observed-versus-population-predicted concentration plot places \(C_{ij}\) on the vertical axis and \(\widehat C^{\mathrm{pop}}_{ij}\) on the horizontal axis, with one point for each observed concentration. The identity line \(y=x\) marks exact agreement between the observed and predicted values. A point above this line has an observed concentration greater than its population prediction; a point below it has an observed concentration less than its population prediction.

The usual target is that observations be distributed around the identity line without a systematic trend or important subgroup displacement. This is a calibration target, not a claim that each observed concentration should equal its prediction. Measurement error and between-subject variability produce scatter even when a model is adequate.

Calibration, Scatter, and the Identity Line

Calibration means that predictions agree, on average, with the outcomes they are intended to predict. For a predictor representing the conditional mean concentration on the plotted scale, a simple calibration target is \(E(C_{ij}\mid \widehat C^{\mathrm{pop}}_{ij}=c)\approx c\). The identity line visualizes this target. In a real plot, the conditional average is not known; the point cloud provides only a noisy diagnostic of it.

A vertical departure from the identity line can be expressed as a signed difference:

$$ d_{ij}=C_{ij}-\widehat C^{\mathrm{pop}}_{ij}. $$

Positive differences indicate underprediction and negative differences indicate overprediction. Their size must be interpreted in concentration units. A difference of \(0.5\ \mathrm{mg/L}\) may be substantial near a predicted concentration of \(1\ \mathrm{mg/L}\), but relatively small near \(50\ \mathrm{mg/L}\). Ratios can complement differences when proportional error is relevant, but they do not replace choosing a scale suited to the model and clinical question.

The cloud’s width also matters. Points widely dispersed around the identity line indicate that population predictions alone do not closely describe individual observations. That can reflect genuine between-subject variability, residual error, model misspecification, or sparse information about the subjects. A narrow cloud is not automatically evidence of a good model: the observations may cover only a limited concentration range, or a systematic offset may be present.

The identity line is a reference, not a fitted regression line. Fitting a line to the same points can summarize an apparent intercept and slope, but does not turn the plot into a validation study. Repeated observations from one subject are also not independent replicates. A dense cluster from a few subjects can dominate the visual impression, so inspect subject and group representation rather than treating every point as independent evidence.

Worked Example: Reading Differences Across the Concentration Range

Suppose a hypothetical dataset has four observations, with concentrations in \(\mathrm{mg/L}\). The population predictions are \(0.9,\ 1.8,\ 3.6,\) and \(7.2\), and the observed values are \(1.0,\ 2.0,\ 4.0,\) and \(8.0\). Calculate the signed differences and compare the observations with the identity line.

For the first observation, \(d=1.0-0.9=0.1\ \mathrm{mg/L}\). For the second, \(d=2.0-1.8=0.2\ \mathrm{mg/L}\). The remaining differences are \(4.0-3.6=0.4\ \mathrm{mg/L}\) and \(8.0-7.2=0.8\ \mathrm{mg/L}\).

All four points are above the identity line, so each observed concentration exceeds its population prediction. The absolute difference increases with the predicted concentration. However, the ratios are \(1.0/0.9=1.11\), \(2.0/1.8=1.11\), \(4.0/3.6=1.11\), and \(8.0/7.2=1.11\), to three significant figures. The observations are therefore about 11% higher than the predictions throughout this small example. This consistent proportional departure is more apparent from ratios or a proportional-error interpretation than from the absolute differences alone. Four points are far too few to establish a general pattern.

Why Groupwise Bias Can Be Hidden

A satisfactory-looking overall cloud can conceal systematic errors in a clinically relevant subgroup. Groups may be defined by categorical covariates, such as formulation or renal-function category, or by other scientifically justified characteristics. The plot can be inspected with group symbols or colors, but the group definitions should be chosen for a reason; arbitrary subdivision can make chance fluctuations look meaningful.

A useful basic check is to compute signed differences within each group. A group mean near zero is not proof of adequate calibration: positive and negative differences can cancel, and a small group can give an unstable estimate. Examine the individual points, the range of predictions, and the number of subjects represented alongside any summary. When group sizes or prediction ranges differ substantially, an overall mean difference can be particularly misleading.

Worked Example: Opposing Bias in Two Groups

Consider hypothetical observations from two groups, with concentrations in \(\mathrm{mg/L}\). Each group has population predictions \(2,\ 2,\ 4,\) and \(4\). Group A has observations \(2.2,\ 2.1,\ 4.3,\) and \(4.0\); Group B has observations \(1.6,\ 1.8,\ 3.2,\) and \(3.6\). Calculate the group mean signed differences.

For Group A, the differences are \(0.2,\ 0.1,\ 0.3,\) and \(0.0\ \mathrm{mg/L}\). Their mean is \((0.2+0.1+0.3+0.0)/4=0.15\ \mathrm{mg/L}\). The positive mean indicates slight underprediction in Group A.

For Group B, the differences are \(-0.4,\ -0.2,\ -0.8,\) and \(-0.4\ \mathrm{mg/L}\). Their mean is \((-0.4-0.2-0.8-0.4)/4=-0.45\ \mathrm{mg/L}\). The negative mean indicates overprediction in Group B. Across all eight observations, the mean difference is \((0.6-1.8)/8=-0.15\ \mathrm{mg/L}\).

The pooled mean suggests only modest overall overprediction, but it obscures the opposing group patterns. A plotted display with group symbols would make that distinction visible. These summaries describe this hypothetical sample; they do not establish that the population model has a stable group effect or identify its cause.

A Justification for the Identity-Line Target

The identity line follows directly from the definition of calibration. Suppose the population prediction \(P\) represents the conditional mean of the observed concentration \(C\) on the chosen scale. Then, under calibration, \(E(C\mid P)=P\). Subtracting \(P\) from both sides gives the expected signed difference at that prediction:

$$ E(C-P\mid P)=E(C\mid P)-P=P-P=0. $$

Thus, under this condition, the average signed difference should be zero at each prediction value, corresponding to observations centered on \(y=x\). The argument is conditional: it relies on the prediction representing the conditional mean on the plotted scale. For example, a model for log concentration does not automatically imply that an exponentiated prediction is the mean concentration; transformation and the error distribution matter. The plotted scale and the meaning of the model prediction must therefore be consistent.

This justification concerns average calibration, not individual agreement. Even if \(E(C\mid P)=P\), the conditional variance \(Var(C\mid P)\) can be large. The identity-line target alone says nothing about whether the observed scatter is clinically acceptable.

Worked Example: High Correlation Despite Miscalibration

For a small hypothetical dataset, suppose the population predictions are \(1,\ 2,\ 4,\) and \(8\ \mathrm{mg/L}\), while the observations are \(2,\ 3,\ 5,\) and \(9\ \mathrm{mg/L}\). Explain why a strong linear relationship would not establish calibration.

For every observation, \(C=P+1\ \mathrm{mg/L}\). The signed differences are therefore \(1,\ 1,\ 1,\) and \(1\ \mathrm{mg/L}\), and every point lies one unit above the identity line. The points fall exactly on the line \(C=P+1\), so their correlation is \(1\), but their predictions are systematically low by \(1\ \mathrm{mg/L}\). The average signed difference is \(1\ \mathrm{mg/L}\), not zero.

This example shows why correlation or a visually straight trend is not the calibration criterion. Agreement with the identity line is the relevant reference. It also illustrates a scale issue: a constant difference may have a large relative impact at the lowest predicted concentrations and a smaller one at the highest.

Assumptions, Diagnostics, and Limits

Before interpreting the plot, verify that each plotted pair refers to the same subject, observation time, dose history, and concentration units. Predictions must be generated from the fitted population model using the intended population-prediction convention. Data errors, mismatched records, or predictions calculated for the wrong time after dose can create apparent model failure that is not a model property.

The display is most informative when the dataset covers a meaningful range of predicted concentrations and includes adequate representation of subjects and relevant groups. Sparse sampling, a narrow concentration range, and few observations in a subgroup limit what can be seen. Repeated measurements from a subject share information and often have correlated errors; consequently, a visually dense plot should not be read as if it contained the same amount of independent evidence as an equally sized sample of unrelated subjects.

Warning signs include a consistent displacement above or below the identity line, curvature, changing scatter with concentration, or distinct group clusters. These patterns are prompts for investigation, not automatic diagnoses of a particular model defect. A trend may reflect a structural model issue, an observation-model issue, data handling, or a population not represented adequately by the fitted model. Use the plot alongside residual diagnostics, including the weighted-residual material, and check whether the pattern is plausible in the study context.

The main alternative displays answer different questions. Residual-versus-prediction plots show signed or standardized departures more directly and can make trends in residual spread easier to see. Group-stratified displays emphasize differences between selected groups. Observed-versus-individual-predicted plots use a different prediction target, so they should not be treated as interchangeable with this population-prediction plot. No one display summarizes all aspects of fit.

Most importantly, a plot using observations that contributed to model fitting is an in-sample diagnostic. The model may appear to agree because its parameters were estimated using those data, and sparse or noisy observations may conceal important shortcomings. Even an apparently good in-sample plot does not demonstrate accurate predictions for new subjects, new regimens, or settings outside the observed covariate and concentration ranges. Predictive performance requires evidence appropriate to that question, not visual agreement alone.

Common Pitfall

  • Wrong reasoning: “The points have a strong upward trend and a high correlation, so the population predictions are accurate.”
  • Correction: A high correlation describes how strongly the two variables move together, not whether their values agree. Check displacement from the identity line, the signed differences, concentration-dependent patterns, and groupwise behavior. Then remember that an in-sample diagnostic is not external predictive validation.
  • Wrong reasoning: “The overall mean difference is close to zero, so no group is biased.”
  • Correction: Opposing group differences can cancel. Inspect relevant groups and their representation before drawing an overall conclusion.
Key takeaway: Plot observed concentrations against population predictions and use the identity line to assess calibration. Evaluate scatter and relevant group patterns separately, and treat visual agreement as a diagnostic rather than proof of predictive performance.

Check Your Understanding

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

  1. A point has an observed concentration of \(3.2\ \mathrm{mg/L}\) and a population prediction of \(2.8\ \mathrm{mg/L}\). Where does it lie relative to the identity line, and what is its signed difference?
  2. Two subgroups have mean signed differences of \(0.6\ \mathrm{mg/L}\) and \(-0.5\ \mathrm{mg/L}\). Why might a pooled mean near zero be misleading?
  3. Why does a correlation of \(1\) not guarantee that population predictions are calibrated?
  4. Name one reason why many plotted points may provide less evidence than the point count suggests.

Solutions. (1) It lies above the identity line, with \(d=3.2-2.8=0.4\ \mathrm{mg/L}\), indicating underprediction for that observation. (2) The group biases oppose one another and can cancel in the pooled summary; inspect the subgroup patterns and their sample sizes. (3) A perfect linear relationship can still have a nonzero intercept or an incorrect slope relative to \(y=x\); calibration requires agreement with the identity line, not merely a straight trend. (4) Multiple observations may come from the same subject and are not independent; a few subjects with many measurements can dominate the display.

Clinical Biostats