A relaxation of the gamma frailty (Burr) model

Frailty models are used in univariate data to account for individual heterogeneity. In the popular gamma frailty model the marginal hazard has the form of a Burr model. Although the Burr model is very useful and can offer insight on the data, it is far from perfect. The estimation of the covariate effects is linked to the baseline hazard and this makes the model coefficients hard to interpret. At the same time, the frailties are assumed constant over time, while biological reasoning in some cases may indicate that frailties may be time dependent. In this paper we present a relaxation of the Burr model which is based on loosening the link between the estimation of the covariate effects and the baseline hazard. This can be achieved by replacing the cumulative baseline hazard in the Burr model by a set of time functions, and the frailty variance by a vector of coefficients directly estimated from the data using a partial likelihood. We illustrate the similarities of the model with the Burr model and a further extension of the latter, a model with an autoregressive stochastic process for the frailty. We compare the models on simulated data sets with constant and time-dependent frailties and show how the relaxed Burr models performs on two different real data sets. We show that the relaxed Burr model serves as a good approximation to the Burr model when the frailty is constant, and furthermore it gives better results when the frailty is time dependent.