Tutorials › Model-Informed Drug Development (MIDD) › Binary Responder Endpoints and Threshold Selection

Model-Based Endpoint and Sample Size Planning · Tutorial 315 of 596

Binary Responder Endpoints and Threshold Selection

Learn to connect a clinically justified responder threshold to event probabilities, treatment contrasts, and defensible power and sample-size calculations.

Advanced 13 min read

What You'll Learn

  • Define a responder event from a continuous outcome and a clinically justified threshold.
  • Examine how threshold changes affect response probabilities and treatment contrasts.
  • Calculate approximate equal-allocation sample size for a risk-difference contrast.
  • Plan sample size using risk-ratio or odds-ratio scales and interpret each in context.
  • Check assumptions and recognize when binary-endpoint approximations may be unreliable.

From a Clinical Outcome to a Responder Event

A responder endpoint turns an outcome into a binary event: each participant either meets a prespecified response criterion or does not. This can make a result easier to communicate, but the event is not defined by the data alone. Its threshold determines which outcomes count as responses, which in turn determines the event probabilities, treatment contrast, power, and required sample size.

The threshold should represent a clinically meaningful outcome, not merely a value that happens to produce a favorable analysis. For planning, the responder definition, assessment time, and handling of the underlying outcome must be specified in terms consistent with the estimand. The earlier material on estimands and model-based endpoint contrasts is relevant here: a carefully computed contrast is not interpretable unless it represents the intended treatment effect.

Definition: Let \(X\) be an outcome oriented so that larger values indicate greater benefit, and let \(c\) be a prespecified threshold. The responder indicator is \(Y=I(X\ge c)\), equal to 1 when the threshold is met and 0 otherwise. In treatment group \(a\), the response probability is \(p_a(c)=P(X\ge c\mid A=a)\). The threshold and assessment rule define the binary endpoint.

For an outcome where smaller values are better, the event would instead be defined with the corresponding inequality, such as \(Y=I(X\le c)\). The direction must be explicit. Response probabilities depend on the threshold: raising a threshold for a benefit-oriented outcome cannot increase the response probability in either group. The treatment contrast may change as well, and need not move in the same direction as either group’s response probability.

Threshold Choice and Its Consequences

A threshold is clinically justified when it distinguishes a meaningful outcome from one that does not meet the intended response criterion. Evidence for that choice may come from clinical interpretation of the outcome scale, external evidence, or a prespecified clinical rationale. Statistical convenience alone—such as choosing the threshold that maximizes an observed treatment contrast—does not establish clinical meaning and risks an optimistic estimate of benefit.

The following hypothetical example illustrates the arithmetic. The normal distributions are a planning assumption for the continuous outcome, not a claim that all clinical outcomes follow a normal distribution.

Worked Example: How the Threshold Changes Response Rates

Suppose a continuous improvement score \(X\), in points, is approximately normally distributed at the assessment time. In a hypothetical control group, \(X\sim N(4,10^2)\); in a hypothetical treatment group, \(X\sim N(10,10^2)\). A clinical rationale supports defining response as an improvement of at least 8 points.

For a normal variable, standardize the threshold using \(Z=(X-\mu)/\sigma\), where \(\sigma=10\). Thus, the control response probability is \(P(X\ge8)=P(Z\ge(8-4)/10)=P(Z\ge0.4)\approx0.3446\). The treatment probability is \(P(X\ge8)=P(Z\ge(8-10)/10)=P(Z\ge-0.2)\approx0.5793\).

The risk difference, treatment minus control, is \(0.5793-0.3446=0.2347\), or 23.47 percentage points. The risk ratio is \(0.5793/0.3446\approx1.681\). The odds ratio is \(\{0.5793/(1-0.5793)\}/\{0.3446/(1-0.3446)\}\approx2.619\). These are three different summaries of the same pair of response probabilities.

