The Additive Genetic Gamma Frailty Model

In this paper the additive genetic gamma frailty model is defined. Individual frailties are correlated as a result of an additive genetic model. An algorithm to construct additive genetic gamma frailties for any pedigree is given so that the variance–covariance structure among individual frailties equals the numerator relationship matrix times a variance. The EM algorithm can be used to estimate the parameters in the model. Calculations are similar using the EM algorithm in the shared frailty model, however the E step is not correspondingly simple. This is illustrated re‐analysing data, analysed by the shared frailty model in Nielsen et al. (1992), from the Danish adoptive register. Goodness of fit of the additive genetic gamma frailty model can be tested after analysing data with the correlated frailty model. Doing so, a “defect” in the often used and otherwise well behaving likelihood was found