我正在尝试用 R 语言在 lavaan 中指定一个奇怪的模型。该模型如下所示：
我的规范尝试如下所示。我发现难以实现的是将观察到的变量的唯一误差固定为唯一项目的两个相关性的总和。
例如项目 y*1,2 与 y*1,3 和 y*2,3 共变，它的错误应该是 cov y*1,3 + cov y*1,3。
如何在下面的 lavaan 语法中显式修复项目错误以等于这些协方差的总和？
cfa_model_spesification<-'
C=~ #C4_12*i10i11+C4_13*i10i12+
#C5_12*i13i14+C5_13*i13i15+
#C6_12*i17i18+C6_13*i17i19+
C1_12*i1i2+C1_13*i1i3+
C2_12*i4i5+C2_13*i4i6+
C3_12*i7i8+C3_13*i7i9
R=~ #R4_23*i10i11+R4_12*i11i12+
#R5_23*i13i14+R5_12*i14i16+
#R6_23*i17i18+R6_12*i18i19+
R1_12*i1i2+R1_23*i2i3+
R2_12*i4i5+R2_23*i5i6+
R3_12*i7i8+R3_23*i8i9
O=~ #O4_13*i10i12+O4_23*i11i12+
#O5_13*i13i15+O5_23*i14i16+
#O6_13*i17i19+O6_23*i18i19+
O1_13*i1i3+O1_23*i2i3+
O2_13*i4i6+O2_23*i5i6+
O3_13*i7i9+O3_23*i8i9
O~~1*O
C~~1*C
R~~1*R
O~~C+R
R~~C
R1_23==-R1_12
R2_23==-R2_12
R3_23==-R3_12
R1_23>0
R2_23>0
R3_23>0
# R1_12<0
# R2_12<0
# R3_12<0
O1_13<0
O1_23<0
O2_13<0
O2_23<0
O3_13<0
O3_23<0
i1i2~~i1i3
i1i2~~i2i3
i1i3~~i2i3
i4i5~~i4i6
i4i5~~i5i6
i4i6~~i5i6
i7i8~~i7i9
i7i8~~i8i9
i7i9~~i8i9
i1i2~~1*i1i2
i4i5~~1*i4i5
i7i8~~1*i7i8
# i1i3~~equal("i1i3~~i1i2+i1i3~~i2i3")*i1i3
# i2i3~~equal("i2i3~~i1i2+i2i3~~i1i3")*i2i3
# i4i6~~equal("i4i6~~i4i5+i4i6~~i5i6")*i4i6
# i5i6~~equal("i5i6~~i4i5+i5i6~~i4i6")*i5i6
# i7i9~~equal("i7i9~~i7i8+i7i9~~i8i9")*i7i9
# i8i9~~equal("i8i9~~i7i8+i8i9~~i7i9")*i8i9
'
mplus 中 this 的语法如下所示
TITLE:
Example
DATA:
FILE IS triplets.dat;
VARIABLE:
NAMES=i1i2 i1i3 i2i3 i4i5 i4i6 i5i6 i7i8 i7i9 i8i9 i10i11 i10i12 i11i12;
CATEGORICAL=i1i2-i11i12;
ANALYSIS:
ESTIMATOR=ulsmv;
PARAMETERIZATION=THETA;
MODEL:
Trait1 BY
i1i2*1 i1i3*1 (L1)
i4i5*-1 i4i6*-1 (L4)
i7i8*1 i7i9*1 (L7)
i10i11*1 i10i12*1 (L10);
Trait2 BY
i1i2*-1 (L2_n)
i2i3*1 (L2)
i4i5*-1 (L5_n)
i5i6*1 (L5)
i7i8*-1 (L8_n)
i8i9*1 (L8)
i10i11*1 (L11_n)
i11i12*-1 (L11);
Trait3 BY
i1i3*-1 i2i3*-1 (L3_n)
i4i6*-1 i5i6*-1 (L6_n)
i7i9*1 i8i9*1 (L9_n)
i10i12*-1 i11i12*-1 (L12_n);
Trait1-Trait3@1
Trait1 WITH Trait2*-0.4 Trait3*0;
Trait2 WITH Trait3*0.3;
i1i2*2 (e1e2);
i1i3*2 (e1e3);
i2i3*2 (e2e3);
i4i5*2 (e4e5);
i4i6*2 (e4e6);
i5i6*2 (e5e6);
i7i8*2 (e7e8);
i7i9*2 (e7e9);
i8i9*2 (e8e9);
i10i11*2 (e10e11);
i10i12*2 (e10e12);
i11i12*2 (e11e12);
i1i2 WITH i1i3*1 (e1);
i1i2 WITH i2i3*-1 (e2_n);
i1i3 WITH i2i3*1 (e3);
i4i5 WITH i4i6*1 (e4);
i4i5 WITH i5i6*-1 (e5_n);
i4i6 WITH i5i6*1 (e6);
i7i8 WITH i7i9*1 (e7);
i7i8 WITH i8i9*-1 (e8_n);
i7i9 WITH i8i9*1 (e9);
i10i11 WITH i10i12*1 (e10);
i10i11 WITH i11i12*-1 (e11_n);
i10i12 WITH i11i12*1 (e12);
MODEL CONSTRAINT:
L2_n=-L2;
L5_n=-L5;
L8_n=-L8;
L11_n=-L11;
e1e2=e1-e2_n;
e1e3=e1+e3;
e2e3=-e2_n+e3;
e4e5=e4-e5_n;
e4e6=e4+e6;
e5e6=-e5_n+e6;
e7e8=e7-e8_n;
e7e9=e7+e9;
e8e9=-e8_n+e9;
e10e11=e10-e11_n;
e10e12=e10+e12;
e11e12=-e11_n+e12;
e1=1;
e4=1;
e7=1;
e10=1;