Tutorials › Adaptive Clinical Trial Design › Identifiability Problems in Two-Parameter CRM

Model-Based Dose Finding: CRM and Extensions · Tutorial 577 of 799

Identifiability Problems in Two-Parameter CRM

Learn to diagnose when sparse dose-finding data cannot separate the two logistic parameters, and to assess how that uncertainty affects posterior toxicity estimates and dose choices.

Advanced 11 min read

What You'll Learn

  • Distinguish structural nonidentifiability from weak identification in a two-parameter logistic CRM.
  • Derive why observations at only one dose identify a linear predictor, not separate intercept and slope values.
  • Use a likelihood calculation to see how the observed dose locations inform the posterior.
  • Calculate posterior mean toxicity probabilities and apply the target-dose rule.
  • Diagnose parameter and recommendation instability without treating a prior-driven result as data-driven.
  • Compare a two-parameter model with a fixed-intercept CRM and recognize each approach’s limitations.

Why Two Parameters Can Be Hard to Learn

A two-parameter logistic CRM can flexibly represent dose-toxicity curves, but that flexibility creates a practical challenge: early in a trial, there may be too little information to determine both the curve’s level and its steepness. Several combinations of parameters can then fit the observed data similarly while implying different toxicity probabilities at untried doses. This is an identifiability problem, not merely a computational inconvenience. It can make a posterior dose recommendation depend substantially on the prior, even when the posterior calculation is performed correctly.

The issue is especially clear when all participants have been treated at one dose. Their outcomes can inform the toxicity probability at that dose, but the data alone do not generally determine how the curve changes on either side of it. A dose recommendation may nevertheless be produced by the prespecified CRM rule; interpreting that recommendation requires attending to the uncertainty behind it.

Model and Meaning of Identifiability

Let \(x_j\) be the fixed numerical score assigned to dose level \(j\), with larger scores representing higher doses. For participant \(i\), let \(Y_i=1\) indicate a dose-limiting toxicity (DLT) and \(Y_i=0\) otherwise. The two-parameter logistic CRM specifies

$$ P(Y_i=1\mid x_i,a,b)=p(x_i;a,b) =\operatorname{expit}(a+bx_i) =\frac{\exp(a+bx_i)}{1+\exp(a+bx_i)}. $$

Here \(a\) is the intercept, \(b\) is the slope, and the usual monotone dose-toxicity specification constrains \(b>0\). Conditional on the assigned doses and model parameters, the DLT outcomes are assumed independent Bernoulli observations with the probabilities specified above. The dose scores are fixed and known, DLT status is completely observed for each participant, and no unmodelled time trend or other covariate is included in this model. Bayesian inference additionally requires a specified prior on \((a,b)\), supported on the allowed parameter space.

Definition: A model parameter is structurally identifiable from a design if distinct allowed parameter values imply distinct distributions for the observations under that design. Parameters are weakly identified in a particular data set when the likelihood provides little information to distinguish among values, even though the design could in principle distinguish them with more or better-placed observations.

For this model, each observed dose informs the linear predictor \(a+bx\) at that dose. Identifying both \(a\) and \(b\) requires information about the curve at different dose scores. A prior can regularize the posterior when the likelihood is weak, but it cannot make the observed data contain information they do not contain. In that situation, it is important to separate what the data support from what follows from the prior.

Why One Observed Dose Creates a Ridge

Suppose \(n\) participants are all treated at the same score \(x_*\), and \(y\) experience a DLT. Apart from the binomial coefficient, their likelihood is

$$ L(a,b)\propto \{\operatorname{expit}(a+bx_*)\}^{y} \{1-\operatorname{expit}(a+bx_*)\}^{n-y}. $$

The likelihood depends on \(a\) and \(b\) only through \(a+bx_*\). Thus, any two parameter pairs satisfying \(a+bx_*=\eta_*\) have identical likelihood for these data. They form a line in the parameter plane, restricted only by the allowed parameter space. This line is often called a likelihood ridge.

The same point appears in the Fisher information. If the probability at \(x_*\) is \(p_*\), the information from these \(n\) observations is proportional to

$$ n p_*(1-p_*) \begin{pmatrix} 1 & x_*\\ x_* & x_*^2 \end{pmatrix}. $$

The determinant of this matrix is zero: its second row is \(x_*\) times its first row. There is therefore no information in this one-dose design to estimate both parameters separately. With observations at two distinct dose scores, the information matrix can have full rank. But full rank is not a guarantee of precision: few participants, few DLTs, or dose scores that are close together can still leave a broad, strongly correlated posterior for \(a\) and \(b\).

