Я пытаюсь вычислить 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)