The use of Gaussian quadrature for estimation in frailty proportional hazards models
- 1 October 2007
- journal article
- Published by Wiley in Statistics in Medicine
- Vol. 27 (14), 2665-2683
- https://doi.org/10.1002/sim.3077
Abstract
In this paper, we propose a novel Gaussian quadrature estimation method in various frailty proportional hazards models. We approximate the unspecified baseline hazard by a piecewise constant one, resulting in a parametric model that can be fitted conveniently by Gaussian quadrature tools in standard software such as SAS Proc NLMIXED. We first apply our method to simple frailty models for correlated survival data (e.g. recurrent or clustered failure times), then to joint frailty models for correlated failure times with informative dropout or a dependent terminal event such as death. Simulation studies show that our method compares favorably with the well‐received penalized partial likelihood method and the Monte Carlo EM (MCEM) method, for both normal and Gamma frailty models. We apply our method to three real data examples: (1) the time to blindness of both eyes in a diabetic retinopathy study, (2) the joint analysis of recurrent opportunistic diseases in the presence of death for HIV‐infected patients, and (3) the joint modeling of local, distant tumor recurrences and patients survival in a soft tissue sarcoma study. The proposed method greatly simplifies the implementation of the (joint) frailty models and makes them much more accessible to general statistical practitioners. Copyright © 2007 John Wiley & Sons, Ltd.Keywords
This publication has 48 references indexed in Scilit:
- Frailty Survival Model Analysis of the National Deceased Donor Kidney Transplant Dataset Using Poisson Variance StructuresJournal of the American Statistical Association, 2005
- Shared parameter models for the joint analysis of longitudinal data and event timesStatistics in Medicine, 2005
- Separate and Joint Modeling of Longitudinal and Event Time Data Using Standard Computer PackagesThe American Statistician, 2004
- Evolution of Recurrent Asthma Event Rate Over Time in Frailty ModelsJournal of the Royal Statistical Society Series C: Applied Statistics, 2003
- Penalized Survival Models and FrailtyJournal of Computational and Graphical Statistics, 2003
- Analyzing Recurrent Event Data With Informative CensoringJournal of the American Statistical Association, 2001
- Panel Data with Survival: Hospitalization of HIV-Positive PatientsJournal of the American Statistical Association, 1998
- Considerations in choice of a clinical endpoint for AIDS clinical trialsStatistics in Medicine, 1994
- Estimating a Multivariate Proportional Hazards Model for Clustered Data Using the EM Algorithm, with an Application to Child Survival in GuatemalaJournal of the American Statistical Association, 1992
- Regression Analysis of Multivariate Incomplete Failure Time Data by Modeling Marginal DistributionsJournal of the American Statistical Association, 1989