Tutorials › Model-Informed Drug Development (MIDD) › Choosing Regimens Across Population Variability

Dose Optimization and Titration · Tutorial 260 of 596

Choosing Regimens Across Population Variability

Learn how exposure variability and explicit patient or subgroup weights can change which regimen is preferred, even when average exposures appear similar.

Advanced 12 min read

What You'll Learn

  • Describe predicted exposure distributions and calculate the probability of falling within a target range.
  • Distinguish a population’s subgroup proportions from additional decision weights assigned to patients or groups.
  • Compare regimens by weighted target attainment and by expected utility for patient outcomes.
  • Identify assumptions and diagnostics that affect population-level regimen comparisons.
  • Explain why mean exposure or unweighted attainment alone can conceal clinically important differences.

Why Compare Regimens Across Patients?

A regimen that produces a desirable exposure for the “average” patient may miss the target for many individuals. Clearance, body size, organ function, adherence, and other sources of variability spread predicted exposures across a population. The distribution matters: two regimens with the same mean exposure can have different proportions of patients below, within, or above a clinically relevant range.

A population comparison also requires a clear answer to whose outcomes count, and how much? A population-weighted average gives each patient the same weight, so groups contribute according to their prevalence. A decision-maker may instead assign explicit weights to groups or outcomes to reflect a stated clinical priority. Those weights should be visible in the analysis rather than hidden in a choice of summary statistic.

Definition: For a regimen \(r\), let \(C_r\) denote a patient’s predicted exposure under that regimen, with the distribution representing variability across the target population. For an exposure target \([L,U]\), target attainment is \(P(L\leq C_r\leq U)\). A population regimen comparison summarizes predicted exposures or outcomes across patients, using stated population proportions and, when appropriate, explicit decision weights.

From an Exposure Distribution to Target Attainment

Suppose predicted exposure is lognormally distributed under a regimen: \(\ln C_r\sim N(\mu_r,\sigma_r^2)\). Here \(\mu_r\) and \(\sigma_r\) are the mean and standard deviation of log exposure, respectively. The median exposure is \(e^{\mu_r}\), and the geometric standard deviation is \(e^{\sigma_r}\). This is a modeling assumption, not a universal property of exposure data.

For a target range with \(0<L<U\), standardizing the log exposure gives the target-attainment probability. The calculation uses the standard normal cumulative distribution function \(\Phi\).

$$ P(L\leq C_r\leq U) = \Phi\!\left(\frac{\ln U-\mu_r}{\sigma_r}\right) - \Phi\!\left(\frac{\ln L-\mu_r}{\sigma_r}\right). $$

The two standardized limits locate the target boundaries in the distribution of log exposure. The difference between their cumulative probabilities is the fraction of the modeled population predicted to fall inside the target range. It is not the probability that a particular patient will attain the target after that patient’s exposure is known.

Worked Example: Target Attainment Under Two Regimens

Consider two hypothetical once-daily regimens and an exposure target of 60–120 exposure units. Assume the predicted exposure is lognormal for each regimen. Regimen A has median exposure 80 units and geometric standard deviation 1.5; regimen B has median exposure 100 units and geometric standard deviation 1.8. Which has higher predicted target attainment?

For A, \(\mu_A=\ln(80)\) and \(\sigma_A=\ln(1.5)\). Its standardized lower and upper limits are:

$$ z_L=\frac{\ln(60/80)}{\ln(1.5)}\approx-0.7095, \qquad z_U=\frac{\ln(120/80)}{\ln(1.5)}=1.000. $$

Therefore,

$$ P(60\leq C_A\leq120) = \Phi(1.000)-\Phi(-0.7095) \approx 0.84134-0.23900 =0.60234. $$

Regimen A has approximately 60.2% target attainment in the modeled population.

For B, \(\mu_B=\ln(100)\) and \(\sigma_B=\ln(1.8)\). The standardized limits are:

$$ z_L=\frac{\ln(60/100)}{\ln(1.8)}\approx-0.8691, \qquad z_U=\frac{\ln(120/100)}{\ln(1.8)}\approx0.3102. $$

Thus,

$$ P(60\leq C_B\leq120) \approx 0.62179-0.19241 =0.42938. $$

Regimen B has approximately 42.9% target attainment. Under the stated target and distributional assumptions, A has higher attainment by about 17.3 percentage points. This comparison alone does not establish that A is clinically preferable: the target’s relevance and the consequences of exposures outside it must also be justified.

