Comparing Kaplan–Meier curves with delayed treatment effects: applications in immunotherapy trials

Deconstructing the Kaplan-Meier curve: Quantification of treatment effect using the treatment effect process

In studies of survival and its association with treatment and other prognostic variables, elapsed time alone will often show itself to be among the strongest, if not the strongest, of the predictor variables. Kaplan-Meier curves will show the overall survival of each group and the general differences between groups due to treatment. However, the time-dependent nature of treatment effects is not always immediately transparent from these curves. More sophisticated tools are needed to spotlight the treatment effects. An important tool in this context is the treatment effect process. This tool can be potent in revealing the complex myriad of ways in which treatment can affect survival time. We look at a recently published study in which the outcome was relapse-free survival, and we illustrate how the use of the treatment effect process can provide a much deeper understanding of the relationship between time and treatment in this trial.

Testing for Differences in Survival When Treatment Effects Are Persistent, Decaying, or Delayed

A statistical test for the presence of treatment effects on survival will be based on a null hypothesis (absence of effects) and an alternative (presence of effects). The null is very simply expressed. The most common alternative, also simply expressed, is that of proportional hazards. For this situation, not only do we have a very powerful test in the log-rank test but also the outcome is readily interpreted. However, many modern treatments fall outside this relatively straightforward paradigm and, as such, have attracted attention from statisticians eager to do their best to avoid losing power as well as to maintain interpretability when the alternative hypothesis is less simple. Examples include trials where the treatment effect decays with time, immunotherapy trials where treatment effects may be slow to manifest themselves as well as the so-called crossing hazards problem. We review some of the solutions that have been proposed to deal with these issues. We pay particular attention to the integrated log-rank test and how it can be combined with the log-rank test itself to obtain powerful tests for these more complex situations.

Explained variation under the additive hazards model

We study explained variation under the additive hazards regression model for right-censored data. We consider different approaches for developing such a measure, and focus on one that estimates the proportion of variation in the failure time explained by the covariates. We study the properties of the measure both analytically, and through extensive simulations. We apply the measure to a well-known survival dataset as well as the linked surveillance, epidemiology, and end results-Medicare database for prediction of mortality in early stage prostate cancer patients using high-dimensional claims codes.