This article proposes a robust method of statistical inference for the Cox'
s proportional hazards model with frailties. We use the Metropolis-Hastings
algorithm and the bootstrap method. We present a computationally efficient
algorithm with a customized data structure to implement this method and de
monstrate this technique with real data.