访问R中列表中的双元素的问题

问题描述 投票:0回答:1

我已经对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。

我在下面添加了数据图片:

https://ibb.co/4fMGv8R

有人对此有解决方案吗?

r bootstrapping dimensions
1个回答
0
投票

似乎您需要apply系列功能。您的数据不可复制,因此我无法确认是否可以,但是如果您这样做:

result <- sapply(JA_lc_fitM_boot1[["bootParameters"]], function(var) var[["bx"]][[70]])

您应该得到想要的东西。

© www.soinside.com 2019 - 2024. All rights reserved.