Generalized gamma frailty model

In this article, we present a frailty model using the generalized gamma distribution as the frailty distribution. It is a power generalization of the popular gamma frailty model. It also includes other frailty models such as the lognormal and Weibull frailty models as special cases. The flexibility of this frailty distribution makes it possible to detect a complex frailty distribution structure which may otherwise be missed. Due to the intractable integrals in the likelihood function and its derivatives, we propose to approximate the integrals either by Monte Carlo simulation or by a quadrature method and then determine the maximum likelihood estimates of the parameters in the model. We explore the properties of the proposed frailty model and the computation method through a simulation study. The study shows that the proposed model can potentially reduce errors in the estimation, and that it provides a viable alternative for correlated data. The merits of proposed model are demonstrated in analysing the effects of sublingual nitroglycerin and oral isosorbide dinitrate on angina pectoris of coronary heart disease patients based on the data set in Danahy et al. (sustained hemodynamic and antianginal effect of high dose oral isosorbide dinitrate. Circulation 1977; 55:381–387). Copyright © 2005 John Wiley & Sons, Ltd.

[1]  M. Greenwood,et al.  An Inquiry into the Nature of Frequency Distributions Representative of Multiple Happenings with Particular Reference to the Occurrence of Multiple Attacks of Disease or of Repeated Accidents , 1920 .

[2]  E. Stacy.  A Generalization of the Gamma Distribution , 1962 .

[3]  R. Prentice.  A LOG GAMMA MODEL AND ITS MAXIMUM LIKELIHOOD ESTIMATION , 1974 .

[4]  W. Aronow,et al.  Sustained Hemodynamic and Antianginal Effect of High Dose Oral Isosorbide Dinitrate , 1977, Circulation.

[5]  D. Rubin,et al.  Maximum likelihood from incomplete data via the EM - algorithm plus discussions on the paper , 1977 .

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

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

[8]  P. Hougaard.  A class of multivanate failure time distributions , 1986 .

[9]  P. Hougaard.  Survival models for heterogeneous populations derived from stable distributions , 1986 .

[10]  McGilchrist Ca,et al.  Regression with frailty in survival analysis. , 1991 .

[11]  K. Yamaguchi.  Accelerated Failure-Time Regression Models with a Regression Model of Surviving Fraction: An Application to the Analysis of “Permanent Employment” in Japan , 1992 .

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

[13]  Odd O. Aalen,et al.  Modelling Heterogeneity in Survival Analysis by the Compound Poisson Distribution , 1992 .

[14]  J P Klein,et al.  Semiparametric estimation of random effects using the Cox model based on the EM algorithm. , 1992, Biometrics.

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

[16]  M. Tanner.  Tools for statistical inference: methods for the exploration of posterior distributions and likeliho , 1994 .

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

[18]  Andrew Pickles,et al.  Generalizations and applications of frailty models for survival and event data , 1994, Statistical methods in medical research.

[19]  N. L. Johnson,et al.  Continuous Univariate Distributions. , 1995 .

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

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

[22]  Ross Ihaka,et al.  Gentleman R: R: A language for data analysis and graphics , 1996 .

[23]  A Marginal Likelihood Approach to Estimation in Frailty Models , 1997 .

[24]  J W Denham,et al.  A generalized F mixture model for cure rate estimation. , 1998, Statistics in medicine.

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

[26]  J. Lawless.  Statistical Models and Methods for Lifetime Data , 2002 .

[27]  J. Kalbfleisch,et al.  The Statistical Analysis of Failure Time Data: Kalbfleisch/The Statistical , 2002 .

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