!job for Falconer GxE model

! Analyse the simulated data
G1: model parameters
Data Calc NGroups=3
Begin Matrices;
X Lower 2 2 Fr   ! genetic structure
Y Lower 2 2 Fr   ! common environmental structure
Z Lower 2 2 Fr   ! specific environmental structure
End Matrices;

!Fix X 1 2 1
Fix X 1 2 2
Start .3 X 1 1 1 X 1 2 2 
Start .1 Y 1 1 1 Y 1 2 2 
Start .5 Z 1 1 1 Z 1 2 2 
Begin Algebra;
A= X*X' ;
C= Y*Y' ;
E= Z*Z' ;
k=\stnd(X*X');	! calculate the genetic correlation
l=\stnd(C);		! common E correlation
m=\stnd(E);		! unique E correlation
End Algebra;
End

G2: MZ twin pairs
Data NInput_vars=4 NObservations=100
CMatrix Full File=mzsim.cov
Matrices= Group 1
Covariances A+C+E | A+C _
            A+C   | A+C+E /
Option RSiduals
End

G3: DZ twin pairs
Data NInput_vars=4 NObservations=100
CMatrix Full File=dzsim.cov
Matrices= Group 1
H Full 1 1
Covariances A+C+E  | H@A+C _
            H@A+C  | A+C+E /
Matrix H .5
End