Semiparametric estimation of random effects using the Cox model based on the EM algorithm

Consider a survival experiment where individuals within a certain subset of the population share a common, unobservable, random frailty. Such a frailty could be an unobservable genetic or early environmental effect if individuals were in sibling groups or an environmental effect if individuals were grouped by households. Suppose that if the frailty, omega, is known, the Cox proportional hazards model for the observable covariates is valid with the consequence of the random effect being a multiplicative factor on the hazard rate. Assuming tht the random frailties follow a gamma distribution, estimates of the fixed and random effects are obtained by using an EM algorithm based on a profile likelihood construction. The method developed is applied to the Framingham Heart Study to examine the risks of smoking and cholesterol levels, adjusting for potential random effects.