If the threshold were raised to 12 points, the probabilities would be \(P(X\ge12\mid control)=P(Z\ge0.8)\approx0.2119\) and \(P(X\ge12\mid treatment)=P(Z\ge0.2)\approx0.4207\). The risk difference would then be \(0.2088\), rather than \(0.2347\). This threshold comparison shows why the response definition must be chosen for clinical reasons and established before examining treatment outcomes.

Power and Sample Size on the Risk-Difference Scale

Let \(p_1\) and \(p_0\) be the planned response probabilities in the treatment and control groups. The risk difference is \(\Delta=p_1-p_0\), measured in probability units; multiplying by 100 expresses it in percentage points. For two independent groups with equal sample size \(n\), a common approximate calculation uses a two-sided test of \(H_0:p_1-p_0=0\), significance level \(\alpha\), and target power \(1-\beta\).

The approximation below uses a normal critical value, the pooled response probability for the null standard error, and the two separate probabilities for the standard error under the planned alternative. It assumes independent participant outcomes, fixed group allocation, and sufficiently large expected responder and nonresponder counts for the normal approximation to be reasonable.

Planning formula: With equal allocation, put \(\bar p=(p_1+p_0)/2\) and \(\Delta=p_1-p_0\). An approximate sample size per group is
$$ n\approx \frac{\left[ z_{1-\alpha/2}\sqrt{2\bar p(1-\bar p)} +z_{1-\beta}\sqrt{p_1(1-p_1)+p_0(1-p_0)} \right]^2}{\Delta^2}. $$
Here \(z_q\) is the \(q\)-quantile of the standard normal distribution. The result is rounded up to a whole participant per group.

This formula is a planning approximation, not a guarantee of the power of every particular analysis procedure. It should be paired with an analysis plan that tests the same contrast and uses the same endpoint definition. For small samples or probabilities near 0 or 1, use a method suited to the planned test or confirm performance by simulation rather than relying uncritically on a normal approximation.

Worked Example: Sample Size for a Risk Difference

Suppose a hypothetical trial plans for a control response probability \(p_0=0.40\) and a treatment response probability \(p_1=0.55\). The clinically relevant difference is therefore \(\Delta=0.15\), or 15 percentage points. Plan for two-sided \(\alpha=0.05\) and 80% power, with equal allocation.

The pooled planning probability is \(\bar p=(0.40+0.55)/2=0.475\). Use \(z_{0.975}=1.96\) and \(z_{0.80}=0.8416\). Substitution gives

$$ n\approx \frac{\left[ 1.96\sqrt{2(0.475)(0.525)} +0.8416\sqrt{(0.55)(0.45)+(0.40)(0.60)} \right]^2}{(0.55-0.40)^2}. $$

The first square root is \(\sqrt{0.49875}\approx0.7062\), and the second is \(\sqrt{0.4875}\approx0.6982\). Thus,

$$ n\approx\frac{(1.96(0.7062)+0.8416(0.6982))^2}{0.0225} =\frac{(1.3842+0.5876)^2}{0.0225} \approx\frac{3.888}{0.0225} \approx172.8. $$

Round up to 173 participants per group, or 346 total, before any additional allowance needed for the trial’s operational assumptions. Under the stated planning probabilities and approximation, this is the estimated sample size for 80% power to detect a 15-percentage-point difference with a two-sided 5% test. It is not a claim that the realized power will be exactly 80%.

Planning on Relative Scales

A risk ratio compares response probabilities multiplicatively: \(RR=p_1/p_0\). An odds ratio compares odds, where the odds in a group are \(p/(1-p)\). These measures are not interchangeable with a risk difference. For example, the same relative effect can correspond to different absolute differences at different baseline probabilities. In addition, odds ratios can differ from risk ratios even when the outcome is not rare; choosing an odds-ratio scale is not simply a relabeling of the risk ratio.

For planning a risk ratio with a large-sample Wald test of \(H_0:RR=1\), the approximate variance of the estimated log risk ratio in two independent groups is

