Why dose-range boundaries need special attention
A BLRM can calculate a posterior toxicity probability at any positive dose. That mathematical ability does not mean the calculation is equally well supported everywhere. At the lowest or highest candidate dose, there may be no candidate dose on one side to help contextualize a recommendation. More importantly, the candidate dose may lie outside the range actually studied. In that case, the model is extrapolating: its estimate depends on the assumed dose-response form and the posterior for its parameters, rather than on observations at nearby doses.
The distinction matters in dose escalation. A boundary dose can be the next dose considered, the dose currently recommended, or simply the endpoint of the protocol’s candidate grid. These are not interchangeable. The lowest candidate dose need not be the lowest observed dose, and the highest candidate dose need not be supported by observations at or near it. A boundary location alone does not determine whether a dose is safe, toxic, admissible, or appropriate to recommend. Those judgments depend on the posterior evidence and the protocol’s decision rules.
The BLRM at a boundary
For a binary dose-limiting toxicity (DLT) outcome, the two-parameter BLRM represents the probability of a DLT at dose \(d\) as
Here \(d_{\mathrm{ref}}\) is a fixed positive reference dose, \(p(d)\) is the DLT probability at dose \(d\), and \(\alpha\) and \(\beta\) are the intercept and dose slope. The model assumes that, conditional on the dose-specific probabilities, patient DLT outcomes are independent Bernoulli outcomes, and that the log odds change linearly with log dose over the range where the model is used. The model is monotone increasing in dose only when \(\beta>0\); a prior that permits negative values does not guarantee monotonicity in every posterior draw.
Dose-range boundary concerns are therefore about both data support and model assumptions. At a dose beyond the observed range, there is no local observation to check whether the log-linear relationship continues. The posterior still produces an estimate, but its uncertainty reflects uncertainty in the slope as well as the intercept. A strong prior on the slope can make this estimate look precise even when the data contain little information about the extrapolation.
For a given dose, write \(x=\log(d/d_{\mathrm{ref}})\) and \(Z(d)=\alpha+\beta x\). If the posterior for \((\alpha,\beta)\) is approximately bivariate normal, then \(Z(d)\) is normal. Its posterior mean and variance are
The \(x^2\) term shows why uncertainty in the slope can matter more as the dose moves away from the reference dose. The covariance term also matters: depending on the posterior dependence between intercept and slope, it can increase or decrease uncertainty at a particular dose. Thus, posterior variance need not increase steadily with dose. It is not correct to diagnose extrapolation by looking only at the variance or only at the distance from \(d_{\mathrm{ref}}\); compare the dose with the observed range and examine the joint parameter uncertainty.
Using the posterior to assess a boundary dose
The BLRM’s posterior distribution can be used to evaluate uncertainty in toxicity at a boundary. For a fixed overdose threshold \(p_{\mathrm{OD}}\), the event that dose \(d\) exceeds that threshold is equivalent to a threshold event for the linear predictor:
If \(Z(d)\mid\text{data}\) is normal with posterior mean \(m_d\) and standard deviation \(s_d\), then
where \(\Phi\) is the standard normal cumulative distribution function. In practice, this probability is often calculated from posterior draws rather than a normal approximation. Either way, the probability describes posterior uncertainty under the specified model and prior; it is not a direct frequency of future DLTs.
The following calculation also explains why one should inspect the boundary itself, not only a posterior mean toxicity. Applying the logit transformation to a posterior mean of \(\alpha\) and \(\beta\) gives a point estimate, but it does not show how much posterior mass lies above a clinically relevant toxicity threshold.
Worked Example: An Upper Candidate Dose Beyond the Observed Range
Hypothetical data. Six patients were treated at 10 mg, with 1 DLT, and six at 20 mg, with 2 DLTs. The candidate grid includes 10, 20, and 40 mg, so 40 mg is the highest candidate dose and lies above the highest observed dose. Suppose a BLRM fit using \(d_{\mathrm{ref}}=10\) mg gives posterior means \(E(\alpha\mid\text{data})=-2.20\), \(E(\beta\mid\text{data})=1.10\), posterior variances \(0.2025\) for each parameter, and posterior covariance \(-0.0405\). Treat these as illustrative fitted posterior summaries, not as a uniquely determined result from the counts alone.
At 40 mg, \(x=\log(40/10)=\log(4)=1.386\). The posterior mean linear predictor is
Its posterior variance is
so \(s_{40}\approx\sqrt{0.479}=0.692\). The toxicity probability evaluated at the posterior mean linear predictor is \(\operatorname{logit}^{-1}(-0.675)\approx0.337\), or about 33.7%. This is a point calculation, not the posterior mean of \(p(40)\).
For illustration, take an overdose threshold \(p_{\mathrm{OD}}=0.30\). Its logit is \(\log(0.30/0.70)=-0.847\). Under the normal approximation for \(Z(40)\),
The approximate posterior probability of exceeding 30% is 59.8%. This calculation does not, by itself, decide whether 40 mg is admissible or recommended; the protocol’s fixed decision rules must still be applied. It does show that the estimate at 40 mg is a slope-dependent extrapolation with substantial posterior uncertainty, despite the model returning a definite point estimate.
Worked Example: A Lowest Candidate Dose Below the Observed Range
Use the same hypothetical fit, but suppose the candidate grid also includes 5 mg and no patient has received less than 10 mg. The 5 mg dose is now the lowest candidate dose and is an extrapolation below the observed range. With \(d_{\mathrm{ref}}=10\) mg, \(x=\log(5/10)=\log(0.5)=-0.693\). The posterior mean linear predictor is
The corresponding posterior variance is
giving \(s_{5}\approx0.596\). The inverse-logit of the posterior mean linear predictor is \(\operatorname{logit}^{-1}(-2.962)\approx0.049\), or 4.9%. That low point value does not establish that the dose has negligible toxicity. For example, the posterior probability that the toxicity probability exceeds 5% is approximately
because \(\operatorname{logit}(0.05)=-2.944\). Under the stated normal approximation, nearly half the posterior mass is above 5%. This hypothetical calculation illustrates why a low fitted toxicity probability at the lowest dose should not be confused with high certainty about that probability. Whether the dose is useful or appropriate also depends on the trial’s objectives and protocol.
What the examples imply about boundary decisions
The two examples show that extrapolation can occur at either end of a candidate grid. When the model has a positive dose slope, a lower-than-observed dose will generally have a lower fitted toxicity probability, and a higher-than-observed dose a higher one. But the amount of uncertainty in those predictions depends on the posterior for both parameters and their covariance. The posterior may also be sensitive to prior assumptions, particularly when the data cover a narrow dose range or provide little information about the slope.
A short derivation makes the role of distance explicit. Let \(x=\log(d/d_{\mathrm{ref}})\). Under a bivariate normal posterior for \((\alpha,\beta)\), any linear combination of those parameters is normal. Therefore \(Z(d)=\alpha+x\beta\) has the mean and variance given above. If the posterior covariance is zero, the slope-uncertainty contribution to the variance is \(x^2\operatorname{Var}(\beta\mid\text{data})\), which grows with \(|x|\). With nonzero covariance, the additional term \(2x\operatorname{Cov}(\alpha,\beta\mid\text{data})\) can offset or amplify that growth at some doses. This is why neither “the endpoint must be most uncertain” nor “uncertainty always rises toward the highest dose” is a general theorem.
A boundary dose that was itself studied can still be uncertain. Having observations at a dose provides direct information about that dose’s toxicity probability, but the amount of information depends on the number of treated patients and DLT outcomes. Conversely, an unobserved candidate dose may have a narrow posterior if the model and prior impose strong structure, but that precision is not the same as direct empirical support.
Diagnostics and alternatives
Before acting on an endpoint estimate, verify the dose history and candidate grid, then review the posterior in the context of the observed range. Useful checks include:
- Distance from the observed doses: record whether the proposed dose is bracketed by observed doses and how far it lies beyond the nearest observation on a log-dose scale.
- Posterior slope and dependence: examine the posterior distribution of \(\beta\) and its association with \(\alpha\), rather than relying only on a fitted mean curve.
- Posterior toxicity at the endpoint: inspect the distribution or relevant tail probability for \(p(d)\), not just the inverse-logit of the mean linear predictor.
- Prior sensitivity: check whether reasonable changes to the slope prior substantially change the boundary estimate or posterior probability. Large changes signal that the data provide limited information for that extrapolation.
- Model fit across observed doses: compare fitted probabilities with observed cohort counts. A poor fit within the observed range weakens confidence in extending the same curve beyond it.
These checks do not replace the protocol’s admissibility and allocation rules. They help identify when a recommendation is driven by a model’s extension beyond the data rather than by local observations. If evidence is too weak, possible design responses include restricting the candidate range to doses supported by prior or emerging evidence, requiring additional data before moving farther outward, or using scientifically justified dose-response information. Such choices should be prespecified and clinically justified; changing a range after seeing outcomes can alter the operating characteristics of the design.
| Approach | What it contributes | Boundary limitation |
|---|---|---|
| BLRM extrapolation | Uses a continuous parametric dose-response relationship to estimate toxicity at an unobserved dose. | Depends on the log-linear log-odds form and posterior slope information outside the observed range. |
| Restricting the candidate range | Avoids making recommendations at doses beyond the prespecified range. | Does not create information about an endpoint that remains in the range but is poorly observed. |
| Additional observations near the boundary | Can provide more direct information about toxicity at or near the dose of interest. | Requires treating more participants and does not remove uncertainty or guarantee a good model fit. |
| Alternative dose-response structure | Can represent different scientific beliefs about how toxicity changes with dose. | More flexibility may require more data and can make sparse-boundary estimates less stable. |
A common mistake is to say, “The model predicts a low toxicity rate at the lowest dose, so that dose is certainly safe,” or, at the upper endpoint, “The model returns a value, so the dose is supported.” Both statements confuse a model-based estimate with the strength of evidence. The correction is to state whether the dose is inside the observed range, report posterior uncertainty in context, and apply the protocol’s decision rule without treating extrapolation as direct observation.
Boundary analysis has limits. It cannot verify that the assumed dose-response relationship remains valid beyond the data, and posterior calculations are conditional on the model, prior, endpoint definition, and available follow-up. A sensitivity analysis can reveal dependence on assumptions, but it cannot manufacture information at an unobserved dose. For that reason, boundary estimates should be interpreted as conditional predictions and not as empirical guarantees.
Worked Example: Comparing an Observed Dose With an Extrapolated Endpoint
Return to the hypothetical posterior summaries in the first example. Compare 20 mg, the highest observed dose, with 40 mg, the highest candidate dose. At 20 mg, \(x=\log(2)=0.693\), so
The inverse-logit of \(m_{20}\) is about 19.2%. For the same 30% threshold used earlier, the approximate posterior probability of exceeding the threshold is
The analogous probability at 40 mg was approximately 59.8%. These values illustrate how the posterior assessment changes with dose under the fitted model; the 20 mg estimate is informed by observations at 20 mg, whereas 40 mg is extrapolated beyond the observed range. Neither calculation alone establishes the final dose choice. The comparison is a reason to examine the source and robustness of the upper-dose estimate before using it in a protocol decision.
Check Your Understanding
Answer each question, then compare your reasoning with the solutions.
- A candidate dose lies between the lowest and highest observed doses. Is it necessarily well supported? Explain.
- In the BLRM, what part of the posterior variance of the linear predictor reflects uncertainty in the dose slope?
- Why is the inverse-logit of the posterior mean linear predictor not generally the posterior mean toxicity probability?
- What does a large change in an endpoint overdose probability under reasonable slope-prior changes suggest?
Solutions. (1) No. Being within the observed range makes the dose an interpolation point, but sparse observations, uneven dose spacing, or poor model fit can still leave substantial uncertainty. (2) The slope-related term is \(x^2\operatorname{Var}(\beta\mid\text{data})\), together with the covariance contribution \(2x\operatorname{Cov}(\alpha,\beta\mid\text{data})\). (3) The inverse-logit is nonlinear, so transforming the mean of \(Z(d)\) does not generally equal the mean of \(\operatorname{logit}^{-1}\{Z(d)\}\) across the posterior. (4) It suggests that the data provide limited information about extrapolation at that dose and the result depends materially on prior assumptions; it is a diagnostic for caution, not a new decision rule.