Tutorials › Multi-Arm Multi-Stage (MAMS) Designs › Event-Driven Selection Rules

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

Event-Driven Selection Rules

Learn to define event-driven selection looks and calculate the treatment-contrast information available at each look.

Advanced 12 min read

What You'll Learn

  • Define a selection look by the accrued number of intermediate events.
  • Derive the Cox partial-likelihood information contributed by an event.
  • Distinguish raw risk-set proportions from model-implied event probabilities.
  • Calculate cumulative information and uncertainty for treatment-versus-control contrasts.
  • Diagnose when an event target provides little or unstable information.
  • Compare event-count looks with calendar-time and information-based looks.

Why Define a Look by Events?

In a time-to-event MAMS design, an intermediate event can provide information for deciding which experimental arms to retain. A selection look may therefore be triggered when a prespecified number of intermediate events has accrued, rather than at a prespecified calendar date. This makes the rule operationally clear, but the event count alone does not determine how precisely the arms can be compared: the information also depends on which participants remain at risk and on the treatment effects in the model.

Here the comparison is a Cox proportional-hazards comparison of each experimental arm with a shared control, using the intermediate time-to-event outcome. The selection rule can rank arms using their estimated intermediate log hazard ratios. It does not, by itself, establish that an arm improves the final endpoint or specify a calibrated threshold for retaining an arm.

Definition: Let \(D(t)\) be the total number of observed intermediate events across the control and experimental arms by time \(t\). A selection look with event target \(d_l\) occurs at \(\tau_l=\inf\{t:D(t)\ge d_l\}\). The information available at that look is the information in the observed risk sets and events through \(\tau_l\), not simply the event target \(d_l\).

Assume there are \(K\) experimental arms, indexed \(j=1,\ldots,K\), and a shared control indexed \(0\). Use treatment indicators \(x_j\), with the control represented by the zero vector and arm \(j\) represented by a vector whose \(j\)th entry is one. A Cox model specifies the hazard in arm \(j\) as the control hazard multiplied by \(\exp(\beta_j)\), where \(\beta_j\) is the log hazard ratio for the intermediate outcome. The principal assumptions here are proportional hazards over the analysis period, independent participants across individuals, non-informative censoring conditional on the model covariates, and correctly recorded event and risk-set times. For the derivation below, assume distinct event times; tied events require a stated tie method, such as the Breslow or Efron approximation.

Information at an Event

Immediately before event time \(t_i\), let \(Y_j(t_i)\) be the number at risk in arm \(j\). Conditional on one event occurring in that risk set, the Cox model assigns the event to an arm with probability proportional to the number at risk times its relative hazard. Thus the model-implied event probability is

$$ q_{ij}(\boldsymbol{\beta})= \frac{Y_j(t_i)\exp(\beta_j)} {Y_0(t_i)+\sum_{h=1}^{K}Y_h(t_i)\exp(\beta_h)},\qquad j=1,\ldots,K, $$

with \(\beta_0=0\) for control. These \(q_{ij}\) are not generally the raw proportions at risk. The raw proportions \(Y_j(t_i)/\sum_hY_h(t_i)\) equal the event probabilities under the null model \(\boldsymbol{\beta}=0\), but away from the null the relative hazards weight the risk-set counts.

Let \(\mathbf{x}_i\) denote the treatment-indicator vector for the individual who experiences event \(i\). The event's partial log-likelihood contribution, up to terms not involving \(\boldsymbol{\beta}\), is the observed treatment indicator dotted with \(\boldsymbol{\beta}\), minus the log of the weighted risk-set total. Differentiating gives the score and its information contribution:

$$ \ell_i(\boldsymbol{\beta}) =\boldsymbol{\beta}^{T}\mathbf{x}_i- \log\left\{Y_0(t_i)+\sum_{j=1}^{K}Y_j(t_i)\exp(\beta_j)\right\}, \qquad \mathcal{I}_i(\boldsymbol{\beta}) =\operatorname{diag}(\mathbf{q}_{i})-\mathbf{q}_{i}\mathbf{q}_{i}^{T}. $$

Here \(\mathbf{q}_i=(q_{i1},\ldots,q_{iK})^T\) contains the model-implied probabilities for experimental arms; the control probability is omitted from the vector because its treatment-indicator vector is zero. The matrix is the covariance matrix of the treatment-indicator vector under the model-implied event distribution in that risk set. The observed event identity contributes to the score, while the curvature of the partial likelihood is determined by the weighted risk set.

From Event Count to Contrast Precision

At a look, the observed information is the sum of the event contributions over all events observed by that look. In practice, it can be evaluated at the fitted Cox coefficients, or at a prespecified value for design calculations. The estimated covariance matrix of the log hazard-ratio estimates is approximately the inverse information matrix when the model is identifiable and the large-sample approximation is adequate. Shared control appears in the off-diagonal entries: estimates for two treatment-versus-control contrasts are generally correlated.

