R語言模擬人類生活預期壽命動態可視化動畫圖gif

原文鏈接: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年出生的人

 

現在,如果我可以獲得歷史數據集,包括在學校花費的平均時間,退休年齡等,可以將其添加到圖表中。

 

發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章