r - R中的引导覆盖

Question

我会估计引导区间的覆盖率，因为知道真实平均值是895.0385. 我有我的向量b<-c(300,300,200,250,600...)，我做了引导和输出间隔：

mean.fun <- function(dat, idx) mean(dat[idx], na.rm = TRUE)
boot.out <- boot(b, mean.fun, R=999)
boot.ci(boot.out)

但是我如何复制它以获得覆盖概率（它包含真实平均值的多少次）？

Answer 1

不久前，我试图做一些类似的事情。我没有使用引导命令，而是使用了示例命令，但这可能会有所帮助。我也可能 100% 错了，我对 R 还不是很好。

mean.fun <- function(dat, idx) mean(dat[idx], na.rm = TRUE)
bootoutput <- data.frame(
  bootoutput = replicate(10000, boot.ci(boot(b, mean.fun, R=999)))
)