Tutorials › Nonlinear Mixed-Effects Models (NLME) › Prediction Intervals Versus Confidence and Tolerance Intervals

Uncertainty in Individual Predictions · Tutorial 756 of 797

Prediction Intervals Versus Confidence and Tolerance Intervals

Learn to choose and interpret prediction, confidence, and tolerance intervals for different pharmacokinetic questions.

Advanced 12 min read

What You'll Learn

  • Identify whether an interval concerns one future concentration, a population parameter, or a specified fraction of individuals.
  • Relate interval width to the quantity being predicted or estimated.
  • Calculate a classical prediction interval and a confidence interval for a population mean.
  • Interpret a model-based tolerance interval and its confidence level.
  • Recognize assumptions and diagnostics that affect interval validity.

Three Intervals, Three Questions

An interval is useful only when its target is clear. In pharmacokinetics, one question may concern a future concentration for a particular patient; another may concern the population mean concentration; a third may ask whether a dose keeps concentrations within a range for a large fraction of the target population. These are different inferential tasks. Their intervals generally have different widths and cannot be substituted for one another.

For an individual forecast, the conditioning information matters. It could include a planned dose and covariates, or observed concentrations for that individual. A prediction interval describes uncertainty about a future observation under those conditions. A confidence interval describes uncertainty in estimating a population quantity. A tolerance interval describes an interval intended to contain a stated fraction of a population, with a stated confidence level.

Definition: A prediction interval is an interval for a future observation, conditional on the stated information and model. A confidence interval is an interval produced by a procedure that, under repeated sampling and its assumptions, covers a fixed population parameter at the stated confidence level. A tolerance interval is an interval intended to contain at least a specified fraction of a target population, with a stated confidence level.

For example, “95% prediction interval” refers to the probability content of a procedure for future observations, under its model and conditioning conditions. “95% confidence interval” refers to repeated-sample coverage of a parameter. A “95%/95% tolerance interval” aims to contain at least 95% of the target population with 95% confidence. The two percentages in a tolerance statement have distinct roles.

Prediction and Confidence Intervals in a Normal Benchmark

A simple normal model makes the distinction in width transparent. Suppose \(X_1,\ldots,X_n\) are independent observations from \(N(\mu,\sigma^2)\), and \(X_{\text{new}}\) is an independent future observation from the same distribution. Assume normality and common variance, and estimate \(\mu\) and \(\sigma\) by \(\bar{x}\) and \(s\). Under these assumptions, the usual two-sided \(100(1-\alpha)\%\) confidence interval for \(\mu\) is \(\bar{x}\pm t_{1-\alpha/2,n-1}s/\sqrt{n}\). The corresponding prediction interval for \(X_{\text{new}}\) is \(\bar{x}\pm t_{1-\alpha/2,n-1}s\sqrt{1+1/n}\).

The prediction interval is wider because it must account for both uncertainty in the estimated mean and the future observation’s variation around that mean. The confidence interval concerns only the mean parameter.

Justification: The difference \(X_{\text{new}}-\bar{X}\) has mean zero. Independence gives \(\operatorname{Var}(X_{\text{new}}-\bar{X})=\sigma^2+\sigma^2/n=\sigma^2(1+1/n)\). Estimating \(\sigma\) with \(s\) and using the \(t\) reference distribution gives the prediction standard error \(s\sqrt{1+1/n}\). By contrast, \(\operatorname{Var}(\bar{X})=\sigma^2/n\), giving standard error \(s/\sqrt{n}\) for the mean.

These formulas are a benchmark, not a substitute for an NLME forecast. In an NLME model, the target might be a future concentration conditional on an individual’s data, dose, and covariates; the forecast distribution is obtained from the specified model and conditioning information. Its interval is still about a future observation, not automatically about a population parameter or the fraction of the population inside a range.

Worked Example: A Future Concentration Versus the Population Mean

