! Estimate Correlations for DZ Twins by IBD status
! Dutch Adult Twins: Lipid levels

#define nvar 1
#define nvarx2 2
#NGroups 3

G1: DZ IBD2 twins
#include DutchDZ.dat
  Select if pi65cat =2;
  Select
   t1$var
   t2$var ;
 Begin Matrices; 
  M Full nvar nvarx2 Free      ! means
  S Diag nvarx2 nvarx2 Free    ! standard deviations
  R Stnd nvarx2 nvarx2 Free    ! correlations
 End Matrices;
  Matrix M 4 4                 ! starting values 
  Matrix S .5 .5
  Matrix R .5
 Means M;
 Covariance S*R*S';
End

G2: DZ IBD1 twins
#include DutchDZ.dat
  Select if pi65cat =1;
  Select
   t1$var
   t2$var ;
 Begin Matrices;
  M Full nvar nvarx2 Free
  S Diag nvarx2 nvarx2 Free
  R Stnd nvarx2 nvarx2 Free
 End Matrices;
  Matrix M 4 4
  Matrix S .5 .5
  Matrix R .5
 Means M;
 Covariance S*R*S';
End

G3: DZ IBD0 twins
#include DutchDZ.dat
  Select if pi65cat =0;
  Select
   t1$var
   t2$var ;
 Begin Matrices;
  M Full nvar nvarx2 Free
  S Diag nvarx2 nvarx2 Free
  R Stnd nvarx2 nvarx2 Free
 End Matrices;
  Matrix M 4 4
  Matrix S .5 .5
  Matrix R .5
 Means M;
 Covariance S*R*S';
End