!
! Power calcs with MZ's and DZ's and abaw model
!
! md's are means
! nd's are sample sizes for dz
! nm's are sample sizes for mz
!
! For reference, proportions of variance for various values of a (=aw)
!
! CALC>a=.1
!  0.100000000000000
! CALC>d=0
!  0.00000000000000
! CALC>p=.5
!  0.500000000000000
! CALC>q=.5
!  0.500000000000000
! CALC>2*p*q*(a+d*(q-p))**2
!  0.00500000000000000
! CALC>a=.2 0.100000  0.100000
! CALC>2*p*q*(a+d*(q-p))**2
!  0.0200000000000000
! CALC>a=.3
! CALC>2*p*q*(a+d*(q-p))**2
!  0.0450000000000000
! CALC>a=.4
! CALC>2*p*q*(a+d*(q-p))**2
!  0.0800000000000000
! CALC>a=.5
! CALC>2*p*q*(a+d*(q-p))**2
!  0.125000000000000
!

#define md11  0.200000   
#define md21  0.200000 
#define md31  0.200000 
#define md41  0.500000 
#define md51  0.500000 
#define md61  0.500000 
#define md71 -0.200000 
#define md81 -0.200000 
#define md91 -0.200000 

#define md12   0.500000
#define md22   0.200000
#define md32  -0.500000
#define md42   0.500000
#define md52   0.200000
#define md62  -0.500000
#define md72   0.500000
#define md82   0.200000
#define md92  -0.500000

! To obtain sample sizes for the DZ's:
#define nd1   105.469
#define nd2    70.313
#define nd3    11.719
#define nd4    70.313
#define nd5   234.375
#define nd6    70.313
#define nd7    11.719
#define nd8    70.313
#define nd9   105.469

#define nm1   62.500 
#define nm2  125.000 
#define nm3   62.500   

#ngroups 12

Group 1 DZ AAAA
 Data NI=2 NObs=nd1 
 CM 
 1 .3 1
 Means md11 md12
 Begin matrices;
  A Lower nvarx2 nvarx2 Free    ! Cholesky for covariance
  E Full 4 1 
  N Full 1 1  Free              ! Overall mean
  W Full 2 1  Free              ! ab and aw parameters
  V Full 9 1                    ! 
  X Full 9 1
 End Matrices;

 Matrix V  1  .5 0 .5 0 -.5 0 -.5 -1
 Matrix X  0 .5 1 -.5 0 .5 -1 -.5 0
 Matrix E 1 1 1 2
 Matrix A 1 0 1

 Begin Algebra;
  Z = \part((V|X),E)*W | \part((V|-X),E)*W ; ! selects row of V|X and V|-X
 End Algebra;

 Means (N|N) + Z;
 Covariance A*A' ;
End Group

Group 2 DZ AAAa
 Data NI=2 NObs=nd2 
 CM 
 1 .3 1
 Means md21 md22
 Begin matrices = Group 1;
  E Full 4 1 
 End Matrices;

 Matrix E 2 1 2 2

 Begin Algebra;
  Z = \part((V|X),E)*W | \part((V|-X),E)*W ; ! selects row of V|X and V|-X
 End Algebra;

 Means (N|N) + Z;
 Covariance A*A' ;
End Group

Group 3 DZ AAaa
 Data NI=2 NObs=nd3 
 CM 
 1 .3 1
 Means md31 md32
 Begin matrices = Group 1;
  E Full 4 1 
 End Matrices;

 Matrix E 3 1 3 2

 Begin Algebra;
  Z = \part((V|X),E)*W | \part((V|-X),E)*W ; ! selects row of V|X and V|-X
 End Algebra;

 Means (N|N) + Z;
 Covariance A*A' ;
End Group


Group 4 DZ AaAA
 Data NI=2 NObs=nd4 
 CM 
 1 .3 1
 Means md41 md42
 Begin matrices = Group 1;
  E Full 4 1 
 End Matrices;

 Matrix E 4 1 4 2

 Begin Algebra;
  Z = \part((V|X),E)*W | \part((V|-X),E)*W ; ! selects row of V|X and V|-X
 End Algebra;

 Means (N|N) + Z;
 Covariance A*A' ;
