Tutorials › Optimal Design and Model-Informed Dosing › Weak Identification Near Dose-Response Plateaus

Dose Optimization Under Uncertainty · Tutorial 218 of 497

Weak Identification Near Dose-Response Plateaus

Learn why a target close to the Emax plateau can produce a poorly identified dose, and how to diagnose that problem with profiles and target-sensitivity calculations.

Advanced 12 min read

What You'll Learn

  • Derive the target-attaining dose in the Emax model and identify when it becomes highly sensitive to the target.
  • Distinguish a broad profile likelihood interval from a precisely estimated dose.
  • Interpret profile-likelihood cutoffs and report interval endpoints only when supported by profile values.
  • Check how small changes in the target alter the estimated dose and its uncertainty.
  • Compare profile diagnostics with delta-method, bootstrap, and posterior summaries.
  • Recognize assumptions and warning signs when the dose profile is flat or does not close.

Why a Plateau Can Leave the Target Dose Uncertain

In dose optimization, a fitted dose-response curve may appear to approach a clear maximum while leaving the dose needed to attain a high target poorly determined. The reason is geometric: close to the plateau, a large change in dose can produce only a small change in predicted response. Data can therefore support a broad range of doses as plausible solutions to the same target, even when the curve's maximum response seems reasonably well estimated.

The issue is weak identification: the observed data do not sharply distinguish among values of the target-attaining dose. It is not the same as saying that no estimate can be calculated. A point estimate may exist, but small changes in the target or fitted curve can move it substantially, and its uncertainty interval may be wide. This tutorial uses the Emax model and focuses on diagnosing that behavior, rather than choosing a dose under decision constraints.

Definition: For a specified response target, the target-attaining dose is the dose at which the model's predicted mean response equals that target. It is weakly identified when the data provide little information to distinguish that dose from substantially different doses, as reflected, for example, by a broad profile-likelihood interval or strong sensitivity to small target changes.

The Emax Dose and Its Sensitivity Near the Plateau

Consider the increasing Emax model

$$ \mu(d)=E_0+\frac{E_{\max}d}{ED_{50}+d}, $$

where \(d\) is dose, \(E_0\) is the mean response at zero dose, \(E_{\max}>0\) is the additional response at the plateau, and \(ED_{50}>0\) is the dose producing half of that additional response. The asymptotic mean response is \(E_0+E_{\max}\). For this model, assume the mean structure is appropriate over the studied doses and the response observations are modeled with a likelihood whose parameters can be estimated from the data. The equation concerns the mean response; it does not assert that every individual response equals that mean.

Let \(T\) be an absolute target response. Solving \(\mu(d_T)=T\) gives

$$ T-E_0=\frac{E_{\max}d_T}{ED_{50}+d_T} \quad\Longrightarrow\quad d_T=ED_{50}\frac{T-E_0}{E_0+E_{\max}-T}. $$

A finite positive solution requires \(E_0<T<E_0+E_{\max}\). As \(T\) approaches the plateau from below, the denominator approaches zero, so the target dose can grow rapidly. Differentiating with respect to the target gives

$$ \frac{\partial d_T}{\partial T} = \frac{ED_{50}E_{\max}}{(E_0+E_{\max}-T)^2}. $$

Thus, for fixed model parameters, target sensitivity increases sharply as the target approaches the asymptote. In real estimation, the parameters are uncertain too; the formula is a useful diagnostic of the geometry, not a substitute for propagating that uncertainty.

Worked Example: A Small Target Change Near the Plateau

Suppose a hypothetical fitted curve has \(E_0=10\) response units, \(E_{\max}=80\) response units, and \(ED_{50}=40\) mg. Its plateau is \(10+80=90\) response units. Calculate the target doses for \(T=80\) and \(T=82\) response units.

For \(T=80\), the target is below the plateau, so the finite-dose solution exists:

$$ d_{80}=40\frac{80-10}{10+80-80} =40\frac{70}{10} =280\ \text{mg}. $$

For \(T=82\):

$$ d_{82}=40\frac{82-10}{10+80-82} =40\frac{72}{8} =360\ \text{mg}. $$

