Tutorials › Group Sequential Methods and Sample Size Re-Estimation › Robustness to Misspecified Delay and Weight Functions

Non-Proportional Hazards and Delayed Effects · Tutorial 563 of 699

Robustness to Misspecified Delay and Weight Functions

Compare the operating characteristics of prespecified weighted logrank designs across plausible true delay and hazard patterns, including how to assess sequential stopping under misspecification.

Advanced 12 min read

What You'll Learn

  • Express the expected weighted logrank signal and its variance under a specified hazard pattern.
  • Compare power when the true treatment delay differs from the delay assumed in weight selection.
  • Explain how a misspecified weight can affect interim stopping and final power.
  • Distinguish power robustness from Type I error control under the null.
  • Use scenario-based sensitivity analysis to identify when a design is fragile.

Why Misspecification Matters

A weighted logrank design is often chosen because a particular pattern of treatment benefit is considered plausible. For example, a treatment may be expected to have little effect initially and a later benefit. A weight that downweights early events and emphasizes later events can be more sensitive to that pattern than the ordinary logrank test. But the assumed delay and the chosen weight are design assumptions, not guarantees about the true hazard pattern.

If the benefit begins earlier or later than expected, the weighted statistic may give relatively more influence to events carrying little treatment signal and less influence to events carrying more signal. This changes power and can change the chance of crossing an interim efficacy boundary. It does not, by itself, invalidate a test whose weights and boundaries were prespecified and calibrated under the null. The key is to examine operating characteristics under more than one plausible alternative.

Definition: A weighted logrank design is robust to misspecification when its relevant operating characteristics—such as power, probability of early stopping, or expected trial duration—remain acceptable across a prespecified range of plausible true hazard patterns that differ from the pattern used to choose the weights. Robustness is always relative to a stated range of scenarios and stated criteria.

How the True Hazard Pattern Enters the Statistic

Let \(t\) denote event time, and let \(w(t)\) be a fixed, prespecified weight function. For a two-arm comparison, the weighted logrank score accumulates weighted differences between observed and expected treatment-arm events. With a sign convention in which a treatment benefit contributes positively, its null variance is estimated by the corresponding weighted event-information sum. Write the cumulative score and its variance as \(U_w(t)\) and \(V_w(t)\), and the standardized statistic as \(Z_w(t)=U_w(t)/\sqrt{V_w(t)}\).

For interpretation under an alternative, let \(q(t)=-\log\{HR(t)\}\), where \(HR(t)\) is the treatment-to-control hazard ratio. Thus \(q(t)>0\) represents a beneficial treatment effect. Under the usual large-sample, local-alternative approximation, the expected weighted score is approximately the weight- and information-weighted treatment signal. In continuous time, this gives the following useful expression for the standardized mean:

$$ \Delta_w(\tau) \approx \frac{\int_0^\tau w(t)q(t)\,dI(t)} {\sqrt{\int_0^\tau w(t)^2\,dI(t)}} $$

Here \(dI(t)\) denotes the incremental logrank information at time \(t\), and \(\tau\) is the analysis time. This expression is an approximation, not a replacement for the trial’s full design calculations: allocation, censoring, event-time distributions, and the risk sets all affect the information. It nevertheless shows the core robustness issue. The numerator rewards weight placed where the true treatment signal is present, whereas the denominator reflects the variance created by weighting events.

A Short Derivation and the Role of Sequential Monitoring

Divide follow-up into \(J\) time intervals and suppose, for illustration, that the unweighted score increments are approximately independent with variances \(I_j\) and means \(q_j I_j\). Applying a constant weight \(w_j\) within interval \(j\) gives a total weighted score with approximate mean \(\sum_j w_jq_jI_j\) and variance \(\sum_jw_j^2I_j\). Standardizing yields:

$$ E(Z_w) \approx \frac{\sum_{j=1}^{J}w_jq_jI_j} {\sqrt{\sum_{j=1}^{J}w_j^2I_j}}. $$

This derivation makes a useful point about a delayed effect. If \(q_j=0\) early and positive later, a weight that is small early and larger later can increase the signal relative to its standard deviation. If the effect instead begins early, that same weight may discard useful early signal. If the effect begins later than anticipated, weight concentrated in the originally expected benefit period may again be poorly aligned.

