Joint modelling of cause-specific hazard functions with cubic splines: an application to a large series of breast cancer patients

The time of appearance of several kinds of relapses after a therapeutic intervention is of increasing interest in oncology. Typically, in breast cancer patients, events of clinical interest are intra-breast tumor recurrences and distant metastases, which act in a competitive way when considered as first failure. The evaluation of differential effects of clinical and biological variables on each event can improve the knowledge on the course of the disease and the targeting of future therapy. A simple tool for the joint smoothed estimation of cause-specific hazards functions and continuous covariate effects has been developed. Within the framework of generalized linear models with Poisson error, an extension of the piecewise exponential model is proposed, based on grouping follow-up times and continuous covariates. Interpolation of cause-specific hazards is obtained by resorting to cubic splines, which are piecewise polynomials of simple implementation with standard statistical software; their flexibility and smoothness are easily controlled by the number of knots and constraints on polynomial derivatives. The approach was applied to a data set of 2233 breast cancer patients treated with conservative surgery. It allowed modelling time-dependent and cause-specific effects of covariates on the hazard functions.

[1]  A. F. M. Smith,et al.  Automatic Bayesian curve ® tting , 1998 .

[2]  R. Gray.  Hazard Rate Regression Using Ordinary Nonparametric Regression Smoothers , 1996 .

[3]  Alan Agresti,et al.  Categorical Data Analysis , 1991, International Encyclopedia of Statistical Science.

[4]  Murray Aitkin,et al.  A Reanalysis of the Stanford Heart Transplant Data , 1983 .

[5]  N. L. Johnson,et al.  Survival Models and Data Analysis , 1982 .

[6]  R. Lenth,et al.  Spline Estimation of Paths Using Bearings-Only Tracking Data , 1996 .

[7]  Adrian F. M. Smith,et al.  Automatic Bayesian curve fitting , 1998 .

[8]  T. Holford.  Life tables with concomitant information. , 1976, Biometrics.

[9]  John Hinde,et al.  Statistical Modelling in GLIM. , 1990 .

[10]  Nan M. Laird,et al.  Covariance Analysis of Censored Survival Data Using Log-Linear Analysis Techniques , 1981 .

[11]  C. J. Stone,et al.  Hazard Regression , 2022 .

[12]  E Marubini,et al.  Local recurrences and distant metastases after conservative breast cancer treatments: partly independent events. , 1995, Journal of the National Cancer Institute.

[13]  R. Simon,et al.  Flexible regression models with cubic splines. , 1989, Statistics in medicine.

[14]  F E Harrell,et al.  The restricted cubic spline as baseline hazard in the proportional hazards model with step function time-dependent covariables. , 1995, Statistics in medicine.

[15]  M. Larson,et al.  Covariate analysis of competing-risks data with log-linear models. , 1984, Biometrics.

[16]  R. L. Eubank,et al.  Approximate regression models and splines , 1984 .

[17]  Daniel B. Mark,et al.  TUTORIAL IN BIOSTATISTICS MULTIVARIABLE PROGNOSTIC MODELS: ISSUES IN DEVELOPING MODELS, EVALUATING ASSUMPTIONS AND ADEQUACY, AND MEASURING AND REDUCING ERRORS , 1996 .

[18]  Jerome H. Friedman.  Multivariate adaptive regression splines (with discussion) , 1991 .

[19]  Calvin L. Williams,et al.  Modern Applied Statistics with S-Plus , 1997 .

[20]  J. Lindsey.  Relationships Among Sample Size, Model Selection and Likelihood Regions, and Scientifically Important Differences , 1999 .

[21]  H. A. David,et al.  Life Tests under Competing Causes of Failure and the Theory of Competing Risks , 1971 .

[22]  V T Farewell,et al.  The analysis of failure times in the presence of competing risks. , 1978, Biometrics.

[23]  Patricia L. Smith.  Splines as a Useful and Convenient Statistical Tool , 1979 .

[24]  K R Hess,et al.  Hazard function estimators: a simulation study. , 1999, Statistics in medicine.

[25]  M. Friedman.  Piecewise Exponential Models for Survival Data with Covariates , 1982 .

[26]  M Schumacher,et al.  Unbiased assessment of treatment effects on disease recurrence and survival in clinical trials. , 1983, Statistics in medicine.

[27]  Lin-An Chen,et al.  Multivariate regression splines , 1997 .

[28]  T. Holford.  The analysis of rates and of survivorship using log-linear models. , 1980, Biometrics.

[29]  Charles J. Stone,et al.  [Generalized Additive Models]: Comment , 1986 .

[30]  C. Kooperberg,et al.  Hazard regression with interval-censored data. , 1997, Biometrics.

[31]  Mary J. Lindstrom,et al.  Penalized Estimation of Free-Knot Splines , 1999 .

[32]  E Biganzoli,et al.  Modelling cause‐specific hazards with radial basis function artificial neural networks: application to 2233 breast cancer patients , 2001, Statistics in medicine.

[33]  Rosalba Miceli,et al.  Time distribution of the recurrence risk for breast cancer patients undergoing mastectomy: Further support about the concept of tumor dormancy , 2005, Breast Cancer Research and Treatment.

[34]  E Biganzoli,et al.  Feed forward neural networks for the analysis of censored survival data: a partial logistic regression approach. , 1998, Statistics in medicine.

[35]  R. Tibshirani.  Regression Shrinkage and Selection via the Lasso , 1996 .

[36]  M H Pejovic,et al.  A method of analysis taking into account competing events: application to the study of digestive complications following irradiation for cervical cancer. , 1987, Statistics in medicine.

[37]  H. Akaike,et al.  Information Theory and an Extension of the Maximum Likelihood Principle , 1973 .

[38]  J. Friedman.  Multivariate adaptive regression splines , 1990 .

[39]  James F. Watkins,et al.  Analysing Survival Data from Clinical Trials and Observational Studies. , 1995 .