lifetablegen

Сгенерируйте ряд таблицы продолжительности жизни из калиброванной модели выхода из строя

Описание

пример

[qx,lx,dx] = lifetablegen(x,a) генерирует ряд таблицы продолжительности жизни из калиброванной модели выхода из строя.

пример

[qx,lx,dx] = lifetablegen(x,a,lifemodel) генерирует ряд таблицы продолжительности жизни из калиброванной модели выхода из строя использование дополнительного аргумента для lifemodel.

Примеры

свернуть все

Загрузите файл данных таблицы продолжительности жизни.

load us_lifetable_2009

Калибруйте таблицу продолжительности жизни из данных о выживании с помощью heligman-pollard по умолчанию параметрическая модель.

a = lifetablefit(x, lx)
a = 8×3

    0.0005    0.0006    0.0004
    0.0592    0.0819    0.0192
    0.1452    0.1626    0.1048
    0.0007    0.0011    0.0007
    6.2853    6.7641    1.1037
   24.1386   24.2894   53.1848
    0.0000    0.0000    0.0000
    1.0971    1.0987    1.1100

Сгенерируйте ряд таблицы продолжительности жизни из калиброванной модели выхода из строя.

qx = lifetablegen(x,a);
display(qx(1:20,:))
    0.0063    0.0069    0.0057
    0.0005    0.0006    0.0004
    0.0002    0.0003    0.0002
    0.0002    0.0002    0.0002
    0.0001    0.0001    0.0001
    0.0001    0.0001    0.0001
    0.0001    0.0001    0.0001
    0.0001    0.0001    0.0001
    0.0001    0.0001    0.0001
    0.0001    0.0001    0.0001
    0.0001    0.0001    0.0001
    0.0001    0.0001    0.0001
    0.0002    0.0002    0.0001
    0.0002    0.0002    0.0002
    0.0002    0.0003    0.0002
    0.0003    0.0004    0.0002
    0.0004    0.0005    0.0002
    0.0005    0.0006    0.0003
    0.0006    0.0008    0.0003
    0.0007    0.0009    0.0003

Постройте qx ряд и отображение легенда. Серия qx условная вероятность, что человек в возрасте x умрет между возрастом x и следующим возрастом в ряду.

plot(x,log(qx))
legend(series)

Загрузите файл данных таблицы продолжительности жизни.

load us_lifetable_2009

Преобразуйте ряд таблицы продолжительности жизни в таблицы продолжительности жизни с принудительным завершением.

[~, lx] = lifetableconv(x, qx, 'qx');

Калибруйте таблицу продолжительности жизни из данных о выживании с помощью heligman-pollard по умолчанию параметрическая модель.

a = lifetablefit(x, lx)
a = 8×3

    0.0005    0.0006    0.0004
    0.0592    0.0819    0.0192
    0.1452    0.1626    0.1048
    0.0007    0.0011    0.0007
    6.2855    6.7638    1.1158
   24.1384   24.2895   52.6087
    0.0000    0.0000    0.0000
    1.0971    1.0987    1.1099

Сгенерируйте ряд таблицы продолжительности жизни из калиброванной модели выхода из строя.

qx = lifetablegen((0:100), a)
qx = 101×3

    0.0063    0.0069    0.0057
    0.0005    0.0006    0.0004
    0.0002    0.0003    0.0002
    0.0002    0.0002    0.0002
    0.0001    0.0001    0.0001
    0.0001    0.0001    0.0001
    0.0001    0.0001    0.0001
    0.0001    0.0001    0.0001
    0.0001    0.0001    0.0001
    0.0001    0.0001    0.0001
      ⋮

Постройте qx ряд и отображение легенда. Серия qx условная вероятность, что человек в возрасте x умрет между возрастом x и следующим возрастом в ряду.

plot((0:100), log(qx));
legend(series, 'location', 'southeast');
title('Conditional Probability of Dying within One Year of Current Age');
xlabel('Age');
ylabel('Log Probability');

Входные параметры

свернуть все

Увеличение возрастов для необработанных данных в виде N вектор из неотрицательных целочисленных значений. Возрасты должны запуститься в 0 (рождение).

Типы данных: double

Параметры модели для num модели в виде numparam- num матрица, где количество параметров (numparam) зависит от модели, заданной с помощью lifemodel аргумент.

Типы данных: double