$$ \mathcal{I}(\boldsymbol{\beta};\tau_l) =\sum_{i:t_i\le \tau_l} \left\{\operatorname{diag}(\mathbf{q}_{i})-\mathbf{q}_{i}\mathbf{q}_{i}^{T}\right\}, \qquad \widehat{\operatorname{Var}}(\widehat{\boldsymbol{\beta}}) \approx \mathcal{I}(\widehat{\boldsymbol{\beta}};\tau_l)^{-1}. $$

The expression is a local information calculation: for inference it is ordinarily evaluated at the fitted model, while a planning approximation may evaluate it under the null or another specified effect configuration. Under the null, the event probabilities reduce to raw risk-set proportions. If risk-set composition is approximately stable, with null proportions \(p_0,\ldots,p_K\), then \(D\) events give the approximation \(D\{\operatorname{diag}(p_1,\ldots,p_K)-\mathbf{p}\mathbf{p}^T\}\). It is an approximation, not a universal conversion from event count to information.

Worked Example: Correctly Weighting a Risk Set

Suppose a hypothetical risk set just before an intermediate event contains 40 control participants, 30 participants in arm A, and 30 in arm B. For illustration, evaluate the Cox information at hazard ratios 2.0 for A and 1.5 for B. The corresponding log hazard ratios are \(\beta_A=\log(2)\) and \(\beta_B=\log(1.5)\).

The relative-hazard-weighted risk-set totals are 40 for control, \(30(2)=60\) for A, and \(30(1.5)=45\) for B, for a total of 145. Therefore \(q_A=60/145=0.4138\) and \(q_B=45/145=0.3103\). The control event probability is \(40/145=0.2759\).

The treatment-block information for this one event is \(\operatorname{diag}(q_A,q_B)-(q_A,q_B)^T(q_A,q_B)\), giving \[ \mathcal{I}_i= \begin{pmatrix} 0.4138(1-0.4138) & -0.4138(0.3103)\\ -0.4138(0.3103) & 0.3103(1-0.3103) \end{pmatrix} \approx \begin{pmatrix} 0.2426 & -0.1284\\ -0.1284 & 0.2140 \end{pmatrix}. \] The negative off-diagonal term reflects the competition between treatment indicators within a risk set. Using the raw proportions \((0.3,0.3)\) here would instead give diagonal entries 0.2100 and off-diagonal entries \(-0.0900\). That is the null-weighted information, not the information at the stated nonzero hazard ratios.

Worked Information Calculations

Worked Example: A Balanced Null Approximation

Consider two experimental arms and one shared control. Suppose a hypothetical selection look targets 120 events and, throughout the observed event risk sets, the three arms have equal numbers at risk. For a planning approximation under \(H_0:\beta_A=\beta_B=0\), each arm has event probability \(1/3\).

Each event contributes diagonal entries \(1/3-(1/3)^2=2/9\) and off-diagonal entries \(-(1/3)(1/3)=-1/9\). Thus the total information is \[ \mathcal{I}=120 \begin{pmatrix}2/9&-1/9\\-1/9&2/9\end{pmatrix} = \begin{pmatrix}26.67&-13.33\\-13.33&26.67\end{pmatrix}. \] Its inverse is \[ \mathcal{I}^{-1}= \frac{3}{120} \begin{pmatrix}2&1\\1&2\end{pmatrix} = \begin{pmatrix}0.0500&0.0250\\0.0250&0.0500\end{pmatrix}. \] The approximate standard error for each log hazard-ratio estimate is \(\sqrt{0.0500}=0.2236\), and the estimated covariance between the two log hazard-ratio estimates is 0.0250.

Suppose the fitted estimates at this hypothetical look are \(\widehat{\beta}_A=-0.35\) and \(\widehat{\beta}_B=-0.10\). Their hazard-ratio estimates are \(\exp(-0.35)=0.7047\) and \(\exp(-0.10)=0.9048\). If the prespecified selection rule ranks arms by lower estimated hazard ratio, A ranks ahead of B. The standardized estimates are \(-0.35/0.2236=-1.565\) and \(-0.10/0.2236=-0.447\); neither calculation alone supplies a calibrated retention threshold.

Worked Example: Changing Risk Sets at a Fixed Event Target

Suppose 50 hypothetical events have accrued by a look. For 20 events, the risk-set proportions are control 0.50, A 0.25, and B 0.25; for the other 30 events they are control 0.20, A 0.40, and B 0.40. Calculate a null-model information approximation, so model-implied event probabilities equal those proportions.

