Tutorials › Adaptive Clinical Trial Design › Selecting the EWOC Feasibility Quantile

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

Selecting the EWOC Feasibility Quantile

See how the EWOC feasibility quantile turns a chosen safety threshold into an allowable dose bound, including how to handle discrete posterior MTD distributions.

Advanced 10 min read

What You'll Learn

  • Define the EWOC feasibility function using the strict event that the MTD is below a proposed dose.
  • Interpret the feasibility quantile as a control on posterior overdose risk, not as a toxicity probability for an individual.
  • Predict how increasing or decreasing the feasibility threshold changes allowable doses.
  • Calculate feasibility bounds under continuous and discrete posterior MTD distributions.
  • Distinguish the strict-event feasibility function from the usual right-continuous CDF at posterior atoms.
  • Recognize the assumptions and practical limitations involved in choosing a feasibility threshold.

What the Feasibility Quantile Controls

In EWOC, a proposed dose is assessed against uncertainty about the maximum tolerated dose (MTD). If the proposed dose \(d\) exceeds the true MTD \(M\), it is treated as being above the MTD. The EWOC feasibility rule limits the posterior probability of that event. Its threshold, often denoted by \(\alpha\), is therefore a safety-aggressiveness control: smaller values restrict the trial to doses with less posterior probability of exceeding the MTD, while larger values permit more risk under the model.

The threshold does not come from a posterior calculation alone. It is a design choice that must be made in the clinical and operational context. Once the posterior distribution of \(M\) is fixed, changing \(\alpha\) changes which doses satisfy the rule; it does not change that posterior distribution.

Definition: Given the trial data \(D\), define the strict-event feasibility function at dose \(d\) by \(G(d\mid D)=P(M<d\mid D)\). For a chosen feasibility threshold \(0<\alpha<1\), a dose is feasible if \(G(d\mid D)\leq\alpha\). On a bounded clinical dose range \([d_{\min},d_{\max}]\), if \(G(d_{\min}\mid D)\leq\alpha\), the EWOC feasibility bound is \(b_\alpha=\sup\{d\in[d_{\min},d_{\max}]:G(d\mid D)\leq\alpha\}\). If \(G(d_{\min}\mid D)>\alpha\), no dose on the clinical range is feasible, \(b_\alpha\) is undefined, and the protocol must specify an operational response.

The event is written \(M<d\), rather than \(M\leq d\): equality is not classified as a dose above the MTD under this convention. Because \(G(d\mid D)\) is nondecreasing in \(d\), feasible doses form an initial portion of the dose range. When the bound is within that range, doses at or below it satisfy the rule, and doses above it do not. The bound can be understood as a generalized lower-tail quantile of the posterior MTD distribution, with the strict-event convention made explicit.

This is a posterior statement conditional on the data and the model. It is not the probability that a particular patient will experience toxicity, nor is it a frequentist guarantee that the trial will never treat above the true MTD. Its interpretation depends on whether the model’s MTD is a meaningful safety boundary for the trial.

How the Threshold Changes the Bound

The central comparison follows directly from the feasibility condition. If \(\alpha_1\leq\alpha_2\), then every dose feasible at \(\alpha_1\) is also feasible at \(\alpha_2\), because \(G(d\mid D)\leq\alpha_1\) implies \(G(d\mid D)\leq\alpha_2\). Thus the feasible set at the smaller threshold is contained in the set at the larger threshold, and, whenever both bounds are defined, \(b_{\alpha_1}\leq b_{\alpha_2}\). The change need not be smooth: the bound may remain unchanged over a range of thresholds and then jump, particularly when the posterior distribution has point masses.

For a continuous posterior distribution, let \(H(d\mid D)=P(M\leq d\mid D)\) be its usual CDF. With no point masses, \(G(d\mid D)=H(d\mid D)\). If this CDF is strictly increasing across the relevant range and \(\alpha\) is attained at an interior point of the clinical range, the boundary solves \(H(b_\alpha\mid D)=\alpha\). If \(\alpha>H(d_{\max}\mid D)\), every dose on the range is feasible and \(b_\alpha=d_{\max}\); if no dose is feasible, the bound is undefined as specified above. In that case, the feasibility bound is the ordinary \(\alpha\)-quantile of the posterior MTD distribution. This equality is not a safe shortcut for distributions with atoms.

Important distinction: The usual CDF is \(H(d\mid D)=P(M\leq d\mid D)\), while EWOC here uses \(G(d\mid D)=P(M<d\mid D)\). At an atom, \(G(d\mid D)=H(d-\mid D)\), where \(H(d-\mid D)\) is the probability strictly below \(d\). A quantile computed from \(H\) must therefore be checked against \(G\) before it is used as an EWOC boundary.

A short justification explains why the strict-event boundary is well behaved on a bounded dose range. As \(d\) increases, the event \(\{M<d\}\) can only grow, so \(G\) is nondecreasing. Also, as doses approach a value \(d\) from below, the events \(\{M<x\}\) increase to \(\{M<d\}\), so their probabilities approach \(G(d\mid D)\). Consequently, if feasible doses approach a finite supremum within the dose range, the limiting boundary itself remains feasible. This left-continuity is why the strict-event rule can include the boundary even when the posterior has a jump there.