In a group sequential design, let looks occur at cumulative information values \(V_1,\ldots,V_K\). Under the standard canonical joint-normal approximation under the null, the cumulative standardized statistics have correlation \(\operatorname{Corr}(Z_i,Z_j)\approx\sqrt{V_i/V_j}\) for \(i\leq j\). Under an alternative, their means depend on the true hazard pattern and the prespecified weights. Therefore, an overall power calculation must account for the joint distribution and the stopping boundaries; final-look power alone does not determine the probability of stopping early or the total probability of rejecting the null.

Worked Example: A Weight That Helps When the Delay Is Correct

Consider a hypothetical two-period design with 100 units of unweighted information in each period. Suppose there is no treatment effect in period 1 and the log hazard reduction is \(q_2=0.30\) in period 2. Compare equal weights \(w=(1,1)\) with a delayed-effect weight \(w=(0.25,1)\). Use a prespecified one-sided upper-tail test with critical value 1.96 for this illustrative final-look comparison.

For equal weights, the approximate mean standardized statistic is \(0.30(100)/\sqrt{100+100}=30/\sqrt{200}=2.121\). The approximate power is \(\Phi(2.121-1.96)=\Phi(0.161)\approx0.564\).

For the delayed-effect weight, the period 1 signal is zero and the period 2 signal is \(0.30(100)=30\). Its variance is \(0.25^2(100)+1^2(100)=106.25\), so \(E(Z_w)\approx30/\sqrt{106.25}=2.910\). The approximate power is \(\Phi(2.910-1.96)=\Phi(0.950)\approx0.829\).

Under the assumed delayed pattern, the delayed-effect weight improves approximate power in this simplified calculation. This conclusion is conditional on the stated hazard pattern and information approximation; it does not establish that this weight will be best if the true effect starts at a different time.

Worked Sensitivity Analyses

A practical robustness analysis keeps the trial design fixed—weights, information targets, looks, and boundaries—and varies the data-generating hazard pattern. Scenarios should be chosen before examining unblinded interim treatment comparisons. Useful changes include moving the onset of benefit earlier or later, changing the size of the post-onset effect, or allowing a gradual rather than abrupt onset. For each scenario, calculate the operating characteristics using the same design and analysis rule.

Worked Example: Later Benefit Under an Unweighted Analysis

Suppose, hypothetically, that follow-up is divided into three periods with 100 information units in each. The treatment has no effect in the first two periods and \(q_3=0.20\) in the last. For an unweighted final-look statistic, the expected score is \(0+0+0.20(100)=20\), and its variance is \(100+100+100=300\). Thus \(E(Z)=20/\sqrt{300}=1.155\).

For a one-sided upper-tail test with critical value 1.96, approximate power is \(\Phi(1.155-1.96)=\Phi(-0.805)\approx0.2103\). The calculation illustrates how a treatment benefit that occurs only late can yield low power when earlier event information contributes variance but no expected treatment signal. The value 0.2103 is the approximate probability of crossing this single final-look threshold under the stated alternative; it is not the probability that the null hypothesis is true.

Worked Example: Misspecification and an Interim Boundary

Consider a hypothetical two-look design at information values \(V_1=100\) and \(V_2=200\), with a prespecified delayed-effect weight. Suppose the design’s null calibration has established upper boundaries \(b_1=2.80\) and \(b_2=2.00\) for its chosen overall Type I error target. Now consider a true effect that begins only after look 1. At that interim look the treatment signal is zero, so the alternative mean is \(\Delta_1=0\).

The probability of crossing the first boundary under this scenario is approximately \(P(Z_1>2.80)=1-\Phi(2.80)\approx0.0026\). This small probability follows because, before the effect begins, the interim statistic has approximately its null distribution. It is not an estimate of the design’s overall power.

To obtain overall power, one must combine first-look crossing with final-look crossing among trials that did not stop at look 1. Under the canonical approximation, the correlation is \(\rho=\sqrt{100/200}=0.7071\). If the alternative means at the two looks are \(\Delta_1\) and \(\Delta_2\), the overall crossing probability is

$$ P(Z_1>b_1)+ \int_{-\infty}^{b_1} \left[ 1-\Phi\left( \frac{b_2-\Delta_2-\rho(z-\Delta_1)} {\sqrt{1-\rho^2}} \right) \right] \phi(z-\Delta_1)\,dz. $$

