R語言人口期望壽命統(tǒng)計(jì)預(yù)測方法
原文鏈接:http://tecdat.cn/?p=18079
原文出處:拓端數(shù)據(jù)部落公眾號(hào)
?
本文對(duì)人口統(tǒng)計(jì)預(yù)測方法進(jìn)行討論。首先,我們將看到基本的靜態(tài)方法。在使用數(shù)據(jù)集之前,我們使用“標(biāo)準(zhǔn)”?生命表。
download.file(url,"mortal")
tables=readHTMLTable("morta
下面是獲取數(shù)據(jù)集的代碼
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))
可以使用生存函數(shù)來計(jì)算出生時(shí)的期望壽命
sum(TV0$lx)/100000-1
[1] 72.01518
可視化的生存概率
?

?
或死亡概率,即假設(shè)您達(dá)到xx歲,則在某特定年齡xx死亡的概率,也稱為?死亡率
n=nrow(TV0)
px=(TV8$lx[1:(n-1)]-TV8$lx[2:n])/
TV8$lx[1:
?

使用對(duì)數(shù)概率的可視化
plot(x,px,type="l",log="y")
?

?
最后,我們可以計(jì)算死亡年齡的概率密度
?

?
用來計(jì)算期望壽命
sum(x*pbx)
[1] 72.01518
對(duì)于初始數(shù)據(jù),我們可以使用人類死亡率數(shù)據(jù)庫中的表格。
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
?

專欄
精算科學(xué)
關(guān)于結(jié)合數(shù)學(xué)、統(tǒng)計(jì)方法以及程序語言對(duì)經(jīng)濟(jì)活動(dòng)來做風(fēng)險(xiǎn)分析、評(píng)估的見解。
探索專欄??