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

#r

#r

Вопрос:

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

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

  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 между участниками в группе A:

  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

Ответ №1:

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

Комментарии:

1. Это работает, спасибо! Я развернул данные и запустил их так, чтобы x = subset(exampledata, group == "A")$Value и g = subset(exampledata, group == "A")$Measure .