Elementy teorii przeżywalności
Podstawą wszystkich rozważań w teorii ubezpieczeń na życie jest funkcja przeżycia
s(x). Podaje ona prawdopodobieństwo zdarzenia, że noworodek losowo wybrany z danej
populacji dożyje wieku x lat. W literaturze ubezpieczeniowej dla oznaczenia osoby w wie-
ku x lat często stosowany jest zapis symboliczny (x). Przyjmuje się, że model opisujący
czas życie ludzkiego jest pewną ciągłą zmienną losową.
Niech
:
X
Ω →
Ω →
Ω →
Ω →
Ñ oznacza zmienną losową przyjmującą wartości równe czasowi
ż
ycia noworodka losowo wybranego z danej populacji
Ω
. Niech
( )
Pr(
)
F x
X
x
=
≤
=
≤
=
≤
=
≤
,
0
x
≥≥≥≥
oznacza dystrybuantę zmiennej losowej X, czyli prawdopodobieństwo, że zmienna losowa
X przyjmie wartości nie przekraczające x, a zatem prawdopodobieństwo zdarzenia, że
ś
mierć nastąpiła w przedziale [0, x]. Zakładamy, że F(0)=0. Mamy zatem
( )
1
( )
Pr(
)
s x
F x
X
x
= −
=
>
= −
=
>
= −
=
>
= −
=
>
.
Funkcja przeżywalności s(x) jest malejąca.
Przyjmujemy też, że istnieje kres górny trwania życia, tzn. taki wiek w, że
Pr(
X
>
w
) = 0
oraz
Pr(
)
0
X
w
> − ε >
> − ε >
> − ε >
> − ε >
dla każdego
0
>
ε
. (W Polsce przyjmuje się w
= 100 lub 110 lat).
Prawdopodobieństwo warunkowe, że noworodek umrze między wiekiem x a wie-
kiem z, przy założeniu, że przeżył wiek x, obliczamy ze wzoru
( )
( )
( )
( )
Pr(
|
)
1
( )
( )
F z
F x
s x
s z
x
X
z X
x
F x
s x
−
−
−
−
−
−
−
−
<
≤
>
=
=
<
≤
>
=
=
<
≤
>
=
=
<
≤
>
=
=
−−−−
.
Niech T oznacza przyszły czas życia osoby (x), tzn., osoby, która ma x lat. Zatem
T = X – x jest też zmienną losową przyjmującą wartości w przedziale [0, w – x]. Niech
( )
Pr( ( )
)
G t
T x
t
=
≤
=
≤
=
≤
=
≤
,
0
t
≥≥≥≥
będzie dystrybuantą zmiennej losowej T(x). W dalszej części wykładu dla uproszczenia
zapisu będziemy zakładać, że kres górny trwania życia jest równy +
∞
.
Zmienne losowe X oraz T są to ciągłe zmienne losowe
W rachunku aktuarialnym przyjmuje się następujące oznaczenia:
t
x
p
−
prawdopodobieństwo, że osoba obecnie w wieku x przeżyje następnych t lat, tzn.
dożyje wieku x+t;
t
x
q
−
prawdopodobieństwo, że osoba obecnie w wieku x nie przeżyje kolejnych t lat,
tzn. nie dożyje wieku x+t.
Oczywiście zachodzi równość
1
t
x
t
x
p
q
+
=
+
=
+
=
+
=
.
Umowa
. Jeśli t = 1, to przyjmujemy:
1
x
x
p
p
====
oraz
1
x
x
q
q
====
.
Z powyższego widać, że
t
x
q
jest dystrybuantą zmiennej losowej T wyznaczającej
dalszy czas trwania życia osoby (x), tzn. zachodzą równości:
Pr( ( )
)
t
x
q
T x
t
=
≤
=
≤
=
≤
=
≤
,
1
Pr( ( )
)
t
x
t
x
p
q
T x
t
= −
=
>
= −
=
>
= −
=
>
= −
=
>
.
Ponadto mamy
Pr(
|
)
Pr( ( )
)
x
X
z X
x
T x
z
x
<
≤
>
=
≤ −
<
≤
>
=
≤ −
<
≤
>
=
≤ −
<
≤
>
=
≤ −
.
Przyjmijmy jeszcze jedno oznaczenie:
|
Pr(
( )
)
t u
x
q
t
T x
t
u
=
<
≤ +
=
<
≤ +
=
<
≤ +
=
<
≤ +
.
Oznacza ono prawdopodobieństwo, że x–latek przeżyje jeszcze t lat, a następnie umrze w
przeciągu czasu u.
Twierdzenie. Zachodzą następujące równości
(a)
(
)
( )
t
x
s x
t
p
s x
++++
====
,
(b)
( )
(
)
( )
t
x
s x
s x
t
q
s x
−
+
−
+
−
+
−
+
====
,
(c)
|
t u
x
t
x
u
x t
q
p
q
++++
=
⋅
=
⋅
=
⋅
=
⋅
,
(d)
t u
x
t
x
u
x t
p
p
p
+
+
+
+
+
+
+
+
=
⋅
=
⋅
=
⋅
=
⋅
,
(e)
1
0
(
)
n
n
x
k
x
x k
k
q
p
q
−−−−
++++
====
=
⋅
=
⋅
=
⋅
=
⋅
∑
∑
∑
∑
, gdzie
0
1
x
p
====
.
Równości te wynikają z definicji prawdopodobieństwa warunkowego oraz ze wzoru
Pr(
|
)
Pr( ( )
)
x
X
x
k X
x
T x
k
<
≤ +
>
=
≤
<
≤ +
>
=
≤
<
≤ +
>
=
≤
<
≤ +
>
=
≤
.
Wprowadźmy teraz model, w którym czas dalszego trwania życia jest całkowity.
Przez K(x) oznaczmy zmienną losową, która opisuje liczbę pełnych lat, jakie pozostały do
przeżycia osobie będącej obecnie w wieku x.
Rozkład prawdopodobieństwa zmiennej losowej K(x) określony jest następująco:
1
( )
Pr( ( )
)
Pr(
( )
1)
k
x
k
x
g k
K x
k
k
T x
k
p
p
++++
=
=
=
≤
< + =
−
=
=
=
≤
< + =
−
=
=
=
≤
< + =
−
=
=
=
≤
< + =
−
|
k
x
x k
k
x k
p
q
q
+
+
+
+
+
+
+
+
=
⋅
=
=
⋅
=
=
⋅
=
=
⋅
=
.
Dystrybuanta zmiennej losowej K(x) jest funkcją schodkową postaci
1|
|
0
( )
...
k
n
x
x
x
k
x
n
G k
q
q
q
q
====
=
=
+
+ +
=
=
=
+
+ +
=
=
=
+
+ +
=
=
=
+
+ +
=
∑
∑
∑
∑
2
1
1
1
(1
)
(
)
...
(
)
1
x
x
x
k
x
k
x
k
x
k
x
p
p
p
p
p
p
q
+
+
+
+
+
+
+
+
+
+
+
+
−
+
−
+ +
−
= −
=
−
+
−
+ +
−
= −
=
−
+
−
+ +
−
= −
=
−
+
−
+ +
−
= −
=
Intensywność wymierania (natężenie umieralności)
Chcielibyśmy czasem mieć możliwość oceny prawdopodobieństwa zgonu nie w
pewnym przedziale czasu, ale lokalnie w danym momencie t. Prawdopodobieństwo
x
x
p
jest równe 0. Musimy więc rozważać niezerowe przedziały i dokonać przejścia graniczne-
go, czyli innymi słowy posłużyć się pojęciem pochodnej.
Definicja. Intensywnością wymierania nazywamy następująco określoną funkcję
0
Pr(
)
lim
Pr(
)
x
x
x
X
x
x
x
X
x
∆ →
∆ →
∆ →
∆ →
<
< + ∆
<
< + ∆
<
< + ∆
<
< + ∆
µ =
µ =
µ =
µ =
∆
>
∆
>
∆
>
∆
>
,
0
x
≥≥≥≥
.
Jeśli dystrybuanta F(x) ma pochodną – oznaczmy ją jako f(x) – to funkcję intensyw-
ności wymierania definiuje się następująco:
( )
( )
1
( )
( )
x
f x
f x
F x
s x
µ =
=
µ =
=
µ =
=
µ =
=
−−−−
.
Funkcja ta nosi też miano funkcji hazardu. Związek między funkcją intensywności
wymierania
x
µµµµ
i funkcją przeżycia
( )
s x
wyraża się wzorem
( )
ln ( )
( )
x
s x
d
s x
s x
dx
′′′′
µ = −
= −
µ = −
= −
µ = −
= −
µ = −
= −
.
Korzystając z tego wzoru można prawdopodobieństwa
t
x
p
i
t
x
q
wyrazić za pomocą
funkcji natężenia umieralności
( )
x
µµµµ
. Całkując powyższą równość w przedziale (x, x + t)
otrzymujemy
(
)
ln
ln
( )
x t
y
t
x
x
s x
t
dy
p
s x
++++
++++
µ
= −
= −
µ
= −
= −
µ
= −
= −
µ
= −
= −
∫∫∫∫
.
Stąd
exp
x t
t
x
y
x
p
dy
++++
=
− µ
=
− µ
=
− µ
=
− µ
∫∫∫∫
.
Można też pokazać, że
0
( )
x t
t
t
x
y
x
x y
x
q
f y dy
p
dy
++++
++++
=
=
µ
=
=
µ
=
=
µ
=
=
µ
∫
∫
∫
∫
∫
∫
∫
∫
.
Omówione wyżej charakterystyki trwania życia mogą być wyznaczone tylko wów-
czas, gdy znany jest rozkład trwania życia. Nie ma jednego ustalonego rozkładu przyszłe-
go życia. Każda populacja osób w wieku x lat może mieć inny rozkład przyszłego czasu
ż
ycia.
Tablice trwania życia są konstrukcją teoretyczną umożliwiającą prowadzenie szcze-
gółowej analizy procesu wymierania badanej populacji. Podstawę do budowy każdej tabli-
cy życia stanowią dokładne informacje o:
−
strukturze ludności według płci i roczników urodzenia (lub o strukturze według miej-
sca zamieszkania, lub grup społecznych),
−
zgonów według tych cech.
Konieczność posiadania tych informacji sprawia, że TTś sporządza się na ogół dla
okresów bliskich spisom ludności.
Podstawą szacowania parametrów prawdopodobieństwa zgonu były między innymi
dane o zgonach. Przy budowie TTś przyjmuje się, że zgony mają miejsce na początku ro-
ku. Z tego wynika, że na podstawie TTś możemy określić rozkład zmiennej losowej sko-
kowej K(x) związanej z dalszym trwaniem życia osoby w wieku x. Chcąc uzyskać rozkład
ciągłej zmiennej losowej K(x), musimy uwzględnić fakt, że zgony mogą zdarzyć się w
każdym momencie. W tym celu przyjmuje się jedno z następujących założeń dotyczących
wymieralności między całkowitymi liczbami lat x i x + 1.
1.Jednostajna umieralność w ciągu roku (UDD) (ang. uniform distribution of deaths).
Zakłada się, że rozkład zgonów miedzy całkowitymi liczbami lat jest równomierny:
(
)
(1
) ( )
(
1)
s x
t
t s x
t s x
+ = −
+ ⋅
+
+ = −
+ ⋅
+
+ = −
+ ⋅
+
+ = −
+ ⋅
+
,
0
1
t
≤ ≤
≤ ≤
≤ ≤
≤ ≤
oraz
0,1, 2,... .
x
====
Z jednostajności rozkładu zgonów w ciągu roku wynika liniowość prawdopodobieństwa
t
x
q
względem t w przedziale [0, 1), czyli
t
x
x
q
t q
= ⋅
= ⋅
= ⋅
= ⋅
.
2.Stała intensywność wymieralności. Zakłada się, że funkcja
x t
++++
µµµµ
ma dla każdego
(0,1)
t
∈
∈
∈
∈
stałą wartość równą:
ln
x
p
µ = −
µ = −
µ = −
µ = −
. Przy tym założeniu mamy
(
)
( )
t
s x
t
s x
e
−µ
−µ
−µ
−µ
+ =
⋅
+ =
⋅
+ =
⋅
+ =
⋅
3.Założenie Balducciego. Założenie to jest określone wzorem
1
(1
)
t
x
x
q
t q
−−−−
= −
= −
= −
= −
.
Idea tego założenia polega na liniowej interpolacji odwrotności funkcji przeżycia:
1
1
1
(1
)
(
)
( )
(
1)
t
t
s x
t
s x
s x
= −
+
= −
+
= −
+
= −
+
+
+
+
+
+
+
+
+
, 0
1
t
≤ ≤
≤ ≤
≤ ≤
≤ ≤
.