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

#NGroups 5
#define nvar2 2

Title 1: MZf data
 #include ozbmi.dat
 Select if zyg =1
 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;
 Options RSiduals
End

Title 2: DZf data
 #include ozbmi.dat
 Select if zyg =3
 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;
 Options RSiduals
End

Title 3: MZm data
 #include ozbmi.dat
 Select if zyg =2
 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;
 Options RSiduals
End

Title 4: DZm data
 #include ozbmi.dat
 Select if zyg =4
 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;
 Options RSiduals
End

Title 5: DZfm data
 #include ozbmi.dat
 Select if zyg =5
 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;
 Options RSiduals
End