Tutorials › Adaptive Clinical Trial Design › Covariate-Conditional Dose Adaptation

Adaptive Designs With Pharmacometric Models · Tutorial 514 of 799

Covariate-Conditional Dose Adaptation

Learn how to translate covariate-dependent pharmacometric predictions into a prespecified, clinically interpretable dose rule.

Advanced 11 min read

What You'll Learn

  • Specify a conditional exposure model and the assumptions needed to use it for dose selection.
  • Calculate the probability that exposure falls within a target window at a candidate dose.
  • Translate body size, organ function, or genotype into covariate-conditional predictions.
  • Combine target attainment, predicted response, and safety constraints in a dose rule.
  • Identify model diagnostics, implementation risks, and limitations of covariate-based dosing.

Why Make Dose Decisions Conditional on Covariates?

A dose that gives suitable exposure for a typical patient may underexpose some patients and overexpose others. Body size, organ function, and genotype can help explain these differences. In an adaptive trial, the practical aim is not to create a separate informal decision for every participant. It is to define a prespecified rule that maps measured covariates to a dose using predictions from a pharmacometric model.

The earlier material on selecting dose regimens by predicted exposure considers exposure-based choices. Here the choice is conditional on patient characteristics: the same candidate dose may be suitable for one covariate profile and unsuitable for another. The rule should use covariates measured before the dose decision, a clinically justified target, and predictions that account for relevant variability.

Definition: A covariate-conditional dose adaptation rule is a prespecified function that maps a patient’s observed covariates \(\mathbf{x}\), together with the available trial information, to a dose or regimen \(D(\mathbf{x})\). The function is constructed from a model for exposure and, when relevant, response or toxicity. It must specify the candidate doses, decision criteria, and what to do if no candidate satisfies those criteria.

Conditional Exposure and Target Attainment

Let \(A\) denote an exposure measure, such as area under the concentration-time curve (AUC), and let \(D\) be dose. For illustration, suppose that, conditional on covariates \(\mathbf{x}\), exposure follows a lognormal model:

$$ \log A \mid D,\mathbf{x} \sim N\!\left(\log m(D,\mathbf{x}),\sigma^2\right). $$

Here \(m(D,\mathbf{x})\) is the conditional median exposure and \(\sigma\) is the standard deviation on the log scale. This model requires positive exposure. It also assumes that the specified dose-covariate relationship is adequate for the population and dose range in which the rule will be used. In practice, predictive variability should reflect the individual-to-individual variability relevant to the decision. Parameter uncertainty can also be included in predictions; examples below use a fixed fitted model to keep the calculations transparent.

Suppose the desired exposure window is from \(L\) to \(U\), where \(0<L<U\). The target-attainment probability at dose \(D\) for covariates \(\mathbf{x}\) is the conditional probability that exposure falls inside that window:

$$ T(D,\mathbf{x}) = P(L \leq A \leq U \mid D,\mathbf{x}) = \Phi\!\left(\frac{\log U-\log m(D,\mathbf{x})}{\sigma}\right) - \Phi\!\left(\frac{\log L-\log m(D,\mathbf{x})}{\sigma}\right). $$

A rule might require \(T(D,\mathbf{x})\geq q\), for a prespecified minimum probability \(q\). This is a model-based prediction for a patient profile, not a guarantee that an individual patient’s exposure will be in range.

Worked Example: Calculating Target Attainment for a Reference Profile

Consider hypothetical data summarized by a fitted model with \(\sigma=0.20\). The target exposure window is 80 to 120 mg·h/L. For a particular covariate profile, the predicted median exposure at 75 mg is 97.5 mg·h/L, and exposure is assumed to scale proportionally with dose. Calculate target attainment at 75 mg and compare it with a threshold \(q=0.65\).

At 75 mg, substitute the median and log-scale standard deviation into the target-attainment formula:

$$ T(75,\mathbf{x}) = \Phi\!\left(\frac{\log(120/97.5)}{0.20}\right) - \Phi\!\left(\frac{\log(80/97.5)}{0.20}\right) = \Phi(1.0382)-\Phi(-0.9891) \approx 0.8504-0.1613=0.6891. $$

Thus the predicted probability of exposure within the target window is about 0.6891, or 68.91%, for this profile at 75 mg. It exceeds the prespecified threshold of 0.65.

Under dose-proportional exposure, the predicted median is \(97.5(D/75)\) mg·h/L. At 50 mg, it is 65 mg·h/L, giving

$$ T(50,\mathbf{x}) = \Phi\!\left(\frac{\log(120/65)}{0.20}\right) - \Phi\!\left(\frac{\log(80/65)}{0.20}\right) \approx 0.9989-0.8504=0.1485. $$

