Force of mortality

In actuarial science, force of mortality represents the instantaneous rate of mortality at a certain age measured on an annualized basis. It is identical in concept to failure rate, also called hazard function, in reliability theory.

Motivation and definition

In a life table, we consider the probability of a person dying from age x to x + 1, called qx. In the continuous case, we could also consider the conditional probability of a person who has attained age (x) dying between ages x and x + Δx, which is

P_{x}(\Delta x)=P(x<X<x+\Delta\;x\mid\;X>x)=\frac{F_X(x+\Delta\;x)-F_X(x)}{(1-F_X(x))}

where FX(x) is the cumulative distribution function of the continuous age-at-death random variable, X. As Δx tends to zero, so does this probability in the continuous case. The approximate force of mortality is this probability divided by Δx. If we let Δx tend to zero, we get the function for force of mortality, denoted by \mu(x):

\mu\,(x)=\frac{F'_X(x)}{1-F_X(x)}

Since fX(x)=F 'X(x) is the probability density function of X, and S(x) = 1 - FX(x) is the survival function, the force of mortality can also be expressed variously as:

\mu\,(x)=\frac{f_X(x)}{1-F_X(x)}=-\frac{S'(x)}{S(x)}=-{\frac{d}{dx}}\ln[S(x)].

To understand conceptually how the force of mortality operates within a population, consider that the ages, x, where the probability density function fX(x) is zero, there is no chance of dying. Thus the force of mortality at these ages is zero.

The force of mortality \mu(x) can be interpreted as the conditional density of failure at age x, while f(x) is the unconditional density of failure at age x.[1] The unconditional density of failure at age x is the product of the probability of survival to age x, and the conditional density of failure at age x, given survival to age x.

This is expressed in symbols as

\,\mu(x)S(x) = f_X(x)

or equivalently

\mu(x) = \frac{f_X(x)}{S(x)}.

In many instances, it is also desirable to determine the survival probability function when the force of mortality is known. To do this, integrate the force of mortality over the interval x to x + t

 \int_{x}^{x+t} -\frac{d}{dy} \ln[S(y)]\, dy .

By the fundamental theorem of calculus, this is simply

 \ln[S(x + t)] - \ln[S(x)],

and taking the exponent to the base e results in

 \frac{S(x+t)}{S(x)} = S_x(t).

Therefore, the survival probability of an individual of age x is written in terms of the force of mortality as

 S_x(t) = e^{-\int_x^{x+t}\mu(y)\, dy\,}.

Examples

This example is taken from.[2] A survival model follows Makeham's law if the force of mortality is

\mu(y) = A + Bc^{y} \quad\text{for } y \geqslant 0.

Using the last formula, we have

\int_x^{x+t} A + Bc^{y} dy = At + B(c^{x+t}-c^{x})/\ln[c].

Then

S_x(t) = e^{-(At + B(c^{x+t}-c^{x})/\ln[c])} = e^{-At} g^{c^{x} (c^{t}-1)}

where g=e^{-B/\ln[c]}.

See also

References

  1. R. Cunningham, T. Herzog, R. London (2008). Models for Quantifying Risk, 3rd Edition, Actex.
  2. Dickson, David C.M., Cambridge (2009). Actuarial Mathematics for Life Contingent Risks, First Edition, Cambridge University Press.

External links