Tutorials › Model-Informed Drug Development (MIDD) › Stress-Testing Dose Recommendations Against Assumptions

Dose Justification Across Development Phases · Tutorial 253 of 596

Stress-Testing Dose Recommendations Against Assumptions

Use structured sensitivity analyses to determine whether a dose recommendation is robust or depends on a particular modeling or data assumption.

Advanced 12 min read

What You'll Learn

  • Distinguish a robustness check from choosing the model that produces the preferred dose.
  • Recalculate dose predictions under linear and saturating dose-response forms.
  • Assess how endpoint definitions can change the apparent dose ranking.
  • Compare missing-data assumptions and identify when they reverse a ranking.
  • Quantify the influence of individual studies on pooled dose-effect estimates.
  • Separate stable conclusions from decisions that depend on an assumption.

Why Stress-Test a Dose Recommendation?

A dose recommendation is supported by evidence interpreted through assumptions: a chosen dose-response model, an endpoint definition, a rule for handling missing outcomes, and decisions about which studies to include. A recommendation can appear clear under the primary analysis yet change when one of these assumptions changes. Stress-testing makes that dependence visible before a recommendation is treated as robust.

The purpose is not to try many analyses and select the one that gives the preferred answer. Instead, define a small set of scientifically plausible alternatives, recalculate the relevant dose comparisons under each, and report whether the decision changes. The earlier material on dose-response interpretation and dose selection supplies the clinical context: statistical differences alone do not establish that one dose is meaningfully better or acceptably safe.

Definition: A dose stress test is a prespecified sensitivity analysis that changes one material modeling or data assumption at a time, or examines a coherent alternative set of assumptions, and then recalculates the dose ranking or decision. A result is decision-sensitive when a plausible alternative changes the selected dose, the ordering of doses, or whether a stated decision criterion is met.

A useful starting point is to state the decision rule in clinical terms. For example, a dose might need to exceed a minimum efficacy threshold while remaining below a safety limit. Keep that rule fixed while testing assumptions; otherwise, a changed result could be due to the rule rather than the evidence model. If the rule itself is uncertain, that is a separate sensitivity question and should be identified as such.

Changing the Dose-Response Model

A model form determines how observed dose-response information is translated into predictions at candidate doses. A linear model assumes a constant change in response per unit increase in dose over the modeled range. A saturating model allows the response to approach a maximum, so that the incremental gain tends to shrink as dose rises. Neither form is automatically correct; the data, scientific plausibility, and range of doses matter.

Definition: In a simple linear model, the predicted efficacy response at dose \(d\) is \(E(d)=\beta_0+\beta_1d\). In a simple Emax model, it is \(E(d)=E_0+E_{\max}d/(ED_{50}+d)\), where \(E_0\) is the baseline response, \(E_{\max}\) is the maximum modeled increase, and \(ED_{50}\) is the dose producing half that increase. These formulas describe mean responses under the specified model; their predictions are only as credible as the model assumptions and supporting dose range.

A model comparison should examine more than which fitted curve looks smoother. Check whether each form adequately describes the observed pattern, whether its parameters are estimated with useful precision, and whether the decision depends on predictions outside the observed dose range. A flexible model can follow noise; a simple model can miss curvature. When the ranking changes, the appropriate conclusion is that the recommendation is model-sensitive—not that whichever model predicts the desired dose should be preferred.

Worked Example: Linear Versus Saturating Predictions

Suppose a hypothetical development program considers doses of 10, 20, and 40 mg. Efficacy is measured as improvement in a symptom score, in points, with larger values better. For illustration, the fitted linear model is \(E(d)=2+0.22d\), while a plausible saturating model is \(E(d)=2+12d/(20+d)\). The decision rule requires predicted improvement of at least 7 points and a serious adverse-event risk no greater than 10%. Hypothetical risks are 2%, 5%, and 12% at 10, 20, and 40 mg, respectively.

Under the linear model, predictions are \(E(10)=2+0.22(10)=4.20\) points, \(E(20)=2+0.22(20)=6.40\) points, and \(E(40)=2+0.22(40)=10.8\) points. Only 40 mg reaches the 7-point efficacy threshold, but its 12% event risk exceeds the 10% safety limit. No dose meets both criteria under this model.

Under the saturating model, \(E(10)=2+12(10)/(20+10)=6.00\) points, \(E(20)=2+12(20)/(20+20)=8.00\) points, and \(E(40)=2+12(40)/(20+40)=10.0\) points. The 20-mg dose reaches the efficacy threshold and has a 5% event risk. The 40-mg dose remains outside the safety limit. Thus, 20 mg meets both criteria under the saturating model, while no dose does under the linear model.

