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...
Uloženo v:
| Vydáno v: | Computational statistics & data analysis Ročník 40; číslo 1; s. 173 - 187 |
|---|---|
| Hlavní autoři: | , |
| Médium: | Journal Article |
| Jazyk: | angličtina |
| Vydáno: |
Amsterdam
Elsevier B.V
28.07.2002
Elsevier Science Elsevier |
| Edice: | Computational Statistics & Data Analysis |
| Témata: | |
| ISSN: | 0167-9473, 1872-7352 |
| On-line přístup: | Získat plný text |
| Tagy: |
Přidat tag
Žádné tagy, Buďte první, kdo vytvoří štítek k tomuto záznamu!
|
| Shrnutí: | 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 |