Suppose, hypothetically, that 36 independent concentrations at a specified time and under a common regimen are modeled as normal observations. The sample mean is \(10.0\) mg/L and the sample standard deviation is \(2.00\) mg/L. Assume the normal model and common variance are adequate, and use a 95% interval with \(df=35\) and \(t_{0.975,35}\approx2.030\).

For one future concentration from the same distribution, the prediction interval is

$$ 10.0 \pm 2.030(2.00)\sqrt{1+\frac{1}{36}} \approx 10.0\pm4.116 =[5.884,\ 14.116]\ \text{mg/L}. $$

The interval is for one future concentration under the assumed common distribution. It is not an interval for the population mean. The model assumptions are explicit: independent observations, normality, a shared mean and variance, and a future observation from that same distribution.

For comparison, the 95% confidence interval for the population mean is

$$ 10.0 \pm 2.030\frac{2.00}{\sqrt{36}} =10.0\pm0.6767 =[9.323,\ 10.677]\ \text{mg/L}. $$

The confidence interval is narrower because it estimates the mean, rather than covering a future individual observation. In repeated samples under the assumptions, 95% of intervals produced by this confidence procedure would cover the fixed population mean; this does not mean there is a 95% probability that this particular fixed mean lies in the calculated interval.

Tolerance Intervals and Population Coverage

A tolerance interval addresses a population-level coverage question. For example, a clinician or development team may ask whether a concentration range contains at least 95% of concentrations for a defined target population at a specified time and regimen. The target population must be stated: it may refer to individuals with specified covariates, or to a broader covariate distribution. Changing that target changes the population coverage being assessed.

In a fitted NLME model, the distribution relevant to tolerance coverage includes the model’s variation across individuals and the observation process if the target is measured concentrations. The exact distribution depends on the question. A range for latent individual typical concentrations is not necessarily a range for observed concentrations, which also reflect residual error. Likewise, a population-wide tolerance interval is not an interval for a particular patient’s next observation.

One practical model-based approach uses a parametric bootstrap. Repeatedly generate data under the fitted model, refit the model, and estimate the target population fraction inside a candidate interval for each refit. The lower tail of those estimated fractions quantifies uncertainty in coverage due to sampling and fitting. This is an approximate tolerance procedure: its repeated-sample interpretation depends on the model, bootstrap design, and accuracy of the coverage calculation.

Worked Example: Checking a 95%/95% Tolerance Interval

Consider a hypothetical NLME analysis for observed concentration at a fixed time under a specified regimen. The target is the population with the stated covariate distribution. A proposed interval is \([4.0,\ 16.0]\) mg/L. Assume the fitted model is suitable for generating data and describing this target population.

To assess the interval, the analyst conducts 1,000 parametric bootstrap replicates. In each replicate, the model is refitted and 100,000 simulated individuals are used to estimate the fraction whose concentration falls between 4.0 and 16.0 mg/L. Suppose the empirical fifth percentile of those 1,000 estimated fractions is 0.952.

The lower 5th percentile is 0.952, which exceeds the target content of 0.95. On this bootstrap assessment, the interval is a candidate 95%/95% tolerance interval: the procedure gives approximately 95% confidence that the interval covers at least 95% of the target population. This conclusion is approximate, not a guarantee; it relies on the fitted model and the bootstrap approximation. Increasing the number of simulated individuals reduces Monte Carlo noise in each replicate’s estimated fraction, but does not repair a misspecified model.

A Confidence Interval for a Population Parameter

A confidence interval can target a parameter other than the mean, such as a fixed-effect coefficient or a variance parameter. Its interpretation remains parameter-focused: its repeated-sample coverage concerns that parameter. For an estimated population mean, the normal-theory interval above uses the sampling uncertainty in \(\bar{x}\); it does not describe the spread of individuals around the mean.

Worked Example: A 97.5% Confidence Interval for Mean Concentration

