Bias correction in the hierarchical likelihood approach to the analysis of multivariate survival data.

A simple and easy-to-implement bias correction method for the H-likelihood estimators under a shared frailty model, which incorporates complex dependence structure within clusters and is also applicable to the penalized likelihood estimators.