exponenta event banner

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.2843    6.7637    1.1038
   24.1387   24.2895   53.1783
    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)

Figure contains an axes. The axes contains 3 objects of type line. These objects represent All, Male, Female.

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

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.2847    6.7634    1.1036
   24.1385   24.2898   53.1918
    0.0000    0.0000    0.0000
    1.0971    1.0987    1.1100

Создайте серию таблиц жизни из калиброванной модели смертности.

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');

Figure contains an axes. The axes with title Conditional Probability of Dying within One Year of Current Age contains 3 objects of type line. These objects represent All, Male, Female.

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

свернуть все

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

    h(x) = A exp(Bx) + C
    для возрастов x0, с параметрами A, B, C0.

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

    h(x) = A exp(Bx) + C + D exp(-Ex)
    для возрастов 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] Ариас, Э. «Таблицы жизни Соединенных Штатов». Национальные отчеты по статистике естественного движения населения, Министерство здравоохранения и социальных служб США. Том 62, № 7, 2009.

[2] Кэрриер, Ф. «Параметрические модели для жизненных таблиц». Операции Общества актуариев. Том 44, 1992, стр. 77-99.

[3] Гомперц, Б. «О характере функции, выражающей закон о смертности человека, и о новом способе определения ценности жизненных непредвиденных обстоятельств». Философские сделки Королевского общества. Т. 115, 1825, стр. 513-582.

[4] Хелигмен, L. M. A. и Дж. Х. Поллард. «Возрастная структура смертности.» Журнал Института актуариев том 107, стр. 1, 1980, стр. 49-80.

[5] Makeham, W. M. «О законе смертности и построении таблиц аннуитета». Журнал Института актуариев т. 8, 1860, с. 301-310.

[6] Силер, В. «Модель конкурирующего риска для смертности животных». Экология том 60, стр. 750-757, 1979.

[7] Сайлер, W. «Параметры смертности в популяциях человека с широко изменяющейся продолжительностью жизни». Статистика в медицине Том 2, 1983, стр. 373-380.

Представлен в R2015a