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