本文对人口统计预测方法进行讨论。首先,我们将看到基本的静态方法。在使用数据集之前,我们使用“标准” 生命表。
- 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:
 

使用对数概率的可视化
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
 

    









