vlambda博客
学习文章列表

原文链接：http://tecdat.cn/?p=13552

这周，我在http://waitbutwhy.com/上发现了一张图片  ，它代表了典型的人类生活，

我觉得很有趣。但是第一个问题是，我不理解下面的年龄限制：90年，不是平均寿命。第二个问题是，当您查看图片时，它不能像看起来的那样静态。0岁时的预期寿命与30岁或50岁时的预期寿命不同。因此，我尝试使用预期寿命表制作动画图表。这是一个代码，用于生成不同时期的生活表（我在这里区分是男性还是女性）

L2 <- lifetable(france.fcast)ex2=L2$exL1=lifetable(fr.mort,series="female")ex1=L1$exexF=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年出生的人

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