! Calculate power for univariate ACE model
! 	30% additive genetic
! 	20% common environment
!    	50% nonshared environment

#NGroups 3

G1: model parameters
Calculation
Begin Matrices;
	A full 1 1 fixed
	C full 1 1 fixed
	E full 1 1 fixed
End Matrices;

! Specify A, C and E directly
Matrix A 0.3
Matrix C 0.2	
Matrix E 0.5

End

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

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


! Second part of script
! Fit the wrong model to the simulated data 
! to calculate power

#NGroups 3
G1 : model parameters
Calculation
Begin Matrices;
	X lower 1 1 free
	Y lower 1 1 fixed
	Z lower 1 1 free
End Matrices;

Begin Algebra;
	A = X*X' ;
	C = Y*Y' ; 
	E = Z*Z' ; 
End Algebra;
End

G2 : MZ twins
Data NInput_vars=2 NObservations=350
CMatrix Full File=mzsim.cov
Matrices= Group 1
Covariances A+C+E | A+C   _
            A+C   | A+C+E /
Option RSiduals
End

G3 : DZ twins
Data NInput_vars=2 NObservations=350
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 /
Matix H 0.5
Option RSiduals

! Power for alpha = 0.05 and 1 df
Option Power= 0.05,1 
End