R中的引导覆盖

【问题标题】：Bootstrap coverage in RR中的引导覆盖
【发布时间】：2012-12-05 11:29:04
【问题描述】：

知道真实平均值为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)

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

【问题讨论】：

查看?replicate函数。
B

标签： r statistics-bootstrap

【解决方案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)))
)

【讨论】：

是的，我误会了，当我看到时，我找不到通过编辑删除它的方法。 =/