619/FirstModelSyntax.R at master · joyceylee/619 · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
```{r}
#First model (stability model) for 619 project

load("H:\\BSF Processing Data\\Datasets\\Dataset for R\\Data Processing\\semdataCtl\\Everything4.RData")

library(lavaan)
Model1_2<-'

###Measurement model
#Coparenting self-report by father
CopW2.d =~ CO1_A.d.W2 + CO1_B.d.W2 + CO1_C.d.W2 + CO1_D.d.W2 + CO1_E.d.W2 + CO1_F.d.W2 + CO1_G.d.W2 + CO1_H.d.W2 + CO1_I.d.W2 + CO1_J.d.W2
CopW3.d =~ CO1_A.d.W3 + CO1_B.d.W3 + CO1_C.d.W3 + CO1_D.d.W3 + CO1_E.d.W3 + CO1_F.d.W3 + CO1_G.d.W3 + CO1_H.d.W3 + CO1_I.d.W3 + CO1_J.d.W3

#Coparenting self-report by mother
CopW2.m =~ CO1_A.m.W2 + CO1_B.m.W2 + CO1_C.m.W2 + CO1_m.d.W2 + CO1_E.m.W2 + CO1_F.m.W2 + CO1_G.m.W2 + CO1_H.m.W2 + CO1_I.m.W2 + CO1_J.m.W2
CopW3.m =~ CO1_A.m.W3 + CO1_B.m.W3 + CO1_C.m.W3 + CO1_m.d.W3 + CO1_E.m.W3 + CO1_F.m.W3 + CO1_G.m.W3 + CO1_H.m.W3 + CO1_I.m.W3 + CO1_J.m.W3

#Second order coparenting latent variables
CopW2 =~ CopW2.d + CopW2.m
CopW3 =~ CopW3.d + CopW3.m

#Father engagement in caregiving
CareW2.d =~ CO3_F.W2 + CO3_G.W2 + CO3_H.W2
CareW3.d =~ CO3_F.W3 + CO3_G.W3 + CO3_H.W3

#Father engagement in play
PlayW2.d =~ CO3_B.W2 + CO3_C.W2 + CO3_D.W2 + CO3_E.W2
PlayW3.d =~ CO3_B.W3 + CO3_C.W3 + CO3_D.W3 + CO3_E.W3

##Correlated residuals for repeated measures
CO1_A.d.W2 ~~ CO1_A.d.W3
CO1_B.d.W2 ~~ CO1_B.d.W3
CO1_C.d.W2 ~~ CO1_C.d.W3
CO1_D.d.W2 ~~ CO1_D.d.W3
CO1_E.d.W2 ~~ CO1_E.d.W3
CO1_F.d.W2 ~~ CO1_F.d.W3
CO1_G.d.W2 ~~ CO1_G.d.W3
CO1_H.d.W2 ~~ CO1_H.d.W3
CO1_I.d.W2 ~~ CO1_I.d.W3
CO1_J.d.W2 ~~ CO1_AJ.d.W3

CO1_A.m.W2 ~~ CO1_A.m.W3
CO1_B.m.W2 ~~ CO1_B.m.W3
CO1_C.m.W2 ~~ CO1_C.m.W3
CO1_D.m.W2 ~~ CO1_D.m.W3
CO1_E.m.W2 ~~ CO1_E.m.W3
CO1_F.m.W2 ~~ CO1_F.m.W3
CO1_G.m.W2 ~~ CO1_G.m.W3
CO1_H.m.W2 ~~ CO1_H.m.W3
CO1_I.m.W2 ~~ CO1_I.m.W3
CO1_J.m.W2 ~~ CO1_AJ.m.W3

CO3_F.W2 ~~ CO3_F.W3
CO3_G.W2 ~~ CO3_G.W3
CO3_H.W2 ~~ CO3_H.W3

CO3_B.W2 ~~ CO3_B.W3
CO3_C.W2 ~~ CO3_C.W3
CO3_D.W2 ~~ CO3_D.W3
CO3_E.W2 ~~ CO3_E.W3

##Suggestions from modification indices
CO3_C.W2 ~~ CO3_D.W2
CO3_C.W3 ~~ CO3_D.W3
PlayW2.d ~~ CareW2.d
CopW2.d ~~ CopW3.d
CopW2.m ~~ CopW3.m

###Structural model (controls upfront)
CopW2 ~ BlackOnly + HispanicCombo + OtherCombo + Work + BSF + HSDiploma + EduOther + MPF + IPVW2.d2 + IPVW2.m2
CareW2.d ~ BlackOnly + HispanicCombo + OtherCombo + Work + BSF + HSDiploma + EduOther + MPF + IPVW2.d2 + IPVW2.m2
PlayW2.d ~ BlackOnly + HispanicCombo + OtherCombo + Work + BSF + HSDiploma + EduOther + MPF + IPVW2.d2 + IPVW2.m2

CopW3 ~ BlackOnly + HispanicCombo + OtherCombo + Work + BSF + HSDiploma + EduOther + MPF + IPVW2.d2 + IPVW2.m2 + CopW2
CareW3.d ~ BlackOnly + HispanicCombo + OtherCombo + Work + BSF + HSDiploma + EduOther + MPF + IPVW2.d2 + IPVW2.m2 + CareW2.d
PlayW3.d ~ BlackOnly + HispanicCombo + OtherCombo + Work + BSF + HSDiploma + EduOther + MPF + IPVW2.d2 + IPVW2.m2 + PlayW2.d
'

Model1Fit_2<-sem(Model1_2, data=semdataCtl5.IPV4, estimator="MLR", missing="fiml")
summary(Model1Fit_2, fit.measure=T, standardized=T, rsquare=T)
#residuals(Model1Fit_2, type='cor')
```