Tutorials › Multi-Arm Multi-Stage (MAMS) Designs › Binary Early Response as a Selection Variable

Time-to-Event MAMS With Intermediate Outcomes · Tutorial 307 of 499

Binary Early Response as a Selection Variable

You will construct and compare fixed-time binary-response contrasts, quantify their uncertainty, and distinguish a useful selection ranking from a conclusion of treatment superiority.

Advanced 11 min read

What You'll Learn

  • Define a binary early-response estimand at a common assessment time.
  • Calculate arm-level response rates and treatment-versus-control risk differences.
  • Derive and apply standard errors for response-rate contrasts.
  • Account for the dependence created by a shared control arm.
  • Use intervals to distinguish a point-estimate ranking from a well-resolved ranking.
  • Recognise endpoint, missing-data, and small-sample problems that can undermine selection.

Why a Binary Early Response Can Guide Selection

In a multi-arm, multi-stage trial, final event outcomes may not yet be sufficiently mature when a decision about continued follow-up is needed. A binary response observed at a prespecified assessment time can provide an interim measure for comparing arms. The aim here is specific: estimate arm-level response rates, contrast each treatment with control, and use the estimates and their uncertainty to inform which arms continue. A response contrast can help rank arms, but it is not automatically a substitute for the final endpoint or proof that one treatment is superior.

A meaningful comparison requires a common response definition and assessment time. For example, the endpoint might be whether a participant has achieved a prespecified improvement by week 8. The response rule should specify how to classify values at the threshold, what to do after treatment discontinuation, and how death or other events before week 8 affect classification. Those choices define the outcome; they should not be made after looking at which choice favours a particular arm.

Definition: Let \(t^*\) be a prespecified assessment time and let \(Y_i(t^*)=1\) if participant \(i\) meets the trial's binary response definition by that time, and \(Y_i(t^*)=0\) otherwise. For arm \(j\), \(p_j(t^*)=P\{Y_i(t^*)=1\mid\text{assigned to arm }j\}\) is the response probability at \(t^*\). The treatment-versus-control response-rate contrast is the risk difference \(\Delta_j(t^*)=p_j(t^*)-p_0(t^*)\), where arm 0 is control. Positive values favour treatment when response is desirable.

The target population and analysis set must also be stated. In the usual intention-to-treat interpretation, the probability refers to participants assigned to an arm, regardless of whether they complete treatment. Each participant must have an outcome classified under the prespecified rule. Excluding participants with unavailable assessments can distort the arm-level rates if availability is related to response or treatment. A missing outcome is not automatically a nonresponse: the protocol must specify an appropriate endpoint or missing-data strategy, and the assumptions behind that strategy must be considered.

Estimating a Contrast and Its Uncertainty

Suppose \(n_j\) participants are assigned to arm \(j\), and \(x_j\) of them meet the binary response definition. The observed response rate is \(\hat p_j=x_j/n_j\), and the estimated treatment-versus-control risk difference is \(\hat\Delta_j=\hat p_j-\hat p_0\). Under independent participant outcomes within and between randomized arms, with a common response probability within each arm, the count satisfies \(X_j\sim\operatorname{Binomial}(n_j,p_j)\). The usual estimated standard error for the risk difference is

$$ SE(\hat\Delta_j)= \sqrt{\frac{\hat p_j(1-\hat p_j)}{n_j} +\frac{\hat p_0(1-\hat p_0)}{n_0}}. $$

The variance terms add because the treatment-arm and control-arm samples are independent. A simple approximate 95% confidence interval is \(\hat\Delta_j\pm1.96SE(\hat\Delta_j)\). Its endpoints are in probability-difference units: for example, 0.10 means a difference of 10 percentage points. The interval describes uncertainty in the risk difference under the sampling model; it is not a probability statement about where the fixed parameter lies.

A short derivation explains the standard error. For a Bernoulli outcome with response probability \(p_j\), the sample mean of \(n_j\) independent outcomes has variance \(p_j(1-p_j)/n_j\). Since the treatment and control sample means are independent, the variance of their difference is the sum of their variances:

$$ \operatorname{Var}(\hat p_j-\hat p_0) =\operatorname{Var}(\hat p_j)+\operatorname{Var}(\hat p_0) =\frac{p_j(1-p_j)}{n_j} +\frac{p_0(1-p_0)}{n_0}. $$

Replacing the unknown probabilities with the observed response rates gives the estimated variance in the standard-error formula. This approximation is most comfortable when arm sample sizes are not small and response rates are not close to 0 or 1. With sparse counts or boundary rates, Wald intervals can have poor coverage or even extend beyond the possible risk-difference range from -1 to 1. Score-based or suitable binomial intervals are preferable in those settings.

