В этом примере показано, как использовать базовый рабочий процесс для таблиц жизненного цикла.
Загрузите файл данных таблицы жизненного цикла.
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