Tutorials › Population Pharmacokinetics › Covariate Support and Transportability

Bootstrap and External Validation · Tutorial 432 of 599

Covariate Support and Transportability

You will learn to diagnose univariate and joint covariate support, quantify extrapolation in PK predictions, and distinguish overlap from a change in the external population’s covariate distribution.

Advanced 12 min read

What You'll Learn

  • Define empirical covariate support and distinguish it from population support.
  • Check external covariates against development ranges, categories, and joint combinations.
  • Recognize why a prediction can extrapolate even when each covariate is individually in range.
  • Interpret covariate support gaps as limits on transportability, not automatic proof of model failure.
  • Separate support overlap from differences in covariate distributions within the overlap.
  • Report support assessments with appropriate attention to sparse tails and data quality.

Why Covariate Support Matters

External validation asks how a locked model performs in data collected outside its development setting. A key question is whether the external subjects resemble the development subjects in the covariates used by the model. If an external subject has a body weight, age, renal function, or covariate combination not represented in development data, the model may have to predict beyond the information that informed its fitted relationship. That is extrapolation, and it limits how confidently performance can be transported to that subject or population.

A model can return a numerical prediction for an unusual covariate value without being well supported there. For example, a linear covariate relationship may produce a plausible-looking clearance estimate outside the observed range, but the development data cannot establish that the relationship continues linearly beyond that range. Conversely, a support gap is a warning about evidential limits, not by itself proof that the model’s prediction is wrong.

Definition: The empirical covariate support of a development dataset is the set of covariate values and combinations represented by its observed subjects. An external value is outside empirical support when it lies beyond the observed development data in a relevant covariate dimension or joint combination. This is a sample-based diagnostic: it does not establish the full support of the underlying development population.

The target of this assessment is the covariates that matter to the prediction: those used directly in the model, and those needed to define its intended application. Check units, coding, and measurement definitions first. A value of 90 may mean a body weight of 90 kg in one dataset and 90 lb in another; comparing the numbers without harmonizing units would give a false support assessment. Covariate support is also distinct from the prediction-performance summaries discussed in the earlier material on Prediction Error Metrics for External PK Validation and Subgroup Performance in External Validation.

From Marginal Ranges to Joint Support

For one continuous covariate, a simple initial check is to compare the development and external minima and maxima. A value outside the development range is an obvious extrapolation beyond the observed sample. But being between the minimum and maximum is not a guarantee of strong support: a value may be isolated in a sparse tail, or may occur only in a covariate combination absent from development.

With multiple covariates, assess the joint distribution as well as each marginal distribution. For example, body weight and renal function might each fall within their separate development ranges, while a low weight paired with very low renal function is scarcely or never observed in development. A model that includes both covariates, or their interaction, must predict for that combination nonetheless.

Useful summaries include overlaid histograms or density plots for continuous covariates, tables of category frequencies, scatterplots for pairs of continuous covariates, and counts for clinically meaningful combinations. In higher dimensions, displays are less complete; report the approach and its limitations rather than treating a single distance or range check as a definitive support test. The convex hull—the smallest convex region containing observed covariate vectors—is one geometric summary, but it can classify points as supported even when they are in very sparse regions. No one diagnostic replaces clinical judgment about which covariates and combinations are relevant.

Worked Example: A Prediction Beyond the Observed Age Range

Suppose a hypothetical development dataset contains adults aged 18 to 75 years. The locked clearance model is \( \log(CL_i)=\log(5.00)+0.010(\mathrm{Age}_i-40) \), where clearance \(CL_i\) is in L/h and age is in years. An external subject is 82 years old. The following calculation describes the model’s prediction; it does not establish that the age relationship is valid at 82 years.

At age 82, the linear predictor is \( \log(5.00)+0.010(82-40)=\log(5.00)+0.420 \). Thus, \( CL_i=5.00\exp(0.420)=7.61 \) L/h (rounded to three significant figures).

Age 82 is \(82-75=7\) years above the largest development age, so this prediction is an extrapolation in age. Relative to the model prediction at age 75, the ratio is \( \exp[0.010(82-75)]=\exp(0.070)=1.073 \). The model therefore predicts about 7.3% higher clearance at age 82 than at age 75. That percentage is a consequence of the specified model, not evidence that the relationship continues in this way beyond 75 years.

A suitable assessment would report the number and proportion of external subjects older than 75, their age distribution, and the prediction’s dependence on the assumed continuation of the age effect. If the external dataset contains observed concentrations for these subjects, their performance can be evaluated there; that evaluation still does not establish performance for other unsupported ages or settings.

Why Joint Support Can Fail When Ranges Overlap

