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

#NGroups 2
#define nvar2 2

Title 1: MZ data
 #include ozbmi2.dat
 Select if zyg =1
 Select if agecat =1
 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
End

Title 2: DZ data
 #include ozbmi2.dat
 Select if zyg =3
 Select if agecat =1
 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                   ! continue in multiple fit mode
                                         ! set this model as saturated model
End

Save ozbmiyfeqmv.mxs                     ! save binary file
! 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

Get ozbmiyfeqmv.mxs                      ! load saved binary file
! equate means within zygosity groups
 Equate M 1 1 1 M 1 1 2
 Equate M 2 1 1 M 2 1 2
End

! equate means across zygosity groups
 Equate M 1 1 1 M 1 1 2  M 2 1 1 M 2 1 2
End

! equate variances within zygosity groups
 Equate X 1 1 1 X 1 2 2
 Equate X 2 1 1 X 2 2 2
End

! equate variances across zygosity groups
 Equate X 1 1 1 X 1 2 2  X 2 1 1 X 2 2 2
End