Main Content

寿命表分析的案例研究

打开实时脚本

此示例说明如何使用寿命表的基本工作流。

加载寿命表数据文件。 

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 数列是一个年龄为 x 岁的人将在 x 岁与数列中下一个年龄之间去世的条件概率。 

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');

另请参阅

| |

相关主题