Основы страховой (актуарной) математики - Г.М.Кошкин
.pdfс максимумом в точке x = (n=k)1=(n+1) (см. Задание 2.6.2).
Модель Эрланга
Рассмотрим модель Эрланга 2-го порядка, для которой кривая смертей описывается формулой
f(x) = |
x |
x |
; x 0: |
|
e a |
||
a2 |
В этом случае функция выживания
s(x) = x + a e xa ; a
а интенсивность смертности
x
x = a(x + a):
В заключение отметим, что определенным преимуществом аналити- ческих законов является то, что для них вероятностные характеристики продолжительности жизни можно быстро вычислять по небольшому числу параметров. Это может оказаться важным также и в случаях, когда доступные данные немногочисленны.
2.7 Контрольные задания
Задание 2.7.1
следующиеХарактертаблицы:зависимости функции выживания s(x) от возраста x äàþò
МУЖЧИНЫ (СССР, 1984 1985 г.) [4, с.86]
|
x |
0 |
|
14 |
20 |
|
|
30 |
40 |
50 |
60 |
|||||
|
s1(x) |
1,000 |
|
0,954 |
0,947 |
|
0,922 |
0,878 |
0,795 |
0,651 |
||||||
|
70 |
|
80 |
|
|
|
90 |
|
100 |
|
110 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
0,434 |
|
0,188 |
|
|
|
0,003 |
|
0 |
|
0 |
|
|
|
|
|
ЖЕНЩИНЫ (СССР, 1984 1985 г.) [4, с.86] |
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
||||||
|
x |
0 |
|
14 |
20 |
|
|
30 |
40 |
50 |
60 |
|||||
|
s2(x) |
1,000 |
|
0,964 |
0,961 |
|
0,953 |
0,939 |
0,908 |
0,841 |
||||||
|
70 |
|
80 |
|
90 |
|
100 |
|
110 |
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
||||||||
|
0,700 |
|
0,417 |
|
0,008 |
|
0 |
|
0 |
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
21 |
|
|
|
|
НАСЕЛЕНИЕ США [1, с.7]
|
x |
0 |
10 |
20 |
30 |
|
40 |
50 |
60 |
||||
|
s3(x) |
1,000 |
0,983 |
0,977 |
0,965 |
0,949 |
0,915 |
0,837 |
|||||
|
70 |
|
80 |
|
90 |
|
100 |
|
110 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
0,682 |
|
0,432 |
|
0,142 |
|
0,012 |
|
0 |
|
|
|
|
1) Исходя из смысла функции выживания проанализируйте данные |
|||
таблицы. |
|
|
|
2) Привлекая порядковые статистики, постройте параметрические |
|||
оценки функций выживания для модели де Муавра, воспользовавшись |
|||
приведенными таблицами. |
|
|
|
3) Опишите несколько способов построения параметрических оценок |
|||
функций выживания для модели Вейбулла на основе приведенных та- |
|||
блиц (в модели Вейбулла интенсивность смертности |
x приближается |
||
степенной функцией вида |
|
|
|
kxn). |
|
|
|
Задание 2.7.2 |
|
|
|
В модели Гомпертца интенсивность смертности |
x приближается по- |
||
казательной функцией вида |
|
||
ðû. |
Be x; ãäå > B > 0 некоторые парамет- |
1)Найти точку максимума для кривой смертей в модели Гомпертца.
2)Где можно использовать полученный результат?
3)Найти точку максимума для кривой смертей в модели Вейбулла.
|
Задание 2.7.3 |
|
|||
Пусть имеются две функции выживания: |
|
||||
s1(x) = e x3=12; x 0; |
|
||||
s2(x) = 1 |
x |
|
|
|
|
|
; 0 x < !; > 0: |
|
|||
! |
|
||||
интенсивностиприсмертностиописывает модель де Муавра. |
|||||
Функция1)Найтивыживания s2(x) |
= 1 |
|
|||
функции распределения |
|
|
|
ix; кривые смертей fi(x) è |
|
si(x ; Найтиi = 1; 2вероятности: |
Fi(x); соответствующие функциям выживания |
||||
того, что умрут в промежутке от |
|
||||
2) |
|
|
|
|
10 äî 30 ëåò: |
а) случайно выбранный человек, |
22
б) застрахованный в возрасте 10 ëåò.
Рассмотрим |
Задание 2.7.4 |
|||
f(x) = |
x |
e x=a: |
||
|
||||
|
|
|||
|
|
a2 |
||
2)1) ОпределитеПокажите,чтовидfсоответствующих(x) может рассматриватьсяфункции выживаниякак кривая смертей. |
||||
тенсивности смертности |
|
|
s(x) è èí- |
|
|
x; а также соответствующую среднюю продол- |
жительность жизни
3) Проанализируйтеe : степень соответствия предложенного аналити-
ческого описания реальным данным с помощью таблицы продолжительности жизни (сравните среднее время жизни с точкой максимума кривой смертей).
4) Найдите дисперсию продолжительности жизни.
Задание 2.7.5
Рассмотрим три таблицы значений функций выживания задания 2.18.1. Подсчитайте среднее и дисперсию числа представителей исходлетной. группы в l0 новорожденных, которые умрут в возрасте от 50 äî 70
Задание 2.7.6
Аппроксимируйте выбранную из задания 2.18.1 какую-либо функцию выживания si(x); i = 1; 2; 3; используя модель Вейбулла для целых
n =1)2;Найдите4: параметр
2) Постройте графикиkполученныхспомощью аппроксимацийМНК. при
ясните какая из построенных моделей наилучшая в смыслеnнекоторого= 2; 4; âû-
критерия, выбором которого распорядитесь самостоятельно.
3)4) МожноКакнайтили провестиоптимальноеодновременнозначение nоптимизацию? по параметрам
è |
k |
|
n? |
Задание 2.7.7
Предположим, что кривая смертей описывается формулой
f(x) = ax2 e x=a; x 0:
ни 1) Найдите функцию распределения Fx(t) остаточного времени жиз-
ственно),T (x = X x (X; x момент смерти и возраст индивида, соответ-
Fx(t) = P (T (x t)):
23
2) Покажите, что плотность распределения остаточного времени жиз-
d
íè fx(t) = dtFx(t) является взвешенной суммой экспоненциальной плот-
ности |
1 |
e t=a и эрланговской плотности |
t |
e t=a: |
|
2 |
|||
|
a |
a |
3) Найдите вероятность P (T (x) > t); предел lim P (T (x) > t) и выясните, можно ли применять такую аппроксимациюx!1кривой смертей для больших возрастов x:
24
3ОСТАТОЧНАЯ ПРОДОЛЖИТЕЛЬНОСТЬ ЖИЗНИ
3.1Остаточное время жизни (time-untill-death), его распределение
луЕслиспецификичеловек дожилеедеятельности,до возрастаинтересует,x лет, то страховуювообщеговоря,компанию,уже невегоси-
общая продолжительность жизни X; а остаточное время жизни T (x) = Xтораяx:являетсяНайдем условнойфункциюфункциейраспределенияраспределенияслучайнойвеличинывеличины T (x); êî-
условии, что |
|
|
|
|
|
|
|
|
X x ïðè |
X > x : |
|
|
|
|
|
|
|
||
Fx(t) = P(T (x) t) = P(X x tjX > x) = |
|
||||||||
|
|
j |
|
t |
x |
|
P(X >Tx) |
|
|
= P(X |
|
x + t |
X > x) = |
q |
|
= |
P(X x + t X > x) |
= |
|
|
|
|
|
||||||
|
= |
P(x < X x + t) |
= |
F (x + t) F (x) |
: |
(3.1.1) |
|||
|
|
P(X > x) |
|
|
|
1 F (x) |
|
распределенияЕсли известны табличные значения lx è lx+t; то через них функция Fx(t) выражается следующим образом:
Fx(t) = |
s(x) s(x + t) |
= |
lx=l lx+t=l |
= |
lx lx+t |
: |
s(x) |
|
|||||
|
|
lx=l |
|
lx |
Воспользовавшисьличины (3.1.1), найдем теперь плотность fx(t) случайной ве-
T (x) : |
|
|
|
|
|
|
d |
f(x + t) |
|
|
|
fx(t) = |
|
Fx(t) = |
|
; 0 |
t < 1: |
dt |
1 F (x) |
В качестве примеров рассмотрим законы де Муавра и Мэйкхама. Модель де Муавра. Введем для удобства обозначение
|
|
|
1; |
åñëè x (a; b] |
|
|
|
Ix(a; b] = 0; |
åñëè x 22= (a; b] |
: |
|
||
Пусть распределение продолжительности жизни |
описывается законом |
|||||
де Муавра, для которого плотность распределения и функция выжива- |
||||||
ния, соответственно, задаются формулами: |
|
|
||||
f(x) = |
Ix(0; !) |
; s(x) = Ix(1; !) |
x Ix(0; !) |
: |
||
! |
|
! |
||||
|
|
|
|
25 |
|
|
Òàê êàê 0 < X < !; òî 0 < T (x) < ! x; è
t
Fx(t) = ! xIt(0; ! x) + It[! x; 1);
à äëÿ t 2 (0; ! x)
fx(t) = |
d |
|
t |
|
= |
|
1 |
|
; t 2 (0; ! x]; |
(3.1.2) |
dt |
|
! x |
! x |
накот.е. остаточноена промежуткевремя жизни T (x) также равномерно распределено, одМодель Мэйкхама(0; !. Дляx):этого закона
t = A + B e t; |
s(x) = exp Z0x |
tdt = exp [ Ax B (e x 1)= ] ; |
|||
f(x) = s0 |
(x) = [A + B e x] exp [ Ax B (e x 1)= ] ; |
||||
f(x + t) = [A + B e (x+t)] exp h A(x + t) B (e (x+t) 1)= i: |
|||||
Поэтому плотность распределения остаточного времени жизни |
|||||
|
x + t) |
|
|||
fx(t) = |
f( |
|
|
= [A + B e xe t]exp At B e x(e t 1)= ; |
|
s(x) |
|
|
делениеоткуда следует,Мэйкхамачтосостаточноетакимипараметрами:время жизни T (x) также имеет распре-
Ax = A; Bx = B e x; x = :
3.2 Величины, связанные с T (x): tqx; tx; qx; px; tjuqx; tjqx
Вроятностьактуарной математике вероятность P(T (x) t) и дополнительная ве-
с (3.1.1) определяютсяP(T (x) > t);следующими формулами: q ; в соответствии обозначаемые символами t x è tx
q |
|
= P(T (x) |
|
t) = |
s(x) s(x + t) |
; |
|
|
|
s(x) |
|||||||
t |
x |
|
|
|
|
|||
tx = P(T (x) > t) = 1 tqx = |
s(x + t) |
|||||||
|
|
: |
||||||
s(x) |
|
(3.2.1)
(3.2.2)
â |
промежутке временивыражает вероятность смерти человека возраста |
x |
ëåò |
|
Величина tqx |
|
|
||
доживет до возраста |
(x; x + t]; à tx вероятность, что такой человек |
|||
|
|
x + t: |
|
|
|
|
26 |
|
|
формул t = |
1) |
|
tqx |
|
tx |
|
Ïðè (3.2. |
|
ипередние(3.2.2) получаеминдексы у |
|
è |
|
опускаются, и из общих |
qx = P(T (x) 1) = s(x) s(x + 1) = lx lx+1 ; s(x) lx
x = P(T (x) > 1) = s(x + 1) = lx+1 : s(x) lx
(3.2.3)
(3.2.4)
Переменные qx è x чаще всего используются на практике; величина
qближайшегоравна вероятностигода, а смерти индивидуума возраста x лет в течение
x
ìåðå îäèí ãîä. |
px вероятности, что он проживет еще по крайней |
||||
Учитывая формулу (3.2.3), представим |
|||||
|
|
|
|
|
|
|
|
qx = |
dx |
: |
|
|
|
|
|||
|
|
|
lx |
|
Величина q
äèò â ÒÏÆx является(Приложениетретьей1).основнойТаккак характеристикой, которая вхо-
lx = l s(x); dx l f(x); òî
qx |
f(x) |
= x; т.е. форма кривой qx приближенно совпадает с формой |
|||||||
s(x) |
|||||||||
кривой функции интенсивности смертности. |
|
|
|
||||||
Таким образом, мы выяснили смысл трех основных величин, содер- |
|||||||||
жащихся в ТПЖ: |
|
|
|
||||||
öèåél |
|
|
|
|
связана через мультипликативный множитель |
l |
|
ñ ôóíê- |
|
выживания |
|
|
|||||||
= l |
|
s(x) |
|
|
|
|
|||
x |
|
|
|
|
|
|
|
||
от возраста |
|
|
s(x) и повторяет ее характер изменения в зависимости |
||||||
|
|
|
|
|
|
|
|
|
|
tjuqx = P(t < T (x) t + u); |
(3.2.5) |
мая во внимание соотношения (3.2.1) и (3.2.2), черезq ;основную; или,функциюприни- и она может быть выражена или через функцию t x èëè tx
страховой математики функцию выживания s(x) :
tjuqx = P(t < T (x) t + u) =
27
= P(T (x) t + u) P(T (x) t) = t+uqx tqx; |
(3.2.6) |
|||
tjuqx = P(t < T (x) t + u) = |
|
|||
= P(T (x) > t) P(T (x) > t + u) = tpx t+upx; |
(3.2.7) |
|||
tju |
qx = |
s(x + t) s(x + t + u) |
: |
(3.2.8) |
|
s(x) |
|
нияСноважизни,случайи, как uобычно,= 1 наиболееэтот индексинтересенопускаетсядля. Согласнопрактикиформуламстрахова- (3.2.6), (3.2.7) и (3.2.8) имеем:
q |
|
= |
|
q |
q |
|
= |
p |
x t+1 |
p |
|
= |
s(x + t) s(x + t + 1) |
: (3.2.9) |
|
|
|
|
s(x) |
||||||||||
tj |
x |
|
t+1 |
|
x t |
x |
|
t |
|
x |
|
|
3.3 Среднее остаточное время жизни, его дисперсия. Коэффициент асимметрии и эксцесс
В актуарной математике среднее остаточного времени жизни человека возраста x обозначается
|
|
|
|
|
|
ex= ET (x); |
|
|
|
эту характеристику часто используют при анализе ситуации на рынке |
||||
|
|
|
|
; и, следовательно, сред- |
страховых услуг. Понятно, что ET (0) = EX =e |
|
|||
няя продолжительность жизни |
|
|
||
жизни |
x > 0: |
e больше среднего остаточного времени |
||
Òàêeкакдля любого |
|
|
|
|
x |
|
|
|
воспользовавшисьT (x) являетсятеоремойнеотрицательной2.4.1, получаем случайной величиной, то,
Z 1 Z 1
ex= ET (x) = t dFx(t) = t dP(T (x) t) =
00
|
= Z01 P(T (x) > t) dt = Z0 |
1 tpx dt: |
|
||||||||||
Применив формулу (3.2.2) и преобразовав последний интеграл |
|||||||||||||
1 |
1 |
|
|
1 |
|
|
1 |
1 |
|
||||
Z0 |
tpx dt = |
|
|
Z0 |
s(x + t) dt = |
|
|
|
Zx |
s(u) du; |
|||
s(x) |
s(x) |
||||||||||||
приходим к соотношению |
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
1 |
|
|
|
|
|
|
|
|
|
ex= |
|
Zx |
s(u) du: |
|
|
(3.3.1) |
||||
|
|
|
s(x) |
|
|||||||||
|
|
|
|
|
|
|
28 |
|
|
|
|
|
|
Рассуждая аналогично, для второго начального момента случайной величины T (x) имеем:
|
E[T (x)]2 = 2 Z01 t P(T (x) > t) dt = 2 Z01 t tpx dt = |
|
|
|||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
2 |
|
|
Z0 |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
= |
|
|
|
|
|
|
|
t s(x + t) dt; |
|
|
|
|
|
|
|
|
|
|||||||||||||
откуда |
|
|
|
|
|
|
|
|
|
s(x) |
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
2 |
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
DT (x) = |
|
|
|
|
Z0 |
|
t s(x + t) dt ex |
|
: |
|
|
|
|
|
|||||||||||||||||||||
|
|
|
s(x) |
|
|
|
|
|
|
|||||||||||||||||||||||||||||
Для нахождения коэффициента асимметрии и эксцесса необходимо |
||||||||||||||||||||||||||||||||||||||
определить E[T (x)]3 è E[T (x)]4 : |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
E[T (x)]3 = |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
3 |
|
|
1 |
|
|
|
|
|
|||
= 3 Z0 |
t2P(T (x) > t) dt = 3 Z0 |
t2tpx dt = |
|
|
Z0 |
|
t2s(x + t) dt; |
|||||||||||||||||||||||||||||||
s(x) |
||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
4 |
|
Z0 |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
E[T (x)]4 = |
|
|
|
t3s(x + t) dt: |
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
s(x) |
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
Пример 3.3.1. Найти среднее |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
дратическое отклонение |
|
|
|
|
|
|
|
|
|
|
ex; дисперсию DT (x) и среднее ква- |
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
(T (x)) = |
|
|
|
|
|
|
|
остаточного времени жизни |
||||||||||||||||||||||
T (xÐ) ве моделише н и дее. МуавраДляэтой. моделиpDT (x) |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||
P(T (x) > t) = |
p |
|
= |
|
s(x + t) |
|
= |
|
1 (x + t)=! |
; 0 |
|
t |
|
! |
|
x: |
||||||||||||||||||||||
x |
|
s(x) |
|
|
|
|
||||||||||||||||||||||||||||||||
Таким образом, |
|
t |
|
|
|
|
|
|
|
|
|
|
1 |
|
x=! |
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
e |
= |
! x |
! x t |
dt = |
|
1 |
|
0 |
|
|
u du = |
! x |
: |
|
|
|
(3.3.2) |
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||
x |
Z0 |
|
|
! x |
|
|
|
|
|
|
! x Z! x |
|
|
|
|
2 |
|
|
|
|
|
|||||||||||||||||
Проще формулу (3.3.2) получить из (3.3.1): |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||
|
e |
= |
! |
|
|
|
! |
! t |
dt = |
|
|
1 |
|
|
|
0 |
u du = |
|
! x |
: |
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||
|
x |
|
! x Zx |
! |
|
|
|
|
|
|
! x Z! x |
|
|
2 |
|
|
|
|
||||||||||||||||||||
Также нетрудно показать, что |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
DT (x) = |
(! x)2 |
: |
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
12 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
29 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3.4Смешанное страхование. Частичная остаточная продолжительность жизни
илиРассмотримсмешанноеn-летнеестрахование,страхованиесуть которогона дожитиесостоит(n-yearв том,endowmentчто клиентinsurance)
заключаетлибов моментдоговорсмертистрахованиязастрахованного,на n лет, иесливыплатаона наступилапроизводится:до окон-
чаниялибоn-влетнегоконце периода,
Момент страховойn-летнеговыплатыпериода,выражаетсяесли застрахованныйформулой minостался жив. называется частичной продолжительностью жизни , а соответствую(T (x); n) è-
щее математическое ожидание частичной средней продолжительностью жизни и обозначается
ex:ne:
ex:ne= E min(T (x); n):
Чтобы применить теорему 2.4.1, необходимо найти вероятность P(min(T (x); n) > t): Òàê êàê äëÿ t < n событие min(T (x); n) > t равносильно событию
T (x) > t; òî
P(min(T (x); n) > t) = |
t |
0x; |
åñëè |
|
|
p ; |
åñëè |
Привлекая формулу (3.2.2), получаем:
0 t < n t n:
n |
tpx dt = s(x) Z0 |
n |
s(x + t) dt = s(x) Zx |
x+n |
||
ex:ne= Z0 |
|
s(u) du: |
||||
|
1 |
|
|
1 |
|
|
Для дисперсии частичной продолжительности жизни имеем
n |
t tpx dt (ex:ne)2 = s(x) Z0 |
n |
|
D min(T (x); n) = 2 Z0 |
t s(x + t) dt (ex:ne)2: |
||
|
2 |
|
|
Пример 3.4.1. Найти среднее и дисперсию частичной продолжи- |
|||||||||||||
тельности жизни в модели де Муавра. |
|
|
|
|
|
|
|||||||
Р е ш е н и е. Прежде всего отметим, что x + n < !; поэтому |
|||||||||||||
|
|
|
|
|
x+n |
|
Zx |
x+n |
|
|
|
||
|
ex:ne= s(x) Zx |
s(u) du = 1 x=! |
(1 u=!) du = |
||||||||||
|
|
|
1 |
|
|
|
1 |
|
|
|
|
|
|
= |
1 |
|
x+n(! |
|
u) du = |
2n(! x) n2 |
|
= n |
|
n2 |
: (3.4.1) |
||
! x Zx |
|
|
|
2(! x) |
|||||||||
|
|
|
|
2(! x) |
|
|
|
|
|||||
Из формулы (3.4.1) следует, что |
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
30 |
|
|
|
|
|
|