Tutorials › Physiologically-Based Pharmacokinetics (PBPK) › Setting Predictive Error Criteria Before Seeing Results

Calibration Versus Independent Validation · Tutorial 711 of 799

Setting Predictive Error Criteria Before Seeing Results

Define decision-relevant acceptance rules in advance and assess each PBPK prediction endpoint against its own prespecified error scale and threshold.

Advanced 10 min read

What You'll Learn

  • Explain why acceptable prediction error depends on the intended use and endpoint.
  • Set multiplicative limits for concentrations, AUC, and Cmax and an absolute tolerance for time to peak.
  • Translate fold-error limits into log-ratio limits and explain their symmetry.
  • Apply prespecified rules to hypothetical concentration-time and summary-endpoint data.
  • Distinguish a prediction-acceptance criterion from a statistical hypothesis test.
  • Identify how pooled results, post hoc thresholds, and unsuitable error scales can mislead.

Why Set Prediction Limits in Advance?

A validation result is useful only in relation to a stated standard for acceptable predictive performance. A twofold error might be tolerable for an exploratory estimate of a concentration profile but too large for a decision that depends on a narrow exposure margin. Likewise, a one-hour error in time to peak may be negligible for one sampling schedule and consequential for another. There is no single error limit that is automatically appropriate for every PBPK application.

Setting criteria before inspecting validation outcomes makes the assessment interpretable: the rules represent the intended use rather than a response to the results. It also prevents a model from being labelled acceptable by choosing a forgiving metric or threshold after its errors are known. Criteria should be fixed before unblinding the predictions and observations, and should state the endpoints, error scales, summaries, thresholds, and units used in the decision.

The earlier material on fold error and bias provides useful measures of prediction error. Here the focus is how to choose and apply acceptance criteria for those measures. A criterion is a prespecified decision rule, not a claim that any particular numerical threshold is universally correct.

Define the Endpoint, Scale, and Decision Rule

Definition: A predictive error criterion is a rule, specified before validation outcomes are examined, that identifies which values or summaries of prediction error will count as acceptable for a defined endpoint and intended use. It includes the prediction-observation pairing, error measure, threshold, aggregation rule, and decision rule.

For a positive observed value \(Y\) and its corresponding prediction \(\hat{Y}\), the prediction ratio is \(R=\hat{Y}/Y\). A fold-error bound \(F\), with \(F>1\), classifies a pair as within the bound when \(1/F \leq R \leq F\). This multiplicative scale treats an overprediction and its reciprocal underprediction comparably: a prediction twice the observation has ratio 2, while one half the observation has ratio 0.5.

$$ \text{Within an }F\text{-fold bound} \quad\Longleftrightarrow\quad \frac{1}{F}\leq \frac{\hat{Y}}{Y}\leq F \quad\Longleftrightarrow\quad \left|\log\left(\frac{\hat{Y}}{Y}\right)\right|\leq \log(F) $$

The last form follows by taking logarithms of the positive ratio and expressing the two-sided bounds around zero. It shows why fold errors are symmetric on a log scale, even though the corresponding absolute differences in the original units are not equal. For example, an \(F=2\) bound allows ratios from 0.5 through 2. It does not allow the same number of concentration units above and below the observation.

The criterion must also specify how multiple comparisons are combined. One option is to require a stated proportion of prespecified concentration-time pairs to fall within a fold bound, separately for each validation condition. Another is to require each condition's AUC ratio to fall inside a narrower interval. A pooled pass rate across all studies can hide a poor result in one important condition, so the aggregation rule should match the intended domain of use.

Choose Limits From the Intended Use

An endpoint-specific criterion begins with the decision the model is meant to support. The analyst should ask how large an error could be without changing that decision, which direction of error is most consequential, and how the endpoint is measured. These questions help set a defensible scale and tolerance. They do not justify tuning a rule to make a particular validation data set pass.

Concentrations, AUC, and Cmax are positive quantities and are commonly assessed with ratios or fold-error limits. A concentration-time profile has many prespecified sampling times, so its rule must state whether performance is summarized over all pairs, at each time, or within each condition. AUC and Cmax are single summary endpoints for a dosing scenario and can be assessed with their own ratio intervals. Time to peak is a time-valued endpoint: an absolute difference in hours or minutes is often easier to interpret than a ratio, particularly when the observed time is near zero or limited to scheduled sampling times.

The limits should be justified by the model's intended application, relevant prior knowledge, and the consequences of errors. They should also be specific enough to reproduce the decision. “Predictions should be reasonably close” is not a criterion; “at least 80% of prespecified concentration pairs must be within twofold in each condition” is a rule that can be applied consistently. Neither example threshold is a universal standard.

