使用wilcoxonpairedR测试R中的效应大小

时间:2019-04-15 10:54:36

标签: r

我已经在R中运行了一些简单的配对采样Wilcoxon测试,现在我想测试效果大小。我一直在使用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)

我进行了配对的Wilcoxon检验,以查看A组参与者之间的度量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

1 个答案:

答案 0 :(得分:1)

似乎您正在将错误的参数传递给wilcox.test()wilcoxonPairedR()wilcox.test()的公式接口需要measure1 ~ group之类的公式,而您要提供measure1 ~ measure2。同样,在wilcoxonPairedR()中,您的xg都是度量,而此函数期望x是度量的向量,而g是度量的向量组标签的对应向量。