! G4. Scaler GxE Interaction 
! Australian Depression Data of Concordant and Discordant Married Twins

G1: genetic structure
Data Calc NGroups=9
Matrices
A Lower 1 1 Free
Compute A*A' /
End

G2: dominance structure
Data Calc  
Matrices
D Lower 1 1 Fixed
Compute D*D' /
End

G3: specific environmental structure
Data Calc 
Matrices
E Lower 1 1 Free
Compute E*E' /
End

G4: Concordant single MZ twin pairs
Data NInput_vars=2 NObservations=254
Labels
dep_t1 bmi-t2
CMatrix Symmetric
0.0968 0.0381 0.0896
Matrices
A Symm 1 1 = %E1
D Symm 1 1 = %E2
E Symm 1 1 = %E3
Covariances A+D+E | A+D _
            A+D   | A+D+E /
End 

G5: Concordant single DZ twin pairs
Data NInput_vars=2 NObservations=155
Labels
dep_t1 bmi-t2
CMatrix Symmetric
0.1087 0.0250 0.1182
Matrices
A Symm 1 1 = %E1
D Symm 1 1 = %E2
E Symm 1 1 = %E3
H Full 1 1
Q Full 1 1
Covariances A+D+E   | H@A+Q@D _
            H@A+Q@D | A+D+E /
Matrix H .5
Matrix Q .25
End

G6: Concordant married MZ twin pairs
Data NInput_vars=2 NObservations=177
Labels
dep_t1 bmi-t2
CMatrix Symmetric
0.1002 0.0336 0.0769
Matrices
A Symm 1 1 = %E1
D Symm 1 1 = %E2
E Symm 1 1 = %E3
K Diag 2 2        ! scalar
Covariances K *(A+D+E | A+D _
                A+D   | A+D+E) *K' /
Specify K 4 4
End 

G7: Concordant married DZ twin pairs
Data NInput_vars=2 NObservations=107
Labels
dep_t1 bmi-t2
CMatrix Symmetric
0.0692 0.0076 0.0882
Matrices
A Symm 1 1 = %E1
D Symm 1 1 = %E2
E Symm 1 1 = %E3
H Full 1 1 = H5
Q Full 1 1 = Q5
K Diag 2 2 = K6
Covariances K *(A+D+E   | H@A+Q@D _
                H@A+Q@D | A+D+E ) *K' /
End

G8: Discordant MZ twin pairs
Data NInput_vars=2 NObservations=139
Labels
dep_t1 bmi-t2
CMatrix Symmetric
0.1198 0.0359 0.1009
Matrices
A Symm 1 1 = %E1
D Symm 1 1 = %E2
E Symm 1 1 = %E3
K Diag 2 2
Covariances K *(A+D+E | A+D _
                A+D   | A+D+E  ) *K' / 
Specify K 0 4
Matrix K 1 1
End

G9: Discordant DZ twin pairs
Data NInput_vars=2 NObservations=87
Labels
dep_t1 bmi-t2
CMatrix Symmetric
0.1013 0.0050 0.0694
Matrices
A Symm 1 1 = %E1
D Symm 1 1 = %E2
E Symm 1 1 = %E3
H Full 1 1 = H5
Q Full 1 1 = Q5
K Diag 2 2 = K8
Covariances K *(A+D+E   | H@A+Q@D _
                H@A+Q@D | A+D+E ) *K' /
Start .4 All
Options NDecimals=4
End