Dlaczego błąd „szacunkowa korekta” a „NA” jest generowany z pakietu rozruchowego R podczas obliczania przedziałów ufności przy użyciu metody bca?

Mam wektor liczb, które przesłałem tutaj (... / code / MyData.Rdata) przy użyciu dput.

Chciałbym dostać bca ci, więc napisałem ten kod:

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

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

Ale gdy uruchomię następujące, otrzymuję to:

> 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

Czy ktoś może mi pomóc dowiedzieć się, co powoduje ten błąd? Dziękuję za pomoc!

Jak widać z komunikatu o błędzie, boot.cipołączenia bca.ci. Ponieważ boot.outobiekt nie dostarcza danych L, empiryczne wartości wpływu dla statystyki, którą obliczasz na danych, bca.cipróbuje je obliczyć za pomocą empinffunkcji, a następnie (jak mówi Michael) używa ich do obliczenia stałej przyspieszenia:

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

Ale przy niewielkiej liczbie replikacji empinfczasami zawodzi i zwraca wektor NAwartości. Powoduje to, że nie masz wartości dla L, anie można obliczyć i pojawia się błąd. Jak mówi ocram, zwiększenie liczby replik boostrap to naprawi. RPrawdopodobnie nawet podwojenie do 2000 r.