原文鏈接：http:///?p=17347昨天上午，我們獲得了分娩產(chǎn)婦的平均年齡兩個(gè)圖表，根據(jù)孩子的出生順序排序，區(qū)間是1905-1965年（點(diǎn)擊文末“閱讀原文”獲取完整代碼數(shù)據(jù)）。 
然后是1960-2000年： 
點(diǎn)擊標(biāo)題查閱往期內(nèi)容 這些圖形令人興奮，特別是在過去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)深藍(lán)色的曲線， 
獲取祖母的平均年齡，我們進(jì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模型對年死亡率進(jìn)行建模，并推斷到當(dāng)前世紀(jì)，我們可以推斷出期望剩余壽命。 > 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)在我們也可以對這項(xiàng)快速研究的局限性感到疑惑。特別是，正如有配偶的壽命之間存在很強(qiáng)的相關(guān)性，我們可能會(huì)問，孩子和孫子的出生是否具有對一個(gè)人的剩余生命的影響（或者我們是否可以像這樣假設(shè)獨(dú)立性）。
|
