R語言Lee-Carter模型對年死亡率建模預(yù)測期望壽命
原文鏈接:http://tecdat.cn/?p=17347
原文出處:拓端數(shù)據(jù)部落公眾號
?
昨天上午,我們獲得了分娩產(chǎn)婦的平均年齡兩個圖表,根據(jù)孩子的出生順序排序,區(qū)間是1905-1965年:

然后是1960-2000年:

?
這些圖形令人興奮,特別是在過去30年中觀察到的增長方面,這使我想到了壽命的增長趨勢。我們可以找到其他有趣的數(shù)據(jù)(在這種情況下為平均?出生年齡?)。
> age$Age=as.character(age$AGE)
> age$AGE=as.numeric(substr(age$Age,1,2))+
+ as.numeric(substr(age$Age,4,4))/10
> plot(age$ANNEE+.5,age$AGE,
+ type="l",lwd=2,col="blue")
我們在上面的圖中發(fā)現(xiàn)深藍色的曲線,
?
獲取祖母的平均年齡,我們進一步分析
> tail(age)
AGE ? Age NAIS.MERE NAIS.GRD.MERE age.GRD.MERE
2000 ?2000 30.3 30,3 ? ? 1970.2 ? ? ? 1942.87 ? ? ? ?57.63
2001 ?2001 30.4 30,4 ? ? 1971.1 ? ? ? 1943.80 ? ? ? ?57.70
2002 ?2002 30.4 30,4 ? ? 1972.1 ? ? ? 1944.92 ? ? ? ?57.58
2003 ?2003 30.5 30,5 ? ? 1973.0 ? ? ? 1945.95 ? ? ? ?57.55
2004 ?2004 30.5 30,5 ? ? 1974.0 ? ? ? 1947.05 ? ? ? ?57.45
2005 ?2005 30.6 30,6 ? ? 1974.9 ? ? ? 1948.04 ? ? ? ?57.46
> plot(age$ANNEE+.5,age$age.GRD.MERE,
+ type="l",lwd=2,col="red")
再一次,我們可以形象地看到外婆的出生年齡

?
我們可以通過使用Lee-Carter模型對年死亡率進行建模,并推斷到當(dāng)前世紀,我們可以推斷出期望剩余壽命。
> Deces <- read.table("Dec.txt",header=TRUE)
> Expo ?<- read.table("Expo.txt",header=TRUE,skip=2)
> Deces$Age <- as.nu
> Expo$Age <- as.numeric(as.character(Expo$Age))
> Expo$Age[is.n
Deces$Female/Expo$Female,nL,nC)
> ?POPF <- matrix(Expo$Female,nL,nC)
> ?BASEF <- demogdata(data=MUF, pop=POPF,ages=AGE,
+ years=YEAR, t
> K1 <- LCF$kt
nction(xentier,T){
+ return(ext) }
> EVIE = function(x,T){
+ x1 <- trunc(x)
> tail(age)
AGE ? Age NAIS.MERE NAIS.GRD.MERE age.GRD.MERE ? ? ? EV
2000 30.3 30,3 ? ? 1970.2 ? ? ? 1942.87 ? ? ? ?57.63 29.13876
2001 30.4 30,4 ? ? 1971.1 ? ? ? 1943.80 ? ? ? ?57.70 29.17047
2002 30.4 30,4 ? ? 1972.1 ? ? ? 1944.92 ? ? ? ?57.58 29.39027
2003 30.5 30,5 ? ? 1973.0 ? ? ? 1945.95 ? ? ? ?57.55 29.52041
2004 30.5 30,5 ? ? 1974.0 ? ? ? 1947.05 ? ? ? ?57.45 29.72511
2005 30.6 30,6 ? ? 1974.9 ? ? ? 1948.04 ? ? ? ?57.46 29.80398
換句話說,在最后一行,2005年,一名57.46歲女性的(剩余)期望壽命約為29.80歲。然后,我們不僅可以看到他祖母的平均年齡,還可以看到她的剩余期望壽命,

?
然后我們就可以確定曾祖母的(平均)年齡,
?

以及曾祖母的(剩余)壽命

?
現(xiàn)在我們也可以對這項快速研究的局限性感到疑惑。特別是,正如有配偶的壽命之間存在很強的相關(guān)性,我們可能會問,孩子和孫子的出生是否具有對一個人的剩余生命的影響(或者我們是否可以像這樣假設(shè)獨立性)。

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