我已经在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
答案 0 :(得分:1)
似乎您正在将错误的参数传递给wilcox.test()
和wilcoxonPairedR()
。 wilcox.test()
的公式接口需要measure1 ~ group
之类的公式,而您要提供measure1 ~ measure2
。同样,在wilcoxonPairedR()
中,您的x
和g
都是度量,而此函数期望x
是度量的向量,而g
是度量的向量组标签的对应向量。