A marginal range check can miss unsupported combinations. The following example uses a deliberately simple geometric description to make that distinction exact. Suppose the observed development covariate combinations lie within the triangle whose vertices are (50 kg, 90 units) for body weight and renal function, (90 kg, 90 units), and (90 kg, 50 units). The triangle can be described by the constraints \(50\leq \mathrm{Weight}\leq90\), \(50\leq \mathrm{RenalFunction}\leq90\), and \(\mathrm{Weight}+\mathrm{RenalFunction}\geq140\), with weight in kg and renal function in the stated units.

Worked Example: Both Covariates in Range, but Their Combination Is Not

Consider an external subject with weight 60 kg and renal function 60 units. Each value lies within its development range of 50 to 90. However, \(60+60=120<140\), so the combination is outside the development triangle. It is not represented by any convex combination of the three stated vertices.

Suppose the locked model is \( \log(CL_i)=\log(5.00)+0.010(\mathrm{Weight}_i-70)-0.005(\mathrm{RenalFunction}_i-70) \). For the external subject, substitution gives \( \log(CL_i)=\log(5.00)+0.010(60-70)-0.005(60-70)=\log(5.00)-0.100+0.050=\log(5.00)-0.050 \). Therefore \(CL_i=5.00\exp(-0.050)=4.76\) L/h.

The arithmetic is defined, but the subject’s joint covariate profile lies outside the described development region. The estimated effects of weight and renal function, and any implied behavior of the model for this combination, are not directly supported by development observations. The result should be flagged as joint extrapolation rather than described as an ordinary within-range prediction.

The triangle is only a transparent illustration. Real datasets are noisy clouds, not exact geometric boundaries. A point just outside a convex hull may be close to many observations, while a point inside it may be far from almost all observations. Therefore, a hull check should be paired with examination of local data density and relevant clinical combinations. In high dimensions, apparent support can be especially misleading because the number of possible combinations grows rapidly.

Overlap Is Not the Same as the Same Population

Even when all external covariates fall within development ranges, their distributions can differ. Support overlap means that external covariate values are represented to some degree in development data; it does not mean that the two populations have the same covariate frequencies. A development sample may be mostly younger adults and an external sample mostly older adults, while both cover the same overall age interval.

Worked Example: Complete Range Overlap but a Shift in Age Distribution

Suppose a hypothetical development dataset has 100 subjects: 30 aged 20–39 years, 40 aged 40–59 years, and 30 aged 60–79 years. A hypothetical external dataset has 80 subjects: 8, 24, and 48 in those same age bands, respectively. Assume all recorded external ages fall within the development minimum and maximum.

The development proportions are \(30/100=0.300\), \(40/100=0.400\), and \(30/100=0.300\). The external proportions are \(8/80=0.100\), \(24/80=0.300\), and \(48/80=0.600\). Thus, the oldest band accounts for 60.0% of the external sample but 30.0% of the development sample.

There is no observed age-range violation in this example, but the external population places much more weight on older subjects. Overall external performance may consequently reflect the model’s performance in that age band more strongly than development performance did. Report the shift alongside the performance results and consider the earlier material on subgroup performance when interpreting results by age. The shift alone does not show that the conditional PK model is wrong; it shows that population composition differs.

A useful conceptual distinction is between conditional prediction behavior and the covariate distribution over which performance is summarized. Let \(X\) denote the relevant covariates, let \(Q\) be their distribution in the external population, and let \(L\) be a prediction loss (a measure of prediction error). The external population’s mean loss can be written as

$$ E_Q[L]=\int_{\mathcal{S}_D} E[L\mid X=x]\,dQ(x) +\int_{\mathcal{S}_D^c} E[L\mid X=x]\,dQ(x), $$

where \(\mathcal{S}_D\) denotes the region with development support. The first term averages conditional prediction loss over the part of the external covariate distribution within development support. The second concerns external covariates outside that region. If the external dataset includes outcome data there, it can provide direct evidence about those subjects in that dataset. But development data alone cannot establish the model’s conditional behavior there. If no external observations cover that region either, the corresponding target-population performance is not empirically assessed.

This decomposition also explains why an overall performance estimate is not automatically transportable to every population. Even if conditional prediction behavior is stable, a different distribution \(Q\) changes how much weight each covariate region contributes. If conditional behavior is not stable, overlap alone is insufficient. Transporting predictions therefore requires both adequate support for the target covariates and a defensible expectation that the model’s conditional relationships apply in the target setting.

A Practical Support Assessment

