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

74 篇文章 5 订阅
35 篇文章 0 订阅

最近我们被要求撰写关于期望寿命的研究报告,包括一些图形和统计输出。

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

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

  • 0
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值