For the first group of events, the treatment-block information per event has diagonal entries \(0.25(0.75)=0.1875\) and off-diagonal entries \(-0.25(0.25)=-0.0625\). Its 20-event contribution is \[ 20\begin{pmatrix}0.1875&-0.0625\\-0.0625&0.1875\end{pmatrix} = \begin{pmatrix}3.75&-1.25\\-1.25&3.75\end{pmatrix}. \] For the second group, the corresponding entries are \(0.40(0.60)=0.2400\) and \(-0.40(0.40)=-0.1600\); its contribution is \[ 30\begin{pmatrix}0.2400&-0.1600\\-0.1600&0.2400\end{pmatrix} = \begin{pmatrix}7.20&-4.80\\-4.80&7.20\end{pmatrix}. \] Adding the contributions gives \(\mathcal{I}=\begin{pmatrix}10.95&-6.05\\-6.05&10.95\end{pmatrix}\). The determinant is \(10.95^2-6.05^2=83.30\), so the diagonal of the inverse is \(10.95/83.30=0.1315\) and the off-diagonal is \(6.05/83.30=0.07263\). Each log hazard-ratio estimate consequently has approximate standard error \(\sqrt{0.1315}=0.3626\). The same total of 50 events as in another design could provide different precision if the risk-set composition differed.

What the Event-Driven Rule Does—and Does Not—Guarantee

The information formula justifies why an event-driven rule can be useful: event occurrence contributes directly to the partial likelihood for the hazard contrasts. An event target also gives a reproducible trigger that does not depend on a fixed calendar date. However, an event count does not guarantee a particular amount of information. If one treatment arm has few people at risk, its event-indicator variation may be small; if an arm has no participants at risk, its contrast may not be estimable from that period. The shared-control covariance also means that the precision of the vector of comparisons is not captured by separate standard errors alone.

A practical selection rule should specify the event-count trigger and the comparison statistic in advance. For example, it may rank the fitted log hazard ratios, with more negative values favoring an arm when fewer intermediate events are desirable. The event-driven trigger determines when the statistic is evaluated; it does not determine the ranking direction, a threshold for dropping an arm, or the error properties of the overall selection procedure.

The information calculation follows from the curvature of the Cox partial likelihood. At an event, the negative second derivative is the covariance of the treatment indicators under the weighted risk-set distribution. Summing those curvatures gives the total information, which explains why changing risk-set composition changes information even when the event count is held fixed.

Justification: The partial-likelihood score contribution is \(\mathbf{x}_i-\bar{\mathbf{x}}_i(\boldsymbol{\beta})\), where \(\bar{\mathbf{x}}_i\) is the risk-set mean under weights proportional to \(Y_j(t_i)\exp(\beta_j)\). Differentiating that weighted mean with respect to \(\boldsymbol{\beta}\) gives the weighted covariance of the treatment indicators. Therefore the negative derivative of the score is \(\operatorname{diag}(\mathbf{q}_i)-\mathbf{q}_i\mathbf{q}_i^T\); summing over events yields the stated information matrix.

Diagnostics, Alternatives, and Limitations

At each look, inspect the risk-set sizes by arm, the fitted coefficients, and the eigenvalues or condition of the information matrix. Warning signs include an arm with very few participants at risk, a nearly singular information matrix, estimates with very large standard errors, or substantial changes in coefficients as events accrue. These features indicate that the nominal event target may not have delivered stable comparisons. The model also requires scrutiny: marked departures from proportional hazards or outcome-dependent censoring can undermine the Cox-based comparison, regardless of how many events have accrued.

Three timing choices answer related but distinct operational needs. A calendar-time look is scheduled at a date and may contain widely varying numbers of events. An event-driven look waits for an event target, making the endpoint count explicit while leaving information dependent on the risk sets. An information-driven look targets a prespecified information level, which more directly controls the intended precision but requires monitoring or forecasting information rather than only counting events. None automatically solves model misspecification or selection-threshold calibration.

Tied event times require a consistent partial-likelihood tie method; the distinct-event derivation should not be applied as if tied events were independent sequential failures without qualification. In addition, the inverse-information variance is an asymptotic approximation. Small event counts, sparse arms, heavy censoring, or unstable parameter estimates can make it unreliable. Finally, an event-driven selection rule concerns the intermediate outcome used for the arm decision. Its association with a final endpoint is a separate scientific question.

Common mistake: Treating the observed risk-set proportions as the event probabilities under every treatment effect configuration gives the wrong information away from the null. Use \(q_{ij}\propto Y_j(t_i)\exp(\beta_j)\) and normalize across arms; the raw proportions are appropriate for the null-model calculation only.

Check Your Understanding

Answer each question, then compare with the brief solutions.

  1. In a risk set with 20 control participants and 10 in treatment A, what is the model-implied probability that the next event is in A if the hazard ratio is 2?
  2. For two experimental arms under the null, why can the information matrix have a nonzero off-diagonal entry?
  3. Two selection looks each contain 80 events. Give one reason their information matrices could differ.
  4. What does an event-count trigger specify, and what does it leave unspecified?

Solutions. (1) The weighted totals are 20 for control and \(10(2)=20\) for A, so the probability is \(20/(20+20)=0.5000\). (2) Both treatment indicators are defined relative to the same risk set and shared control; their covariance is generally nonzero, equal to \(-q_Aq_B\) for one event. (3) Their risk-set compositions may differ, or the model-implied event probabilities may differ because fitted treatment effects differ. (4) It specifies when to evaluate the selection statistic, namely once the target number of events is observed; by itself it does not specify the information achieved, the ranking direction, or a calibrated retention threshold.

Clinical Biostats