At 100 mg, the predicted median is 130 mg·h/L and the target probability is approximately \(0.3369\). At 125 mg, the median is 162.5 mg·h/L and the probability is approximately \(0.06457\); at 150 mg, it is approximately \(0.0076\). Among the candidate doses 50, 75, 100, 125, and 150 mg, only 75 mg meets the requirement \(T(D,\mathbf{x})\geq0.65\).

How Covariates Enter the Rule

A useful model makes the covariate relationship explicit. For example, with standardized covariates, a simple dose-proportional model for the conditional median could be

$$ m(D,\mathbf{x}) = D\exp\!\left(\beta_0+\beta_1x_1+\cdots+\beta_px_p\right). $$

The coefficients describe how the predicted exposure changes with covariates while holding dose fixed. A categorical genotype can be represented by indicator variables; a continuous measure such as body size or renal function can be represented continuously. The model form must match the scientific and data context. A linear effect on log exposure is not automatically appropriate for every covariate.

The following short derivation shows why the exposure target can suggest a useful dose before a discrete candidate set is considered. Suppose \(m(D,\mathbf{x})=D\exp(c(\mathbf{x}))\), with fixed \(\sigma\), and let \(y=\log m(D,\mathbf{x})\), \(a=\log L\), and \(b=\log U\). The target probability as a function of \(y\) is

$$ g(y)=\Phi\!\left(\frac{b-y}{\sigma}\right) -\Phi\!\left(\frac{a-y}{\sigma}\right), \qquad g'(y)=\frac{1}{\sigma} \left[ \phi\!\left(\frac{a-y}{\sigma}\right) -\phi\!\left(\frac{b-y}{\sigma}\right) \right]. $$

The standard normal density \(\phi\) is largest at zero and decreases with distance from zero. When \(y<(a+b)/2\), \(y\) is closer to \(a\) than to \(b\), so the first density in brackets is larger and \(g'(y)>0\). When \(y>(a+b)/2\), the reverse holds and \(g'(y)<0\). Therefore target attainment is maximized when \(y=(a+b)/2\): the predicted median is the geometric midpoint \(\sqrt{LU}\). Under this simple model, the dose attaining that median is \(\sqrt{LU}/\exp(c(\mathbf{x}))\). Real trial decisions still need to respect available doses, safety limits, and response goals.

Worked Example: Adjusting for Lower Organ Function

Suppose, hypothetically, that a standardized covariate \(z\) represents organ function, with \(z=0\) for the reference profile and \(z=-1\) for a profile with lower function. A fitted model gives

$$ m(D,z)=97.5\left(\frac{D}{75}\right)\exp(-0.40z) \quad \text{mg·h/L}. $$

For the lower-function profile at 50 mg, the predicted median is

$$ m(50,-1) = 97.5\left(\frac{50}{75}\right)\exp(0.40) = 65(1.4918) \approx 96.97 \text{ mg·h/L}. $$

Using the same target window, log-scale standard deviation, and threshold as before gives

$$ T(50,-1) = \Phi\!\left(\frac{\log(120/96.97)}{0.20}\right) - \Phi\!\left(\frac{\log(80/96.97)}{0.20}\right) \approx 0.8568-0.1682=0.6886. $$

At 75 mg, the median for this profile is approximately \(145.5\) mg·h/L. Its target probability is about 0.1668, below 0.65. For candidate doses 25, 50, and 75 mg, the 25 mg median is approximately 48.48 mg·h/L and its target probability is about 0.006. Thus 50 mg is the only candidate in this set that meets the target-attainment threshold for this profile. The rule recommends a lower dose than for the reference profile because the model predicts higher exposure at a given dose when \(z\) is lower.

Combining Exposure, Response, and Safety

Exposure targets are useful when a clinically justified exposure window is available, but dose choice may also depend on the probability of response and toxicity. For example, a response model might specify \(P(R=1\mid A,\mathbf{x})\), where \(R\) indicates response. The predicted response probability at a dose is obtained by averaging over the conditional exposure distribution:

$$ \pi_R(D,\mathbf{x}) = \int P(R=1\mid A=a,\mathbf{x})f_A(a\mid D,\mathbf{x})\,da. $$

A corresponding model can provide predicted toxicity probability \(\pi_T(D,\mathbf{x})\). A transparent rule can first require target attainment and a safety constraint, then choose among eligible doses to optimize a prespecified benefit measure. The response and toxicity predictions should be conditional on the same covariate profile and relevant dose regimen.

Worked Example: Applying Exposure and Safety Constraints

For a hypothetical reference profile, compare 75 mg and 80 mg. The exposure model above gives \(T(75,\mathbf{x})=0.6891\). At 80 mg the median is \(97.5(80/75)=104\) mg·h/L, so

