bcaメソッドを使用して信頼区間を計算すると、Rブートパッケージから「推定調整 'a'はNA」というエラーが生成されるのはなぜですか？

dputを使用してここにアップロードした数値のベクトル（... / code / MyData.Rdata）があります。

bca ciを取得したいので、次のコードを記述しました。

my.mean <- function(dat, idx){
 return (mean(dat[idx], na.rm = TRUE))
} 

boot.out<-boot(data=my.data, statistic = my.mean, R=1000)

しかし、私が以下を実行すると、これが得られます：

> boot.ci(boot.out)
Error in bca.ci(boot.out, conf, index[1L], L = L, t = t.o, t0 = t0.o,  : 
estimated adjustment 'a' is NA
In addition: Warning message:
In boot.ci(boot.out) : bootstrap variances needed for studentized intervals

誰かがこのエラーの原因を突き止めるのを手伝ってくれる？ご協力ありがとうございました！

エラーメッセージからわかるように、をboot.ci呼び出しますbca.ci。そのためboot.out、オブジェクトが、供給はしないL、あなたはデータ上で計算している統計の経験的影響値をbca.ci使用してそれらを計算しようとするempinf機能を、次に（マイケルが言うように）それが加速定数を計算するためにそれらを使用しています。

L <- empinf(boot.out, index = index, t = t.o, ...)
a <- sum(L^3)/(6 * sum(L^2)^1.5)

ただし、レプリケーションの数が少ないと、empinf失敗して値のベクトルが返されることがありますNA。その結果L、の値aがなく、計算できず、エラーが発生します。ocramが言うように、ブーストラップレプリケーションの数を増やすと、これが修正されます。R2000に2倍にしても、おそらくそれを行うはずです。