Analysing multicentre competing risks data with a mixed proportional hazards model for the subdistribution

In the competing-risks setting, to test the effect of a covariate on the probability of one particular cause of failure, the Fine and Gray model for the subdistribution hazard can be used. However, sometimes, competing risks data cannot be considered as independent because of a clustered design, for instance in registry cohorts or multicentre clinical trials. Frailty models have been shown useful to analyse such clustered data in a classical survival setting, where only one risk acts on the population. Inclusion of random effects in the subdistribution hazard has not been assessed yet. In this work, we propose a frailty model for the subdistribution hazard. This allows first to assess the heterogeneity across clusters, then to incorporate such an effect when testing the effect of a covariate of interest. Based on simulation study, the effect of the presence of heterogeneity on testing for covariate effects was studied. Finally, the model was illustrated on a data set from a registry cohort of patients with acute myeloid leukaemia who underwent bone marrow transplantation.

[1]  J. Kalbfleisch,et al.  The Statistical Analysis of Failure Time Data , 1980 .

[2]  J Crowley,et al.  Estimation of failure probabilities in the presence of competing risks: new representations of old estimators. , 1999, Statistics in medicine.

[3]  O. Aalen,et al.  Heterogeneity in survival analysis. , 1988, Statistics in medicine.

[4]  Robert Gray,et al.  Flexible Methods for Analyzing Survival Data Using Splines, with Applications to Breast Cancer Prognosis , 1992 .

[5]  Aeilko H. Zwinderman,et al.  Centre‐effect on Survival after Bone Marrow Transplantation: Application of Time‐dependent Frailty Models , 2004 .

[6]  R Brookmeyer,et al.  Modelling paired survival data with covariates. , 1989, Biometrics.

[7]  Val Gebski,et al.  On the issue of ‘multiple’ first failures in competing risks analysis , 2002, Statistics in medicine.

[8]  Philip Hougaard,et al.  Analysis of Multivariate Survival Data , 2001 .

[9]  M S Pepe,et al.  Kaplan-Meier, marginal or conditional probability curves in summarizing competing risks failure time data? , 1993, Statistics in medicine.

[10]  J P Klein,et al.  Testing for centre effects in multi-centre survival studies: a Monte Carlo comparison of fixed and random effects tests. , 1999, Statistics in medicine.

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

[12]  Paula R Williamson,et al.  Investigating heterogeneity in an individual patient data meta‐analysis of time to event outcomes , 2005, Statistics in medicine.

[13]  David V Glidden,et al.  Modelling clustered survival data from multicentre clinical trials , 2004, Statistics in medicine.

[14]  Robert Gray,et al.  A Proportional Hazards Model for the Subdistribution of a Competing Risk , 1999 .

[15]  R. Wolfe,et al.  A Frailty Model for Informative Censoring , 2002, Biometrics.

[16]  McGilchrist Ca.  REML estimation for survival models with frailty. , 1993 .

[17]  R. Gray.  A Class of $K$-Sample Tests for Comparing the Cumulative Incidence of a Competing Risk , 1988 .

[18]  T. Yamaguchi,et al.  Investigating centre effects in a multi-centre clinical trial of superficial bladder cancer. , 1999, Statistics in medicine.

[19]  Clelia Di Serio,et al.  The Protective Impact of a Covariate on Competing Failures with an Example from a Bone Marrow Transplantation Study , 1997 .

[20]  J. Sierra,et al.  Effect of centre on outcome of bone-marrow transplantation for acute myeloid leukaemia , 2000, The Lancet.

[21]  P. Grambsch,et al.  Penalized Survival Models and Frailty , 2003 .

[22]  Gerard J. van den Berg,et al.  The identifiability of the mixed proportional hazards competing risks model , 2003 .

[23]  Jianwen Cai,et al.  Some Graphical Displays and Marginal Regression Analyses for Recurrent Failure Times and Time Dependent Covariates , 1993 .