! Estimate thresholds & correlations -  Saturated model
! US MDD data - females                 ordinal data

#NGroups 2
#define nvar2 2

Title 1: MZ data
 Data NInput=3
 Ordinal File=usmdd.ord
 Labels zyg mdd1 mdd2
 Select if zyg =1;
 Select mdd1 mdd2 ;
 Begin Matrices;
  T Full 1 nvar2 Free
  X Stnd nvar2 nvar2 Free
 End Matrices;
  Start 1 T 1 1 - T 1 nvar2
 Threshold T;
 Correlation X;
 Options RSiduals
End

Title 2: DZ data
 Data NInput=3
 Ordinal File=usmdd.ord
 Labels zyg mdd1 mdd2
 Select if zyg =2;
 Select mdd1 mdd2 ;
 Begin Matrices;
  T Full 1 nvar2 Free
  X Stnd nvar2 nvar2 Free 
 End Matrices;
  Start 1 T 1 1 - T 1 nvar2
 Threshold T;
 Correlation X;
 Options RSiduals
 Option Multiple Issat
End

! equal thresholds (means/variances) within zygosity
 Equate T 1 1 1 T 1 1 2
 Equate T 2 1 1 T 2 1 2
End

! equal thresholds (means/variances) across zygosity
 Equate T 1 1 1 T 1 1 2  T 2 1 1 T 2 1 2
End