Population Weights and Decision Weights

When patients belong to distinct groups, estimate or predict attainment within each group before combining results. Let \(p_g\) be the proportion of the target population in group \(g\), and let \(a_{rg}\) be the predicted target-attainment probability for regimen \(r\) in that group. The population-average attainment is \(\sum_g p_g a_{rg}\), where the population proportions sum to one.

Population proportions answer a descriptive question: how would the regimen perform in the specified population mix? Decision weights answer a different question: how much should each group’s result count in a decision? If \(w_g\) are normalized decision weights, with \(\sum_g w_g=1\), a decision-weighted summary is \(\sum_g w_g a_{rg}\). It should not be described as the population-average attainment unless \(w_g=p_g\).

Worked Example: How Group Weights Change the Comparison

In a hypothetical population, 70% of patients are in group 1 and 30% in group 2. Predicted target attainment is 70% and 40% under regimen X, and 60% and 80% under regimen Y, respectively. First compare population-average attainment.

For X, the population-weighted attainment is \(0.70(0.70)+0.30(0.40)=0.49+0.12=0.61\), or 61%. For Y, it is \(0.70(0.60)+0.30(0.80)=0.42+0.24=0.66\), or 66%. In this population mix, Y is higher by 5 percentage points.

Now suppose the decision explicitly gives the two groups equal weight, so \(w_1=w_2=0.50\). The decision-weighted summaries are \(0.50(0.70)+0.50(0.40)=0.55\) for X and \(0.50(0.60)+0.50(0.80)=0.70\) for Y. Y is higher by 15 percentage points under these decision weights.

The second calculation is not a revised estimate of attainment in the stated population. It reflects a different decision rule that places equal weight on the groups despite their unequal prevalence. Reporting both the weights and the interpretation prevents this change in question from being mistaken for a change in the predictions.

Weighting Patient Outcomes, Not Just Target Attainment

A target indicator assigns the same score to every patient within the target and the same score to every patient outside it. That can be too crude if outcomes differ in value or severity. Expected utility provides a general way to summarize predicted outcomes using explicit scores. For each patient or group, define a utility \(u(y)\) for outcome \(y\); the expected utility is the probability-weighted average of those scores. A population-level expected utility then averages across patients using the intended population or decision weights.

$$ EU_r=\sum_g w_g\sum_k P(Y_r=y_k\mid g)\,u(y_k), \qquad \sum_g w_g=1. $$

The scale and direction of utility must be stated. For example, a score of 1 might represent a substantial symptom improvement and 0 no improvement. Such scores are decision inputs, not facts produced by the exposure model. Their clinical rationale and sensitivity should be examined. This framework can summarize a single outcome dimension; it does not by itself resolve trade-offs among multiple objectives.

Worked Example: Comparing Expected Outcome Utility

Suppose two hypothetical regimens are predicted to produce three symptom outcomes: no improvement, moderate improvement, or substantial improvement. Assign utilities 0, 0.5, and 1, respectively. Group 1 makes up 60% of the population and group 2 makes up 40%. The predicted outcome probabilities are as follows.

Regimen and groupNo improvementModerate improvementSubstantial improvement
A, group 10.200.500.30
A, group 20.400.400.20
B, group 10.100.400.50
B, group 20.200.500.30

For A in group 1, expected utility is \(0.20(0)+0.50(0.5)+0.30(1)=0.55\). In group 2 it is \(0.40(0)+0.40(0.5)+0.20(1)=0.40\). Population-weighting these gives \(EU_A=0.60(0.55)+0.40(0.40)=0.49\).

For B in group 1, expected utility is \(0.10(0)+0.40(0.5)+0.50(1)=0.70\). In group 2 it is \(0.20(0)+0.50(0.5)+0.30(1)=0.55\). Therefore, \(EU_B=0.60(0.70)+0.40(0.55)=0.64\). On this stated utility scale and population mix, B has higher expected utility by 0.15 utility units. That conclusion depends on the predicted probabilities and the chosen utilities; it is not a claim that every patient benefits more under B.

Justification, Assumptions, and Diagnostics

The weighted summaries follow from averaging patient-level outcomes. If patient \(i\) has an outcome score \(U_i\), the finite-population average is \(N^{-1}\sum_{i=1}^N U_i\). If patients are represented by groups with proportions \(p_g\), and the model gives the group-specific mean score \(E(U\mid g,r)\), then partitioning the population by group gives:

