! ACE Model for Power Calculations
! Simulated data

! Simulate the data
!    30% additive genetic   (.5477ý=.3)
!    20% common environment (.4472ý=.2)
!    50% random environment (.7071ý=.5)

G1: model parameters
Data Calc NGroups=3
Begin Matrices;
X Lower 1 1 Fixed   ! genetic structure
Y Lower 1 1 Fixed   ! common environmental structure
Z Lower 1 1 Fixed   ! specific environmental structure
End Matrices;
Matrix X .5477 
Matrix Y .4472 
Matrix Z .7071
Begin Algebra;
A= X*X' ;
C= Y*Y' ;
E= Z*Z' ;
End Algebra;
End

G2: MZ twin pairs
Data Calc NInput_vars=2 
Matrices= Group 1
Covariances A+C+E | A+C _
            A+C   | A+C+E /
Options MX%E=mzsim.cov  
End

G3: DZ twin pairs
Data Calc NInput_vars=2 
Matrices= Group 1
H Full 1 1
Covariances A+C+E  | H@A+C _
            H@A+C  | A+C+E /
Matrix H .5
Options MX%E=dzsim.cov
End

! Fit the wrong model to the simulated data

G1: model parameters
Data Calc NGroups=3
Begin Matrices;
X Lower 1 1 Free    ! genetic structure
Y Lower 1 1 Fixed   ! common environmental structure
Z Lower 1 1 Free    ! specific environmental structure
End Matrices;
Begin Algebra;
A= X*X' ;
C= Y*Y' ;
E= Z*Z' ;
End Algebra;
End

G2: MZ twin pairs
Data NInput_vars=2 NObservations=1000
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=2 NObservations=1000
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
Start .5 All
Options RSiduals Power= .05,1     ! .05 sig level & 1 df
End