我已经对Lee Carter模型的2.000个重采样执行了自举,以进行死亡率预测。这个问题并非专门针对死亡率研究,而是在R中更一般的维度。
执行自举后,我得到一个包含2000个元素的列表,每个元素对应于模型的2.000个重新估计。对于每个模型,我的3个变量都有一个估计值:a_x,b_x和k_t。a_x和b_x都是特定年龄的,因此“ x”表示间隔[0:95]中的年龄。
我现在想绘制年龄x = 70的所有b_x值的直方图。
### Performing the bootstrap:
JA_lc_fitM_boot1 <- bootstrap(LCfit_JA_M, nBoot = 2000, type = "semiparametric")
### Plotting the histogram with all b_x for x = 70:
JA_lc_fitM_boot1[["bootParameters"]][1:2000][["bx"]][[70]]
我尝试了多种选择,但无法使其正常工作。触发我的是,我正在使用列表中的一个double。
我在下面添加了数据图片:
有人对此有解决方案吗?
似乎您需要apply
系列功能。您的数据不可复制,因此我无法确认是否可以,但是如果您这样做:
result <- sapply(JA_lc_fitM_boot1[["bootParameters"]], function(var) var[["bx"]][[70]])
您应该得到想要的东西。