REML estimation for survival models with frailty

A method of estimation for generalised mixed models is applied to the estimation of regression parameters in proportional hazards models for failure times when there are repeated observations of failure on each subject. The subject effect is incorporated into the model as a random frailty term. Best linear unbiased predictors are used as an initial step in the computation of maximum likelihood and residual maximum likelihood estimates.