Inference in mixed proportional hazard models with <Emphasis Type="Italic">K</Emphasis> random effects

A general formulation of mixed proportional hazard models with K random effects is provided. It enables to account for a population stratified at K different levels. I then show how to approximate the partial maximum likelihood estimator using an EM algorithm. In a Monte Carlo study, the behavior of the estimator is assessed and I provide an application to the ratification of ILO conventions. Compared to other procedures, the results indicate an important decrease in computing time, as well as improved convergence and stability.