lifetableconv

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

Описание

пример

[qx,lx,dx] = lifetableconv(x0,lx0) преобразует таблицу продолжительности жизни с возрастами x0 и выживание считает lx0 в таблицы продолжительности жизни с завершением.

пример

[qx,lx,dx] = lifetableconv(x0,y0,y0type) преобразует таблицу продолжительности жизни с возрастами x0 и серия y0, заданный дополнительным аргументом y0type, в таблицы продолжительности жизни с завершением.

Примеры

свернуть все

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

load us_lifetable_2009

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

[qx,lx,dx] = lifetableconv(x,lx);
display(qx(1:20,:))
    0.0064    0.0070    0.0057
    0.0004    0.0004    0.0004
    0.0003    0.0003    0.0002
    0.0002    0.0002    0.0002
    0.0002    0.0002    0.0001
    0.0001    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.0002    0.0002    0.0002
    0.0003    0.0004    0.0002
    0.0004    0.0005    0.0002
    0.0005    0.0006    0.0003
    0.0005    0.0007    0.0003
    0.0006    0.0009    0.0004
    0.0007    0.0010    0.0004
display(lx(1:20,:))
   1.0e+05 *

    1.0000    1.0000    1.0000
    0.9936    0.9930    0.9943
    0.9932    0.9926    0.9939
    0.9930    0.9923    0.9937
    0.9927    0.9920    0.9935
    0.9926    0.9919    0.9933
    0.9924    0.9917    0.9932
    0.9923    0.9916    0.9931
    0.9922    0.9914    0.9930
    0.9921    0.9913    0.9929
    0.9920    0.9912    0.9928
    0.9919    0.9911    0.9927
    0.9918    0.9910    0.9926
    0.9917    0.9909    0.9925
    0.9915    0.9907    0.9923
    0.9912    0.9903    0.9921
    0.9908    0.9898    0.9919
    0.9904    0.9892    0.9916
    0.9899    0.9885    0.9913
    0.9892    0.9876    0.9909
display(dx(1:20,:))
  637.2266  698.8750  572.6328
   40.4062   43.9297   36.7188
   27.1875   30.0938   24.1406
   20.7656   23.0781   18.3359
   15.9141   17.2109   14.5625
   14.8672   16.3125   13.3516
   13.3672   14.7891   11.8750
   12.1328   13.3828   10.8203
   10.8125   11.6094    9.9844
    9.4609    9.5781    9.3438
    8.6172    8.1328    9.1172
    9.2656    8.8359    9.7188
   12.5938   13.5078   11.6328
   19.1016   22.9844   15.0234
   27.6719   35.5938   19.3516
   36.6328   48.5703   24.0547
   45.0156   60.7109   28.4844
   53.1406   72.8906   32.2812
   60.8984   85.1172   35.2578
   68.3438   97.2266   37.6875

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

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

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