Worked Example: Same Likelihood, Different Curves

Consider two hypothetical parameter pairs, \((a,b)=(-1.386,0.5)\) and \((a,b)=(-1.386,1.5)\), on dose scores \(x=0,1,2\). Suppose the target toxicity probability is \(\phi=0.25\). At the starting score \(x=0\), both pairs give \(p(0)\approx0.2000\). If five participants receive this dose and one has a DLT, the likelihood under either pair is

$$ {5\choose 1}(0.2)^1(0.8)^4 =5(0.2)(0.4096)=0.4096. $$

The data at \(x=0\) therefore do not distinguish the two pairs. Their predictions away from that dose do differ. For the first pair, \(p(1)=\operatorname{expit}(-0.886)\approx0.2919\) and \(p(2)=\operatorname{expit}(-0.386)\approx0.4047\). For the second, \(p(1)=\operatorname{expit}(0.114)\approx0.5285\) and \(p(2)=\operatorname{expit}(1.614)\approx0.8340\).

The same observed dose and outcome can thus be compatible with quite different risks at higher scores. This is not a contradiction: those higher-dose probabilities were not directly informed by outcomes at \(x=0\). Any posterior preference between the two pairs in this example comes from their prior probabilities, not from a difference in their likelihoods.

How the Posterior Carries Uncertainty Into a Dose Rule

Let \(D\) denote all observed dose assignments and DLT outcomes. The posterior distribution is proportional to the likelihood times the prior. For each candidate dose \(j\), define its posterior mean toxicity probability as \(q_j=E[p(x_j;a,b)\mid D]\). Under the target-dose rule considered here, the recommended dose is the eligible candidate whose \(q_j\) is closest to \(\phi\). The rule selects using a posterior summary; it does not imply that the toxicity probability at that dose is known precisely.

Worked Example: A Recommendation With a Broad Posterior Prediction

Continue the hypothetical one-dose data above, and suppose the prior assigns probability \(0.5\) to each of the two parameter pairs. Because their likelihoods at \(x=0\) are equal, their posterior probabilities remain \(0.5\) and \(0.5\). The posterior mean toxicity probabilities at \(x=0,1,2\) are therefore

$$ q_0=0.2000,\qquad q_1=\frac{0.2919+0.5285}{2}=0.4102,\qquad q_2=\frac{0.4047+0.8340}{2}=0.6194. $$

Their absolute distances from \(\phi=0.25\) are approximately \(0.0500\), \(0.1602\), and \(0.3694\), respectively. The target-dose rule selects \(x=0\), the dose with posterior mean probability closest to the target. Yet the posterior probability at \(x=1\) is split between approximately \(0.2919\) and \(0.5285\), and at \(x=2\) between approximately \(0.4047\) and \(0.8340\). The selected dose is a well-defined result of the rule, but the higher-dose predictions remain uncertain.

Averaging probabilities over the posterior is not the same as averaging the parameters and then applying the logistic function. The dose rule here uses \(E[p(x_j;a,b)\mid D]\), the posterior mean of the toxicity probability. In general, because the logistic function is nonlinear, \(E[\operatorname{expit}(a+bx_j)\mid D]\) need not equal \(\operatorname{expit}(E[a\mid D]+x_jE[b\mid D])\).

What Additional Dose Information Can—and Cannot—Do

Worked Example: A Second Dose Level Adds Information

Suppose a hypothetical trial observes 1 DLT among 5 participants at score \(x=0\), and 3 DLTs among 5 participants at score \(x=2\). The observed proportions are \(0.2\) and \(0.6\). To illustrate the information supplied by distinct dose levels, match the logistic curve to these two group proportions. Their logits are

$$ \operatorname{logit}(0.2)=\log\!\left(\frac{0.2}{0.8}\right) \approx -1.3863,\qquad \operatorname{logit}(0.6)=\log\!\left(\frac{0.6}{0.4}\right) \approx 0.4055. $$

Equating these logits to \(a+bx\) at the two scores gives \(a\approx-1.3863\) at \(x=0\), and \(-1.3863+2b\approx0.4055\) at \(x=2\). Thus

$$ b\approx\frac{0.4055-(-1.3863)}{2}=0.8959. $$

