! Age Correction Model
! Conservatism data from Australian female twins

G1: model parameters
Data Calc NGroups=3
Begin Matrices;
X Lower 1 1 Free    ! genetic structure
Y Lower 1 1 Free    ! common environmental structure
Z Lower 1 1 Free    ! specific environmental structure
S Lower 1 1 Free    ! effect of age on phenotype
V Lower 1 1 Free    ! variance of age
End Matrices;
Begin Algebra;
A= X*X' ;
C= Y*Y' ;
E= Z*Z' ;
G= S*S' ;
End Algebra;
End

G2: female MZ twin pairs
Data NInput_vars=3 NObservations=941
Labels age cons_t1 cons-t2
CMatrix Symmetric File=ozconmzf.cov
Matrices= Group 1
Covariances V*V  | S*V     | S*V    _
            S*V  | A+C+E+G | A+C+G  _
            S*V  | A+C+G   | A+C+E+G /
Option RSidual
End 

G3: female dz twin pairs
Data NInput_vars=3 NObservations=548
Labels age cons_t1 cons-t2
CMatrix Symmetric File=ozcondzf.cov
Matrices= Group 1
H Full 1 1
Covariances V*V  | S*V     | S*V      _
            S*V  | A+C+E+G | H@A+C+G  _
            S*V  | H@A+C+G | A+C+E+G /
Matrix H .5
Start 5 All
Start 15 V 1 1 1
Options NDecimals=4
Option RSidual
End