Worked Example: Two Arms Compared With a Shared Control

Suppose a hypothetical interim assessment at week 8 records responses for 100 participants in each of three arms. Control has 30 responses, treatment A has 45, and treatment B has 40. Assume all participants have a classified outcome, the response rule is identical across arms, and participant outcomes are independent within and between arms.

The observed rates are \(\hat p_0=30/100=0.30\), \(\hat p_A=45/100=0.45\), and \(\hat p_B=40/100=0.40\). Thus the estimated risk differences are \(\hat\Delta_A=0.45-0.30=0.15\), or 15 percentage points, and \(\hat\Delta_B=0.40-0.30=0.10\), or 10 percentage points.

For A, the estimated standard error is

$$ SE(\hat\Delta_A)= \sqrt{\frac{0.45(0.55)}{100}+\frac{0.30(0.70)}{100}} =\sqrt{0.004575}=0.06764. $$

The approximate 95% interval is \(0.15\pm1.96(0.06764)\), which is \(0.15\pm0.1326=(0.0174,0.2826)\). For B,

$$ SE(\hat\Delta_B)= \sqrt{\frac{0.40(0.60)}{100}+\frac{0.30(0.70)}{100}} =\sqrt{0.004500}=0.06708, $$

and the approximate 95% interval is \(0.10\pm1.96(0.06708)=0.10\pm0.1315=(-0.0315,0.2315)\). The point estimates rank A ahead of B, but both intervals are wide, and B's interval includes zero. These are interim selection summaries, not evidence that A is definitively better than B or that either response contrast predicts the final endpoint.

Shared Control and Uncertainty in the Ranking

When several treatments are contrasted with the same control, the estimated contrasts are correlated. In particular, under independent randomized arm samples,

$$ \operatorname{Cov}(\hat\Delta_j,\hat\Delta_k) =\operatorname{Var}(\hat p_0) =\frac{p_0(1-p_0)}{n_0},\qquad j\ne k. $$

Both differences contain the same estimated control rate with a minus sign, so an unusually high or low control rate shifts all treatment-control contrasts together. This dependence matters when comparing contrasts jointly; it is not correct to treat them as independent simply because the treatment arms differ. To compare two treatments directly, the shared control cancels: \(\hat\Delta_j-\hat\Delta_k=\hat p_j-\hat p_k\). Under independent treatment-arm samples, its estimated standard error is based on those two arms alone.

Worked Example: Is the Point-Estimate Ranking Resolved?

Consider hypothetical week-8 data with 120 participants per arm: 48 responses in control, 62 in A, and 55 in B. Suppose the response definition and outcome ascertainment are common across arms, and the binomial assumptions used above are reasonable.

The observed rates are \(0.400\), \(62/120=0.5167\), and \(55/120=0.4583\), respectively. The treatment-control differences are therefore \(0.1167\) for A and \(0.0583\) for B. By point estimate, A ranks ahead of B.

For a direct A-versus-B comparison, the estimated difference is \(0.5167-0.4583=0.0583\). Its standard error is

$$ SE(\hat p_A-\hat p_B)= \sqrt{\frac{0.5167(0.4833)}{120} +\frac{0.4583(0.5417)}{120}} =\sqrt{0.004150}=0.06442. $$

The approximate 95% interval is \(0.0583\pm1.96(0.06442)=0.0583\pm0.1263=(-0.0679,0.1846)\). It includes differences in either direction, so the observed ranking is uncertain. For context, the covariance between the two treatment-control contrasts is estimated as \(0.40(0.60)/120=0.002000\), arising from the shared control. The proper conclusion for selection is not that A has been shown superior, but that A has the larger observed rate while the data do not clearly distinguish the two arms.

An uncertainty-aware selection rule should say in advance how estimates and uncertainty will be used. One transparent approach is to rank by the estimated risk difference, report an interval for each contrast, and flag arms whose rankings are not clearly separated. If continued follow-up is limited to a fixed number of arms, the rule can use a prespecified practical threshold or a defined uncertainty criterion, rather than treating tiny differences between noisy estimates as decisive. The criterion should fit the purpose of selection; an interim ranking is not automatically a confirmatory significance test.

Worked Example: Equal Estimated Response Rates, Different Precision

Suppose a hypothetical assessment yields 24 responses among 80 control participants, 36 among 90 participants in A, and 30 among 75 participants in B. Assume the same fixed assessment time and response definition in every arm, with classified outcomes and independent participant responses.