Worked Example: Writing Criteria for a Hypothetical Use

Suppose a hypothetical PBPK model will be used to compare exposure profiles and summary endpoints for a proposed dose. Before seeing validation results, a team specifies the following illustrative rules. The scenario has concentration samples planned at five fixed times, and time to peak is assessed on an hourly sampling schedule.

Concentrations. In each validation condition, at least 80% of the prespecified matched prediction-observation pairs must have ratios between 0.5 and 2.0. This is a twofold bound because \(1/F=1/2=0.5\) and \(F=2\). Requiring the rule separately by condition prevents good performance in one condition from compensating for poor performance in another.

AUC and Cmax. Each endpoint ratio must lie between 0.80 and 1.25 in each condition. The lower limit is the reciprocal of the upper limit because \(1/1.25=0.80\). Thus the range is symmetric on the log-ratio scale, and its endpoints express the allowed multiplicative discrepancy.

Time to peak. The absolute difference between predicted and observed time to peak must be no greater than one hour, matching one scheduled sampling interval in this hypothetical design. The tolerance is in hours, not a percentage of the observed time.

These rules are internally clear, but their scientific suitability still depends on the intended decision. For another application, narrower or wider limits could be justified. The important point is that the limits and the rule for combining results are set before the validation outcomes are inspected.

Apply Concentration Criteria Without Changing Them

For concentration data, the paired values must refer to the same prespecified time, condition, and comparison. When the criterion uses fold error, the observations must be positive for the ratio to be defined. The analysis plan should identify the time points and conditions that enter the denominator before outcomes are known. Details such as how to handle values outside assay quantification limits require their own prespecified procedures; they should not be resolved by selectively removing inconvenient pairs.

A pass proportion is descriptive: it reports the fraction of the defined pairs meeting the threshold. It is not, by itself, an estimate that future samples will meet the threshold with that probability. In particular, several time points from one profile can be related rather than independent. The criterion should therefore be understood as the agreed performance rule for the specified validation data, not as an independence-based probability statement.

Worked Example: Checking a Hypothetical Concentration Profile

A hypothetical validation condition has five prespecified observed concentrations, in mg/L, and corresponding PBPK predictions. The rule set in advance requires at least 80% of the five pairs to fall within twofold. The pairs are: observed 10 and predicted 6 mg/L; observed 8 and predicted 6.4 mg/L; observed 5 and predicted 6 mg/L; observed 2 and predicted 3.4 mg/L; observed 1 and predicted 2.2 mg/L.

For each pair, calculate \(R=\hat{Y}/Y\). The ratios are \(6/10=0.60\), \(6.4/8=0.80\), \(6/5=1.20\), \(3.4/2=1.70\), and \(2.2/1=2.20\). The prespecified twofold interval is \(0.50 \leq R \leq 2.00\).

The first four ratios are within the interval; the fifth is above 2.00. Thus the pass proportion is \(4/5=0.80\), or 80%. The condition meets the stated concentration criterion exactly. This result does not imply that each concentration is accurate to 20%; the rule was a twofold bound applied to at least 80% of the pairs.

The fifth prediction is an overprediction at the lowest observed concentration. The pass/fail result alone does not explain why it occurred. A plot against time and inspection of the individual ratios can reveal whether errors cluster at a particular part of the profile, but such a pattern should be interpreted as a diagnostic rather than used to revise the acceptance rule after the fact.

Assess Summary Endpoints and Time to Peak

For AUC and Cmax, compute the prediction-to-observation ratio using the same endpoint definition for both values. Apply each endpoint's own prespecified interval. Do not substitute a concentration-profile criterion merely because it is more familiar: integrated exposure and peak concentration answer different questions and can show different errors. If the rule requires each condition to pass, report each condition rather than relying only on an overall average.

For time to peak, let \(\hat{t}_{\max}\) be the predicted time and \(t_{\max}\) the observed time. An absolute-error rule uses \(|\hat{t}_{\max}-t_{\max}|\), with units. The tolerance should reflect the resolution of the sampling schedule and the intended decision. Because observed peak time can depend on when samples are taken, a small difference should not automatically be interpreted as a precise estimate of a biological timing difference.

Worked Example: Applying Three Endpoint-Specific Rules

Consider one hypothetical validation condition with the rules from the earlier design example: AUC and Cmax ratios must each lie between 0.80 and 1.25, and the time-to-peak difference must be at most 1 hour. The observed AUC is 1,000 mg·h/L and the prediction is 920 mg·h/L. The observed Cmax is 10.0 mg/L and the prediction is 13.1 mg/L. The observed time to peak is 2 hours and the predicted time to peak is 3 hours.

