Modeling Spatial Survival Data Using Semiparametric Frailty Models

We propose a new class of semiparametric frailty models for spatially correlated survival data. Specifically, we extend the ordinary frailty models by allowing random effects accommodating spatial correlations to enter into the baseline hazard function multiplicatively. We prove identifiability of the models and give sufficient regularity conditions. We propose drawing inference based on a marginal rank likelihood. No parametric forms of the baseline hazard need to be assumed in this semiparametric approach. Monte Carlo simulations and the Laplace approach are used to tackle the intractable integral in the likelihood function. Different spatial covariance structures are explored in simulations and the proposed methods are applied to the East Boston Asthma Study to detect prognostic factors leading to childhood asthma.