Suppose, hypothetically, that 25 independent concentrations from a normal population have \(\bar{x}=10.0\) mg/L and \(s=2.00\) mg/L. Assume independent sampling and normality, and seek a two-sided 97.5% confidence interval for the population mean. Here \(df=24\), and \(t_{0.9875,24}\approx2.391\).

The standard error is \(SE=s/\sqrt{n}=2.00/\sqrt{25}=0.400\) mg/L. Thus,

$$ 10.0\pm2.391(0.400) =10.0\pm0.9564 =[9.0436,\ 10.9564]\ \text{mg/L}. $$

This interval estimates the population mean concentration. It does not imply that 97.5% of individual concentrations lie between its endpoints; individual values are more variable than the estimated mean.

Choosing, Checking, and Interpreting the Interval

IntervalTargetTypical questionWhat its level refers to
PredictionOne future observationWhere might this patient’s next concentration fall, given the stated information?Coverage of future observations under the predictive model
ConfidenceA fixed population parameterWhat values are compatible with the data for the population mean?Repeated-sample coverage of that parameter
ToleranceA specified fraction of a populationWhat range is expected to contain at least a stated fraction of the target population?Confidence in achieving the stated population coverage

Before reporting an interval, state its target, conditioning information, units, and assumptions. For a prediction interval, check that the model represents the intended future setting and that the forecast distribution accounts for the variability relevant to the observation being predicted. For a confidence interval, check whether its sampling or model-based approximation is appropriate for the parameter and data. For a tolerance interval, check that the target population and coverage definition are unambiguous and that the fitted model reproduces important distributional features, including tails.

Useful diagnostics include comparing simulated and observed concentration distributions, checking residual patterns and extreme values, assessing whether the covariate distribution used for a population target is appropriate, and examining whether bootstrap fits are stable. A tolerance interval can be especially sensitive to misspecified tails: a model that fits the center well may still give unreliable estimates of the fraction outside the interval. A prediction interval may also be poorly calibrated if the individual’s conditioning information or future dosing conditions differ from those represented by the model.

The main alternatives answer different questions. A confidence interval is appropriate for uncertainty in a parameter, not individual spread. A prediction interval is appropriate for a future individual observation, not the fraction of an entire population captured. A tolerance interval directly targets population content, but usually requires more modelling or computation to account for uncertainty in that content. A credible interval for a Bayesian parameter has a posterior probability interpretation conditional on the model and prior; a posterior predictive interval instead concerns future observations. The label alone is not enough: always identify what random quantity the interval concerns.

A common mistake is to call a narrow confidence interval for the mean a “95% range for patients.” That reasoning confuses uncertainty in the estimated mean with variability among individuals. The correction is to choose a prediction interval for one future observation or a tolerance interval for a stated population fraction, depending on the scientific question. Neither interval eliminates model dependence: misspecification, unrepresentative covariates, or an incorrect residual model can make the reported coverage misleading.

Key takeaway: Name the target before calculating or interpreting an interval. Prediction intervals concern future observations, confidence intervals concern population parameters, and tolerance intervals concern population coverage. Their percentages refer to different properties.

Check Your Understanding

Answer each question, then compare with the brief solutions.

  1. A 95% interval is calculated for the population mean concentration. Does it aim to contain 95% of individual concentrations?
  2. In one sentence, distinguish the target of a prediction interval from the target of a tolerance interval.
  3. A 95%/95% tolerance interval is reported for a defined population. What do the two 95% values mean?
  4. Why can a model-based tolerance interval be unreliable even when the fitted model describes the center of the observed concentration distribution well?

Solutions. (1) No. It targets the fixed population mean; individual concentrations generally vary more widely. (2) A prediction interval targets one future observation under stated conditions, whereas a tolerance interval targets a stated fraction of a population. (3) The procedure has 95% confidence that the interval contains at least 95% of the target population. (4) Its population coverage can depend strongly on tail behavior, so a model that misses the tails can estimate the covered fraction poorly.