1
Specify the target covariates.
List the model covariates and any prespecified combinations that matter clinically or structurally. Confirm harmonized units, definitions, coding, and measurement timing before comparing datasets.
2
Compare marginal support.
For continuous covariates, report development and external ranges and inspect distribution plots or quantiles. For categorical covariates, report which levels occur and their counts. Flag external values beyond observed development ranges and levels absent from development.
3
Inspect joint and sparse regions.
Use clinically relevant cross-tabulations and plots for important pairs, and assess whether external subjects occupy combinations or tails with little development representation. State the diagnostic used; do not imply that a geometric boundary proves dense support.
4
Relate support to the prediction.
Identify which model terms determine predictions in flagged regions and describe the assumptions needed to extend them. Distinguish a model-generated value from evidence that the value is reliable.
5
Qualify the transportability claim.
Report support gaps, distribution shifts, and which external subjects contribute to the evaluation. Limit conclusions to the covariate regions and settings actually represented; do not silently generalize to unsupported target groups.

Assumptions, Diagnostics, and Limitations

A support-based transportability argument depends on several assumptions. First, the covariates must be measured comparably across datasets; unit errors or inconsistent definitions invalidate a numerical comparison. Second, the development sample must provide a useful representation of the region claimed as supported. A range based on one unusual subject does not mean that the surrounding tail is well characterized. Third, any claim that a fitted covariate relationship applies externally requires that the relevant conditional PK relationship remains sufficiently stable across settings, given the covariates included in the model. This is a substantive assumption, not a consequence of overlap.

Warning signs include external values beyond development minima or maxima, external categories never observed during development, thin development tails, and external combinations absent or rare in development. Also inspect whether a few extreme observations determine the apparent range, whether missingness patterns differ, and whether the covariates were recorded at comparable times. These checks diagnose limitations in the support argument; they do not replace assessment of observed prediction performance.

Support diagnostics have limits. Sample extrema are unstable summaries of population boundaries; small development samples can make genuine population support look absent, while a single extreme observation can make the observed range look broader than it is. Joint-support diagnostics become harder as the number of covariates increases. A convex hull, standardized distance, or density estimate is a summary that depends on choices and data size, not a universal pass/fail test. Report practical details, such as flagged counts and the covariate regions involved, rather than relying on an unexplained label of “overlap.”

The main alternative to a support-focused assessment is to summarize external prediction performance alone. That directly answers how predictions performed for the external subjects observed, but it can obscure where those subjects lie relative to development data and cannot establish performance for unobserved target regions. Conversely, a covariate-overlap analysis describes the populations but does not itself measure prediction error. These approaches answer complementary questions and should not be treated as substitutes. Re-estimating or revising the model to address a support gap changes the evaluation question; support assessment of a locked model should make the gap visible rather than conceal it through model updating.

Common Mistakes

  • “Every external value is between the development minimum and maximum, so there is no extrapolation.” This checks only marginal ranges. Joint combinations can be unsupported, and values in sparse tails may have weak empirical support.
  • “An unsupported prediction is necessarily wrong.” A prediction can happen to be accurate outside observed development support. The correction is that its accuracy is not justified by development observations there; direct external data may evaluate it for the observed subjects.
  • “The two samples overlap, so their performance should transport unchanged.” Overlap does not imply identical covariate distributions or stable conditional relationships. Both composition and model behavior matter.
  • “The model returned a finite value, so it is supported.” Numerical calculability is not empirical support. State which covariates or combinations are extrapolated and which assumptions extend the model there.
Key takeaway: Assess support in the relevant joint covariate space, not only by separate ranges. Support gaps identify where development data provide limited evidence for PK predictions; distribution shifts within overlap also affect how external performance should be interpreted.

Check Your Understanding

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

  1. Development ages range from 18 to 70 years. An external subject is 74 years old. Is this a marginal support gap, and what does that imply?
  2. Development body weights range from 45 to 100 kg and renal function ranges from 30 to 120 units. An external subject has weight 50 kg and renal function 35 units. Why can marginal ranges alone not establish support?
  3. Development and external datasets have the same age range, but older adults make up a larger proportion of the external dataset. Does this establish extrapolation? What should be reported?
  4. A model gives predictions for an external covariate combination not observed in development. Does that prediction establish transportability to all subjects with that combination?

Solutions. (1) Yes. Since \(74>70\), the subject is beyond the observed development age range by \(74-70=4\) years. The age-related prediction requires extrapolation; this is a limitation on development-based evidence, not proof that the prediction is inaccurate. Report how many external subjects are affected and the model relationship being extended. (2) Both values are within their separate ranges, but that does not show the pair occurs in development. Compare joint combinations or relevant multivariate support, and inspect how densely nearby development subjects are represented. (3) No marginal age extrapolation is established if the observed ranges overlap, but the covariate distribution has shifted. Report the age-band frequencies or another distribution summary and interpret overall performance in light of the greater external weight on older subjects. (4) No. A computed prediction is not evidence that it generalizes to all subjects with that combination. It may be evaluated for external subjects with observed outcomes, but broader transportability also requires relevant support and a defensible assumption that conditional PK relationships remain applicable.

Clinical Biostats