Parametric accelerated failure time models with random effects and an application to kidney transplant survival

Accelerated failure time models with a shared random component are described, and are used to evaluate the effect of explanatory factors and different transplant centres on survival times following kidney transplantation. Different combinations of the distribution of the random effects and baseline hazard function are considered and the fit of such models to the transplant data is critically assessed. A mixture model that combines short- and long-term components of a hazard function is then developed, which provides a more flexible model for the hazard function. The model can incorporate different explanatory variables and random effects in each component. The model is straightforward to fit using standard statistical software, and is shown to be a good fit to the transplant data.

[1]  P. Hougaard,et al.  Fundamentals of Survival Data , 1999, Biometrics.

[2]  M. Cowles.  Modelling Survival Data in Medical Research (2nd ed.) (Book) , 2004 .

[3]  S. Self,et al.  Semiparametric Likelihood Estimation in the Clayton–Oakes Failure Time Model , 1999 .

[4]  O. Aalen.  Effects of frailty in survival analysis , 1994, Statistical methods in medical research.

[5]  P. Royston,et al.  Flexible parametric proportional‐hazards and proportional‐odds models for censored survival data, with application to prognostic modelling and estimation of treatment effects , 2002, Statistics in medicine.

[6]  J. Copas,et al.  Estimating the Risk of Reoffending by using Exponential Mixture Models , 1997 .

[7]  K Borch-Johnsen,et al.  Heterogeneity models of disease susceptibility, with application to diabetic nephropathy. , 1994, Biometrics.

[8]  Richard D. Gill,et al.  A counting process approach to maximum likelihood estimation in frailty models , 1992 .

[9]  T. Louis,et al.  Survival analysis using a scale change random effects model , 1995 .

[10]  D J Sargent,et al.  A general framework for random effects survival analysis in the Cox proportional hazards setting. , 1998, Biometrics.

[11]  Joel L. Horowitz,et al.  Semiparametric Estimation of a Proportional Hazard Model with Unobserved Heterogeneity , 1999 .

[12]  J. Klein,et al.  Modeling Random Effects for Censored Data by a Multivariate Normal Regression Model , 1999, Biometrics.

[13]  N. Keiding,et al.  The role of frailty models and accelerated failure time models in describing heterogeneity due to omitted covariates. , 1997, Statistics in medicine.

[14]  S. Gore,et al.  RENAL TRANSPLANT REJECTION: TRANSIENT IMMUNODOMINANCE OF HLA MISMATCHES , 1990, Transplantation.

[15]  S M Gore,et al.  High energy neutron treatment for pelvic cancers: study stopped because of increased mortality. , 1991, BMJ.

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

[17]  R Crouchley,et al.  A comparison of frailty models for multivariate survival data. , 1995, Statistics in medicine.

[18]  Bani K. Mallick,et al.  Hierarchical Generalized Linear Models and Frailty Models with Bayesian Nonparametric Mixing , 1997 .

[19]  S. Pocock,et al.  Regression Models and Non‐Proportional Hazards in the Analysis of Breast Cancer Survival , 1984 .

[20]  W Pan,et al.  Using Frailties in the Accelerated Failure Time Model , 2001, Lifetime data analysis.

[21]  James J. Heckman,et al.  Advances in Econometrics: The identification problem in econometric models for duration data , 1983 .

[22]  Vic Barnett,et al.  Outliers in Statistical Data , 1980 .

[23]  Jerzy Neyman,et al.  OUTLIER PRONENESS OF PHENOMENA AND OF RELATED DISTRIBUTIONS , 1971 .