The rates are \(\hat p_0=24/80=0.30\), \(\hat p_A=36/90=0.40\), and \(\hat p_B=30/75=0.40\). Both treatment arms have an estimated risk difference of \(0.10\), or 10 percentage points. For A,

$$ SE(\hat\Delta_A)= \sqrt{\frac{0.40(0.60)}{90}+\frac{0.30(0.70)}{80}} =\sqrt{0.0052916667}=0.07274, $$

giving the approximate interval \(0.10\pm1.96(0.07275)=(-0.0426,0.2426)\). For B,

$$ SE(\hat\Delta_B)= \sqrt{\frac{0.40(0.60)}{75}+\frac{0.30(0.70)}{80}} =\sqrt{0.005825}=0.07632, $$

giving \(0.10\pm1.96(0.07632)=(-0.0496,0.2496)\). The estimated rates tie, while A's interval is slightly narrower because its treatment arm is larger. The small precision difference alone does not establish that A is better; it is useful context if a prespecified selection rule must choose between otherwise similar candidates.

Assumptions, Diagnostics, and Alternatives

Before using these summaries for selection, check the design and data that make them meaningful. The assessment time should be the same across arms, and participants should have had a comparable opportunity to reach it. Confirm that the response threshold and rules for discontinuation and events before assessment were applied consistently. Check counts and denominators by arm, the extent and pattern of unavailable outcomes, and whether any data-quality issue could have changed classification. Also check whether the independent-binomial approximation is appropriate: clustering, repeated contributions from the same participant, or other dependence requires a variance method that reflects that structure.

A raw response rate among only those assessed is not necessarily an estimate of the randomized-arm response probability. If the proportion assessed differs by arm, or assessment availability is related to response, the observed proportions may be misleading. Presenting missing participants as nonresponders can be defensible only when that is part of the defined endpoint; it is not a neutral default. The selection decision should acknowledge the analysis strategy and its assumptions rather than hiding the missingness problem in a denominator.

The risk difference is directly interpretable in percentage points, but it is not the only possible effect scale. Risk ratios and odds ratios answer different questions; when treatment arms share a control, risk difference, risk ratio, and odds ratio preserve the same treatment-arm ranking, but rankings can differ when comparing effects across settings or strata with different baseline probabilities. The chosen scale should match the selection objective and be prespecified. For small samples or rates near the boundaries, score-based intervals or other suitable binomial methods are alternatives to the Wald interval. In a design with several treatment comparisons, the contrasts also share control information; methods that jointly account for that dependence are preferable to calculations that assume independent comparisons.

A further limitation is selection on noisy estimates. When the arm with the largest observed response is chosen, random variation tends to make that arm's observed rate look unusually favourable. The introductory material on post-selection estimation addresses the consequences of selecting on observed effects; here the practical safeguard is to preserve uncertainty in the selection summary and avoid presenting the winning point estimate as an unbiased measure of its effect. Finally, a binary response at \(t^*\) is only an intermediate selection measure. Its association with, or ability to predict, a later event outcome must not be assumed from the response-rate comparison alone.

Key takeaway: Define response at a common fixed time, estimate each arm's response rate and treatment-control risk difference, and report uncertainty alongside any ranking. Shared control makes treatment-control contrasts dependent, and a point-estimate leader is not necessarily a clearly better treatment.

Check Your Understanding

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

  1. In a hypothetical trial, control has 18 responses among 60 participants and treatment A has 27 among 60. Calculate the response-rate difference and its approximate standard error.
  2. If a control arm has 100 participants and an observed response rate of 0.20, what is the estimated covariance between two treatment-control response-rate differences using that control?
  3. Why can response rates measured at different assessment times give a misleading treatment ranking?
  4. If the largest observed treatment-control contrast has a confidence interval that includes zero, what can and cannot be concluded for interim selection?

Solutions. 1. The rates are \(18/60=0.30\) and \(27/60=0.45\), so the estimated difference is \(0.45-0.30=0.15\). Its standard error is \(\sqrt{0.45(0.55)/60+0.30(0.70)/60}=\sqrt{0.007625}=0.08732\). 2. The estimated covariance is \(0.20(0.80)/100=0.001600\); it is positive because both contrasts use the same estimated control rate. 3. Participants observed for different durations have different opportunities to meet the response definition, so the rates do not describe the same fixed-time outcome. 4. The arm has the largest observed contrast, but the interval indicates substantial uncertainty and does not establish a positive treatment-control difference at the stated confidence level. A selection rule may still retain it for further follow-up, but the ranking should not be described as proof of superiority.

Clinical Biostats