load us_lifetable_2009

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

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((0:120),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

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

[~,~,dx] = lifetableconv((0:120),qx,'qx');
display(dx(1:20,:))
  630.9935  686.9471  571.6100
   48.7919   55.1033   40.9870
   24.8017   26.3778   23.6166
   17.0833   17.5877   17.0317
   13.6183   13.8188   13.6142
   11.8664   12.0076   11.6314
   10.9785   11.1573   10.4905
   10.5999   10.8605    9.9488
   10.5760   10.9396    9.8952
   10.8792   11.3613   10.2718
   11.6084   12.2508   11.0419
   12.9918   13.9270   12.1764
   15.3470   16.8833   13.6482
   18.9922   21.6788   15.4301
   24.1370   28.7659   17.4943
   30.7981   38.3208   19.8133
   38.7691   50.1484   22.3602
   47.6516   63.6906   25.1097
   56.9291   78.1271   28.0384
   66.0577   92.5264   31.1256

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

plot((0:119),dx(1:end-1,:));
legend(series, 'location', 'northwest');
title('\bfLife Table Yearly Decrements');
xlabel('Age');
ylabel('Number Dying within a Given Year');

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

свернуть все

Увеличение возрастов для необработанных данных в виде неотрицательных целочисленных значений в N0 вектор.

Вектор из возрастов x должен содержать неотрицательные целочисленные значения. Если входной ряд является дискретной функцией выживания lx, затем стартовый возраст должен только быть неотрицательным. В противном случае стартовым возрастом должен быть 0.

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

Набор num стандартизированный ряд оставшегося в живых в виде N0- num матрица. Вход lx0 ряд является количеством людей, живых в возрасте x, учитывая 100 000 живых при рождении. Значения 0 или NaN во входной таблице lx0 проигнорированы.

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

Набор num ряд таблицы продолжительности жизни, который будет преобразован в виде N0- num матрица. y0 по умолчанию рядом является lx0.

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

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

  • 'qx' — Вход является таблицей дискретных опасностей (qx).

  • 'lx' — Вход является таблицей дискретных количеств выживания (lx).

  • 'dx' — Вход является таблицей дискретных декрементов (dx).

Принимая во внимание, что выходные ряды обеспечили завершение, входной ряд (y0) может иметь один из нескольких типов завершения:

  • Естественное завершение заканчивается последнему человеку так, чтобы lx переходит к 0qx переходит к 1, и dx переходит к 0. Для получения дополнительной информации смотрите Естественное Завершение.

  • Усеченное завершение останавливается в терминальном возрасте так, чтобы lx положительно, qx  меньше 1, и dx положительно. Любыми возрастами после терминального возраста является NaN значения. Для получения дополнительной информации смотрите Усеченное Завершение.

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

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

свернуть все

Дискретная функция опасности, возвращенная как N0- num матрица с принудительным завершением. Для получения дополнительной информации смотрите Принудительное Завершение.

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

Дискретная функция выживания, возвращенная как N0- num матрица с принудительным завершением. Для получения дополнительной информации смотрите Принудительное Завершение.

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

Дискретная функция декрементов, возвращенная как N0- num матрица с принудительным завершением. Для получения дополнительной информации смотрите Принудительное Завершение.

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

Больше о

свернуть все

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

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

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

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

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

Естественное завершение

Прежде 1970, таблицы продолжительности жизни часто публиковались с данными, которые включали все возрасты, для которых люди, сопоставленные с данным рядом, были все еще живы. Таблицы в этой форме имеют "естественное" завершение. В естественном завершении последняя строка таблицы продолжительности жизни для каждого ряда считает смертельные случаи или вероятности смертельных случаев последнего остающегося человека в соответствующем возрасте. Таблицы в этой форме могут быть проблематичными из-за гранулярности данных и того, что группы рядов могут завершить работу в отличных возрастах. Естественное завершение проиллюстрировано в следующей выборке последних нескольких лет таблицы продолжительности жизни.

Эта форма для таблиц продолжительности жизни излагает много проблем, которые идут вне очевидных статистических проблем. Во-первых, таблица l x слева завершает работу в 108 лет, 109, 109, и 113 для четырех рядов в таблице. Технически, числами после этих возрастов является 0, но может также быть NaN значения, потому что никакой человек не жив после этих возрастов завершения. Во-вторых, вероятности q, который x справа отключает с колеблющимися вероятностями, которые идут от 0 к 1 в некоторых случаях. В этом случае, однако, все вероятности 1qx вероятности (различающиеся принудительные вероятности завершения). Можно утверждать, что вероятностями после возрастов завершения может быть 1 (любой живой в этом возрасте, как ожидают, умрет в следующем году), 0 (возраст находится вне поддержки вероятностного распределения), или NaN значения.

Усеченное завершение

Усеченное завершение происходит при усечении таблиц продолжительности жизни в произвольном возрасте. Например, от 1970–1990, таблицы продолжительности жизни Соединенных Штатов, усеченные в 85 лет. Этот формат проблематичен, потому что вероятности таблицы продолжительности жизни должны или отключить с вероятностью 1 (принудительное завершение) или данные об отбрасывании, которые превышают завершающий работу возраст. Эта выборка последних нескольких лет таблицы продолжительности жизни иллюстрирует усеченное завершение. Необработанные данные для этой таблицы являются серией l x. Серия q x выведена из этого ряда.

Этот формат таблицы продолжительности жизни создает проблемы для завершения, потому что, например, более чем 27% населения для четвертой серии l x все еще живы в 85 лет. Утверждать, что вероятность смерти за все возрасты после возраста 85 является 100%-й силой быть верной, но неинформативна. Несмотря на статистические проблемы, однако, таблицы в этой форме заполнены принудительным завершением.

Ссылки

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

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