R中的Bootstrap覆盖

如果知道真实平均值是

895.0385

，我会估计平均值的引导间隔覆盖率。

---

我有我的向量

b我刚才想做一些类似的事情。我没有使用boot命令，我使用了示例命令，但这可能会有所帮助。我也可能是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)))
)

mean.fun查看
？replicate
函数。B是的，我搞错了，当我看到时，我找不到通过编辑删除它的方法=/
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)))
)