This hypothetical calculation does not establish which model is correct. It shows that the decision depends on dose-response shape. Before relying on the 20-mg conclusion, examine fit, uncertainty, parameter plausibility, and whether the available data support the predicted curvature.

Endpoint Definitions Can Change the Ranking

A dose ranking is always a ranking for a particular endpoint and analysis population. A continuous change score and a responder endpoint can order doses differently because they answer different questions. A responder definition also depends on a threshold: changing the threshold can change which patients count as responders, even when the underlying measurements are unchanged. Therefore, endpoint sensitivity should compare clinically meaningful, defensible definitions, not thresholds chosen after seeing which one favors a dose.

Worked Example: Continuous Outcome Versus Responder Rate

Consider hypothetical results at two doses. At 10 mg, the mean symptom improvement is 6.0 points and 58 of 100 patients meet a prespecified responder definition. At 20 mg, the mean improvement is 7.0 points and 56 of 100 meet that definition. Suppose the responder definition is an improvement of at least 5 points.

On the continuous endpoint, the 20-mg dose ranks higher: \(7.0-6.0=1.0\) point greater mean improvement. On the responder endpoint, the rates are \(58/100=0.580\), or 58%, at 10 mg and \(56/100=0.560\), or 56%, at 20 mg. The 10-mg dose ranks higher by 2 percentage points on this endpoint.

The results are not arithmetically inconsistent: a higher mean can coexist with a lower proportion above a threshold. They describe different features of the response distribution. Report the endpoint definition and clinical interpretation alongside each ranking; do not describe one endpoint as confirming the other unless the evidence supports that conclusion.

Missing Outcomes and Study Influence

Missing-data assumptions matter when the probability of an unobserved outcome is related to treatment, observed information, or the unobserved outcome itself. A complete-case analysis describes patients with observed outcomes and may not represent all randomized patients. A nonresponder-imputation analysis treats missing outcomes as failures for a responder endpoint. A missing-at-random analysis instead relies on the assumption that, conditional on included observed information, missingness does not depend on the unobserved outcome. These assumptions are not interchangeable, and none can be verified from the observed data alone.

Worked Example: Missing-Data Assumption Reverses a Ranking

In a hypothetical responder analysis, 100 patients are assigned to each dose. At 10 mg, 54 of 90 patients with observed outcomes respond; 10 outcomes are missing. At 20 mg, 42 of 60 patients with observed outcomes respond; 40 outcomes are missing.

With nonresponder imputation, missing patients are counted as nonresponders. The estimated rates among all randomized patients are \(54/100=54\%\) at 10 mg and \(42/100=42\%\) at 20 mg, so 10 mg ranks higher. Under a simplified missing-at-random calculation that imputes each arm’s observed response rate to its missing patients, the rates are \(54/90=60\%\) and \(42/60=70\%\), respectively. The implied all-randomized response counts are \(54+10(0.60)=60\) and \(42+40(0.70)=70\); the corresponding rates are 60% and 70%. Under this assumption, 20 mg ranks higher.

The ranking changes because the assumptions assign different outcomes to missing patients, especially in the 20-mg group. The calculation is a sensitivity illustration, not proof that either assumption is appropriate. Investigators should examine why outcomes are missing, use observed patient information where justified, and show how plausible alternatives affect the decision.

A study can also exert disproportionate influence on a pooled estimate. For independent study estimates with known standard errors, the inverse-variance weighted mean is the estimate-weighted average using weights equal to the inverse variances. The following derivation explains why a study with a small standard error receives greater weight.

$$ \hat{\theta}=\frac{\sum_{i=1}^{k} w_i\hat{\theta}_i}{\sum_{i=1}^{k}w_i}, \qquad w_i=\frac{1}{SE_i^2} $$

To justify the formula, assume \(\hat{\theta}_i\) are independent normally distributed estimates of a common effect \(\theta\), with variances \(SE_i^2\). Their log likelihood, apart from constants, is \(-\frac{1}{2}\sum_i(\hat{\theta}_i-\theta)^2/SE_i^2\). Differentiating with respect to \(\theta\) and setting the derivative to zero gives \(\sum_i(\hat{\theta}_i-\theta)/SE_i^2=0\). Rearranging yields the weighted-mean formula above. The common-effect and independence assumptions are important: if studies estimate different underlying effects, this calculation may not be an adequate pooling model.

Worked Example: Excluding an Influential Study

