Life Insurance: 4.1 Survival Distribution and Life Tables
Life Insurance: 4.1 Survival Distribution and Life Tables
Life Insurance
4.1 Survival Distribution And Life Tables
Introduction
• X, Age-at-death
• T (x), time-until-death
• Life Table
– Engineers use life tables to study the reliability of complex mechanical and
electronic systems.
– Biostatistician use life tables to compare the effectiveness of alternative
treatments of serious disease.
– Demographers use life tables as tools in population projections.
– In this text, life tables are used to build models for insurance systems
designed to assist individuals facing uncertainty about the time of
their death.
1
Probability for the Age-at-Death
1 The Survival Function
s(x) has traditionally been used as a starting point for further development
in Actuarial science and demography. In statistics and probability, the d.f.
usually plays this role.
2
Time-until-Death for a Person Age x
The conditional probability that a newborn will die between the ages x and
z, given survival to age x, is
• xp0 = s(x), x ≥ 0.
s(x + t) − s(x + t + u)
t|u qx =
s(x)
s(x + t) s(x + t) − s(x + t + u)
=
s(x) s(x + t)
= tpx · uqx+t
4
Curtate-future-life time K(x) A discrete random variable associated with
the future lifetime is the number of future years completed by (x) prior to death.
k
X
FK(x)(y) = h| qx = k+1qx,
h=0
y ≥ 0 and k is the greatest integer in y
5
Force of Mortality
FX (x + ∆x) − FX (x) fX (x)∆x
Pr(x < X ≤ x + ∆x|X > x) = =
1 − FX (x) 1 − FX (x)
fX (x) −s0(x)
µ(x)=
ˆ =
1 − FX (x) s(x)
In reliability theory, µ(x) is called failure rate or hazard rate or, hazard
rate function.
6
Force of mortality can be used to specify the distribution of X.
Z n
n px = exp[− µ(x + s)ds].
0
Z n
n p0 = exp[− µ(s)ds].
0
Z x
FX (x) = 1 − s(x) = 1 − exp[− u(s)ds]
0
Z x
0
FX (x) = fX (s) = exp[− µ(s)ds]µ(x) = xp0µ(x).
0
FT (x)(t) and fT (x)(t) denote, respectively, the d.f. and p.d.f. of T (x), the
future lifetime of (x). Then we have
d d s(x + t)
fT (x)(t) = q
t x = 1 −
dt dt s(x)
0
s(x + t) s (x + t)
= −
s(x) s(x + t)
= t px µ(x + t), t ≥ 0 7
8
Example Ā refer to the complement of the event A within the sample space
and Pr(Ā) 6= 0, the following expression and identity in probability theory:
Rewrite this identity in actuarial notation for the events A = [T (x) ≤ t] and
B = [t < T (x) ≤ 1], 0 < t < 1.
Solution: Pr(A ∪ B) becomes Pr[T (x) ≤ 1] = qx0 Pr(A) is tqx, and Pr[B|Ā]
is 1−tqx+t, hence
qx = tqx + tpx1−tqx+t.
9
Life Tables
Ration of Life Table Functions to the Survival Function
s(x + t) s(x + 1)
t qx = 1 − ⇒ qx = 1 −
s(x) s(x)
A group of l0 = 100, 000 newborns. Let L(x) denote the group’s number of
survivors to age x. We index these lives by j = 1, 2, . . . , l0 and observe that
l0
X
L(x) = Ij
i=0
l0
X
lx = E[L(x)] = = l0s(x).
j=1
10
n Dx denotes the number of deaths between ages x and x + n from among
the initial l0 lives.
1 dlx 1 ds(x)
− =− = µ(x) and − dlx = lxµ(x)dx.
lx dx s(x) dx
Since
lxµ(x) = l0xp0µ(x) = l0fX (x)
We note further that Z x
lx = l0 exp[− µ(y)dy],
0
Z x+n
lx+n = lx exp[− µ(y)dy],
x
Z x+n
lx − lx+n = ly µ(y)dy
11
x
Life Table Example
• The 1979-81 U.S. Life Table was not construct by observing 100,000 newborns
until the last survival died.
• In using the random survivorship group concept with this table, we must
make the assumption that the probabilities derived from the table will be
appropriate for the lifetimes of those who belong to the survivorship group.
• Most studies of human mortality for insurance purpose use the representation
s(x) = lx/l0. Since 100, 000s(x) is displayed only for integer value of x,
there is a need to interpolate in evaluating s(x) for noninteger values.
12
Example On the basis of The 1979-81 U.S. Life Table, evaluate the probability
that (20) will
a. Live to 100
b. Die before 70
Soltion:
s(70)−s(20) 68,248
b. s(20) = 1 − ll70
20
= 1 − 97,741 = 0.3017
13
14
15
16
• The function lxµ(x) is proportional to the p.d.f of the age-at-death of a
newborn. Since lxµ(x) is the expected density of death at age x, under the
random survivorship group idea, the graph of lxµ(x) is called the curve of
deaths.
• There is a local minimum of lxµ(x) at about age 10. The mode of the
distribution of deaths — the age at which the maximum of the curve of
death occurs — is around age 80.
d2
d d d
lxµ(x) = − lx = − 2 lx .
dx dx dx dx
17
2. Life Insurance
Introduction In this section we develop models for life insurances designed to
reduce the financial impact of the random event of untimely death.
18
Insurances Payable at the Moment of Death
• T = T (x), The elapsed time from policy issue to the death of the issued is
the insured’s future-lifetime random variable.
The first step in our analysis of a life insurance will be to define bt and vt. The
next step is to determine some characteristics of the probability distribution of
Z that are consequences of an assumed distribution for T , and we work through
these steps for several conventional insurances.
19
Level Benefit Insurance
n-year term life insurance A payment only if the insured dies within the n-year
term of an insurance commencing at issue.
vT , T ≤ n
1, t ≤ n t
bt = vt = v , t ≥ 0, Z=
0, t > n, 0, T > n,
20
Expected value of Z
Z ∞ Z n
Ā1x:n| = E[Z] = E[zT ] = ztfT (t)dt = v ttpxµx(t)dt.
0 0
Rule of Moments
Z n Z n
E[Z j ] = (v t)j tpxµx(t)dt = e−(δj)ttpxµx(t)dt
0 0
Hence
E[Z j ]@δt = E[Z]@jδt
and
V ar(Z) = 2Ā1x:n| − (Ā1x:n|)2
where 2Ā1x:n| is the actuarial present value for an n-year term insurance for a
unit amount calculated at force of interest 2δ.
21
Whole life insurance provides for a payment following the death of the insured
at any time in the future. If the payment is to be a unit amount at the moment
of death (x), then
bt = 1, vt = v t, Z = v T , t ≥ 0, T ≥ 0.
22
Example 1. The p.d.f of the future lifetime T , for (x) is assumed to be
1/80, 0 ≤ t ≤ 80
fT (t) =
0, elsewhere.
b. The variance
23
24
Example 2. For the assumption in Example 1, determine
a. Z’s d.f.
b. Z’s p.d.f
• Is age x
0.04
E[Z 2] = 102 2Āx = 100 = 25.
0.04 + 2(0.06)
and Var(Z) = 9. Hence
E[S] = 100(4) = 400, Var(S) = 100(9) = 900.
Observations:
• The 49.35 difference between this initial fund of 449.35 and the expectation
of the present value of all payments, 400, is the risk loading.
• This example used the individual risk model and a normal approximation to
the probability distribution S.
• A graph of the amount in the fund during the first 2 years for a payout
pattern when one death occurs at each time 1/8, 7/8, 9/8,13/8 and 15/8,
and two deaths occur at time 10/8 is shown in Figure 1. Between the
benefit payments, represented by the discontinuities, are exponential arcs
representing the growth of the fund at δ = 0.06.
• There are infinitely many payout patterns, each with its own graph. Both the
number of claims and times of those claims affect the fund. For example, had
the seven claims all occurred within the first instant, instead of the payout
pattern of the Figure, the fund would have dropped immediately to 379.35
27
and grown to 427.72 by the end of the second year.
28
Endowment Insurance
An n-year pure endowment provides for a payment at the end of the n years
if and only if the insured survives at least n years from the time of policy issue.
0, t ≤ n 0, T ≤n
bt = vt = v n, t ≥ 0, Z=
1, t > n, v n, T > n,
29
N-year endowment insurance provides for an amount to be payable either
following the death of the insured or upon the survival of the insured to the end
of the n-year term, whichever occurs first. If the insurance is for a unit amount
and the death benefit is payable at the moment of death, then
t
vT , T ≤ n
v , t≤n
bt = 1, t ≥ 0, vt = Z=
v n, t > n, v n, T > n,
and
Āx:n| = Ā1x:n| + A1x:n|
Var(Z) = 2Ā1x:n| − (Ā1x:n|)2 + 2A1x:n| − (A1x:n|)2 − 2Ā1x:n|A1x:n|.
30
Deferred Insurance
M-year deferred insurance provides for a benefit following the death of the
insured only if the insured dies at least m years following policy issue.
vT , T > m
1, t > m t
bt = vt = v , t ≥ 0, Z=
0, t ≤ m, 0, T ≤ n,
31
Example 4. Consider a 5-year deferred whole life insurance payable at the
moment of the death of (x). The individual is subject to a constant force of
mortality µ = 0.04. For the distribution of the present value of the benefit
payment, at δ = 0.10:
• Calculate the expectation
• Calculate the variance
• Display the distribution function
• Calculate the median ξZ0.5.
Observations:
1. The largest value of Z with nonzero probability density in this example is
e−0.1(5) = 0.6065, corresponding to T = 5.
2. The distribution of Z in this example is highly skewed to the right. While
its total mass is in the interval [0, 0.6065] and its mean is 0.1419, its median
is only 0.0573. This skewness in the direction of large positive values is
characteristic of many claim distributions in all field of insurance.
32
33
34
3. Life Annuities
Introduction
35
Life annuities play a major role in life insurance operations
• Some type of life insurance carry this concept even further and instead of
featuring a lump sum payable on death, provide stated forms of income
benefits. Thus for example, there may be monthly income payable to a
surviving spouse or to a retired insured.
36
Life annuities also have a role in disability and worker’s compensation
insurances
37
Continous Life Annuities
39
A relationship, familiar from compound interest theory
Variance of āT |
T
Var(v T ) 2Āx − (Āx)2
1−v
Var(āT |) = Var = 2
= 2
.
δ δ δ
a. āx = E[āT |]
b. Var(āT |)
40
N-year temporary life annuity of 1 per year, payable continuously while (x)
survives during the next n year, is
āT |, 0 ≤ T < n
Y =
ān|, T ≥ n
(
1−v T 1 − Āx:n|
āT | = δ , 0≤T <n 1−Z
Y = 1−v n
=⇒E[Y ] = āx:n| = E =
ān| = δ , T ≥n δ δ
2
VarZ Āx:n| − (Āx:n|)2 2 2 2
Var(Y ) = 2
= 2
= (āx:n| − āx:n| ) − (āx:n| ) .
δ δ δ
41
N- year deferred whole life annuity
Z ∞
n| āx = E[Y ] = v nāt−n|tpxµ(x + t)dt
n
Z ∞
= v n n px ās|spx+nµ(x + n + s)ds = nExāx+n.
0
Z ∞
Var(Y ) = v 2n(āt−n|)2tpxµ(x + t)dt − (n|āx)2
n
2 2n
= v npx(āx+n − 2āx+n) − (n|āx)2
δ
42
N-year certain and life annuity This is a whole life annuity with a guarantee
of payments for the first n years.
ān|, 0 ≤ T < n
Y =
āT |, T ≥ n
Z n Z ∞
āx:n| = E[Y ] = ān|tpxµ(x + t)dt +āt|tpxµ(x + t)dt
0 n
Z ∞ Z ∞
= nqxān| + āt|tpxµ(x + t)dt = ān| + v ttpxdt.
n n
Y is the sum of a constant ān| and the random variable for the n-year deferred
annuity. Hence
Further more, since Var(Y − ān|) = Var(Y ), the variance for the n-year certain
and life annuity is the same as that of the n-year deferred annuity given above.
43
d
An Expression of dx āx
Z ∞ Z ∞
d ∂
ā = vt t px dt = v ttpx[µ(x) − µ(x + t)]dt
dx 0 ∂x 0
= µ(x)āx − Āx = µ(x)āx − (1 − δāx)
= [µ(x) + δ]āx − 1.
∂
a. ∂x āx:n| = [µ(x) + δ]āx:n| − (1 − nEx).
∂
b. ∂n n| āx = −v nnpx
44
45
46
Discrete Life Annuities
Whole life annuity-due: This annuity pays a unit amount at the beginning of
each year that the annuitant (x) survives. The present value random variable is
denoted by
Y = äK+1|
where the random variable K is the curtate-future lifetime of (x).
The actuarial present value of the annuity:
∞
X
äx = E[Y ] = E[äK+1|] = äk+1|k pxqx+k
k=0
∞
X X∞
= 1+ v k+1k+1px = v k k px .
k=0 k=0
The last term of the equation above is the current payment form of the actuarial
present value for a whole life-annuity-due where the k px term is the probability
of a payment of size 1 being made at time k.
47
We can also obtian in succession
K+1
1−v 1 − Ax
äx = E =
d d
and
äx = ä∞| − ä∞|Ax, 1 = däx + Ax.
and the variance formula is
K+1
Var(v K+1) 2Ax − (Ax)2
1−v
Var(äK+1|) = Var = = .
d d2 d2
48
N-year temporary life annuity-due of 1 per year is
äK+1|, 0 ≤ K < n
Y =
än|, K≥n
n−1
X n−1
X
äx:n| = E[Y ] = äk+1|k pxqx+k + än|npx = v k k px .
k=0 k=0
v K+1, 0 ≤ K < n
Z=
v n, K≥n
is the present -value random variable for a unit of endowment insurance, payable
at the end of the year or at maturity, we have
1 − E[Z] 1 − Ax:n|
äx:n| = = =⇒1 = däx:n| + Ax:n|. 49
d d
To calculate the variance, we can use
∞
X
E[Y ] = n|äx = nExäx+n = äx − äx:n| = v k k px .
k=n
2 2n
Var(Y ) = v npx(äx+n − 2äx+n) + n|2äx − (n|ä)2.
d
50
N-year certain and life annuity-due
än|, 0≤K<n
Y =
äK+1|, K ≥ n
∞
X ∞
X
äx:n| = E[Y ] = än| nqx+ äk+1| k px qx+k = än|+ v k k px = än|+äx−äx:n|.
k=n k=n
The procedures used above for annuities-due can be adapted for annuities
immediate where payments are made at the ends of the payment period.
For instance, for a whole life annuity-immediate, the present-value random
variable is Y = aK|. Then
∞ ∞
X X 1 − (1 + i)Ax
ax = E[Y ] = k px qx+k ak| = v k k px = .
i
k=0 k=1
51
52
4. Benefit Premiums
Introduction
• benefit
• and margins for profit and for offsetting possible unfavorable experience.
53
Example 1. An insurer is planning to issue a policy to a life age 0 whose
curtate-future-lifetime, K is governed by the p.f.
k| q0 = 0.2, k = 0, 1, 2, 3, 4.
The policy will pay 1 unit at the end of the year of death in exchange for
the payment of a premium P at the beginning of each year, provided the life
survives. Find the annual premium, P , as determined by:
a. Principle I: P will be the least annual premium such that the insurer has
probability of a positive financial loss of at most 0.25.
b. Principle II: P will be the annual premium such that the insurer, using a
utility of wealth function u(x) = x will be indifferent between accepting and
not accepting the risk.
c. Principle III: P will be the annual premium such that the insurer, using
a utility of wealth function u(x) = −e−0.1x, will be indifferent between
accepting and not accepting the risk.
For all three parts assume the insurer will use an annual effective interest rate
of i = 0.06. 54
55
• Percentile premiums Although the principle is attractive on surface, it is
easy to show that it can lead to quite unsatisfactory premiums.
• Equivalence principle It has many application. If we define the insurer’s
loss, L, as the random variable of the present value of benefit to be paid
by insurer less the annuity of premiums to be paid by the insured. The
requirement of this principle is that
E[L] = 0.
Benefit premiums
E[present value of benefits] = E[present value of benefit premiums]
Single benefit premium When the equivalence principle is used to determine
a single premium at policy issue for a life insurance or a life annuity, the
premium is equal to the actuarial present value of the benefit payments.
• Exponential premiums These premiums are nonproportional in the sense
that the premium for the policy with a benefit level of 10 is more than 10
times the premium for a policy with a benefit level of 1. This is consistent
for a risk averse utility function.
56
Fully Continuous Premiums
• Fully continuous level annual benefit premium for a unit whole life insurance
payable immediately on the death of (x).
l(t) = v t − P̄ āt|,
the present value of the loss to the issuer if death occurs at time t.
– l(0) = 1.
– l(t) → −P̄ /δ as t → ∞.
– if l(t0) = 0, death before t0 results in a positive loss, whereas death after
t0 produces a negative loss.
57
Consider the loss random variable
L = l(T ) = v T − P̄ āT |.
Solution: These assumptions yields āx = 10, Āx = 0.4, and 2Āx = 0.25.
Then we obtain
Āx
P̄ (Āx) = = 0.04,
āx
and
0.25 − 0.16
Var(L) = 2
= 0.25.
(0.6)
60
• Under the constant force of mortality assumption,
µ(µ + δ)−1
P̄ (Āx) = −1
= µ,
(µ + δ)
which does not depend on the force of interest interest or the age at issue.
• For a variety of fully continuous life insurances, the formulas for annual
premiums can be determined by
E[bT vT ]
bT vT − P̄ Y = Z − P̄ Y and E[bT vT − P̄ Y ] = 0 =⇒ P̄ = .
E[Y ]
• For an n-year deferred whole life annuity of 1 per year payable continuously.
In this case bT vT = 0, T ≤ n and bT vT = āT −n|v n, T > n Then
Example 4. Find the 25th percentile premium for an insured age 55 for the
following plans of insurances:
a. 20-year endowment
b. 20-year term
c. 10-year term.
Assume a fully continuous basis with a force of interest, δ = 0.06 and mortality
following the illustrative life Table.
62
63
Example 5. This example builds on Example 4.except that T (55) has a De
Movivre distribution with p.d.f
For the three loss variables, display the d.f. of L and determine the parameter
P̄ as the smallest non-negative number such that Pr(L > 0) ≤ 0.25.
64
65
Fully Discrete Premiums
• Px, the level annual benefit premium for a unit whole life insurance.
• the absence of (Āx) means that the insurance is payable at the end of policy
year of death.
• The loss of this insurance is
L = v K+1 − PxäK+1|, k = 0, 1, 2, . . . .
which yields Ax
Px =
äx
2
and Ax − (Ax)2
Var(L) = .
(däx)2
66
67
Example 6. If k|qx = c(0.96)k+1, k = 0, 1, 2, . . . where c = 0.04/0.96 and
i = 0.06, calculate Px and Var(L).
E[bK+1vK+1]
E[bK+1vK+1 − P Y ] = 0 =⇒ P =
E[Y ]
68
Example 8. Consider a 10,000 fully discrete whole life insurance. Let π
denote an annual premium for this policy and L(π) denote the loss-at issue
random variable for one such policy on the basis of the Illustrative Life Table,
6% interest and issue age 35
a. Determine the premium, πa, such that the distribution of L(πa) has mean 0.
Calculate the variance of L(πa).
b. Approximate the smallest non-negative premium, πb, such that the probability
is less than 0.5 that the loss L(πb) is positive. Find the variance of L(πb).
c. Determine the premium, πc, such that the probability of a positive total loss
on 100 such independent policies is 0.05 by the normal approximation.
69