When both bounds are defined: $$ \alpha_1\leq\alpha_2 \quad\Longrightarrow\quad \{d:G(d\mid D)\leq\alpha_1\} \subseteq \{d:G(d\mid D)\leq\alpha_2\} \quad\Longrightarrow\quad b_{\alpha_1}\leq b_{\alpha_2}. $$

Worked Example: A Continuous Posterior MTD Distribution

Suppose, hypothetically, that the posterior distribution of the MTD is uniform from 20 to 80 mg/m². This is a simplified model for illustrating the threshold; assume the clinical dose range includes those values. For doses between 20 and 80 mg/m², the strict-event probability is \(G(d\mid D)=(d-20)/60\). Below or at 20 it is zero, and above 80 it is one.

For \(\alpha=0.25\), solve the boundary equation:

$$ \frac{b_{0.25}-20}{60}=0.25 \quad\Longrightarrow\quad b_{0.25}=20+60(0.25)=35\ \text{mg/m}^2. $$

At 35 mg/m², \(G(35\mid D)=(35-20)/60=15/60=0.25\), so the boundary satisfies the rule. A dose of 50 mg/m² has \(G(50\mid D)=30/60=0.50\), so it is not feasible at this threshold.

If the threshold is raised to \(\alpha=0.50\), the same calculation gives \(b_{0.50}=20+60(0.50)=50\) mg/m². At \(\alpha=0.75\), it gives \(b_{0.75}=20+60(0.75)=65\) mg/m². Thus, with the posterior held fixed, increasing the threshold allows a higher boundary: the posterior probability of the proposed dose exceeding the MTD may be larger.

Posterior Atoms and Quantile Conventions

A posterior distribution may have point masses, for example when model structure or posterior computation leaves probability concentrated on a finite set of MTD values. At an atom, the usual inverse-CDF quantile and the EWOC feasibility boundary can differ. This does not invalidate the rule; it means that feasibility must be checked using the event specified by the rule.

Worked Example: Checking Feasibility at Posterior Atoms

Suppose, hypothetically, that the posterior MTD equals 10 mg/m² with probability 0.20, 20 mg/m² with probability 0.50, and 30 mg/m² with probability 0.30. For doses \(d\) in successive intervals, the strict-event function is:

$$ G(d\mid D)= \begin{cases} 0, & d\leq 10,\\ 0.20, & 10<d\leq 20,\\ 0.70, & 20<d\leq 30,\\ 1, & d>30. \end{cases} $$

With \(\alpha=0.20\), every dose up to and including 20 mg/m² is feasible, since \(G(20\mid D)=P(M<20\mid D)=0.20\). Any dose just above 20 mg/m² has \(G(d\mid D)=0.70\), so the feasibility bound is 20 mg/m². The mass at 20 is excluded from \(G(20\mid D)\), but included as soon as the proposed dose is greater than 20.

Now use \(\alpha=0.50\). The bound is still 20 mg/m²: the strict-event probability is 0.20 at 20, but jumps to 0.70 immediately above 20. In contrast, the usual right-continuous CDF has \(H(10\mid D)=0.20\) and \(H(20\mid D)=0.70\), so its conventional 0.50 quantile is 20 mg/m². In this case that value is also the feasibility bound, but that agreement does not hold in every case.

For example, return to \(\alpha=0.20\). The usual 0.20 quantile is 10 mg/m², while the strict-event EWOC bound is 20 mg/m². The relevant check is \(G(20\mid D)=0.20\), not \(H(20\mid D)=0.70\). Using the usual CDF without attending to the strict event can therefore give a different answer.

The example also shows why a reported quantile should include its convention. One common generalized inverse of the usual CDF is \(Q_\alpha^H=\inf\{d:H(d\mid D)\geq\alpha\}\). An EWOC bound, by contrast, is defined through the largest dose satisfying \(G(d\mid D)\leq\alpha\). Under a continuous posterior, these descriptions coincide when the equation \(H(d\mid D)=\alpha\) has a unique solution and that solution lies in the clinical range. If \(H\) is flat at the threshold, they can differ; with atoms, calculate or inspect the strict-event probability directly.

Choosing a Threshold in Practice

The posterior distribution determines what each candidate threshold permits, but it does not by itself determine which threshold is appropriate. Choosing \(\alpha\) requires a protocol-level judgment about the acceptable posterior risk of dosing above the MTD. A smaller \(\alpha\) is more conservative under this particular rule; a larger \(\alpha\) is more permissive. Neither direction is universally correct without considering the trial’s clinical context and objectives.

A disciplined choice separates two tasks. First, specify the acceptable risk threshold before interpreting the resulting bound as a recommendation. Second, calculate the bound from the current posterior and apply any other protocol requirements separately. Those may include dose availability, starting-dose rules, or restrictions on escalation. The feasibility threshold alone does not settle every operational decision.

Worked Example: Translating a Protocol Risk Limit