Suppose three hypothetical studies estimate the efficacy contrast for 20 mg versus a reference dose as 8 points (standard error 2), 2 points (standard error 2), and −1 point (standard error 1). For comparison, suppose a separate 40-mg analysis gives an estimate of 3 points. Use inverse-variance weights for the 20-mg estimates.

The weights are \(1/2^2=0.25\), \(1/2^2=0.25\), and \(1/1^2=1.00\). The pooled 20-mg estimate is \((0.25(8)+0.25(2)+1.00(-1))/(0.25+0.25+1.00)=(2+0.5-1)/1.5=1.00\) point. Compared with the 3-point estimate for 40 mg, the pooled point estimates rank 40 mg higher.

Now exclude the third study as a sensitivity analysis. The remaining pooled estimate is \((0.25(8)+0.25(2))/(0.25+0.25)=2.5/0.5=5.00\) points. It exceeds the 3-point estimate for 40 mg, so the point-estimate ranking reverses. This exclusion does not establish that the third study is invalid: the reason for exclusion must be scientifically or methodologically defensible. It shows that the ranking depends materially on that study under this calculation.

Interpreting Results and Diagnosing Fragility

For every stress test, record the baseline assumption, the alternative, the recalculated result, and whether the decision changes. Also check whether the comparison is being made on the same endpoint, dose range, population, and decision criteria. A changed estimate is not necessarily a changed decision; conversely, a small change can matter if it moves a dose across a clinical or safety threshold.

1
Specify the decision.
State the candidate doses, endpoint, clinical thresholds, and safety constraints before examining alternative analyses.
2
Choose plausible alternatives.
Test model forms, endpoint definitions, missing-data assumptions, and study exclusions that have a scientific or methodological rationale.
3
Recalculate consistently.
Keep the target population and decision rule fixed where possible, and report estimates, units, and any changed assumptions.
4
Classify the result.
State whether the dose order, threshold status, or selected dose changes, then explain which assumption drives that change.

Useful diagnostics include residual patterns and poor fit for a dose-response model; sparse data near a predicted threshold; excessive sensitivity to a responder cutoff; substantial or imbalanced missingness; and large changes when a study is removed. Also inspect whether predictions rely on extrapolation, whether one study has unusually high weight, and whether its results differ for a known clinical or design reason. These signals do not automatically invalidate an analysis, but they identify where a recommendation needs qualification or further evidence.

Alternative approaches answer related but different questions. A primary analysis estimates effects under its chosen assumptions; a sensitivity analysis checks how conclusions vary under alternatives. Model averaging can represent uncertainty across several dose-response forms, but requires a defensible set of models and a way to combine them. A random-effects meta-analysis allows study effects to vary, unlike the common-effect calculation illustrated above, but it does not make influential studies or missing-data assumptions irrelevant. These methods complement, rather than replace, transparent stress tests.

Stress tests have limits. A stable ranking across a few alternatives does not prove that all important assumptions are correct. Conversely, instability may reflect sparse evidence or endpoints that capture distinct clinical features rather than a defective analysis. Avoid selecting assumptions after observing which dose they favor, treating a nonsignificant contrast as proof of equality, or interpreting a sensitivity analysis as independent confirmation. The practical outcome is a calibrated statement: which dose is favored, under which assumptions, and how easily that conclusion changes.

Common pitfall: “The dose remains preferred in the primary analysis, so the recommendation is robust.” This reasoning ignores the possibility that the primary model, endpoint, or missing-data assumption drives the result. Robustness requires recalculating the decision under plausible alternatives and showing that the relevant ranking or criteria remain stable.

Check Your Understanding

Answer each question, then compare with the brief solutions below.

  1. A linear model predicts 5 points of improvement at 10 mg and 9 points at 30 mg. What is its predicted improvement at 20 mg, and what assumption makes this interpolation valid?
  2. At one dose, 48 of 80 patients with observed outcomes respond, and 20 of 100 randomized patients have missing outcomes. What are the observed response rate and the rate under nonresponder imputation?
  3. Why can a dose rank higher on mean change but lower on responder rate?
  4. In the pooled-study example, why did excluding the third study increase the 20-mg estimate?

Solutions: (1) The prediction is halfway between the two fitted values: \(5+(20-10)(9-5)/(30-10)=7\) points. This requires the linear model’s constant slope to be reasonable over that interval. (2) The observed rate is \(48/80=60\%\); nonresponder imputation gives \(48/100=48\%\). (3) The mean uses the magnitude of all measured changes, whereas the responder rate counts how many cross a threshold; the distributions can differ in spread or location. (4) The third study had weight 1.00 and an estimate of −1 point, pulling the weighted average down more than either study with weight 0.25.