Основы страховой (актуарной) математики - Г.М.Кошкин
.pdf1)
ex:ne< n;
2)åñëè n = ! è x = 0; òî e0:!e=e = !=2;
Также3) еслинетрудноx = ! n;получитьтоe выражение= n=2: для дисперсии:
(! n):ne
D min(T (x); n) = |
n3 |
|
n4 |
|
|||
|
|
|
|
: |
(3.4.2) |
||
3(! x) |
4(! x)2 |
||||||
Из формулы (3.4.2) следует, что при n = ! è x = 0 |
|
||||||
!2 |
|
|
|
||||
D min(T (x); !) = |
|
|
; |
|
|
||
|
|
|
|
12
т.е. получаем известный результат для дисперсии равномерного закона распределенияРассмотримнадляинтервале (0; !):
! = 90 лет два конкретных случая при n = 5 è nжительностей= 10 лет. Значенияжизни,среднихвычисленныеи дисперсийпо формуламчастичных(3.4.средних1)и (3.4.продол2),для-
ряда возрастов и n = 5; 10 сведены в следующую таблицу:
x |
0 |
10 |
20 |
30 |
40 |
50 |
|
|
|
|
|
|
|
ex:5e |
4,853 |
4,844 |
4,812 |
4,792 |
4,750 |
4,688 |
D min(T (x); 5) |
0,444 |
0,496 |
0,563 |
0,651 |
0,771 |
0,943 |
|
|
|
|
|
|
|
ex:10e |
9,445 |
9,375 |
9,286 |
9,167 |
9,000 |
8,750 |
D min(T (x); 10) |
2,777 |
3,836 |
4,252 |
4,861 |
5,666 |
6,771 |
60 |
70 |
80 |
85 |
90 |
4,583 |
4,375 |
3,750 |
2,500 |
|
1,215 |
1,632 |
3,604 |
2,080 |
|
8,333 |
7,500 |
5,000 |
|
|
8,333 |
13,542 |
8,333 |
|
|
3.5Округленное остаточное время жизни, его распределение, среднее и дисперсия
Врассматриваютактуарной математикеего целую частьнаряду с остаточным временем жизни T (x)
ленной остаточной продолжительностьюK(x) = [T (xжизни)]; которую(curtateназывают-future-округlifetime)- .
Это связано со следующими причинами:
1) человек обычно считает свой возраст в целых годах;
2) договоры страхования жизни, как правило, заключаются на целое число лет;
31
3) в ТПЖ приводятся данные для возрастов в целых годах. образом,Если Tслучайная(x) = 10 летвеличина9 месяцев = 10; 75 ëåò, òî K(x) = 10 лет. Таким
личиной, принимающей целыеKзначения(x) является.Как известно,дискретнойисчерпывающейслучайной ве-
характеристикой такой случайной величины является набор вероятностей
Понятно, что
Òàê êàê T (x) непрерывная случайная величина, то
P(T (x) = k) = P(T (x) = k + 1) = 0;
P(K(x) = k) = P(k < T (x) k + 1) =
= s(x + k) s(x + k + 1) = kpx k+1px: s(x)
и округленногоТеперь, учитывая,временичтожизниX = T (0); можно определить распределение
K(0) = [X] :
P(K(0) = k) = |
s(k) s(k + 1) |
= s(k) |
|
s(k + 1) = |
lk lk+1 |
= |
dk |
: |
s(0) |
|
|
||||||
|
|
|
l |
l |
Íî òàê êàê dx l f(x); òî
P(K(0) = k) f(k); |
(3.5.1) |
стигде приближенногоf(x) плотностьравенстваслучайной(3.5.1),величинывообще говоря,X; причемправильнеев правойписатьча-
f(kПоэтому) 1 ãîä; такравенствокак в его3.левой5.1подразумевает,части находитсячто безразмерная величина.
P(K(0) = k) f(k) 1ãîä;
откуда видно, что кривая смертей тесно связана с распределением округленного времени жизни.
остаточнойСреднее случайнойпродолжительностьювеличины K(жизниx) называется,обозначаетсясредней округленной
ex = E K(x);
32
и согласно определению математического ожидания для дискретной случайной величины
1
X
ex = k P(K(x) = k):
k=1
|
Òàê êàê |
s(x + k) s(x + k + 1) |
|
|
P(K(x) = k) = |
||
è |
s(x) |
||
|
1
X
k[s(x + k) s(x + k + 1)] = 1 s(x + 1) + 2 s(x + 2) +
k=1
1
X
1 s(x + 2) 2 s(x + 3) = s(x + k);
òî |
|
|
|
|
k=1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
1 |
|
|
|
|
|
|
X |
|
|
ex = |
s(x) |
|
s(x + k): |
|
Аналогично находится |
|
|
|
k=1 |
момент: |
|
|||||
второй начальный |
|||||
1 |
|
|
|
||
|
|
X |
k2 P(K(x) = k) = |
||
E [K(x)]2 = |
k=0
= 1 s(x)
1
( 1 |
1 |
) |
XX
|
k2 s(x + k) k2 s(x + k + 1) = |
k=0 |
k=0 |
( 1
= 1 Xk2 s(x + k) s(x)
k=0
1)
X |
|
X |
|
|
||
(k + 1)2 s(x + k + 1) + (2k + 1) s(x + k + 1) = |
||||||
k=0 |
|
k=0 |
|
|
||
1 |
1 |
2 |
|
1 |
||
|
|
X |
|
|
X |
|
= |
s(x) |
(2k 1) s(x + k) = |
s(x) |
|
k s(x + k) ex; |
|
òàê êàê |
|
k=1 |
|
|
k=1 |
|
|
|
|
|
|||
1 |
|
1 |
|
|
|
1 |
X |
X |
|
|
|
X |
|
k2 s(x + k) (k + 1)2 s(x + k + 1) + |
(2k + 1) s(x + k + 1) = |
|||||
k=0 |
k=0 |
|
|
|
k=0 |
|
|
|
|
33 |
|
|
|
=12 s(x + 1) + 22 s(x + 2) + 32 s(x + 3) +
12 s(x + 1) 22 s(x + 2) 32 s(x + 3)
1
X
+ 1 s(x + 1) + 3 s(x + 2) + 5 s(x + 3) + = (2k 1)s(x + k):
Теперь легко получаем дисперсию |
k=1 |
||
|
|||
|
|
|
|
|
2 |
1 |
|
|
|
|
X |
|
D K(x) = E [K(x)]2 ex2 = |
s(x) |
k s(x + k) ex ex2 : |
|
|
|
k=1 |
Все числовые характеристики округленного времени жизни выражаются через функции выживания, поэтому, в силу равенства (2.1.2), их также можно находить и по данным ТПЖ. В частности,
|
|
|
|
|
|
|
|
|
1 |
1 |
|
|
1 |
|
1 |
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
X |
|
|
|
|
|
|
X |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
e = l |
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
ex = lx k=1 lx+k; |
|
|
|
k=1 lk; |
|
|
|
||||||||||||
|
|
2 |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
2 |
1 |
|||||||||
|
|
|
|
|
|
|
|
X |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
X |
||
E [K(x)]2 = lx |
k lx+k ex; |
|
E [K(0)]2 |
= l |
|
|||||||||||||||||||||
k=1 |
|
|
|
|
k lk e : |
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
k=1 |
||
Пример 3.5.1. По данным ТПЖ (Приложение 1) отдельно для |
||||||||||||||||||||||||||
мужчин и женщин найдите: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
2)1) дисперсиисредниеокругленныеокругленныхпродолжительностипродолжительностейжизнижизниe89; e88; e84; |
||||||||||||||||||||||||||
Р е ш е н и е. Для мужчин: |
|
|
|
|
|
|
|
|
|
|
|
|
|
DK(89); DK(88): |
||||||||||||
e |
|
= |
290 |
|
= 0; 2; |
|
DK(89) = |
2 290 |
|
|
|
0; 2 |
|
|
|
(0; 2)2 = 0; 16; |
||||||||||
89 |
|
|
|
1449 |
|
|||||||||||||||||||||
|
1449 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
e88 = |
1449 + 290 |
|
= 0; 48; |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
3623 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
DK(88) = |
2(1 1449 + 2 290) |
|
0; 48 |
|
(0; 48)2 = 0; 41; |
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
3623 |
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
e84 = |
9063 + 7546 + 6037 + 3623 + 1449 + 290 |
= 2; 6: |
||||||||||||||||||||||
|
|
|
||||||||||||||||||||||||
Для женщин: |
|
|
|
10735 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
e |
|
= |
806 |
= 0; 2; |
|
DK(89) = |
2 806 |
|
|
|
0; 2 |
|
|
|
(0; 2)2 = 0; 16; |
|||||||||||
89 |
|
|
4030 |
|
||||||||||||||||||||||
|
4030 |
|
|
|
|
|
|
|
|
|
|
|
|
|
34
e88 = 4030 + 806 = 0; 48; 10075
DK(88) = 2(1 4030 + 2 806) 0; 48 (0; 48)2 = 0; 41; 10075
e84 = 24265 + 20988 + 16791 + 10075 + 4030 + 806 = 2; 75: 27665
35
3.6 Контрольные задания
Задание 3.6.1
Предположим, что кривая смертей описывается формулой
|
x |
f(x) = |
a2 e x=a; x 0: |
ни 1) Найдите функцию распределения Fx(t) остаточного времени жиз- T2)(xПокажите,) = X x:что плотность распределения остаточного времени жиз-
d
íè fx(t) = dtFx(t) является взвешенной суммой экспоненциальной плот- ности a1 e t=a и эрланговской плотности at2 e t=a:
3) Найдите вероятность P (T (x) > t); предел lim P (T (x) > t) и выясните, можно ли применять такую аппроксимациюx!1кривой смертей для больших возрастов x:
Задание 3.6.2
Используя ТПЖ Приложения 1 оцените вероятности того, что индивидуум (21):
1) доживет до 70, 80, 90 лет;
2) умрет до 70, 80, 90 лет; 3) умрет от 60 до 70, от 70 до 80, от 80 до 90 лет.
Задание 3.6.3
Раскройте смысл следующих обозначений актуарной математики:
q60 q80 p21; 5p21; q21; 5q21; 1jq21; 3jq25; 3j4q29; q20 ; q20 :
Используя ТПЖ Приложения 1 оцените приведенные выше величины.
Доказать, что |
Задание 3.6.4 |
|
tjuqx = tpx uqx+t: |
||
|
36
Задание 3.6.5
Используя выбранную из Задания 2.18.1 какую-либо таблицу для функции выживания определить вероятность того, что остаточное время70 дожизни80 лет(20). лежит в промежутке от 30 до 40 лет, от 40 до 50 лет, от
Задание 3.6.6
Найдите дисперсию частичной продолжительности жизни в модели де Муавра, постройте ее графики в зависимости от возраста страхующегося для n = 5 è n = 10 ëåò.
|
|
|
|
Задание 3.6.7 |
|
Пусть функции выживания задаются формуламè: |
|||||
s1 |
(x) = a |
e x=a; |
x 0; s2(x) = r1 110; 0 x 110: |
||
|
|
x + a |
|
x |
Найдите:
1) среднюю продолжительность жизни
e
2) полную вероятную продолжительность ;жизни
e
3) дисперсию среднего остаточного времени жизниx;
DT (x);
4)частичную среднюю продолжительность жизни
5)дисперсию частичной продолжительности жизниex:ne;
Dfmin(ИзобразитеT (x); n)графическиg. полученные зависимости, проведите сравнительный анализ.
Задание 3.6.8
По данным ТПЖ (Приложение 1) отдельно для мужчин и женщин найдите:
1) средние округленные продолжительности жизни e ; e10; e21; e40; e70;2) дисперсии округленных продолжительностей жизни
DK(0); DK(21);
DKОбсудит(70); DK(84):
37
4 ДРОБНАЯ ПРОДОЛЖИТЕЛЬНОСТЬ ЖИЗНИ
4.1Сплайновые аппроксимации для дробных возрастов (fractional ages)
(в годах),Реальнаячтостатистикаобусловленодоступнакак определеннымиобычно толькотрадициямидля целыхизначенийудобствомx
сбора статистических данных, так и формой их представления в ТПЖ, большинствогде аргументыклиентовx; как правило,в свой деньпринимаютрождениязначенияв страховую0; 1; 2компанию; : : : : Так какне
приходят, то для работы с конкретными индивидуумами нужно уметь находить приближения некоторых вероятностных характеристик для дробныхРассмотреннаявозрастовзадачапоих известныместьтипичнаязначениямзадачадляинтерполяции,целых x: причем
такможнокакограничитьсядругие величиныее решениемможно выразитьтолькодлячерезфункции выживания s(x), В актуарной математике эту задачу решают сs(помощьюx). òàê íàçû-
ваемых сплайнов. Рассмотрим три постулата, дающие различные приближения:
1. равномерное распределение смертей,
2. постоянная интенсивность смертности,
3. предположение Балдуччи (Balducci).
Равномерное распределение смертей
В этом случае функция выживания интерполируется линейной функ- s x) = a +èçbÒÏÆ),x ïðè nсоставляемx n +уравнения1. Òàê êàê s(n) è s(n + 1)
известныцией вида(например,n n
an + bnn = s(n)
an + bn(n + 1) = s(n + 1)
и определяем неизвестные |
an |
è |
bn (вычитаем из второго уравнения пер- |
|
вое уравнение): |
|
|
||
|
bn = s(n + 1) s(n); |
|||
|
an = s(n)(n + 1) s(n + 1)n: |
|||
Отметим,Следовательно,что |
на отрезке |
|
|
|
bn < 0: |
|
|
|
руется линейным сплайном видаn x n + 1 функция s(x) аппроксими-
s(x) = (n + 1 x)s(n) + (x n)s(n + 1); |
(4.1.1) |
38 |
|
nностиx n + 1: Отсюда для кривой смертей f(x) и интенсивности смерт-x получаем соответственно:
|
0 |
|
|
|
|
|
|
||
f(x) = s (x) = s(n) s(n + 1); n < x < n + 1; |
|
|
|||||||
x = |
f(x) |
= |
|
s(n) s(n 1) |
|
= |
|
|
|
|
|
|
|
|
|
|
|||
|
s(x) (n + 1 |
x)s(n) + (x n)s(n + 1) |
|
|
|||||
= |
|
|
|
s(n) s(n + 1) |
|
: |
(4.1.2) |
||
Таким образом,(n + |
1)s(n) ns(n + 1) x[s(n) s( |
n + 1)] |
|
|
|||||
n < x < n + 1: |
|
f(x) = bn = s(n) s(n + 1); |
|
|
|
|
|||
|
|
|
|
|
s(n) s(n + 1) |
|
|||
С помощью ранее введенной величины qn = |
|
||||||||
ной вероятности того, что человек в возрасте |
s(n) |
, ðàâ- |
ближайшего года, преобразуем формулу (4.1.2) кnболеелетумретудобномув течениивиду:
x = |
s(n) s(n + 1) |
= |
|
qn |
; |
s(n) + (x n)(s(n + 1) s(n)) |
|
(x n)qn |
|||
|
1 |
|
интенсивностиn < x < n + 1:смертностиВидим, что такое приближение влечет возрастание
|
|
|
(4.1.3) |
x = |
qn |
||
|
|
||
|
1 (x n)qn |
|
|
между узлами интерполяции (n < x < n+1), а плотность распределения
f(x) = |
const не меняется, причем в целочисленных точках |
f(x) |
è |
x |
íå |
|||||||||
определены. |
|
|
|
|
|
|
|
|
|
|
|
|
||
Заметим, что qn > s(n) s(n + 1), òàê êàê s(n) < 1: |
|
|
|
|
|
|||||||||
Постоянная интенсивность смертности |
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
ane bnx. Ïðè ýòîì |
x |
n + 1 убывающей |
||||||||
показательнойПриблизимфункциейфункцию s(x) на отрезке n |
|
|
||||||||||||
|
|
|
|
|
уравнения принимают вид |
|||||||||
|
|
|
|
ane bnn = s(n) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ane bn(n+1) = s(n + 1); |
|
|
|
|
|
|
|
|
|
|
и, разделив второе из них на первое, получаем |
|
|
|
|
|
|
|
|
|
|||||
|
|
n + 1) |
|
|
|
s(n + 1) |
|
n |
|
|
|
|||
|
bn = ln |
s( |
|
; an = s(n)ebnn = s(n) |
|
|
; |
|
|
|
||||
|
s(n) |
s(n) |
|
|
|
|||||||||
|
|
|
|
39 |
|
|
|
|
|
|
|
|
|
|
ò.å. |
b |
|
= |
|
ln p |
; a = s(n)p n; |
|
|
n |
||||||
|
|
|
n |
n |
n |
гдекрайнейp естьмеревероятностьеще один годтого,. что человек в возрасте n лет проживет по
n
В этом случае
s(x) = a |
e bnx = s(n)p ne(ln pn)x = s(n)px n; n |
|
x |
|
n + 1; |
|
n |
n |
n |
|
|
и основные вероятностные характеристики приближенно выражаются следующим образом:
|
|
|
s(x) = s(n)px n; |
n |
|
x |
|
n + 1; |
|
|
|
|
|||||
|
|
|
|
|
n |
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
f(x) |
|
|
|
|
|
f(x) = |
|
s0(x) = |
|
s(n)px n ln p |
; |
|
= |
= |
|
ln p |
; |
||||||
|
|
|
|||||||||||||||
|
|
|
n |
n |
|
|
|
x |
|
s(x) |
|
n |
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x = const.
ln pn > s(n) ln pn, òàê êàê s(n) < 1.
Предположение Балдуччи
В этом случае линейной функцией вместо s(x) интерполируется
s 1(x): Формально заменяя в (4.1.1) s(x) дим к соотношению
|
1 |
|
= |
n + 1 x |
+ |
x n |
; |
|
s(x) |
s(n) |
s(n + 1) |
||||||
|
|
|
||||||
откуда |
|
|
|
|
|
|
|
íà
1
s(x), немедленно прихо-
n x n + 1;
s(x) = |
|
|
|
|
s(n)s(n + 1) |
|
|
= |
|
|
s(n + 1) |
|
= |
||||
(n + 1 x)s(n + 1) + (x n)s(n) |
|
(n + 1 x)pn + x n |
|||||||||||||||
|
|
|
|
|
|
s(n + 1) |
|
|
|
|
|
|
|
|
|||
|
|
|
= |
|
; n x n + 1; |
|
|
||||||||||
|
|
|
pn + (x n)qn |
|
|
||||||||||||
f(x) = |
|
s0 |
(x) = |
s(n)s(n + 1)(s(n) s(n + 1)) |
= |
|
|||||||||||
|
|
|
|
|
|
(s(n + 1)(n + 1 |
|
x) + (x |
|
n)s(n))2 |
|
|
|||||
|
|
|
|
|
|
s(n + 1)qn |
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
= |
|
|
; n < x < n + 1; |
|
|
||||||||||
|
|
|
(pn + (x n)qn)2 |
|
|
40