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

#NGroups 2
#define nvar2 2

Title 1: MZ 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: DZ 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