Asymptotic and small sample statistical properties of random frailty variance estimates for shared gamma frailty models
Vu, H.T.V., Segal, M.R., Knuiman, M.W. and James, I.R. (2001) Asymptotic and small sample statistical properties of random frailty variance estimates for shared gamma frailty models. Communications in Statistics - Simulation and Computation, 30 (3). pp. 581-595.
*Subscription may be required
This paper concerns maximum likelihood estimation for the semiparametric shared gamma frailty model; that is the Cox proportional hazards model with the hazard function multiplied by a gamma random variable with mean 1 and variance θ. A hybrid ML-EM algorithm is applied to 26400 simulated samples of 400 to 8000 observations with Weibull hazards. The hybrid algorithm is much faster than the standard EM algorithm, faster than standard direct maximum likelihood (ML, Newton Raphson) for large samples, and gives almost identical results to the penalised likelihood method in S-PLUS 2000. When the true value θ0 of θ is zero, the estimates of θ are asymptotically distributed as a 50-50 mixture between a point mass at zero and a normal random variable on the positive axis. When θ0 > 0> the asymptotic distribution is normal. However, for small samples, simulations suggest that the estimates of θ 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 light of this, p-values and confidence intervals need to be adjusted accordingly. We indicate an approximate method for carrying out the adjustment.
|Publication Type:||Journal Article|
|Murdoch Affiliation:||School of Mathematical and Physical Sciences|
|Publisher:||Taylor and Francis|
|Item Control Page|