Journal of Probability and Statistics

Volume 2018, Article ID 2378798, 5 pages

https://doi.org/10.1155/2018/2378798

## Frailty in Survival Analysis of Widowhood Mortality

Department of Mathematics and Statistics, UiT The Arctic University of Norway, Norway

Correspondence should be addressed to Elinor Ytterstad; on.tiu@datsretty.ronile

Received 20 March 2018; Revised 19 June 2018; Accepted 17 September 2018; Published 2 October 2018

Academic Editor: Dejian Lai

Copyright © 2018 Elinor Ytterstad. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

Heterogeneity between individuals has attracted attention in the literature of survival analysis for several decades. Widowed individuals also differ; some are more frail than others and thereby have a higher risk of dying. The traditional hazard rate in a survival model is a measure of population risk and does not provide direct information on the unobservable individual risk. A frailty model is developed and applied on a large Norwegian data set of 452 788 widowed individuals. The model seemed to fit the data well, for both widowers and widows in all age groups. The random frailty term in the model is significant, meaning that widowed persons may have individual hazard rates.

#### 1. Introduction

In survival analysis of time to death, or to any event, the hazard rate is a measure of the risk in the population at large. However, some individuals are more frail and thus have a higher risk than others. Individuals with a high risk are more likely to die early and the remaining individuals are mainly those with lower risk. This will give a "pulled down" population hazard rate as time passes [1, chapter 6]. It is therefore important to keep in mind that the observable population risk is not a measure of the individual risk. Some of the individual variation can be expressed by covariates and then included in the model. The unobservable differences between individuals are usually considered as random variation, but in survival analysis, observing over time, this individual variation cannot be ignored.

Several studies of mortality after the death of a spouse, summarized in a meta-analysis [2], demonstrated that the population hazard rate is substantially elevated immediately after the death of the spouse, before gradually declining. This is also the case after controlling for age, sex, and other covariates. Thus, the hazard rate for the widowed population is observed as being "pulled down" by time in widowhood. Many of the studies find this "pulled down" pattern to fade out at a substantial higher level than that of a married population [3–7].

The unobservable individual hazard rate can be modeled as a random variable from which the population hazard is derived. Such models are called frailty models [8], first suggested by Vaupel, Manton and Stellard [9] in mortality studies, and by Lancaster [10] in a study of unemployment duration. The mathematical properties of frailty models are thoroughly discussed in [11].

The aim of this study is to develop a frailty model to fit the characteristic widowhood mortality. Properties are discussed, and the model is applied on a large Norwegian data set.

#### 2. A Frailty Model

A commonly used and simple model for the unobserved heterogeneity between individuals is the individual proportional hazard rate: , where is a basic rate and a gamma distributed random variable with expectation . It is shown by [11] that this frailty effect will "pull down" the population hazard rate over time. Widowed mortality is previously shown to be elevated immediately after the death of the spouse and thereafter decrease, but still remaining higher than the married mortality several years into widowhood [7]. Assuming constant individual hazard ratio (HR) compared to married, but with the characteristic "pulled down" shaped population HR, our proposed individual hazard rate iswhere , , and is a known mortality rate for a married population of same sex and age as the widowed.

This will give us a constant individual HRwith expected value .

The individual survival function is , where is the known cumulative mortality rate for married population.

The population survival function is found by integrating over the distribution of , after replacing the individual married rates by an average rate for the group at risk The expectation on the right hand side is the Laplace transform of , and for a gamma distributed with and , the population survival function becomes

The population hazard rate is then derived as and the population HR can be expressed asa decreasing function of time, starting at when and declining to when , which is the long run population excess mortality. The decline towards is steeper for larger than for smaller .

The random HR in (2) is a constant function of time for each widowed person, with expected value . The population HR in (6) decreases towards over time . We have chosen a least squares’ (LS) approach to demonstrate estimation of , , and . The minimized sum of squares is , with the Nelson-Aalen survival estimator and the survival function in (4). Bootstrapping and simulation results are also included. The population HR in (6) is then estimated using these LS estimates. Comparing other estimation methods is beyond the scope of this paper.

Many of the common frailty models in the literature have a multiplicative random element in the hazard rate, and a variety of estimation procedures are available in standard statistical software [12]. However, to the author’s knowledge, there is no software available for our model (1).

#### 3. Data Example

The population of married individuals in Norway on January 1, 1975, were followed until the end of 2006 for changes in marital status and date of death. Data were provided by Statistics Norway and are described in detail elsewhere [7]. Individuals dying on the same day as their spouse are not recorded as a widower or widow. During the follow-up period, 452 807 individuals became widowed at the age of 55 to 94 years. We further excluded 19 widowed individuals dying within three weeks after their spouse and of a similar external cause, assuming their death was caused by the same accident as their spouse. The remaining sample containing 452 788 widowers and widows is used here for estimation of the frailty model.

The parameters , and in model (4) are estimated separately for widowers and widows in four age groups. The total numbers in each group are presented in Table 1. Time is days after spousal death. The LS estimations are performed in the statistical computing language R [13], and the R-code is available in Supplementary Table 2. The cumulative mortality rates for married are calculated from life tables given in [7].