! Estimate means and variances - Saturated model
! OZ BMI data - young females

#NGroups 2
#define nvar2 2

Title 1: MZ data
 #include ozbmimz.dat
 Select bmi1 bmi2 ;
 Begin Matrices = Group 1;
  M Full 1 nvar2 Free
  X Symm nvar2 nvar2 Free
 End Matrices;
  Start 20 M 1 1 - M 1 nvar2
  Start 1 X 1 1 X 2 2
 Means M;
 Covariance X;
 Option RSiduals
End

Title 2: DZ data
 #include ozbmidz.dat
 Select bmi1 bmi2 ;
 Begin Matrices;
  M Full 1 nvar2 Free
  X Symm nvar2 nvar2 Free
 End Matrices;
  Start 20 M 1 1 - M 1 nvar2
  Start 1 X 1 1 X 2 2
 Means M;
 Covariance X;
 Option RSiduals
 Option Multiple Issat
End

! equate means and variances
 Equate M 1 1 1 M 1 1 2  M 2 1 1 M 2 1 2
 Equate X 1 1 1 X 1 2 2  X 2 1 1 X 2 2 2
End