! C5. Age-correction Model 
! Australian Conservatism Data of Female Twins

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

G2: common environmental structure
Data Calc 
Matrices
C Lower 1 1 Free
Compute C*C' /
End

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

G4: age
Data Calc 
Matrices
Y Lower 1 1 Free   ! effect of age on phenotype
Compute Y*Y' /
End

G5: female MZ twin pairs
Data NInput_vars=3 NObservations=941
Labels
age cons_t1 cons-t2
CMatrix Symmetric
179.522
 63.463 143.361
 66.167  97.159 141.530
Matrices
A Symm 1 1 = %E1  ! expected matrix of group 1
C Symm 1 1 = %E2  ! expected matrix of group 2
E Symm 1 1 = %E3  ! ..
Y Symm 1 1 = %E4
S Lower 1 1 =Y4
V Lower 1 1 Free  ! standard deviation of age
Covariances V*V  | S*V     | S*V    _
            S*V  | A+C+E+Y | A+C+Y  _
            S*V  | A+C+Y   | A+C+E+Y /
End 

G6: female dz twin pairs
Data NInput_vars=3 NObservations=548
Labels
age cons_t1 cons-t2
CMatrix Symmetric
168.343
 66.365 141.840
 76.316  88.083 157.763
Matrices
A Symm 1 1 = %E1
C Symm 1 1 = %E2
E Symm 1 1 = %E3
D Symm 1 1 = %E4
H Full 1 1
Y Symm 1 1 = %E4
S Lower 1 1 =Y4
V Lower 1 1 =V5
Covariances V*V  | S*V     | S*V    _
            S*V  | A+C+E+Y | H@A+C+Y  _
            S*V  | H@A+C+Y | A+C+E+Y /
Matrix H .5
Start 5 All
Start 15 V 5 1 1
Options NDecimals=4
End