R语言人口期望寿命统计预测方法

原文链接：http://tecdat.cn/?p=18079

本文对人口统计预测方法进行讨论。首先，我们将看到基本的静态方法。在使用数据集之前，我们使用“标准” 生命表。

download.file(url,"mortal")
tables=readHTMLTable("morta

下面是获取数据集的代码

tables\[\[2\]\]
a1=as.numeric(as.character(TV8\[,1\]))
a2=as.numeric(as.cha
TV0=data.frame(x=c(a1,a2),lx=as.numeric(c(b1,b2))

可以使用生存函数来计算出生时的期望寿命

sum(TV0$lx)/100000-1
\[1\] 72.01518

可视化的生存概率

或死亡概率，即假设您达到xx岁，则在某特定年龄xx死亡的概率，也称为 死亡率

n=nrow(TV0)
px=(TV8$lx\[1:(n-1)\]-TV8$lx\[2:n\])/
TV8$lx\[1:

---

点击标题查阅往期内容

R语言Lee-Carter模型对年死亡率建模预测预期寿命

左右滑动查看更多

01

02

03

04

使用对数概率的可视化

plot(x,px,type="l",log="y")

最后，我们可以计算死亡年龄的概率密度

用来计算期望寿命

sum(x*pbx)
\[1\] 72.01518

对于初始数据，我们可以使用人类死亡率数据库中的表格。

download.file(url,"E.txt")

以下代码可用于读取这些文件。

read.table("C.txt",skip = 3,header=TRUE)

      Year  Age Female  Male  Total
22195 2015  101 242.66 25.86 323.82
22196 2015  104 132.95 16.39 213.34
22197 2015  101 101.87 9.50 145.37
22198 2015  103  57.27  4.07  64.34
22199 2015  109  31.93  2.59  32.52
22200 2015 110+  33.03  1.61  33.64

点击文末“阅读原文”

获取全文完整资料。