Ошибка вычисления доверительного интервала начальной загрузки с использованием метода BCa ('w' бесконечно) - PullRequest
1 голос
/ 17 июня 2020

Я пытаюсь вычислить CI с начальной загрузкой, используя метод BCa для моей оценки W Кендалла, но получаю следующее сообщение об ошибке:

Error in bca.ci(boot.out, conf, index[1L], L = L, t = t.o, t0 = t0.o,  : 
  estimated adjustment 'w' is infinite

Мой фактический набор данных составляет 4908 строк. Все остальные методы (norm, per c, et c.) Работают с этим набором данных, и у меня нет пропущенных значений. Я также безуспешно пытался с 10000 повторений. Аналогичный вопрос был задан здесь , но, к сожалению, ответов не было. Интересно, знал ли кто-нибудь, что не так?

Вот тестовый набор данных с моими кодами:

library(rcompanion)

testdata = structure(list(var1 = structure(c(3, 1, 1, 1, 3, 0, 3, 
                                             0, 1, 0, 0, 1, 0, 0, 
                                             0, 0, 0, 1, 1, 0), label = "Variable 1", class = c("labelled", "numeric")), 
                          var2 = structure(c(1, 0, 0, 0, 3, 0, 3, 
                                             1, 1, 0, 1, 0, 0, 0, 
                                             1, 1, 0, 0, 1, 0), label = "Variable 2", class = c("labelled","numeric"))), 
                   row.names = c(NA, -20L), class = c("tbl_df", "tbl", "data.frame"))

testci = rcompanion::kendallW(testdata, correct = TRUE, ci=TRUE, conf=0.95, type="bca", R=1000) 

...