R语言模拟人类生活预期寿命动态可视化动画图gif

这周,我在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年出生的人

现在,如果我可以获得历史数据集,包括在学校花费的平均时间,退休年龄等,可以将其添加到图表中。


可下载资源

关于作者

Kaizong Ye拓端研究室(TRL)的研究员。

本文借鉴了作者最近为《R语言数据分析挖掘必知必会 》课堂做的准备。

​非常感谢您阅读本文,如需帮助请联系我们!

 
QQ在线咨询
售前咨询热线
15121130882
售后咨询热线
0571-63341498

关注有关新文章的微信公众号


永远不要错过任何见解。当新文章发表时,我们会通过微信公众号向您推送。

技术干货

最新洞察

This will close in 0 seconds