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

私は少し前にこのようなことを少しやろうとしていました。サンプルコマンドを使用した 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)))
)