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年出生的人
现在,如果我可以获得历史数据集,包括在学校花费的平均时间,退休年龄等,可以将其添加到图表中。