!
! 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.500000   
#define md21  0.500000 
#define md31  0.500000 
#define md41  0.200000 
#define md51  0.200000 
#define md61  0.200000 
#define md71 -0.500000 
#define md81 -0.500000 
#define md91 -0.500000 

#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:
!
! CALC>p=.1
!  0.100000000000000
! CALC>q=.9
!  0.900000000000000
! CALC>p**4+p**3*q+p*p*q*q/4.
!  0.00302500000000000
! CALC>750*(p**4+p**3*q+p*p*q*q/4.)
!  2.26875000000000
! CALC>750*(p**3+p**2*q**2/2)
!  3.78750000000000
! CALC>750*(p**2*q**2/4)
!  1.51875000000000
! CALC>750*(p**3*q+p**2*q**2/2)
!  3.71250000000000
! CALC>750*(p**3*q+3*p**2*q**2+p*q**3)
!  73.5750000000000
! CALC>750*(p**2*q**2/2+p*q**3)
!  57.7125000000000
! CALC>750*(p*p*q*q/4)
!  1.51875000000000
! CALC>750*(p**2*q**2/2+p*q**3)
!  57.7125000000000
! CALC>750*(p**2*q**2/4+p*q**3+q**4)
!  548.268750000000
!

#define nd1   2.26875
#define nd2    3.7875
#define nd3    1.51875
#define nd4    3.7125
#define nd5   73.575
#define nd6    57.7125
#define nd7    1.51875
#define nd8    57.7125
#define nd9   548.26875

! and for the MZ's
! CALC>250*p*p
!  2.50000000000000
! CALC>250*2*p*q
!  45.0000000000000
! CALC>250*q*q
!  202.500000000000
!
#define nm1   2.500 
#define nm2  45.000 
#define nm3   202.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