Skew-t fits to mortality data--can a Gaussian-related distribution replace the Gompertz-Makeham as the basis for mortality studies?

J Gerontol A Biol Sci Med Sci. 2013 Aug;68(8):903-13. doi: 10.1093/gerona/gls239. Epub 2012 Dec 10.

Abstract

Gompertz-related distributions have dominated mortality studies for 187 years. However, nonrelated distributions also fit well to mortality data. These compete with the Gompertz and Gompertz-Makeham data when applied to data with varying extents of truncation, with no consensus as to preference. In contrast, Gaussian-related distributions are rarely applied, despite the fact that Lexis in 1879 suggested that the normal distribution itself fits well to the right of the mode. Study aims were therefore to compare skew-t fits to Human Mortality Database data, with Gompertz-nested distributions, by implementing maximum likelihood estimation functions (mle2, R package bbmle; coding given). Results showed skew-t fits obtained lower Bayesian information criterion values than Gompertz-nested distributions, applied to low-mortality country data, including 1711 and 1810 cohorts. As Gaussian-related distributions have now been found to have almost universal application to error theory, one conclusion could be that a Gaussian-related distribution might replace Gompertz-related distributions as the basis for mortality studies.

Keywords: Graduation by mathematical formula; Maximum likelihood estimation; Parametric fits; bbmle; mle2.

Publication types

  • Research Support, Non-U.S. Gov't

MeSH terms

  • Adult
  • Aged
  • Aged, 80 and over
  • Bayes Theorem
  • Cohort Studies
  • Databases, Factual / statistics & numerical data
  • Female
  • Humans
  • Life Tables
  • Likelihood Functions
  • Male
  • Middle Aged
  • Models, Statistical*
  • Mortality*
  • Normal Distribution