From Dose-Ranging Data to a Phase II Dose Choice
A phase II dose-ranging trial should inform more than whether a drug has activity. It should help identify a dose, or a plausible range, whose expected benefit is large enough to matter without choosing a dose solely because it produced the largest observed mean. A model can connect the randomized dose groups to exposure and response, interpolate between studied groups, and quantify uncertainty. The model does not replace the observed group results: a useful dose decision checks whether its predictions are consistent with those results and makes the target and uncertainty criteria explicit.
The case below is entirely hypothetical. The endpoint is change from baseline on a symptom scale at week 12, with larger values representing greater improvement. Effects are reported as placebo-adjusted mean improvements, in scale points. A clinically relevant target has been set at an improvement of at least 8 points. The analysis is about selecting a dose based on efficacy and exposure evidence; it does not establish that a dose is acceptably safe or appropriate as a final regimen.
Hypothetical Trial Data and Exposure-Response Model
Suppose participants were randomized equally to placebo or one of four active doses. The endpoint is measured at baseline and week 12; the table shows the placebo-adjusted week-12 mean improvement. Each group has 40 participants. The exposure column gives the group's mean steady-state average concentration, \(C_{\mathrm{avg}}\), in mg/L. The dose-to-exposure relationship in this hypothetical trial is approximately \(C_{\mathrm{avg}}=0.20 \times \text{dose}\), with dose in mg.
| Assigned dose (mg) | Mean \(C_{\mathrm{avg}}\) (mg/L) | Observed mean improvement (points) | SE of observed mean (points) |
|---|---|---|---|
| Placebo | 0 | 0.0 | 0.45 |
| 5 | 1 | 3.8 | 0.52 |
| 10 | 2 | 5.7 | 0.58 |
| 20 | 4 | 7.9 | 0.62 |
| 40 | 8 | 9.0 | 0.70 |
For illustration, an exposure-response analysis fit a saturable \(E_{\max}\) model to the individual-level efficacy and exposure data. The fitted parameter estimates are \(E_{\max}=12\) points and \(EC_{50}=2\) mg/L. This model describes an increasing response that approaches a plateau as exposure increases. It is being used here as a compact description of the hypothetical data, not as a universal biological law.
Here, \(E(C)\) is the model-predicted placebo-adjusted mean improvement at exposure \(C\). \(E_{\max}\) is the limiting response as exposure becomes very large, and \(EC_{50}\) is the exposure at which the model predicts half of \(E_{\max}\). For this model to be appropriate, response must be adequately represented by a monotone saturating curve over the observed range, and the exposure and endpoint measurements must be suitable for the stated interpretation. A fitted curve also assumes that important differences between participants or groups have not made the common exposure-response relationship misleading.
Worked Example: Predicting Responses at the Studied Exposures
At the 20 mg dose, the mean exposure is 4 mg/L. Substituting the fitted estimates into the model gives
The model predicts an 8.00-point mean improvement at that exposure. At 40 mg, the mean exposure is 8 mg/L:
Doubling exposure from 4 to 8 mg/L raises the predicted response by only 1.60 points, because the curve is approaching its plateau. These predictions are for mean exposure values. If exposure varies substantially within a dose group, a more appropriate group prediction averages \(E(C)\) over that group's exposure distribution; in general, that average need not equal \(E\) evaluated at the group's mean exposure.
Finding the Exposure Associated With the Target
A model can translate a response target into an exposure target. To find the exposure at which the fitted curve reaches 8 points, solve \(E(C)=8\). For a target \(T\) below \(E_{\max}\), rearranging the model gives \(C=T EC_{50}/(E_{\max}-T)\). A finite positive solution requires \(0<T<E_{\max}\); a target at or above the plateau cannot be reached at finite exposure under this model.
Given the hypothetical mean exposure relationship \(C_{\mathrm{avg}}=0.20 \times \text{dose}\), a mean exposure of 4 mg/L corresponds to \(4/0.20=20\) mg. This is a model-based mapping, not a guarantee that every participant assigned 20 mg will have an exposure of 4 mg/L or an 8-point response. Exposure varies among participants, and the dose-to-exposure relationship is estimated rather than exact.
Worked Example: Comparing Model Predictions With Observed Groups
At 10 mg, the observed mean improvement is 5.7 points and the mean exposure is 2 mg/L. The model prediction is
The observed-minus-predicted difference is \(5.7-6.0=-0.30\) points. At 20 mg, the prediction is 8.00 points and the observed mean is 7.9, giving \(7.9-8.0=-0.10\) points. At 40 mg, the prediction is 9.60 points while the observed mean is 9.0, giving \(9.0-9.6=-0.60\) points.
These differences are small relative to the response scale and do not show a dramatic mismatch in this compact example. That does not prove the model is correct. In a real analysis, inspect the pattern across all dose groups and participants, and check whether the fitted curve misses systematic features such as a downturn at high exposure, a delayed response, or distinct patterns among relevant participant subgroups. The observed means are also uncertain; their SEs describe sampling uncertainty in those means, not uncertainty in the model's fitted predictions.
Using Uncertainty to Screen Candidate Doses
A dose should not be declared adequate just because its point prediction reaches the target. The model-based uncertainty for each predicted mean should also be considered. Suppose the fitted model's covariance calculation gives prediction SEs of 0.55, 0.62, and 0.78 points at 10, 20, and 40 mg, respectively. These are SEs for predicted mean responses, including parameter uncertainty; they are not the observed-group SEs in the data table. For a prediction \(E(C)\), the delta-method approximation uses the gradient of the model with respect to its fitted parameters and their covariance matrix \(V\):
The gradient for parameters \((E_{\max},EC_{50})\) is \(\left(C/(EC_{50}+C),\ -E_{\max}C/(EC_{50}+C)^2\right)\). In practice, the covariance calculation must reflect the fitted model and the data; an SE should not be inferred from the curve's shape alone.
If three candidate doses are screened and the aim is to make a one-sided 90% family-wise lower-bound statement, a Bonferroni adjustment allocates a tail probability of \(0.10/3\) to each dose. Under a normal approximation, the one-sided critical quantile is \(z_{1-0.10/3}\approx1.834\). For each candidate, calculate the lower bound as prediction minus 1.834 times its prediction SE. Bonferroni's union-bound argument does not require the three predictions to be independent: the probability that any of the three bounds fails to cover its corresponding mean is no greater than the sum of the three individual failure probabilities, which is at most \(3(0.10/3)=0.10\).
Worked Example: Applying Simultaneous Lower Bounds
The three candidate model predictions are 6.00 points at 10 mg, 8.00 points at 20 mg, and 9.60 points at 40 mg. Their supplied prediction SEs are 0.55, 0.62, and 0.78 points. The Bonferroni one-sided 90% family-wise lower bounds are:
Rounded to two decimal places, the bounds are 4.99, 6.86, and 8.17 points. Of these three candidates, only 40 mg has a lower bound above the prespecified 8-point target. Under the model and normal-approximation assumptions, the Bonferroni procedure provides at least 90% simultaneous coverage for the three corresponding true mean predictions. This is not a 90% probability that a particular parameter or mean lies in its interval; it is a repeated-sampling coverage statement about the procedure.
This criterion is deliberately conservative. The model's point prediction at 20 mg meets the target, but its lower bound does not. If the decision rule requires the simultaneous lower bound to clear the target, the hypothetical data support 40 mg as the efficacy candidate among these doses. A different prespecified rule might allow selection based on a point estimate or a less conservative uncertainty criterion; the resulting choice should be described accordingly.
Assumptions, Diagnostics, and Alternative Summaries
The dose recommendation is only as credible as the analysis behind it. Randomization supports comparisons of assigned-dose groups, provided allocation was implemented appropriately and endpoint data are handled without bias. Exposure is measured after randomization, however, and participants with different exposures may also differ in characteristics related to response. An exposure-response association alone is therefore not automatically a causal effect of changing exposure. The randomized dose-group results remain an important check on the model-based exposure interpretation.
The \(E_{\max}\) fit assumes an increasing saturating relationship with no downturn in the studied range. It also assumes that the chosen week-12 endpoint adequately captures the relevant response and that the model accounts for the exposure variation relevant to prediction. If efficacy changes over time, a week-12 summary can conceal an earlier benefit, delayed effect, or loss of response. The analysis should use a consistent endpoint definition and account for missing observations and departures from the planned treatment or measurement schedule.
Diagnostics should ask whether the model is failing in ways that matter to the dose decision. Plot observed group means and uncertainty alongside the predictions; inspect individual-level residual patterns against exposure and time; and look for curvature, systematic overprediction at high exposure, influential observations, or subgroup patterns. Compare alternative plausible curve shapes when justified by the data. A pattern of increasingly negative observed-minus-predicted differences at high exposure, for example, could suggest that the assumed curve or its exposure distribution is inadequate. A visually close fit to five group means is not by itself sufficient evidence, especially when the model was fitted to those same data.
| Approach | What it contributes | Main limitation for dose selection |
|---|---|---|
| Observed randomized dose-group means | Directly summarize outcomes by assigned dose without requiring an exposure-response curve. | Only studied doses are directly represented; group means may be imprecise and do not alone interpolate between doses. |
| Exposure-response model | Connects individual exposure to response and can predict at exposures between observed values. | Predictions depend on model form, exposure measurement, and assumptions about the exposure-response relationship. |
| Model-based predictions checked against randomized groups | Combines a structured exposure-response description with a check against observed treatment-group behavior. | Agreement in one trial does not remove uncertainty, confounding concerns for exposure, or the need for suitable diagnostics. |
The strongest practical interpretation uses these approaches together rather than treating them as competitors. The observed randomized means show what happened at the tested doses; the exposure-response model helps explain the pattern and estimate response at relevant exposures; and uncertainty bounds make the decision rule visible. If these sources disagree materially, investigate the disagreement rather than selecting whichever result supports a preferred dose.
Common Pitfall and Decision Limits
A common mistake is to select the dose with the largest observed mean and call it optimal. The highest mean can arise partly from sampling variability, and a modest additional response at the highest dose may not justify increased exposure. The correction is to define a clinically meaningful target and an uncertainty criterion in advance, compare observed results with model predictions, and state clearly what the evidence does and does not establish. Here, 40 mg is the only candidate whose simultaneous lower bound clears the target; this supports an efficacy-based candidate under the stated rule, not a complete clinical recommendation.
The target itself is a decision input, not something the fitted curve determines. A different target, an individual-responder criterion, or a different uncertainty requirement could change the selected dose. In addition, efficacy alone cannot establish benefit-risk balance, and a mean response target does not guarantee that an individual patient benefits. These limitations matter when translating a phase II result into later development decisions.
Define the endpoint, time point, response direction, clinically relevant threshold, and uncertainty rule before comparing candidate doses.
State the model assumptions, compare predictions with randomized dose-group results, and investigate systematic residual patterns.
Use uncertainty bounds as prespecified, and describe the result as an efficacy-based candidate rather than a safety or final regimen conclusion.
Check Your Understanding
Answer each question using the hypothetical case and the stated model.
- What response does the model predict at \(C=1\) mg/L, and how does it compare with the observed mean at 5 mg?
- What exposure corresponds to a target of 6 points? Under the stated dose-to-exposure relationship, what dose corresponds to that exposure?
- Why is the model prediction at mean exposure not necessarily the same as the mean model-predicted response for a dose group?
- For the 20 mg candidate, calculate the Bonferroni lower bound using its prediction SE of 0.62 points and explain whether it clears the 8-point target.
- Why does an exposure-response association not, by itself, prove that increasing exposure causes greater efficacy?
Solutions. (1) \(E(1)=12(1)/(2+1)=4.00\) points, compared with the observed 3.8 points at 5 mg; the observed-minus-predicted difference is \(-0.20\) points. (2) Solve \(6=12C/(2+C)\): \(12+6C=12C\), so \(C=2\) mg/L; at 0.20 mg/L per mg, this corresponds to \(2/0.20=10\) mg. (3) The response curve is nonlinear, so in general the average of \(E(C)\) over individual exposures differs from \(E\) at the average exposure. (4) \(8.00-1.834(0.62)=8.00-1.13708=6.86292\), or 6.86 points; it does not clear the 8-point target. (5) Exposure is measured after randomization and may be associated with participant characteristics that also affect response, so the association may not represent the causal effect of changing exposure.