frailtypack: A computer program for the analysis of correlated failure time data using penalized likelihood estimation

Correlated survival outcomes occur quite frequently in the biomedical research. Available software is limited, particularly if we wish to obtain smoothed estimate of the baseline hazard function in the context of random effects model for correlated data. The main objective of this paper is to describe an R package called frailtypack that can be used for estimating the parameters in a shared gamma frailty model with possibly right-censored, left-truncated stratified survival data using penalized likelihood estimation. Time-dependent structure for the explanatory variables and/or extension of the Cox regression model to recurrent events are also allowed. This program can also be used simply to obtain directly a smooth estimate of the baseline hazard function. To illustrate the program we used two data sets, one with clustered survival times, the other one with recurrent events, i.e., the rehospitalizations of patients diagnosed with colorectal cancer. We show how to fit the model with recurrent events and time-dependent covariates using Andersen-Gill approach.

[1]  J. Ramsay.  Monotone Regression Splines in Action , 1988 .

[2]  N Mantel,et al.  Mantel-Haenszel analyses of litter-matched time-to-response data, with modifications for recovery of interlitter information. , 1977, Cancer research.

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

[4]  D. Marquardt.  An Algorithm for Least-Squares Estimation of Nonlinear Parameters , 1963 .

[5]  P. Grambsch,et al.  Modeling Survival Data: Extending the Cox Model , 2000 .

[6]  Peter J Green,et al.  Lead discussion on 'Varying coefficient models' by Hastie & Tibshirani , 1993 .

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

[8]  D. Commenges,et al.  Maximum Penalized Likelihood Estimation in a Gamma-Frailty Model , 2003, Lifetime data analysis.

[9]  D. Clayton.  A model for association in bivariate life tables and its application in epidemiological studies of familial tendency in chronic disease incidence , 1978 .

[10]  Luc Duchateau,et al.  Evolution of recurrent asthma event rate over time in frailty models , 2003 .

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

[12]  Patrick J. Kelly,et al.  A Review of Software Packages for Analyzing Correlated Survival Data , 2004 .

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

[14]  K. Liang,et al.  Asymptotic Properties of Maximum Likelihood Estimators and Likelihood Ratio Tests under Nonstandard Conditions , 1987 .

[15]  R. Tibshirani,et al.  Varying‐Coefficient Models , 1993 .

[16]  L. J. Wei,et al.  Regression analysis of multivariate incomplete failure time data by modeling marginal distributions , 1989 .

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

[18]  A Alioum,et al.  PHMPL: a computer program for hazard estimation using a penalized likelihood method with interval-censored and left-truncated data. , 1999, Computer methods and programs in biomedicine.

[19]  P. Hougaard,et al.  Frailty models for survival data , 1995, Lifetime data analysis.

[20]  J. Klein,et al.  Estimating Random Effects In the Framingham Heart Study , 1992 .

[21]  Esteve Fernandez,et al.  Sex differences in hospital readmission among colorectal cancer patients , 2005, Journal of Epidemiology and Community Health.