End Group

Group 5 DZ AaAa
 Data NI=2 NObs=nd5 
 CM 
 1 .3 1
 Means md51 md52
 Begin matrices = Group 1;
  E Full 4 1 
 End Matrices;

 Matrix E 5 1 5 2

 Begin Algebra;
  Z = \part((V|X),E)*W | \part((V|-X),E)*W ; ! selects row of V|X and V|-X
 End Algebra;

 Means (N|N) + Z;
 Covariance A*A' ;
End Group

Group 6 DZ Aaaa
 Data NI=2 NObs=nd6 
 CM 
 1 .3 1
 Means md61 md62
 Begin matrices = Group 1;
  E Full 4 1 
 End Matrices;

 Matrix E 6 1 6 2

 Begin Algebra;
  Z = \part((V|X),E)*W | \part((V|-X),E)*W ; ! selects row of V|X and V|-X
 End Algebra;

 Means (N|N) + Z;
 Covariance A*A' ;
End Group

Group 7 DZ aaAA
 Data NI=2 NObs=nd7 
 CM 
 1 .3 1
 Means md71 md72
 Begin matrices = Group 1;
  E Full 4 1 
 End Matrices;

 Matrix E 7 1 7 2

 Begin Algebra;
  Z = \part((V|X),E)*W | \part((V|-X),E)*W ; ! selects row of V|X and V|-X
 End Algebra;

 Means (N|N) + Z;
 Covariance A*A' ;
End Group

Group 8 DZ aaAa
 Data NI=2 NObs=nd8 
 CM 
 1 .3 1
 Means md81 md82
 Begin matrices = Group 1;
  E Full 4 1 
 End Matrices;

 Matrix E 8 1 8 2

 Begin Algebra;
  Z = \part((V|X),E)*W | \part((V|-X),E)*W ; ! selects row of V|X and V|-X
 End Algebra;

 Means (N|N) + Z;
 Covariance A*A' ;
End Group

Group 9 DZ aaaa
 Data NI=2 NObs=nd9 
 CM 
 1 .3 1
 Means md91 md92
 Begin matrices = Group 1;
  E Full 4 1 
 End Matrices;

 Matrix E 9 1 9 2

 Begin Algebra;
  Z = \part((V|X),E)*W | \part((V|-X),E)*W ; ! selects row of V|X and V|-X
 End Algebra;

 Means (N|N) + Z;
 Covariance A*A' ;
End 

Group 10 MZ AAAA
 Data NI=2 NObs=nm1 
 CM 
 1 .3 1
 Means md11 md12
 Begin matrices = Group 1;
  E Full 4 1 
 End Matrices;

 Matrix E 1 1 1 2

 Begin Algebra;
  Z = \part((V|X),E)*W | \part((V|-X),E)*W ; ! selects row of V|X and V|-X
 End Algebra;

 Means (N|N) + Z;
 Covariance A*A' ;
End Group

Group 11 MZ AaAa
 Data NI=2 NObs=nm2 
 CM 
 1 .3 1
 Means md51 md52
 Begin matrices = Group 1;
  E Full 4 1 
 End Matrices;

 Matrix E 5 1 5 2

 Begin Algebra;
  Z = \part((V|X),E)*W | \part((V|-X),E)*W ; ! selects row of V|X and V|-X
 End Algebra;

 Means (N|N) + Z;
 Covariance A*A' ;
End Group

Group 12 MZ aaaa
 Data NI=2 NObs=nm3 
 CM 
 1 .3 1
 Means md91 md92
 Begin matrices = Group 1;
  E Full 4 1 
 End Matrices;

 Matrix E 9 1 9 2

 Begin Algebra;
  Z = \part((V|X),E)*W | \part((V|-X),E)*W ; ! selects row of V|X and V|-X
 End Algebra;

 Means (N|N) + Z;
 Covariance A*A' ;
 Option Multiple
End 

Drop W 1 2 1
Option power=.05,1
End

Drop W 1 1 1
End
!
!Specify 1 W 99 99
!End