(Необязательно) Параметрический тип модели выхода из строя в виде вектора символов с одним из следующих значений:

  • 'heligman-pollard' — Модель Хелигмен-Полларда с восемью параметрами (версия 1), заданная в терминах дискретной функции опасности:

    q(x)1q(x)=A(x+B)C+Dexp(E(logxF)2)+GHX

    для ages x0, параметрами A, B, C, D, E, F, G, H0.

  • 'heligman-pollard-2' — Модель Хелигмен-Полларда с восемью параметрами (версия 2), заданная в терминах дискретной функции опасности:

    q(x)1q(x)=A(x+B)C+Dexp(E(logxF)2)+GHX1+GHX

    для ages x0, параметрами A, B, C, D, E, F, G, H0.

  • 'heligman-pollard-3' — Модель Хелигмен-Полларда с восемью параметрами (версия 3), заданная в терминах дискретной функции опасности:

    q(x)=A(x+B)C+Dexp(E(logxF)2)+GHX

    для ages x0, параметрами A, B, C, D, E, F, G, H0.

  • 'gompertz' — Модель Two-parameter Gompertz, заданная в терминах непрерывной функции опасности:

    h(x) = A exp(Bx)
    для ages x0, параметрами A, B0.

  • 'makeham' — Модель Gompertz-Makeham с тремя параметрами, заданная в терминах непрерывной функции опасности:

    h(x) = A exp(Bx) + C
    для ages x0, с parametersA, B, C0.

  • 'siler' — Модель Siler с пятью параметрами, заданная в терминах непрерывной функции опасности:

    h(x) = A exp(Bx) + C + D exp(-Ex)
    для ages x0, параметрами A, B, C, D, E0.

Типы данных: char

Выходные аргументы

свернуть все

Условные вероятности смерти за N возрасты и num ряд, возвращенный как N- num матрица. Серия qx условная вероятность, что человек в возрасте x умрет между возрастом x и следующим возрастом в ряду. Для последнего возраста, qx представляет вероятности или значит все возрасты после последнего возраста.

Последняя строка N- num выведите для qx значения для всех возрастов на или после последнего возраста в x (из-за Принудительного Завершения). Поэтому последняя строка qx содержит 1 (100%-я вероятность смерти на или после последнего возраста).

Выживание значит N возрасты и num ряд, возвращенный как N- num матрица. Серия lx количество людей, живых в возрасте x, учитывая 100 000 живых при рождении.

Декрементные счета для N возрасты и num ряд, возвращенный как N- num матрица. Серия dx количество людей из 100 000 живых при рождении, которые умирают между возрастом x и следующий возраст в ряду. Для последнего возраста, dx представляйте вероятности или счета для всех возрастов после последнего возраста.

Последняя строка N- num выведите для dx значения для всех возрастов на или после последнего возраста в x (из-за Принудительного Завершения). Поэтому последняя строка dx содержит остающееся количество 100 000 человек, живых при рождении, кто не умер последним возрастом.

Больше о

свернуть все

Принудительное завершение

Большинство современных таблиц продолжительности жизни “обеспечило” завершение. Принудительное завершение означает, что последняя строка таблицы продолжительности жизни запрашивает всех людей с возрастами на или после последнего возраста в таблице продолжительности жизни.

Эта выборка иллюстрирует обеспеченное завершение.

В этом случае последняя строка таблицы продолжительности жизни запрашивает всех людей в возрасте 100 или более старый. А именно, вероятности qx 1qx целую вечность меньше чем 100 и, технически, ∞qx для возраста 100.

Принудительное завершение имеет терминальные значения возраста, которые применяются ко всем возрастам после терминального возраста так, чтобы lx положительно, qx 1, и dx положительно. Возрастами после терминального возраста является NaN значения, несмотря на то, что lx и dx может быть 0 и qx может быть 1 для входного ряда. Принудительное завершение инициировано естественно завершающим работу рядом, последним возрастом в усеченном ряду или первым NaN значение в ряду.

Ссылки

[1] Арии, E. “Таблицы продолжительности жизни Соединенных Штатов”. Национальные жизненные отчеты статистики, американское министерство здравоохранения и социального обеспечения. Издание 62, № 7, 2009.

[2] Carriere, F. “Параметрические Модели для Таблиц продолжительности жизни”. Транзакции Общества Актуариев. Издание 44, 1992, стр 77–99.

[3] Gompertz, B. “По Природе Функции, Выразительной из Закона Человеческого Выхода из строя, и на Новом Режиме Определения Значения Жизненных Непредвиденных обстоятельств”. Философские Транзакции Королевского общества. Издание 115, 1825, стр 513–582.

[4] Хелигмен, L. M. A. и Дж. Х. Поллард. “Шаблон Возраста Выхода из строя”. Журнал Института Издания 107 Актуариев, Pt. 1, 1980, стр 49–80.

[5] Makeham, W. M. “На Законе Выхода из строя и Конструкции Аннуитетных Таблиц”. Журнал Института Издания 8, 1860 Актуариев, стр 301–310.

[6] Siler, W. “Модель Конкурирующего Риска для Выхода из строя Животных”. Издание 60 экологии, стр 750–757, 1979.

[7] Siler, W. “Параметры Выхода из строя в Народонаселении с Широко Различными Продолжительностями жизни”. Статистика в Издании 2, 1983 Медицины, стр 373–380.

Представленный в R2015a
Для просмотра документации необходимо авторизоваться на сайте