我正在处理调查数据并尝试使用 SEM 对其进行R
分析lavaan
数据无缺失值,共1723例。
我的模型有 2 个因子和两个序数内生变量。
模型如下所示:
#model1
model1 <- '
# measurement model
F1=~ Q15+Q16+Q17
F2=~ Q18+Q19+Q20+Q22+Q23
# regressions
V1~ V2+ V3+ V4
HKER ~ V2+ V3+ V4 + V5
F2~ Ordin1+ V5
F1~ V5
Ordin2~ parscor + Ordin1+ F2 + F1
# residual correlations
Q15 ~~ Q16+Q17
Q16 ~~ Q17
Q18 ~~ Q19+Q20+Q22+Q23
Q19 ~~ Q20+Q22+Q23
Q20 ~~ Q22+Q23
Q22 ~~ Q23
V1~~ V2 + V3+ V4
V2~~ V3 + V4
V3 ~~ V4
'
所有变量,无论是外生变量还是内生变量,都是序数变量。
该程序显示两个警告:
- 优化器 (NLMINB) 声称模型收敛,但并非所有梯度元素都(接近)为零;优化器可能没有找到局部解;使用 check.gradient = FALSE 跳过此检查。
- 构造 W 矩阵的麻烦;对 nlminb 中的 A11 子矩阵误差使用广义逆(start = start.x,objective = objective_function,gradient = GRADIENT,:
我按照以下说明进行操作:
http://lavaan.ugent.be/tutorial/cat.html
并将 Ordin1 和 Ordin2 声明为序数变量,ordered
但它不起作用。并且所有估计的标准误差都丢失了。
有谁知道这里发生了什么?