$$ T(80,\mathbf{x}) = \Phi\!\left(\frac{\log(120/104)}{0.20}\right) - \Phi\!\left(\frac{\log(80/104)}{0.20}\right) \approx 0.7628-0.0948=0.6680. $$

Both doses meet the target threshold of 0.65. Suppose a hypothetical fitted response model, averaged over exposure, predicts response probabilities of 0.55 at 75 mg and 0.60 at 80 mg. Suppose the corresponding toxicity probabilities are 0.04 and 0.07. The protocol requires toxicity probability no greater than 0.10 and defines a benefit score \(U=100\pi_R-250\pi_T\). Both doses pass the safety constraint. Their scores are

$$ U(75,\mathbf{x})=100(0.55)-250(0.04)=55-10=45, \qquad U(80,\mathbf{x})=100(0.60)-250(0.07)=60-17.5=42.5. $$

Under this specifically stated rule, 75 mg is selected: it meets the exposure and safety requirements and has the higher benefit score. The weights in this score are hypothetical and are not universal clinical values; in a trial they require clinical justification and prespecification.

Building, Checking, and Interpreting the Rule

A practical decision algorithm should be fully specified before its use. One orderly approach is:

1
Define the decision inputs.
Specify the covariates, how they are measured and coded, the candidate doses, and the exposure and outcome models.
2
Calculate conditional predictions.
For each relevant covariate profile and candidate dose, calculate target attainment and any response or toxicity probabilities required by the rule.
3
Apply constraints and choose.
Retain doses satisfying prespecified efficacy and safety criteria, apply the stated selection objective, and specify the fallback if no dose qualifies.
4
Check performance before use.
Evaluate the complete rule over plausible covariate distributions and model conditions, and confirm that the proposed adaptations can be implemented consistently.

Diagnostics should assess whether predicted exposure agrees with observed exposure across dose levels and covariate ranges, not only overall. Check residual patterns, influential observations, adequate representation of important covariate values, and whether the dose range used for the rule is supported by the data. If predictions are systematically too high for one genotype or too low at one end of organ function, a pooled average fit can conceal a consequential problem. Sparse data at extreme covariate values make conditional predictions particularly uncertain.

The main alternatives answer different questions. A fixed dose rule is simpler but does not respond to predicted covariate-related differences. A stratified rule assigns categories, which can be operationally clear but loses information when a continuous covariate is divided into broad groups. A continuous conditional model can preserve that information, but depends more directly on the correctness and support of its functional form. None of these approaches removes the need to justify the exposure target or assess the whole decision rule.

Important limitations include measurement error in covariates, omitted predictors, changes in adherence or concomitant treatment, and uncertainty in model parameters. A model fitted over a narrow range should not be used to make confident predictions far outside that range. If no candidate dose meets the criteria, the rule should not silently select the “closest” dose: it should specify a safe fallback, such as withholding assignment pending review or using a prespecified conservative regimen. Operational instructions must also make clear which measurements are available at the decision time.

Common Mistake

A tempting but incorrect argument is: “The average predicted exposure is in range, so this dose is suitable for every covariate group.” An average can conceal underexposure in one profile and overexposure in another. The correction is to evaluate the conditional predictive distribution and decision criteria for the covariate profiles that matter. Conversely, a covariate association alone does not establish that a different dose improves clinical outcomes; the exposure-response and safety assumptions still need support.

Key takeaway: A covariate-conditional dose rule is a model-based mapping from measured patient characteristics to a protocol-specified dose. Its credibility depends on conditional predictions, justified clinical criteria, supported covariate ranges, and a complete fallback rule.

Check Your Understanding

Answer each question, then compare with the solutions below.

  1. With a lognormal exposure model, what does \(m(D,\mathbf{x})\) represent: the mean, median, or a guaranteed exposure?
  2. For a target window from 80 to 120 mg·h/L and log-scale standard deviation 0.20, a dose gives a predicted median of 100 mg·h/L. Write the target-attainment probability formula for this case.
  3. Why can a model using broad organ-function categories produce a different rule from one using a continuous organ-function measure?
  4. What should a protocol specify if no candidate dose meets the target and safety criteria?

Solutions. (1) It is the conditional median exposure; it is neither a guarantee nor, generally, the arithmetic mean. (2) The probability is \(\Phi(\log(120/100)/0.20)-\Phi(\log(80/100)/0.20)\), or approximately \(\Phi(0.9116)-\Phi(-1.1157)\). (3) Categorizing a continuous covariate discards within-category differences and can create abrupt boundaries, while a continuous model uses the measured value but relies on its chosen functional form. (4) It should state an explicit fallback, such as withholding assignment for review or using a prespecified conservative regimen; the algorithm should not improvise a dose.

Clinical Biostats