The Gompertz–Makeham law of mortality is a mathematical model for the age pattern of death rates. It expresses the force of mortality (instantaneous hazard of death) as the sum of two components: an age-dependent term, the Gompertz function, that increases approximately exponentially with age, and an approximately age-independent background term known as the Makeham term. In populations where deaths from external causes are rare, the background component is often small and mortality can be well described by the simpler Gompertz law of mortality.
For adult humans, the Gompertz component captures the empirical regularity that the individual risk of death rises steeply with age. After early adulthood, many populations are reasonably described by an almost exponential increase in hazard, with the probability of dying in a given year approximately doubling every eight years. The model is not intended to describe the elevated mortality of infancy and early childhood, and at the very highest ages estimates of the mortality trajectory are sensitive to data quality and modelling assumptions. Some demographic studies report a slowing of the increase in death rates or a plateau among the oldest-old, while others find that cleaned data remain compatible with Gompertz–Makeham-type behaviour over a wide age range.
Because of its simple functional form and interpretable parameters, the Gompertz–Makeham law is widely used in actuarial science for constructing life tables and pricing life insurance and pension products, in demography and gerontology for modelling adult mortality, and in reliability theory of ageing and longevity as a parametric survival distribution for biological and technical systems.
Contents
History
In 1825 the English actuary Benjamin Gompertz published an analysis of human mortality in which he proposed that the age-specific force of mortality increases approximately exponentially with age. This "law of human mortality" was soon adopted in actuarial work on life table construction and annuity valuation, and later became known simply as the Gompertz law of mortality.
In 1860 the actuary William Makeham generalised Gompertz's formulation by adding an approximately age-independent term to represent deaths from external causes that affect all ages in a population. The resulting Gompertz–Makeham law expressed the force of mortality as the sum of an exponentially increasing intrinsic component and a constant background component, and it quickly became a standard model in actuarial science and demography.
During the twentieth century the Gompertz–Makeham model was revisited in the emerging fields of biodemography and the reliability theory of ageing and longevity, where it was used to describe regularities in adult human mortality and to compare patterns of ageing across species and populations.
Mathematical formulation
In continuous-time form, the Gompertz–Makeham law specifies the hazard function (or force of mortality) at age
x
{\displaystyle x}
as
μ
(
x
)
=
α
e
β
x
+
λ
,
x
≥
0
,
α
>
0
,
β
>
0
,
λ
≥
0.
{\displaystyle \mu (x)=\alpha e^{\beta x}+\lambda ,\qquad x\geq 0,\ \alpha >0,\ \beta >0,\ \lambda \geq 0.}
The term
Moments and life expectancy
Let
X
{\displaystyle X}
denote a non-negative random variable with Gompertz–Makeham hazard
μ
(
x
)
=
α
e
β
x
+
λ
{\displaystyle \mu (x)=\alpha e^{\beta x}+\lambda }
and survival function
S
(
x
)
{\displaystyle S(x)}
as above. The mean lifetime (life expectancy at birth) can be written as
E
[
X
]
=
∫
0
∞
S
Rate of increase and mortality doubling time
The Gompertz component of the hazard,
μ
G
(
x
)
=
α
e
β
x
{\displaystyle \mu _{\text{G}}(x)=\alpha e^{\beta x}}
, has a logarithm that is linear in age:
log
μ
G
(
x
)
=
log
α
+
β
x
.
{\displaystyle \log \mu _{\text{G}}(x)=\log \alpha +\beta x.}
The constant slope
β
{\displaystyle \beta }
Discrete-time approximation and annual death probabilities
In life-table work, mortality is often described in terms of the one-year probability of death between ages
x
{\displaystyle x}
and
x
+
1
{\displaystyle x+1}
, usually denoted
q
x
{\displaystyle q_{x}}
. For a continuous-time model with hazard
μ
(
t
)
{\displaystyle \mu (t)}
, the corresponding one-year death probability can be written as
q
x
=
1
−
exp
(
−
∫
x
x
Alternative parametrisations
Several alternative parametrisations of the Gompertz–Makeham law are used in the literature. One simple reparameterisation is to take a reference age
x
0
{\displaystyle x_{0}}
and define the background-free Gompertz hazard at that age,
h
0
=
μ
G
(
x
0
)
=
α
e
β
x
0
,
{\displaystyle h_{0}=\mu _{\text{G}}(x_{0})=\alpha e^{\beta x_{0}},}
so that the full hazard can be written as
μ
(
x
)
=
λ
+
h
Parameter estimation
Parameters of the Gompertz–Makeham model are typically estimated from observed lifetimes or grouped mortality data. For independent lifetimes
T
1
,
…
,
T
n
{\displaystyle T_{1},\ldots ,T_{n}}
with right-censoring indicators
δ
i
{\displaystyle \delta _{i}}
, the log-likelihood under a continuous-time hazard
μ
(
x
;
θ
)
{\displaystyle \mu (x;\theta )}
and survival function
S
(
x
;
θ
)
{\displaystyle S(x;\theta )}
can be written in the standard form
Illustrative example
Suppose lifetimes follow a Gompertz–Makeham law with hazard
μ
(
x
)
=
α
e
β
x
+
λ
{\displaystyle \mu (x)=\alpha e^{\beta x}+\lambda }
where
α
=
0.00003
,
β
=
0.095
,
λ
=
0.0004.
{\displaystyle \alpha =0.00003,\qquad \beta =0.095,\qquad \lambda =0.0004.}
The cumulative hazard, and survival functions are
H
(
x
)
=
λ
Empirical patterns
Studies of low-mortality populations typically find that the force of mortality is well approximated by a Gompertz–Makeham curve from roughly ages 40 to 90, with smaller deviations at younger and older ages. The Makeham term improves the fit at younger adult ages by capturing age-unrelated risks, while the Gompertz term captures the increase in intrinsic mortality with age.
Across many countries and time periods, estimates of the Gompertz–Makeham parameters imply that the individual risk of death in adulthood increases steeply with age. A commonly cited rule of thumb is that adult human mortality rates approximately double every eight years of age, although the exact rate of increase varies somewhat between populations, sexes and historical periods. In semi-logarithmic plots of the hazard function this appears as an almost straight line over a wide adult age range, which is often used as a visual check of Gompertz–Makeham behaviour.
Historical analyses indicate that much of the decline in human mortality over the twentieth century was initially driven by reductions in the age-independent (Makeham) component, due to improvements in sanitation, control of infectious diseases and safer environments. In terms of the Gompertz–Makeham parameters this appears as a marked fall in the background risk parameter with relatively stable age-dependent parameters, leading to more people surviving to old age and a more “rectangular” survival curve. Later in the century, gradual improvements at older ages affected the age-dependent component in many countries, contributing to longer remaining lifetimes at retirement ages.
At very young ages the pattern of mortality is different: mortality is relatively high in infancy, falls through childhood and only later begins the roughly exponential increase that the Gompertz–Makeham law describes well. The model therefore does not account for the specific causes of elevated infant and early-childhood mortality, which are usually handled by additional components or separate parametric models in life-table construction.
The behaviour of human mortality at the highest ages remains an area of active research. Some biodemographic studies report a slowing of the increase in death rates or the emergence of a late-life mortality plateau among the oldest-old, which would represent a deviation from a purely Gompertz–Makeham trajectory. Other work argues that much of the apparent late-life mortality deceleration can be explained by data-quality problems such as age misreporting and cohort mixing, and finds that after careful data cleaning the observed trajectories remain close to Gompertz–Makeham forms well into the tenth decade of life. As a result, the extent and universality of late-life deviations from the Gompertz–Makeham law in humans is still debated, and applications that focus on extreme old age often employ extended or alternative models.