Volume 31 - Article 19 | Pages 553-592 Author has provided data and code for replicating results

A matrix approach to the statistics of longevity in heterogeneous frailty models

By Hal Caswell

Print this page  Facebook  Twitter


Date received:25 Jan 2014
Date published:02 Sep 2014
Word count:6125
Keywords:frailty, heterogeneity, individual stochasticity, longevity, Markov chains, mortality
Additional files:readme.31-19 (text file, 730 Byte)
 demographic-research.31-19 (zip file, 3 kB)


Background: The gamma-Gompertz model is a fixed frailty model in which baseline mortality increases exponentially with age, frailty has a proportional effect on mortality, and frailty at birth follows a gamma distribution. Mortality selects against the more frail, so the marginal mortality rate decelerates, eventually reaching an asymptote. The gamma-Gompertz is one of a wider class of frailty models, characterized by the choice of baseline mortality, effects of frailty, distributions of frailty, and assumptions about the dynamics of frailty.

Objective: To develop a matrix model to compute all the statistical properties of longevity from the gamma-Gompertz and related models.

Methods: I use the vec-permutation matrix formulation to develop a model in which individuals are jointly classified by age and frailty. The matrix is used to project the age and frailty dynamics of a cohort and the fundamental matrix is used to obtain the statistics of longevity.

Results: The model permits calculation of the mean, variance, coefficient of variation, skewness and all moments of longevity, the marginal mortality and survivorship functions, the dynamics of the frailty distribution, and other quantities. The matrix formulation extends naturally to other frailty models. I apply the analysis to the gamma-Gompertz model (for humans and laboratory animals), the gamma-Makeham model, and the gamma-Siler model, and to a hypothetical dynamic frailty model characterized by diffusion of frailty with reflecting boundaries. The matrix model permits partitioning the variance in longevity into components due to heterogeneity and to individual stochasticity. In several published human data sets, heterogeneity accounts for less than 10% of the variance in longevity. In laboratory populations of five invertebrate animal species, heterogeneity accounts for 46% to 83% of the total variance in longevity.

Author's Affiliation

Hal Caswell - Universiteit van Amsterdam, Netherlands [Email]

Other articles by the same author/authors in Demographic Research

» The sensitivity analysis of population projections
Volume 33 - Article 28

» Lifetime reproduction and the second demographic transition: Stochasticity and individual variation
Volume 33 - Article 20

» Demography and the statistics of lifetime economic transfers under individual stochasticity
Volume 32 - Article 19

» Why do lifespan variability trends for the young and old diverge? A perturbation analysis
Volume 30 - Article 48

» Reproductive value, the stable stage distribution, and the sensitivity of the population growth rate to changes in vital rates
Volume 23 - Article 19

» Perturbation analysis of nonlinear matrix population models
Volume 18 - Article 3

Most recent similar articles in Demographic Research

» Attrition in heterogeneous cohorts
Volume 23 - Article 26    | Keywords: frailty, heterogeneity, mortality

» The relative tail of longevity and the mean remaining lifetime
Volume 14 - Article 7    | Keywords: frailty, heterogeneity, mortality

» Demography and the statistics of lifetime economic transfers under individual stochasticity
Volume 32 - Article 19    | Keywords: individual stochasticity, Markov chains

» Sex ratios at sexual maturity and longevity: Evidence from Swedish register data
Volume 29 - Article 31    | Keywords: longevity, mortality

» The mystery of Japan's missing centenarians explained
Volume 26 - Article 11    | Keywords: longevity, mortality


»Volume 31





Similar Articles



Jump to Article

Volume Page
Volume Article ID