原文鏈接:http://tecdat.cn/?p=13552
這周,我在http://waitbutwhy.com/上發現了一張圖片 ,它代表了典型的人類生活,
我覺得很有趣。但是第一個問題是,我不理解下面的限制:90年,不是平均壽命。第二個問題是,當您查看圖片時,它不能像看起來的那樣靜態。0歲時的預期壽命與30歲或50歲時的預期壽命不同。因此,我嘗試使用預期壽命表製作動畫圖表。這是一個代碼,用於生成不同時期的生活表(我在這裏區分是男性還是女性)
L2 <- lifetable(france.fcast)
ex2=L2$ex
L1=lifetable(fr.mort,series="female")
ex1=L1$ex
exF=cbind(ex1,ex2)
france.LC1 <- lca(fr.mort,adjust="e0",series="male",years=c(1900,2100))
根據這些生命表,我們可以得出某個特定年份(例如1950年)出生的人在各個年齡段(例如50、51、52等)的剩餘預期壽命。根據這些預期的剩餘壽命,我們可以繪製
plot(va%%52,va%/%52,cex=.6,pch=15,col=c("light yellow","pink","white")[1+
(va>=age*52)*1+(va>(age+F[age+1])*52)*1],ylim=c(100,0),axes=FALSE,xlab="Week",
ylab="Age",main=paste("Woman, born on ",yearborn,
", age ",age,sep=""))
例如,如果我們想要上面的圖,對於1980年出生的30歲的某人,我們使用
picture(1980,30)
現在,如果我們運行代碼以獲取動畫gif,就可以爲1950年出生的人獲得
對於2000年出生的人
現在,如果我可以獲得歷史數據集,包括在學校花費的平均時間,退休年齡等,可以將其添加到圖表中。