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

我之前尝试过这样的事情。我没有使用我使用sample命令的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)))
)

R中的Bootstrap覆盖范围

1 个答案: