Почему при расчете доверительных интервалов с использованием метода bca генерируется ошибка «оценочная корректировка« a »является NA» из загрузочного пакета R?

У меня есть вектор чисел, который я загрузил здесь (... / code / MyData.Rdata), используя dput.

Я хотел бы получить 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вы не можете рассчитать, и вы получите свою ошибку. Как говорит Окрам, увеличение количества репликаций Boostrap это исправит. Даже удвоение Rдо 2000 года, вероятно, должно сделать это.