Lavaan - CFA - 分类变量 - 最后一个阈值很奇怪

Question

我想在lavaan与R进行多组CFA。

我有几个分类变量，一些变量包含11个类别。所以这些变量将有10个阈值。在下面的结果中，您可以看到第10个阈值小于第9个阈值，即它不在折痕顺序中。

具有11个类别的若干变量具有相同的问题。

Question:

为什么阈值会失真？

R-code:

model2<-'range = ~ NA*gvjbevn + gvhlthc + gvslvol + gvslvue + gvcldcr + gvpdlwk

goals = ~ NA*sbprvpv + sbeqsoc + sbcwkfm

range~~1*range

goals~~1*goals

gvhlthc ~~ gvslvol

gvcldcr ~~ gvpdlwk

'

cfa.model2<-cfa(model2, ordered=varcat, estimator="WLSMV",data=sub)

summary(cfa.model2,fit.measures=TRUE,standardized=TRUE, modindices=TRUE)

Answer 1

好吧，由于没有足够的声誉，似乎我无法添加评论，所以我只能回复一个答案，虽然这不是一个正确的答案（它肯定不会解决你的问题，但我希望它指向右边方向）。

为了使您的示例具有可重现性，您应该为社区提供适合模型的数据。

另一方面，我猜您的问题必须与该类别的性质有关：您的第11类可能并不意味着与该项目达成“最高级别的协议”，或者响应类别不是从1开始订购的到11，或类似的东西。鉴于其余阈值似乎准确地表示连续的，单调递增的比例，并且同样的问题恰好发生在不同变量的同一类别中（至少是你所展示的两个），必须有一些回应这些项目的规模。

总之，它似乎更多地是解释模型参数而不是统计问题的问题。

Answer 2

阈值的标签分配按字母顺序排序，也就是c('t1','t10','t2','t3'....)，但summary()将其“正确”排序。

您可以尝试添加其他因素来检查您的比例是否与以下内容相对应：

c('t1','t10','t11','t12',...,'t2','t3'....)

除了了解哪一行是你的每个因素之外，你不能做多少事情。