$$ E(U\mid r) = \sum_g P(g)E(U\mid g,r) = \sum_g p_g E(U\mid g,r). $$

The same reasoning applies to target attainment by setting \(U_i=1\) when exposure is in the target and \(U_i=0\) otherwise. Replacing \(p_g\) with decision weights \(w_g\) is a separate normative choice, not a consequence of the population-average identity.

A defensible comparison depends on several conditions. The target population and its group proportions should be relevant to the decision. Exposure predictions should represent the intended regimen, dosing conditions, and sources of variability; the exposure-response relation and target range should be clinically justified. If a lognormal model is used, its fit should be checked rather than assumed. The comparison should also make clear whether uncertainty in model parameters and population proportions is represented or whether predictions are conditional on fixed estimates.

Useful diagnostics include comparing predicted exposure distributions with observed data, checking whether residuals or subgroup patterns suggest a poor distributional fit, and plotting the full exposure distributions alongside target boundaries. Compare predicted and observed target-attainment summaries where suitable data exist. Inspect subgroup-specific results: an overall mean can hide a poorly served group. Repeat the calculation under plausible alternative population proportions, target ranges, and utility scores. Large changes in the ranking signal that the decision depends on those choices and should be reported.

1
Define the decision population.
Specify the patient groups, their prevalence, and whether the goal is population-average performance or a differently weighted decision summary.
2
Predict distributions and outcomes.
For each regimen, obtain the exposure distribution and, where needed, group-specific probabilities of clinically meaningful outcomes.
3
Apply stated targets and weights.
Calculate group-specific attainment or expected utility, then combine using the declared population or decision weights.
4
Check sensitivity and interpret in context.
Assess model fit and variation in weights, targets, and utilities; state what the comparison does and does not establish.

Alternatives, Limitations, and a Common Pitfall

Comparing mean or median exposure is simpler and can be useful for describing central tendency, but it does not show how much of the population lies within a target range. Target attainment is directly interpretable when the range is clinically justified, but it reduces a continuous exposure distribution to an inside-or-outside indicator. Expected utility retains distinctions among outcomes if meaningful scores are available, but its result depends on those scores. Reporting the distribution, subgroup results, and a weighted summary together is often more informative than relying on any one measure.

Model-based predictions are not direct observations for every patient. They can be sensitive to distributional assumptions, covariate representation, adherence assumptions, and uncertainty in estimated parameters. A high average attainment can coexist with low attainment in a subgroup, while a favorable decision-weighted score can reflect a deliberate priority rather than the population’s actual composition. These summaries support regimen comparison; they do not establish causal benefit beyond the evidence and assumptions underlying the models.

  • Pitfall: choosing the regimen with the closest mean exposure to the target midpoint. That reasoning ignores distribution width and tails. A regimen centered at the midpoint may still have low attainment if between-patient variability is large. Calculate target attainment from the full predicted distribution.
  • Pitfall: calling a decision-weighted result “population average.” If weights differ from population proportions, the summary answers a different question. Label the weights and state whose outcomes receive extra emphasis.
  • Pitfall: treating utility scores as model outputs. The exposure or outcome model predicts probabilities; decision-makers supply the utility scale. Justify the scores and test whether the regimen ranking changes under plausible alternatives.
Key takeaway: A population regimen comparison should connect the predicted exposure distribution to clinically meaningful attainment or outcomes, then combine patient-level predictions using transparent population proportions or explicitly justified decision weights.

Check Your Understanding

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

  1. Why can two regimens with the same median exposure have different target-attainment probabilities?
  2. A population has 80% group 1 and 20% group 2. A regimen’s group-specific attainments are 0.50 and 0.90. Calculate population-average attainment.
  3. Using the same group attainments, calculate the decision-weighted summary if both groups receive equal weight. Is it the population average?
  4. What does an expected utility of 0.64 mean in the worked example, and what does it not mean?

Solutions. (1) Their distribution widths or shapes may differ, so different proportions can fall between the target boundaries despite equal medians. (2) \(0.80(0.50)+0.20(0.90)=0.40+0.18=0.58\), or 58%. (3) \(0.50(0.50)+0.50(0.90)=0.70\), or 70%. It is not the population average because the equal weights differ from the population proportions. (4) It is the population-weighted mean of the stated utility scores under regimen B and the hypothetical outcome probabilities. It is not a 64% probability of benefit, nor does it mean every patient has utility 0.64.