$$ \operatorname{Var}(\log\widehat{RR})\approx \frac{1-p_1}{n p_1}+\frac{1-p_0}{n p_0} =\frac{1/p_1+1/p_0-2}{n}. $$

The expression follows from the approximate variance of a sample proportion, \(p(1-p)/n\), and the delta method for the logarithm: the derivative of \(\log(p)\) is \(1/p\), so the approximate variance of \(\log(\hat p)\) is \(\{p(1-p)/n\}/p^2=(1-p)/(np)\). Independent groups contribute additive variances. Under the planned alternative, the log risk ratio is \(\log(RR)\); dividing its squared value by the variance gives the approximate squared signal-to-noise ratio. Setting this ratio to \((z_{1-\alpha/2}+z_{1-\beta})^2\) yields the planning expression below.

Planning formula: For equal group sizes and a two-sided Wald test of \(RR=1\),
$$ n\approx \frac{(z_{1-\alpha/2}+z_{1-\beta})^2 \left(1/p_1+1/p_0-2\right)} {\{\log(p_1/p_0)\}^2}. $$
This is a large-sample approximation and requires positive response probabilities in both groups. It should not be used without scrutiny when the planned event counts are sparse.

Worked Example: Sample Size for a Risk Ratio

In a hypothetical study, assume \(p_0=0.20\) and \(p_1=0.30\), so the target \(RR=1.50\). Plan for two-sided \(\alpha=0.05\) and 80% power. The variance coefficient is \(1/0.30+1/0.20-2=3.3333+5-2=6.3333\), and \(\log(1.50)\approx0.4055\).

Using \(z_{0.975}=1.96\) and \(z_{0.80}=0.8416\), the sum of critical values is approximately 2.802. The calculation is

$$ n\approx \frac{(2.802)^2(6.3333)}{\{\log(1.50)\}^2} \approx\frac{49.72}{0.1644} \approx302.4. $$

The required approximate sample size is 303 participants per group, or 606 total. The target is a 50% relative increase in response probability, from 20% to 30%; the corresponding absolute difference is 10 percentage points. The calculation assumes independent binary outcomes and an adequate large-sample approximation for the log risk ratio.

For an odds ratio \(\theta\), a planned control probability \(p_0\) and target odds ratio imply the treatment probability

$$ p_1=\frac{\theta p_0}{1-p_0+\theta p_0}. $$

For a large-sample Wald calculation of \(\log(\widehat{\theta})\), its approximate variance with equal group sizes is \(\{1/[p_1(1-p_1)]+1/[p_0(1-p_0)]\}/n\). The resulting approximate sample size is the sum of these variance coefficients multiplied by \((z_{1-\alpha/2}+z_{1-\beta})^2\), then divided by \(\{\log(\theta)\}^2\). This relies on adequate counts in all four cells of the two-by-two table.

Worked Example: Translating an Odds Ratio Into a Probability

Suppose a hypothetical planning target is an odds ratio of \(\theta=1.80\) with control response probability \(p_0=0.20\). First convert the odds ratio to the implied treatment probability:

$$ p_1=\frac{(1.80)(0.20)}{1-0.20+(1.80)(0.20)} =\frac{0.36}{1.16} \approx0.3103. $$

The corresponding risk difference is about \(0.3103-0.20=0.1103\), or 11.03 percentage points. For two-sided \(\alpha=0.05\) and 80% power, the variance coefficient is \(1/[0.3103(0.6897)]+1/[0.20(0.80)]\approx4.672+6.250=10.922\). Therefore,

$$ n\approx \frac{(1.96+0.8416)^2(10.922)}{\{\log(1.80)\}^2} \approx\frac{85.72}{0.3455} \approx248.1. $$

Round up to 249 participants per group. This calculation targets an odds ratio of 1.80—not a risk ratio of 1.80—and its treatment response probability is derived from the specified control probability. The large-sample Wald approximation should be checked if the actual planned event counts are small.

