Тестирование величины эффекта в R с использованием wilcoxonpairedR

Question

Я провел несколько простых тестов Вилкоксона с парными выборками в R, и теперь я хочу проверить размеры эффекта. Я использовал wilcoxonPairedR из пакета effsize, но почему-то не могу заставить его работать с моими данными; Я получаю сообщение об ошибке Error in validObject(.Object) : invalid class “SymmetryProblem” object: FALSE.

Вот некоторые примеры данных (обратите внимание, что в группе B есть NA, которые я отфильтрую для теста ниже:

 exampledata <- tribble(~subject, ~group, ~Measure1, ~Measure2,
     "1", "A", .8, .23,
     "2", "B", NA, .79,
     "3", "A", .6, .28,
     "4", "B", NA, .18,
     "5", "A", .86, .90,
     "6", "A", .34, .23,
     "7", "B", NA, .12,
     "8", "B", NA, .27)

Я запускаю парный тест Вилкоксона, чтобы увидеть, отличается ли Мера 1 от Меры 2 между участниками в группе А:

 wilcox.test(Measure1 ~ Measure2, subset(exampledata, group =="A"), paired=T)

Затем я проверяю размер эффекта:

 wilcoxonPairedR(x = subset(exampledata, group == "A")$Measure1,
                 g = subset(exampledata, group == "A")$Measure2)

Который генерирует сообщение об ошибке. Когда я фильтрую данные вручную, а затем тестирую с использованием отфильтрованных подмножеств данных, это также не работает:

 exampledataA <- exampledata %>% filter(group == "A")

 wilcoxonPairedR(x = exampledataA$Measure1,
                 g = exampledataA$Measure2)

Error in validObject(.Object) : invalid class “SymmetryProblem” object: FALSE

Answer 1

Похоже, что вы передаете неверные аргументы wilcox.test() и wilcoxonPairedR().Интерфейс формул wilcox.test() ожидает формулу, такую ​​как measure1 ~ group, тогда как вы предоставляете measure1 ~ measure2.Аналогично, в wilcoxonPairedR() ваши x и g являются мерами, тогда как эта функция ожидает, что x будет вектором мер, а g будет соответствующим вектором групповых меток.