我正在尝试使用R的vcdExtra软件包中的Cochran-Mantel-Haenszel测试在两个变量(两级分类和七级序数)之间找到关联。我是该测试的新手,我'已经看到许多人将它用于两个两个级别的变量,但是我想知道使用带有七个级别的序数变量的测试是否仍然有效。
我首先尝试使用Wilcoxon检验,但由于我的序数变量是一个因素而不是数字,所以无法使用。
以下是Wilcoxon检验的输出
wilcox.test(testdata$freq_use ~ testdata$gender)
Error in wilcox.test.default(x = c(7L, 4L, 3L, 5L, 1L, 7L, 1L, 7L, 3L, :
'x' must be numeric
CMH测试的输出有效,并且代码运行正常(如下)。我只是想知道鉴于我的序数变量具有多个级别,因此测试本身是否有效。
library(vcdExtra)
CMH.test(table(testdata$gender,testdata$freq_use)
AltHypothesis<chr> Chisq<dbl> Df<dbl>
cor Nonzero correlation 83.32795 1
rmeans Row mean scores differ 114.66140 6
cmeans Col mean scores differ 83.32795 1
general General association 114.66140 6