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