A hybrid ML-EM algorithm for calculation of maximum likelihood estimates in semiparametric shared frailty models

This paper describes a generalised hybrid ML-EM algorithm for the calculation of maximum likelihood estimates in semiparametric shared frailty models, the Cox proportional hazard models with hazard functions multiplied by a (parametric) frailty random variable. This hybrid method is much faster than...

Ausführliche Beschreibung

Gespeichert in:
Bibliographische Detailangaben
Veröffentlicht in:Computational statistics & data analysis Jg. 40; H. 1; S. 173 - 187
Hauptverfasser: Vu, Hien T.V., Knuiman, Matthew W.
Format: Journal Article
Sprache:Englisch
Veröffentlicht: Amsterdam Elsevier B.V 28.07.2002
Elsevier Science
Elsevier
Schriftenreihe:Computational Statistics & Data Analysis
Schlagworte:
ISSN:0167-9473, 1872-7352
Online-Zugang:Volltext
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
Beschreibung
Zusammenfassung:This paper describes a generalised hybrid ML-EM algorithm for the calculation of maximum likelihood estimates in semiparametric shared frailty models, the Cox proportional hazard models with hazard functions multiplied by a (parametric) frailty random variable. This hybrid method is much faster than the standard EM method and faster than the standard direct maximum likelihood method (ML, Newton–Raphson) for large samples. We have previously applied this method to semiparametric shared gamma frailty models, and verified by simulations the asymptotic and small sample statistical properties of the frailty variance estimates. Let θ 0 be the true value of the frailty variance parameter. Then the asymptotic distribution is normal for θ 0>0 while it is a 50–50 mixture between a point mass at zero and a normal random variable on the positive axis for θ 0=0. For small samples, simulations suggest that the frailty variance estimates are approximately distributed as an x−(100− x)% mixture, 0⩽ x⩽50, between a point mass at zero and a normal random variable on the positive axis even for θ 0>0. In this paper, we apply this method and verify by simulations these statistical results for semiparametric shared log-normal frailty models. We also apply the semiparametric shared gamma and log-normal frailty models to Busselton Health Study coronary heart disease data.
ISSN:0167-9473
1872-7352
DOI:10.1016/S0167-9473(01)00099-8