Confidence Bands for the Difference of Two Survival Curves Under Proportional Hazards Model

A common approach to testing for differences between the survival rates of two therapies is to use a proportional hazards regression model which allows for an adjustment of the two survival functions for any imbalance in prognostic factors in the comparison. When the relative risk of one treatment to the other is not constant over time the question of which therapy has a survival advantage is difficult to determine from the Cox model. An alternative approach to this problem is to plot the difference between the two predicted survival functions with a confidence band that provides information about when these two treatments differ. Such a band will depend on the covariate values of a given patient. In this paper we show how to construct a confidence band for the difference of two survival functions based on the proportional hazards model. A simulation approach is used to generate the bands. This approach is used to compare the survival probabilities of chemotherapy and allogeneic bone marrow transplants for chronic leukemia.

[1]  D. Harrington,et al.  Counting Processes and Survival Analysis , 1991 .

[2]  Simultaneous M-estimator of the common location and the scale-ratio in the two-sample problem , 1982 .

[3]  H. Heimpel,et al.  Survival with bone marrow transplantation versus hydroxyurea or interferon for chronic myelogenous leukemia. The German CML Study Group. , 1998, Blood.

[4]  Laurence L. George,et al.  The Statistical Analysis of Failure Time Data , 2003, Technometrics.

[5]  David R. Cox,et al.  Regression models and life tables (with discussion , 1972 .

[6]  O. Aalen.  A linear regression model for the analysis of life times. , 1989, Statistics in medicine.

[7]  Lee-Jen Wei,et al.  Confidence bands for survival curves under the proportional , 1994 .

[8]  D. Dabrowska,et al.  Methods for comparing cumulative hazard functions in a semi-proportional hazard model. , 1992, Statistics in medicine.

[9]  R. Gill,et al.  Cox's regression model for counting processes: a large sample study : (preprint) , 1982 .

[10]  L. J. Wei,et al.  Testing Goodness of Fit for Proportional Hazards Model with Censored Observations , 1984 .

[11]  Danyu Lin,et al.  Goodness-of-Fit Analysis for the Cox Regression Model Based on a Class of Parameter Estimators , 1991 .

[12]  Michael Parzen,et al.  Simultaneous Confidence Intervals for the Difference of Two Survival Functions , 1997 .

[13]  D. Cox.  Regression Models and Life-Tables , 1972 .

[14]  Mei-Jie Zhang,et al.  Statistical Challenges in Comparing Chemotherapy and Bone Marrow Transplantation as a Treatment for Leukemia , 1996 .

[15]  Gilbert MacKenzie,et al.  The Statistical Analysis of Failure Time Data , 1982 .

[16]  Elja Arjas,et al.  A Graphical Method for Assessing Goodness of Fit in Cox's Proportional Hazards Model , 1988 .

[17]  Kjell A. Doksum,et al.  Graphical comparison of cumulative hazards for two populations , 1989 .

[18]  Zhiliang Ying,et al.  Estimating a distribution function with truncated and censored data , 1991 .

[19]  M. Woodroofe.  Estimating a Distribution Function with Truncated Data , 1985 .

[20]  E. Kaplan,et al.  Nonparametric Estimation from Incomplete Observations , 1958 .

[21]  N. Breslow,et al.  Analysis of Survival Data under the Proportional Hazards Model , 1975 .

[22]  Per Kragh Andersen,et al.  Testing Goodness of Fit of Cox's Regression and Life Model , 1982 .

[23]  J. Klein,et al.  Statistical Models Based On Counting Process , 1994 .