! Estimate Correlations for MZ & DZ Twins
! Dutch Adult Twins: Lipid levels

#define nvar 1
#define nvarx2 2
#NGroups 2

G1: MZ twins
#include DutchMZ.dat
  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 .6
 Means M;
 Covariance S*R*S';
End

G2: DZ twins
#include DutchDZ.dat
  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';
 Option Multiple Issat
End

! Test for familial resemblance
! Set MZ & DZ correlations equal
 Equate R 1 2 1 R 2 2 1
Exit