! Cholesky Model
! US BMI data 

#NGroups 4
#define nvar 5
#define nvarx2 10

G1: female model parameters
Calculation
 Begin Matrices;
  X Lower nvar nvar Free    ! genetic cholesky
  Y Lower nvar nvar Free    ! common env cholesky
  W Lower nvar nvar Fixed   ! dominance cholesky
  Z Lower nvar nvar Free    ! specific env cholesky
  I Iden 2 2
  H Full 1 1
  Q Full 1 1
 End Matrices;
  Matrix H .5
  Matrix Q .25
  Start .4 All
  Start .5 Z 1 1 1 Z 1 2 2 Z 1 3 3 Z 1 4 4 Z 1 5 5 
 Begin Algebra;
  A= X*X' ;
  C= Y*Y' ;
  D= W*W' ;
  E= Z*Z' ;
 End Algebra;
End

G2: MZ twin pairs
#include bmi5mzf.dat
 Begin Matrices =Group 1
  M Full 1 nvarx2 Free
 End Matrices;
  Start 10 M 1 1 to M 1 nvarx2
 Means M ;
 Covariances  (A+C+D+E | A+C+D _
               A+C+D   | A+C+D+E) ;
Option No_Output
End 

G3: DZ twin pairs
#include bmi5dzf.dat
 Begin Matrices =Group 1
  M Full 1 nvarx2 Free
 End Matrices;
  Start 10 M 1 1 to M 1 nvarx2
 Means M ;
 Covariances  (A+C+D+E   | H@A+C+Q@D _
               H@A+C+Q@D | A+C+D+E  ) ;
Option No_Output
End

G4
Calculation
 Begin Matrices = Group 1
  I Iden nvar nvar
 End Matrices;
 Begin Algebra;
  V= A+C+D+E ;
  S= \sqrt(I.(V)) ;
  M= ((S~*X)_ (S~*Y)_ (S~*W)_ (S~*Z));
 End Algebra;
 Option NDecimals=3
End