Nonparametric Comparison of Survival Distributions

Abstract

Testing the equivalence of two groups is a familiar problem in statistics. Typically we are interested in testing a null hypothesis that two population means are equal versus an alternative that the means are not equal (for a two-sided test) or that the mean for an experimental treatment is greater than that for a standard treatment (one-sided test). We compute a test statistic from the observed data, and reject the null hypothesis if the test statistic exceeds a particular constant. The significance level of the test is the probability that we reject the null hypothesis when the null hypothesis is in fact true. A widely known test is the two-sample Students t-test for continuous observations, which requires the assumption that the observations are normally distributed. If the normal distribution assumption is in doubt, a rank-based test called the Mann-Whitney test may be used, which gives valid test results without making parametric assumptions. With survival data, if we are willing to assume that the data follow a particular parametric distribution, we can use likelihood theory to construct a test for equivalence of the two distributions, as we shall see in Chap. 10However, as we have discussed in the previous chapters, survival data from biomedical experiments or clinical trials generally doesn't lend itself to analysis by parametric methods. Thus, we shall construct nonparametric tests of equivalence of two survival functions, H0:{\thinspace}S1(t){\thinspace}={\thinspace}S0(t). Typically, S1 and S0 will represent the survival distributions for, respectively, an experimental and a control therapy. Now, a statistical hypothesis test (in the classical hypothesis testing framework) also requires us to specify an alternative hypothesis, and one might at first try to specify a one-sided alternative HA:{\thinspace}S1(t){\thinspace}>{\thinspace}S0(t) or two-sided alternative HA:{\thinspace}S1(t) ≠ S0(t). Unfortunately, things aren't so simple in survival analysis, since the alternative can take a wide range of forms. What if the survival distributions are similar for some values of t and differ for others? What if the survival distributions cross? How do we want our test statistic to behave under these different scenarios? One solution is to consider what is called a Lehman alternative, HA:S1(t)=S0(t)$ψ${\$}{\$}H{\_}{\{}A{\}}: S{\_}{\{}1{\}}(t) = {\backslash}left [S{\_}{\{}0{\}}(t){\backslash}right ]^{\{}{\backslash}psi {\}}{\$}{\$}. Equivalently, we can view Lehman alternatives in terms of proportional hazards as h1(t){\thinspace}={\thinspace}$ψ$ h0(t). Either way we would construct a one sided test as H0:{\thinspace}$ψ${\thinspace}={\thinspace}1 versus HA:{\thinspace}$ψ${\thinspace}