Diagnostics, Alternatives, and Limitations

Before accepting a calculation, check that the endpoint, analysis contrast, and planning assumptions align. The threshold and assessment time should be clinically defensible and fixed in advance. The planned probabilities should be plausible for the intended population and treatment context. Check expected responders and nonresponders in each group; for instance, the risk-ratio and odds-ratio approximations can become unstable when a probability is near zero or one. Also verify that the sample-size formula matches the planned test, sidedness, allocation ratio, and significance level.

The normal approximations are less dependable with sparse data, severe imbalance, or very high or low probabilities. Alternatives include score-based procedures for comparing two proportions, exact procedures for sparse tables, and regression models when covariate adjustment is planned. They may have different finite-sample operating characteristics, so sample size should be based on the intended analysis rather than chosen from a formula solely because it is convenient. For complex designs or uncertain response probabilities, simulation can estimate power directly by generating binary outcomes under the planned design and applying the planned analysis to each simulated trial.

A binary endpoint also discards distinctions among participants on the same side of the threshold. Someone just below the response cutoff is classified with someone whose outcome is far worse, while two outcomes just on opposite sides of the cutoff are treated as different response categories. If the full outcome scale is clinically informative, a continuous-endpoint analysis may retain more information; the appropriate choice depends on the estimand, interpretation, and design. Threshold-based results can be sensitive to the cutoff, so sensitivity analyses using clinically plausible thresholds can help describe that dependence, but they do not justify selecting the most favorable threshold after observing results.

  • Pitfall: “The risk ratio is 1.5, so the treatment adds 50 percentage points to the response rate.” This confuses a relative increase with an absolute difference. If the control probability is 0.20, an RR of 1.5 implies a treatment probability of 0.30 and a 10-percentage-point difference.
  • Correction: State the scale and its units: a risk difference is in probability units or percentage points, while a risk ratio and odds ratio are dimensionless relative measures. Translate to probabilities where possible for clinical interpretation.
  • Pitfall: Choosing the threshold that yields the smallest planned sample size or the largest observed treatment contrast.
  • Correction: Justify the cutoff clinically before using trial outcomes to compare groups; assess how plausible alternative cutoffs affect planning without treating them as outcome-driven replacements.
Key takeaway: A responder sample-size calculation begins with a clinically meaningful, prespecified event definition. The planned event probabilities and effect scale then determine the contrast and its uncertainty; assumptions and approximation quality must be checked against the intended analysis.

Check Your Understanding

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

  1. For a benefit-oriented score, what does \(Y=I(X\ge c)\) record, and why should \(c\) be justified before examining treatment outcomes?
  2. If \(p_0=0.25\) and \(p_1=0.40\), calculate the risk difference and risk ratio. State what each measure means.
  3. Why can a risk ratio of 1.50 not be interpreted as a 50-percentage-point increase?
  4. Name two conditions that can make a large-sample Wald sample-size approximation unreliable.
  5. Why is a threshold chosen to maximize the observed treatment contrast not, by itself, a clinically justified threshold?

Solutions. (1) The indicator equals 1 for a participant whose outcome is at least \(c\), and 0 otherwise. The cutoff determines the event probabilities and therefore the clinical meaning and statistical contrast; outcome-driven selection can make the apparent effect optimistic. (2) The risk difference is \(0.40-0.25=0.15\), or 15 percentage points. The risk ratio is \(0.40/0.25=1.60\), meaning the treatment response probability is 1.60 times the control probability. (3) A ratio is relative, not an absolute difference: at a control probability of 0.25, an RR of 1.50 implies a treatment probability of \(1.50(0.25)=0.375\), a difference of 0.125, or 12.5 percentage points. (4) Examples include sparse expected event counts, probabilities near 0 or 1, and severe group imbalance. (5) Maximizing the observed contrast is a statistical selection rule, not evidence that the cutoff distinguishes clinically meaningful outcomes; it can also capitalize on random variation.