A target increase of only 2 response units changes the point estimate by \(360-280=80\) mg. The local derivative at \(T=80\) is \(40(80)/(10^2)=32\) mg per response unit; the finite change is larger because sensitivity increases as the target moves closer to the 90-unit plateau. This calculation signals potential instability, but by itself it does not quantify sampling uncertainty.

Using a Profile to Diagnose Identification

A profile likelihood examines how well different candidate target doses fit the data while allowing the other model parameters to adjust. For each candidate \(d_T\), impose the constraint \(\mu(d_T)=T\), maximize the likelihood over the remaining free parameters, and compare that constrained maximum with the unconstrained maximum. A profile that rises slowly as the candidate dose moves away from its optimum indicates that the data do not strongly distinguish those doses.

For a one-dimensional target dose, a commonly used approximate 95% likelihood-ratio confidence set contains values satisfying a profile deviance no greater than 3.84. This cutoff uses the approximate chi-square reference with one degree of freedom and relies on regularity conditions, including adequate identifiability and an interior parameter value. Near a plateau those conditions may be questionable, so the cutoff is a diagnostic convention, not a guarantee of exact 95% coverage. The profile should be inspected, and its results interpreted alongside the design and model fit.

Worked Example: Estimating Profile-Interval Endpoints From Supported Values

In a hypothetical study, a profile is calculated for the dose attaining a target of 80 response units. The table gives the profile deviance at selected doses. Use the 3.84 cutoff and linear interpolation between the two listed values surrounding each crossing.

Dose (mg)Profile deviance
355.0
453.0
601.4
1000.5
1801.0
2602.8
3003.4
3304.2

The lower crossing lies between 35 mg, where the deviance is 5.0, and 45 mg, where it is 3.0. Linear interpolation gives

$$ 35+\frac{5.0-3.84}{5.0-3.0}(45-35) =35+\frac{1.16}{2.0}(10) =40.8\ \text{mg}. $$

The upper crossing lies between 300 mg, at 3.4, and 330 mg, at 4.2:

$$ 300+\frac{3.84-3.4}{4.2-3.4}(330-300) =300+\frac{0.44}{0.8}(30) =316.5\ \text{mg}. $$

Under the stated interpolation rule, the approximate profile interval is therefore 40.8 to 316.5 mg, or about 41 to 317 mg. It is broad: doses spanning roughly 276 mg remain inside the cutoff. This supports a diagnosis of weak identification. The endpoints are approximate because they depend on local interpolation between evaluated doses; a finer profile grid or direct numerical root-finding would refine them. Without profile values near the crossings, one could report only which evaluated doses fall inside or outside the cutoff, not invented endpoint estimates.

Why the Profile Can Be Flat

Near a plateau, dose-response predictions change slowly with dose. Consequently, shifting a candidate target dose can often be offset by changes in curve parameters, such as \(E_{\max}\) or \(ED_{50}\), with little loss of likelihood. This produces a shallow profile or a ridge of parameter combinations that fit nearly as well. The profile is useful precisely because it allows nuisance parameters to readjust instead of pretending they are known.

A short justification follows from the target-dose equation. Write \(A=E_0+E_{\max}\). Then \(d_T=ED_{50}(T-E_0)/(A-T)\). If \(A-T\) is small, a small change in \(A\) has a large effect on the denominator. In fact,

$$ \frac{\partial d_T}{\partial A} =-ED_{50}\frac{T-E_0}{(A-T)^2}. $$

The magnitude also grows as \(A-T\) shrinks. Thus uncertainty in the plateau estimate can create substantial uncertainty in the target dose, even if other fitted quantities appear stable. This is one reason to inspect the dose profile rather than infer precision from the point estimate alone.

Worked Example: Two Plausible Curves, Different High-Target Doses

Suppose two hypothetical parameter combinations fit the observed responses similarly because the data contain limited information about the plateau. The first is \(E_0=10\), \(E_{\max}=80\), \(ED_{50}=40\) mg; the second is \(E_0=10\), \(E_{\max}=84\), \(ED_{50}=50\) mg. Compare their target doses for 86 response units.

For the first curve, the plateau is 90 response units, and

$$ d_{86}^{(1)} =40\frac{86-10}{10+80-86} =40\frac{76}{4} =760\ \text{mg}. $$

