Why Pair an Early Response With Overall Survival?
A multi-arm, multi-stage design can use a response measured at a fixed time to decide which experimental treatments should continue. The attraction is practical: response may be observed earlier than enough deaths have accrued for an overall-survival comparison. The statistical challenge is that the early outcome is a selection tool, while overall survival remains the final endpoint. A treatment that looks promising on response is not thereby shown to improve survival.
This tutorial considers three randomized arms: experimental treatments A and B, and shared control C. At a prespecified time \(\tau\), the trial applies a response rule to the experimental arms. Arms that meet the rule continue to an overall-survival comparison with control. The rule, the analysis population, the handling of missing response assessments, and the final inferential procedure must be specified before outcomes are examined.
For example, a protocol might define response as a confirmed tumour response by week 12, with death or progression before week 12 classified as nonresponse. That is only one possible definition: the endpoint and the rule for unevaluable assessments must match the clinical question and be fixed in advance. The survival endpoint should also specify its time origin, event definition, and censoring rules.
The usual response-rate calculation treats randomized participants as the analysis population and requires outcome classification to be sufficiently complete and comparable across arms. If response is missing, silently excluding those participants can bias the comparison; a prespecified missing-data rule or sensitivity analysis is needed. For an ordinary logrank analysis of overall survival, participants should be randomized comparably, follow-up and event ascertainment should be reliable, and censoring should be independent of the event time within the analysis groups, at least conditional on variables used in the analysis. A hazard-ratio interpretation also relies on proportional hazards; the logrank test itself can be less informative when treatment effects vary over time.
Specify the Rule Before Seeing the Data
A simple rule compares each experimental arm's observed response rate with the control response rate. For instance, retain arm \(i\) if \(\hat p_i-\hat p_C\geq \delta\), where \(\delta\) is a clinically chosen minimum difference. This is a decision threshold, not automatically a hypothesis test: unless the protocol says otherwise, it does not require a response p-value to cross a significance boundary. A rule that selects the top-ranked experimental arm is another option, but it can retain an arm even when both experimental response rates are poor.
At selection, a dropped arm generally stops recruiting, while already randomized participants continue to be followed as specified. The final comparison should use the participants and control data allowed by the protocol; excluding participants merely because their arm was dropped can distort the planned analysis. Selection also changes the distribution of the final test statistic. The material on familywise error with intermediate-based selection and operating characteristics of intermediate-based rules is directly relevant: final error rates and power must reflect the complete selection-and-testing procedure, not just an isolated survival test.
Worked Example: Apply a Fixed-Time Response Threshold
Suppose a hypothetical trial assesses response at week 12 in 150 participants per arm. There are 66 responders in A, 54 in B, and 45 in C. The prespecified rule retains an experimental arm only if its observed response rate is at least 0.08 higher than the control rate.
The observed rates are
The estimated differences from control are \(0.440-0.300=0.140\) for A and \(0.360-0.300=0.060\) for B. Thus A exceeds the 0.08 rule and continues; B does not and is dropped from further accrual. The decision follows the prespecified rule. It does not establish that A improves overall survival, nor does it establish that B has no benefit.
As a descriptive check, the unpooled standard error for the observed A-versus-C response-rate difference is
The observed difference is about \(0.140/0.05516=2.538\) standard errors from zero under this approximate calculation. That descriptive statistic is not the selection rule here: the protocol specified an 0.08 difference threshold. Substituting a test-based rule after seeing the data would change the design and its operating characteristics.
Continue Follow-Up and Analyze Overall Survival
For a standard two-group logrank comparison of selected arm A with control C, each observed death contributes to the observed event count and to an expected count under the null hypothesis of equal survival distributions. Let \(O_A\) be the observed number of deaths in A, \(E_A\) the number expected in A under the null, and \(V_A\) the variance of the observed-minus-expected count, calculated from the risk sets at event times. The approximate standardized logrank statistic is
With the sign convention above, a negative value means fewer observed deaths in A than expected under the null, consistent with a possible survival advantage for A. Under the null and usual large-sample conditions, \(Z\) is approximately standard normal. This supplies a nominal p-value for the comparison; in a design that selected A using response, the final decision rule may need a selection-adjusted boundary or calibration from the full design.
Worked Example: Calculate the Selected-Arm Logrank Statistic
Continue the hypothetical trial above. At the planned overall-survival analysis, suppose the logrank calculations for A versus C give \(O_A=80\), \(E_A=90\), and \(V_A=45\). These are hypothetical event-analysis summaries; the expected count and variance come from the event-time risk sets, not simply from the total sample sizes.
The logrank statistic is
The approximate two-sided nominal p-value is \(2\Phi(-1.491)\approx 0.136\). At a nominal two-sided 0.05 level, this result is not statistically significant. It is not evidence that the two survival distributions are identical; the data are compatible with a range of effects, and the study may have limited information.
A score-based approximation to the log hazard ratio is \((O_A-E_A)/V_A=-10/45=-0.2222\), giving an approximate hazard ratio of \(\exp(-0.2222)=0.8007\). Its approximate standard error on the log scale is \(1/\sqrt{45}=0.1491\). A corresponding approximate 95% interval is
Under a proportional-hazards interpretation, the point estimate corresponds to an estimated death hazard about 20% lower in A than C over follow-up. The interval includes 1, so it does not rule out no hazard difference at the conventional two-sided 5% level. These score-based values are approximations, not a substitute for the prespecified final model. Moreover, the nominal p-value and interval do not by themselves account for selecting A based on response.
Justification and Design Checks
The logrank standardization follows from the null distribution of its score. At each death time, the risk set determines how many deaths would be expected in A if the hazards were equal. Summing the observed-minus-expected contributions gives \(O_A-E_A\); the corresponding variance contributions sum to \(V_A\). Under the null, this score has mean approximately zero and variance approximately \(V_A\). Standardizing therefore gives an approximately unit-variance statistic:
This approximation is most dependable when there is adequate event information and no severe violation of the assumptions behind the planned analysis. For interpreting the logrank comparison as a constant hazard ratio, proportional hazards is an additional substantive assumption. If survival curves cross or effects emerge only after a delay, a single hazard ratio can be difficult to interpret and the ordinary logrank test can have reduced power.
A useful operational sequence is:
State \(\tau\), what counts as response, and how early death, progression, and missing assessments are handled.
Calculate response rates and compare each experimental arm with control using the prespecified threshold or rule.
Use the designated participants, time origin, event definition, censoring convention, and final analysis time.
Interpret the survival estimate and uncertainty using the inferential procedure planned for the complete selection design.
Threshold Sensitivity, Alternatives, and Failure Modes
A threshold is a design choice with consequences. In the example, if the response-difference threshold were 0.05 instead of 0.08, B's observed difference of 0.060 would also pass. The final stage would then include two experimental-versus-control comparisons, changing resource use and the multiplicity problem. The choice of threshold should be motivated clinically and evaluated by simulation or another valid operating-characteristic calculation under plausible relationships between response and survival. The earlier material on sensitivity analysis for the intermediate-final relationship addresses why that relationship cannot be assumed to be stable or perfectly predictive.
The fixed-difference rule can be compared with two common alternatives. A top-ranked rule always chooses the apparent leader, even if the difference from control is negligible; a minimum-difference rule can drop all experimental arms if none meets its threshold. A confidence-bound or test-based response rule accounts more explicitly for sampling variability, but it is not interchangeable with a point-estimate rule and must be specified and calibrated in advance. No response-based rule eliminates the need to assess overall survival.
For the survival stage, the logrank test is a conventional unadjusted comparison that uses event-time risk sets. A Cox model can additionally provide a hazard-ratio estimate and may adjust for prespecified covariates; under proportional hazards, that ratio has a relatively direct interpretation. Neither approach automatically fixes selection-induced error inflation. A valid final decision procedure must account for the treatment-selection process and any planned multiple comparisons, for example through a design-specific adjustment or calibration. The nominal p-value in the worked example is therefore illustrative, not a claim that a selected-arm test always controls the overall type I error at its nominal level.
Diagnostics should be considered before interpreting the final result. Check response assessment completeness and timing by arm, confirm that the selection rule was applied without post hoc changes, and examine event and censoring follow-up by arm. For survival, inspect whether risk sets become sparse, whether censoring patterns differ materially, and whether the proportional-hazards interpretation is plausible. A discrepancy between response and survival is not a computational error by itself: intermediate and final effects can disagree. Conversely, a favorable response difference cannot be used to declare a survival benefit when the survival analysis does not support that conclusion.
Worked Example: See How the Threshold Changes Selection
Use the same response data: \(\hat p_A=0.440\), \(\hat p_B=0.360\), and \(\hat p_C=0.300\). Compare a minimum difference of 0.08 with a minimum difference of 0.05.
With \(\delta=0.08\), A passes because \(0.440-0.300=0.140\geq0.08\), while B fails because \(0.360-0.300=0.060<0.08\). Only A continues. With \(\delta=0.05\), A still passes and B now passes because \(0.060\geq0.05\). Both experimental arms continue.
Thus a change of 0.03 in the threshold changes the number of final comparisons from one to two. It does not change the observed response data. The design team must evaluate the resulting follow-up, multiplicity, and power before choosing the rule; selecting the threshold after seeing these counts would make the analysis data-dependent in an unplanned way.
Check Your Understanding
Answer the questions, then compare your reasoning with the solutions below.
- In a hypothetical study, \(\hat p_A=0.41\), \(\hat p_B=0.34\), and \(\hat p_C=0.29\). Under a rule requiring an observed difference of at least 0.10, which arms continue?
- For a logrank comparison, \(O_A=42\), \(E_A=50\), and \(V_A=16\). Calculate \(Z\) and describe its direction.
- Why does passing an early response threshold not establish an overall-survival benefit?
- Name two assumptions or diagnostics relevant to interpreting the final survival comparison.
Solutions. (1) A's difference is \(0.41-0.29=0.12\), so A continues. B's is \(0.34-0.29=0.05\), so B does not. (2) \(Z=(42-50)/\sqrt{16}=-8/4=-2.000\); fewer deaths were observed in A than expected under the null, in the direction of a possible survival advantage. A final decision still needs the design's selection-aware procedure. (3) The response endpoint is used to select arms; response and survival effects can differ, and the final survival data must support a survival conclusion. (4) Examples include adequate and comparable event ascertainment, appropriate handling of censoring, sufficient event information, and plausibility of proportional hazards for a constant hazard-ratio interpretation.