我在R中有一个小数据集,其中二元和Likert(0-4)项已合并。将二进制项重新编码为likert,其中1> 1。 0和2> 4.我试图在这些数据上运行polychoric correlation(r package:psych)和ordinal alpha,但是,我收到错误,因为我的二进制项本质上不是连续的。
我正在寻找一种方法来告诉R我的变量的可能值范围是0-4,尽管事实上我的所有受访者都没有为这些项选择值1-3,所以我可以运行序数alpha数据集。我们的PI已经考虑过将二进制项分成它们自己的分量表是否更好;这对我们的目的来说是不够的。
感谢您提供的任何见解!
编辑:包括示例。
DATA:
q1 <- c(4,2,0,1,3)
q2 <- c(1,4,2,3,0)
q3 <- c(0,4,4,0,4)
q4 <- c(0,4,4,0,4)
q5 <- c(4,0,4,4,4)
data <- data.frame(q1,q2,q3,q4,q5)
PACKAGES:
库(迷幻) 库(GPArotation)
POLYCOR:
polychoric(数据)
POLYCOR OUTPUT:
通话:polychoric(x =数据)
多重相关
q1 q2 q3 q4 q5
q1 1.0
q2 -0.4 1.0
q3 NA NA 1.0
q4 NA NA NA 1.0
q5 NA NA NA NA 1.0
与tau 0 1 2 3 q1 -0.84 -0.25 0.25 0.84 q2 -0.84 -0.25 0.25 0.84 q3 -0.25 -0.25 -0.25 -0.25 q4 -0.25 -0.25 -0.25 -0.25 q5 -0.84 -0.84 -0.84 -0.84 有50个或更多警告(使用警告()查看前50个)
警告讯息: 1:优化(polyF,interval = c(-1,1),rc = taux,cc = tauy,...: NA / Inf由最大正值代替(所有50都相同)。
CREATE POLYCOR MATRIX:
datmtx&lt; - polychoric(data)
Polycor矩阵结果:
警告讯息: 1:优化(polyF,interval = c(-1,1),rc = taux,cc = tauy,...: NA / Inf替换为最大正值
ALPHA:
α(datmtx $ RHO)
具有缺失值的可能变量是q3 q4 q5
校长错误(x):
对不起,相关矩阵中的缺失值(NAs)不允许我继续。
请删除这些变量,然后重试。