For the second, the plateau is 94 response units, and

$$ d_{86}^{(2)} =50\frac{86-10}{10+84-86} =50\frac{76}{8} =475\ \text{mg}. $$

Both curves have a finite solution because 86 is below both plateaus, yet the target doses differ by \(760-475=285\) mg. This calculation does not establish that either parameter combination is supported: that judgment requires the data and likelihood. It demonstrates why uncertainty in the plateau and curve scale matters especially for a high target.

Diagnostics, Alternatives, and Limitations

Assess the profile over a sufficiently wide dose grid and inspect its shape, not only its minimum. Warning signs include a broad region under the likelihood-ratio cutoff, a profile that remains nearly flat at the edge of the scanned range, marked asymmetry, or an interval endpoint that changes substantially when the grid is refined. A profile that does not close on one side does not provide a finite endpoint from that scan; it is not evidence that the endpoint equals the last dose evaluated. Check that the model fits the observed dose-response pattern and that the target is below plausible plateau values. A poor structural fit can make a profile misleading as a measure of uncertainty within the chosen model.

The delta method, bootstrap distributions, and posterior distributions are complementary approaches covered in the earlier material on target-dose uncertainty. The delta method summarizes local uncertainty and can be unreliable when the target-dose transformation is strongly nonlinear. A bootstrap can show skewness or long tails, but its results depend on the resampling design and may be unstable when target estimates are frequently extreme or undefined. A posterior distribution incorporates prior information and can be informative when likelihood information is weak, but its concentration then depends in part on the prior. A profile makes the likelihood support across candidate doses visible; it does not include prior information and does not, by itself, diagnose model misspecification.

No single diagnostic repairs weak identification. More observations near the plateau may help estimate the asymptote, while observations across informative parts of the curve may help separate its parameters. However, adding doses or selecting a study design is a separate design question. This tutorial's central limitation is that its diagnostics are conditional on the chosen Emax structure and likelihood. A narrow profile cannot rule out a wrong model, and a wide profile should not be converted into a falsely precise point recommendation.

Common Mistake

Common pitfall: Treating a calculable target-dose estimate as evidence that the dose is precisely identified. For example, substituting fitted parameters into the Emax formula gives one number, but that arithmetic does not show how much the data distinguish that dose from alternatives. Inspect the profile interval and target sensitivity; report only endpoint estimates supported by profile values or an explicitly stated interpolation procedure.

Also avoid interpreting a dose beyond the largest studied dose as though the data directly observed the response there. Whether and how to address such extrapolation is a distinct issue; a profile diagnoses likelihood support under the fitted model, not the clinical credibility of every dose in its range.

Check Your Understanding

Try each question, then compare your reasoning with the solutions.

  1. For \(E_0=5\), \(E_{\max}=60\), \(ED_{50}=30\) mg, and target \(T=60\), calculate the target dose and state whether it is below the plateau.
  2. For a profile cutoff of 3.84, a dose with deviance 3.2 is inside or outside the approximate profile confidence set? What does that comparison mean?
  3. A profile has deviance 4.6 at 20 mg and 3.0 at 30 mg. Using linear interpolation, estimate the lower crossing of the 3.84 cutoff.
  4. Why can a small increase in a high target produce a large change in the target dose near the Emax plateau?
  5. Name one limitation of a likelihood profile as a diagnosis of target-dose uncertainty.

Solutions. (1) The plateau is \(5+60=65\) response units, so 60 is below it. The dose is \(30(60-5)/(5+60-60)=30(55)/5=330\) mg. (2) Since \(3.2<3.84\), the dose is inside the approximate set; its constrained fit has not deteriorated enough to cross the stated likelihood-ratio cutoff. (3) The crossing is \(20+(4.6-3.84)/(4.6-3.0)(30-20)=20+(0.76/1.6)(10)=24.75\) mg, approximately 24.8 mg under linear interpolation. (4) The derivative \(ED_{50}E_{\max}/(E_0+E_{\max}-T)^2\) grows as the target approaches the plateau, so the dose changes increasingly fast with the target. (5) The profile is conditional on the specified model and likelihood; it cannot establish that the model is correctly specified.

Clinical Biostats