Suppose a hypothetical protocol specifies that a proposed dose must have no more than a 0.30 posterior probability of being above the MTD. The current posterior MTD distribution is uniform from 20 to 80 mg/m², as in the earlier calculation. Here the protocol’s risk limit corresponds to \(\alpha=0.30\); the limit is a design input, not an estimate extracted from the posterior.

The feasibility boundary is:

$$ b_{0.30}=20+60(0.30)=38\ \text{mg/m}^2. $$

At 38 mg/m², \(G(38\mid D)=(38-20)/60=18/60=0.30\), so the dose is at the permitted limit. At 44 mg/m², \(G(44\mid D)=(44-20)/60=24/60=0.40\), which exceeds 0.30; that dose fails this feasibility rule. A separate operational restriction could still rule out 38 mg/m², but it would not change the calculation of this posterior risk bound.

Diagnostics, Alternatives, and Limitations

Before using a computed boundary, check that the posterior MTD distribution and dose scale are being interpreted consistently. The model must provide a meaningful MTD on an ordered dose scale, and the event \(M<d\) must represent the safety concern the trial intends to control. If dose-response behavior is nonmonotone or the model’s MTD is poorly defined, the feasibility calculation may not reflect the clinically relevant risk. The posterior also depends on the model and its assumptions; a numerically precise quantile does not remove model uncertainty.

For posterior draws \(M^{(1)},\ldots,M^{(B)}\), a direct approximation at a candidate dose \(d\) is:

$$ \widehat{G}(d\mid D)=\frac{1}{B}\sum_{b=1}^{B} I\{M^{(b)}<d\}, $$

where \(I\{\cdot\}\) equals 1 when its condition holds and 0 otherwise. This estimate should be checked near the proposed boundary, where Monte Carlo error can affect whether a dose appears to meet the threshold. With discrete posterior draws, use the strict comparison shown: replacing it with \(M^{(b)}\leq d\) changes the calculation when draws equal the candidate dose.

EWOC’s bound is not the same as choosing a dose whose predicted toxicity probability is closest to a target. That kind of decision asks about toxicity at the dose; EWOC as defined here controls posterior probability that the dose lies above the MTD. A fixed dose cap is another alternative, but it does not adapt to the current posterior in the same way. In practice, such restrictions may be combined, with each rule retaining its own interpretation.

Common mistake: “The ordinary \(\alpha\)-quantile is always feasible because its CDF equals \(\alpha\).” This reasoning uses \(H(d\mid D)=P(M\leq d\mid D)\) as though it were the EWOC function. At an atom, EWOC uses \(G(d\mid D)=P(M<d\mid D)=H(d-\mid D)\); evaluate \(G\) at the candidate dose and verify \(G(d\mid D)\leq\alpha\).

A threshold should not be interpreted as a guarantee about realized patient outcomes. It controls a posterior probability under the assumed model and current data. If the model is misspecified, if toxicity is not adequately summarized by the MTD, or if operational rules are applied inconsistently, the nominal feasibility interpretation may not provide the intended protection. Threshold selection is therefore a design judgment, while the posterior calculation is the conditional statistical assessment.

Key takeaway: EWOC feasibility quantiles translate a chosen upper limit on \(P(M<d\mid D)\) into an allowable dose boundary. Raising \(\alpha\) cannot reduce that boundary, but at posterior atoms the boundary must be defined and checked using the strict-event probability, not assumed from the usual CDF quantile.

Check Your Understanding

Answer each question using the strict-event EWOC definition. Solutions follow the questions.

  1. If the feasibility threshold changes from 0.20 to 0.40 while the posterior MTD distribution remains fixed, can the feasible set become smaller? Explain.
  2. Suppose \(M\) is uniform from 10 to 50 mg/m². What is the feasibility bound at \(\alpha=0.25\)?
  3. Suppose \(P(M=5)=0.30\) and \(P(M=15)=0.70\). Calculate \(G(15\mid D)\). Is 15 mg/m² feasible at \(\alpha=0.30\)?
  4. Why can the ordinary CDF quantile fail to equal the EWOC feasibility bound when the posterior has an atom?
  5. Does an EWOC threshold of 0.25 mean that each patient treated at the boundary has a 25% probability of toxicity? Explain.

Solutions. (1) No. Every dose satisfying \(G(d\mid D)\leq0.20\) also satisfies \(G(d\mid D)\leq0.40\), so the feasible set can only expand or stay the same. (2) For a uniform distribution on \([10,50]\), \(G(d\mid D)=(d-10)/40\) within the interval. Thus \(b_{0.25}=10+40(0.25)=20\) mg/m². (3) \(G(15\mid D)=P(M<15\mid D)=0.30\), since only the mass at 5 is strictly below 15. The dose is feasible at \(\alpha=0.30\), including equality in the rule. (4) The ordinary CDF includes mass at the candidate dose, whereas the strict-event function excludes it; they can therefore assign different probabilities at an atom. (5) No. The threshold limits the posterior probability, conditional on the data and model, that the proposed dose exceeds the MTD. It is not an individual patient’s probability of toxicity.

Clinical Biostats