Assessing gamma frailty models for clustered failure time data

Proportional hazards frailty models use a random effect, so called frailty, to construct association for clustered failure time data. It is customary to assume that the random frailty follows a gamma distribution. In this paper, we propose a graphical method for assessing adequacy of the proportional hazards frailty models. In particular, we focus on the assessment of the gamma distribution assumption for the frailties. We calculate the average of the posterior expected frailties at several followup time points and compare it at these time points to 1, the known mean frailty. Large discrepancies indicate lack of fit. To aid in assessing the goodness of fit, we derive and estimate the standard error of the mean of the posterior expected frailties at each time point examined. We give an example to illustrate the proposed methodology and perform sensitivity analysis by simulations.

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

[2]  N. Laird.  Nonparametric Maximum Likelihood Estimation of a Mixing Distribution , 1978 .

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

[4]  T. Louis,et al.  Smoothing the non-parametric estimate of a prior distribution by roughening: A computational study , 1991 .

[5]  D G Clayton,et al.  A Monte Carlo method for Bayesian inference in frailty models. , 1991, Biometrics.

[6]  J. Kalbfleisch,et al.  The effects of mixture distribution misspecification when fitting mixed-effects logistic models , 1992 .

[7]  D. Oakes,et al.  Bivariate survival models induced by frailties , 1989 .

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

[9]  D. Clayton,et al.  Multivariate generalizations of the proportional hazards model , 1985 .

[10]  K G Manton,et al.  Alternative models for the heterogeneity of mortality risks among the aged. , 1986, Journal of the American Statistical Association.

[11]  Thomas A. Louis,et al.  Time-Dependent Association Measures for Bivariate Survival Distributions , 1992 .

[12]  Niels Keiding,et al.  Statistical Models Based on Counting Processes , 1993 .

[13]  D. Pierce.  The Asymptotic Effect of Substituting Estimators for Parameters in Certain Types of Statistics , 1982 .

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

[15]  Guang Guo,et al.  Estimating a Multivariate Proportional Hazards Model for Clustered Data Using the EM Algorithm, with an Application to Child Survival in Guatemala , 1992 .

[16]  T A Louis,et al.  Random effects models with non-parametric priors. , 1992, Statistics in medicine.

[17]  P. Grambsch,et al.  Martingale-based residuals for survival models , 1990 .

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

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

[20]  Larry Lee,et al.  Multivariate distributions having Weibull properties , 1979 .

[21]  Elja Arjas,et al.  A Graphical Method for Assessing Goodness of Fit in Cox's Proportional Hazards Model , 1988 .

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

[23]  Z. Ying,et al.  Checking the Cox model with cumulative sums of martingale-based residuals , 1993 .

[24]  T. Louis,et al.  Inferences on the association parameter in copula models for bivariate survival data. , 1995, Biometrics.

[25]  Brian D. Ripley,et al.  Stochastic Simulation , 2005 .

[26]  J. Heckman,et al.  A Method for Minimizing the Impact of Distributional Assumptions in Econometric Models for Duration Data , 1984 .