В этом примере показано, как использовать базовый рабочий процесс для таблиц жизни.
Загрузите файл данных таблицы жизни.
load us_lifetable_2009Калибровка таблицы жизни из данных о выживании с помощью heligman-pollard по умолчанию параметрическая модель.
a = lifetablefit(x, lx);
Сгенерируйте ряд таблиц жизни из калиброванной модели смертности.
qx = lifetablegen((0:100), a); display(qx(1:40,:))
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
0.0008 0.0011 0.0003
0.0008 0.0012 0.0004
0.0009 0.0013 0.0004
0.0009 0.0014 0.0005
0.0010 0.0014 0.0005
0.0010 0.0015 0.0005
0.0010 0.0015 0.0006
0.0010 0.0015 0.0006
0.0010 0.0015 0.0007
0.0010 0.0014 0.0007
0.0011 0.0014 0.0007
0.0011 0.0014 0.0008
0.0011 0.0014 0.0008
0.0011 0.0014 0.0009
0.0011 0.0014 0.0009
0.0012 0.0015 0.0010
0.0012 0.0015 0.0011
0.0013 0.0016 0.0011
0.0014 0.0017 0.0012
0.0015 0.0018 0.0013
Постройте график qx серия и отображение легенды. Последовательность qx - условная вероятность того, что человек в возрасте умрет между возрастом и следующий возраст в серии
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');
lifetableconv | lifetablefit | lifetablegen