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 rando
m frailty term. Best linear unbiased predictors are used as an initial
step in the computation of maximum likelihood and residual maximum li
kelihood estimates.