R語言模擬人類生活期望壽命動(dòng)態(tài)可視化動(dòng)畫圖gif
原文鏈接:http://tecdat.cn/?p=13552
?
這周,我在http://waitbutwhy.com/上發(fā)現(xiàn)了一張圖片??,它代表了典型的人類生活,

?
我覺得很有趣。但是第一個(gè)問題是,我不理解下面的限制:90年,不是平均壽命。第二個(gè)問題是,當(dāng)您查看圖片時(shí),它不能像看起來的那樣靜態(tài)。0歲時(shí)的期望壽命與30歲或50歲時(shí)的期望壽命不同。因此,我嘗試使用期望壽命表制作動(dòng)畫圖表。這是一個(gè)代碼,用于生成不同時(shí)期的生活表(我在這里區(qū)分是男性還是女性)
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))
根據(jù)這些生命表,我們可以得出某個(gè)特定年份(例如1950年)出生的人在各個(gè)年齡段(例如50、51、52等)的剩余期望壽命。根據(jù)這些期望的剩余壽命,我們可以繪制
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=""))
例如,如果我們想要上面的圖,對(duì)于1980年出生的30歲的某人,我們使用
picture(1980,30)
現(xiàn)在,如果我們運(yùn)行代碼以獲取動(dòng)畫gif,就可以為1950年出生的人獲得

?
對(duì)于2000年出生的人

?
現(xiàn)在,如果我可以獲得歷史數(shù)據(jù)集,包括在學(xué)?;ㄙM(fèi)的平均時(shí)間,退休年齡等,可以將其添加到圖表中。
?