Unlike the one-dose case, the two distinct scores provide two different linear predictors, so both parameters can be determined by this illustrative matching calculation. This is not evidence that the parameter estimates are precise: only five participants were observed at each score. Sampling variation in either group proportion could substantially change the fitted curve. The example shows the difference between being able to distinguish parameters in principle and having enough data to estimate them reliably.

With observations at scores \(x_1\) and \(x_2\), the information matrix is a sum of contributions in the directions \((1,x_1)\) and \((1,x_2)\). If the contributions have positive weights \(w_1\) and \(w_2\), its determinant is \(w_1w_2(x_1-x_2)^2\). Distinct scores make that determinant positive, but a small separation or small weights can make it small. In a trial, the weights depend on the number of observations and the probabilities: information is limited when there are few participants or when observed outcomes are nearly all the same at an informative dose.

Diagnosing Instability and Choosing a Comparison

A useful review of a two-parameter CRM asks not only whether the computation converged, but also whether the available dose-outcome pattern supports the curve used for the recommendation. Examine the posterior distributions jointly and marginally. A pronounced intercept-slope trade-off, broad marginal distributions, or a posterior concentrated near a parameter constraint can indicate weak information. Also compare posterior predictions at each candidate dose and assess whether the recommended dose changes under reasonable prior choices or modest changes in the observed data. Sensitivity is especially informative when the early observations are concentrated at one dose.

These checks are not automatic proof of failure. For example, a correlation between \(a\) and \(b\) is expected when a likelihood ridge runs through the parameter space. It is the resulting uncertainty in clinically relevant predictions and decisions that matters. Report uncertainty in toxicity probabilities, and avoid describing a posterior mean as if it were the risk for every plausible curve. If one dose is clearly closest to target under all plausible posterior curves, a recommendation may be stable despite parameter uncertainty; if plausible curves point to different doses, the rule’s single selected dose can conceal decision uncertainty.

ApproachWhat it doesTrade-off under sparse data
Two-parameter logistic CRMEstimates both the level and steepness of the logistic curve.Flexible, but observations at a narrow set of doses can leave substantial parameter and prediction uncertainty.
Fixed-intercept logistic CRMFixes the intercept and estimates the slope.Reduces the number of unknown parameters, but conclusions depend on the fixed intercept being appropriate.

The fixed-intercept model is a meaningful alternative when its intercept is justified by prior clinical information or a prespecified calibration. It is not a general cure for sparse data: fixing a parameter transfers uncertainty and assumptions into that constraint. Conversely, retaining two parameters does not by itself make the model unsuitable; it makes careful prior specification and uncertainty assessment particularly important.

Common Pitfall

Wrong reasoning: “The posterior dose rule selected a dose, so the data must have identified both logistic parameters.” A rule can return a dose even when the likelihood has a ridge, because the prior and the decision rule still produce a posterior summary and an ordering of candidates.

Correction: Check whether the observed dose levels provide information that distinguishes the parameters, then examine the posterior predictions and recommendation sensitivity. A computed recommendation is a decision under the model and prior; it is not evidence that the curve is precisely identified. When the recommendation is prior-sensitive, state that limitation rather than presenting the selected dose as a data-determined conclusion.

Key takeaway: One observed dose level informs one linear predictor, \(a+bx\), rather than separately identifying \(a\) and \(b\). Distinct dose levels can make both parameters identifiable in principle, but sparse outcomes may still leave posterior toxicity predictions and dose recommendations uncertain.

Check Your Understanding

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

  1. If all participants receive score \(x=1\), which combination of \(a\) and \(b\) determines their likelihood?
  2. Why does observing outcomes at two different dose scores help distinguish the intercept and slope?
  3. In the two-point posterior example, why did the data at \(x=0\) leave the posterior weights unchanged?
  4. Does a stable recommendation under several prior choices prove that both parameters are precisely estimated? Explain.

Solutions. (1) The likelihood depends on \(a+b\), the linear predictor at \(x=1\); parameter pairs with the same sum have equal likelihood. (2) The two scores supply information about two linear predictors, \(a+bx_1\) and \(a+bx_2\); when \(x_1\ne x_2\), these can distinguish the two parameters. (3) Both candidate pairs had the same \(p(0)=0.2000\), so their likelihoods for the observed outcomes were equal and their prior odds were unchanged. (4) No. A dose rule may remain stable even if the parameter posterior is broad; inspect uncertainty in predictions and parameters as well as the recommendation.

Clinical Biostats