Чому помилка "розрахункове коригування" a "є NA" генерується з пакета завантаження R при обчисленні довірчих інтервалів методом bca?

У мене є векторний номер, який я завантажив сюди (... / 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їх не можна обчислити, і ви отримаєте свою помилку. Як каже ocram, збільшення кількості реплікацій boostrap це виправить. Навіть подвоєння Rдо 2000 року повинно це зробити.