【问题标题】:Problems with accessing double elements in a list in R在 R 中访问列表中的双元素的问题
【发布时间】: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


【解决方案1】:

看起来您需要apply 系列函数。您的数据不可重现,因此我无法确认这是否有效,但如果您这样做:

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

你应该得到你想要的。

【讨论】:

    【解决方案2】:

    您可能想看看purrr 包和map 函数系列,或tidyrhoist 函数。

    (如果你想要代码有效,你确实需要提供一些数据!)

    【讨论】:

      猜你喜欢
      • 1970-01-01
      • 2021-12-20
      • 2014-09-24
      • 1970-01-01
      • 2020-05-24
      • 2021-05-17
      • 1970-01-01
      • 1970-01-01
      相关资源
      最近更新 更多