The integral accounts for the conditional distribution of the second-look statistic among trials that have not crossed the first boundary. It should be evaluated for each true hazard scenario using its corresponding \(\Delta_1,\Delta_2\). Reporting only the marginal final-look probability \(P(Z_2>b_2)\) would ignore the interim stopping rule and would not give the group sequential power.

Interpreting Type I Error, Diagnostics, and Alternatives

A crucial distinction is between misspecification under the alternative and misspecification under the null. If a weight function is fixed in advance and the group sequential boundaries are correctly calibrated for that statistic, a different-than-expected treatment effect generally changes power and stopping behavior, not the null Type I error. In contrast, selecting or changing weights using unblinded interim treatment contrasts can change the null distribution of the final procedure. Such adaptation requires an analysis and calibration that account for the selection rule; one cannot assume the original boundaries still control the error rate.

Review the sensitivity analysis for warning signs such as a steep drop in power for a small shift in onset, most of the expected power coming from one narrow scenario, or substantial changes in early stopping probabilities across plausible patterns. Also check that the modeled event information and risk-set behavior are credible for each scenario. A calculation based on identical information accrual across scenarios may be misleading if changing the hazard pattern would materially change event timing or censoring.

The main alternatives are not interchangeable. An ordinary logrank analysis uses constant weights and can be less sensitive to a particular delayed pattern, but it may be less vulnerable to choosing the wrong time emphasis. A prespecified weighted logrank test may gain power for its target pattern, at the cost of potentially losing power when the pattern differs. A design that combines evidence across multiple prespecified weights can reduce dependence on a single weight, but its statistic and sequential boundaries must be calibrated as a combined procedure. Any alternative should be compared using the same plausible scenarios, error target, and stopping rules.

Robustness is not a guarantee against every possible hazard pattern. Scenario analyses depend on assumptions about treatment effects, event accrual, censoring, and the information process. They also do not turn a weighted logrank result into a universal summary of treatment benefit: the statistic answers a design-specific weighted comparison. The appropriate conclusion is therefore conditional—how the planned procedure behaves over the scenarios considered—not that it is optimal for every possible true effect.

1
Specify the design.
Record the weight function, information targets, interim boundaries, and analysis direction before comparing scenarios.
2
Vary the true pattern.
Use plausible earlier, later, weaker, or gradual treatment effects while holding the analysis rule fixed.
3
Recalculate joint operating characteristics.
Assess overall power and interim stopping probabilities with the sequential dependence and boundaries included.
4
Judge sensitivity against design goals.
Decide whether performance across the scenario range is acceptable, and document where the design is fragile.

Common Mistake

A common mistake is to conclude that a prespecified weight is “robust” because the trial controls Type I error under the null. Null calibration and power robustness are different properties. A procedure may control its false-positive rate at the target level and still have poor power, or very different early-stopping probabilities, when the true treatment effect begins at a time other than anticipated. Conversely, poor power under one alternative does not by itself show that null Type I error is inflated.

Key takeaway: Evaluate a weighted group sequential design under multiple plausible true hazard patterns. Keep the prespecified weights and boundaries fixed in that comparison, and report both power and sequential stopping behavior; null error control alone does not establish robustness.

Check Your Understanding

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

  1. If a weight emphasizes late events but the true benefit begins early, which part of the approximate standardized mean may be reduced relative to a well-aligned weight?
  2. In the three-period example, why does the unweighted statistic have variance 300 even though the treatment signal occurs only in the last period?
  3. Does a low interim crossing probability before a delayed effect begins determine overall group sequential power? Why or why not?
  4. What must be reconsidered if weights are selected using unblinded interim treatment comparisons?

Solutions. (1) The numerator, \(\sum_j w_jq_jI_j\), may be reduced because relatively little weight is assigned to periods with positive signal; the variance denominator also depends on the weights. (2) Under the stated approximation, information contributes to the variance in all three periods, including periods with no treatment signal. (3) No. Overall power includes crossing at any look and depends on the joint distribution, boundaries, and the chance of reaching later looks. (4) The null distribution and Type I error calibration may change; the selection rule must be incorporated into a valid analysis and boundary calibration.

Clinical Biostats