Ошибка в упаковке монет или неправильно поставленная проблема?

Question

Я получаю некоторые противоречивые результаты при использовании аргумента весов в пакете монет.В частности, для функций kruskal_test и spearman_test.

С обычными данными все работает нормально и согласуется с kruskal.test в пакете статистики:

> x <- xtabs( ~gear + vs,data=mtcars)
> df <- as.data.frame.table(x)
> kruskal_test(gear ~ as.factor(vs),data=mtcars)

    Asymptotic Kruskal-Wallis Test

data:  gear by as.factor(vs) (0, 1) 
chi-squared = 2.4768, df = 1, p-value = 0.1155
> kruskal.test(gear ~ as.factor(vs),data=mtcars)

    Kruskal-Wallis rank sum test

data:  gear by as.factor(vs) 
Kruskal-Wallis chi-squared = 2.4768, df = 1, p-value = 0.1155

Но, когда одни и те же данные передаются вkruskal_test с частотными весами, я получаю неверный результат.

> kruskal_test(as.numeric(df[[1]]) ~ df[[2]],
+ weights=~as.integer(df[[3]]))

    Asymptotic Kruskal-Wallis Test

data:  as.numeric(df[[1]]) by df[[2]] (0, 1) 
chi-squared = 1.3158, df = 1, p-value = 0.2513

Есть ли проблема с тем, как я настраиваю этот вызов функции?

Answer 1

Это действительно было ошибкой.Торстен ответил, что ранговая трансформация не учитывает веса.Следующий код демонстрирует неквалифицированную версию теста, дающую идентичные результаты:

> oneway_test(as.integer(gear) ~ vs, data = df, weights  = ~ Freq)

    Asymptotic 2-Sample Permutation Test

data:  as.integer(gear) by vs (0, 1) 
Z = -1.1471, p-value = 0.2513
alternative hypothesis: true mu is not equal to 0 

Надеюсь, это будет исправлено в будущем.