【发布时间】:2020-07-19 13:45:54
【问题描述】:
我使用 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]]
我已经尝试了多种选择,但我无法让它发挥作用。 触发我的事情是,我正在使用列表列表中的双精度。
我已经添加了下面数据的图片:
有没有人可以解决这个问题?
【问题讨论】:
-
请使用
dput(yourdata)提供数据的子样本。这大大增加了获得帮助的机会
标签: r dimensions statistics-bootstrap