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