For AUC, the ratio is \(920/1000=0.920\), which is between 0.80 and 1.25. AUC meets its criterion.

For Cmax, the ratio is \(13.1/10.0=1.31\), which exceeds 1.25. Cmax does not meet its criterion. The prediction is 31% above the observed value, while the prespecified upper ratio permits at most 25% above it.

For time to peak, the absolute difference is \(|3-2|=1\) hour. This is equal to, and therefore within, the maximum allowed difference of 1 hour. Time to peak meets its criterion.

The condition meets two of the three endpoint rules but does not meet all of them. If the overall decision rule requires every endpoint to pass, this condition fails that overall rule because Cmax fails. Reporting only the passing AUC and timing results would give an incomplete account.

Justification, Diagnostics, and Alternatives

The log-ratio form provides a short justification for using fold limits when proportional error is the concern. Let \(R=\hat{Y}/Y>0\). The two-sided fold rule \(1/F \leq R \leq F\) becomes \(-\log(F)\leq\log(R)\leq\log(F)\), since the logarithm is increasing and \(\log(1/F)=-\log(F)\). Equivalently, the absolute log ratio cannot exceed \(\log(F)\). This is why reciprocal errors have equal magnitude on the chosen scale: \(|\log(2)|=|\log(1/2)|\). The result justifies the scale's symmetry; it does not establish which value of \(F\) is acceptable.

Diagnostics should check whether the prespecified rule has been applied to the intended data and whether the chosen summary is concealing a failure. Confirm that predictions and observations are paired correctly, units and endpoint definitions agree, and the denominator includes exactly the records defined in advance. Then inspect errors by time, condition, and endpoint. A pooled pass rate can appear satisfactory even when one condition systematically performs poorly; reporting condition-specific results helps expose this. If a diagnostic suggests a problem, report it and investigate under the analysis plan rather than changing the threshold to restore a pass.

Alternative error measures can be useful for different purposes, but they answer different questions. Absolute error in concentration units gives larger weight to differences at high concentrations and is sensitive to the units used. Percentage error and fold error emphasize relative discrepancy, but ratios are unsuitable when the observed value is zero and can be unstable near zero. A single mean error can also conceal large errors in opposite directions that cancel. Choosing a different measure after seeing which one makes the model look best is not a defensible alternative; any additional measures should be specified in advance and reported alongside the primary rule.

A pass/fail criterion also compresses information. A result just inside a limit and one near perfect agreement both count as passes, while a slight miss and a large miss both count as failures. Report the underlying endpoint values and error summaries with the decision so readers can judge the magnitude and pattern of discrepancies. A criterion is a practical standard for a stated use, not proof that the model is universally accurate or a statistical test that the model is correct.

Common Pitfall

Common mistake: After seeing that a Cmax ratio is 1.31, an analyst widens the upper limit from 1.25 to 1.35 and declares success. This changes the decision rule in response to the result, so the apparent pass no longer evaluates a prespecified standard. The correction is to retain and report the original 1.25 limit, state that Cmax failed it, and explain any subsequent model revision or new validation separately.

Another error is to apply the same numerical tolerance to all endpoints simply for convenience. A limit of 1 hour has no meaning for AUC, while a 1.25 ratio may not be appropriate for time to peak. Endpoint-specific criteria preserve interpretable scales and make explicit which predictive performance matters for the intended decision.

Check Your Understanding

Use the stated criteria in each question. The values are hypothetical.

  1. A concentration criterion requires ratios from 0.5 through 2.0 for at least 75% of four prespecified pairs. The ratios are 0.40, 0.75, 1.50, and 2.00. Does the condition pass?
  2. An AUC criterion requires a ratio from 0.80 through 1.25. Observed AUC is 500 mg·h/L and predicted AUC is 625 mg·h/L. Does it pass?
  3. A time-to-peak criterion allows an absolute difference of no more than 0.5 hours. Observed time is 1.0 hour and predicted time is 1.6 hours. Does it pass?
  4. Why is it not valid to widen an error limit after seeing that a validation result falls just outside it?

Solutions. (1) Three of the four ratios are within the interval: 0.75, 1.50, and 2.00. The pass proportion is \(3/4=0.75\), or 75%, so the condition passes exactly. (2) The ratio is \(625/500=1.25\), equal to the allowed upper limit, so AUC passes. (3) The difference is \(|1.6-1.0|=0.6\) hours, greater than 0.5 hours, so time to peak fails. (4) Widening the limit in response to the observed result replaces the prespecified rule with a post hoc rule and biases the assessment toward a pass. The original criterion should be reported and applied unchanged.