! Mdd meta analysis clin samples only

#NGroups 15

G1: Model parameters, male
 Calc 
Begin Matrices;
 H Full 1 1
 Q Full 1 1
 X Lower 1 1 Free
 Y Lower 1 1 Free
 Z Lower 1 1 Free
 W Lower 1 1 
End Matrices;
 Matrix H .5
 Matrix Q .25
 Start .6 All
Begin Algebra;
 A= X*X';
 C= Y*Y';
 E= Z*Z';
 D= W*W';
End Algebra;
End

G2: Model parameters, female
Calc
Begin Matrices;
 H Full 1 1
 Q Full 1 1
 X Lower 1 1 Free
 Y Lower 1 1 Free
 Z Lower 1 1 Free
 W Lower 1 1 
End Matrices;
 Matrix H .5
 Matrix Q .25
 Start .6 All
Begin Algebra;
 A= X*X';
 C= Y*Y';
 E= Z*Z';
 D= W*W';
End Algebra;
End

G3: Constrain male variance
 Constraint NI=1
 Begin Matrices = Group 1;
  U Unit 1 1
 End Matrices;
 Constrain U=A+C+E+D;
End

G4: Constrain female variance
 Constraint NI=1
 Begin Matrices = Group 2;
  U Unit 1 1
 End Matrices;
 Constrain U=A+C+E+D;
 Option Multiple Issat
End

G21:  mcgmzf.dat
#include mcgmzf.dat
Begin Matrices= Group 2;
 T Full 2 1
End Matrices;
Matrix T 0.796055 0.796055
Thresholds T;  
Covariances A+C+D+E | A+C+D _
            A+C+D   | A+C+D+E /
Options RSidual
Option Multiple
End

G22:  mcgdzf.dat
#include mcgdzf.dat
Begin Matrices= Group 2;
 T Full 2 1
End Matrices;
Matrix T 0.796055 0.796055
Thresholds T;
Covariances A+C+D+E   | H@A+C+Q@D _
            H@A+C+Q@D | A+C+D+E /  
Options RSidual
Option Multiple
End

G23:  mcgmzm.dat
#include mcgmzm.dat
Begin Matrices= Group 1;
 T Full 2 1
End Matrices;
Matrix T 1.140687 1.140687
Thresholds T;
Covariances A+C+D+E   | A+C+D _
            A+C+D | A+C+D+E /
Options RSidual
Option Multiple
End

G24:  mcgdzm.dat
#include mcgdzm.dat
Begin Matrices= Group 1;
 T Full 2 1
End Matrices;
Matrix T 1.140687 1.140687
Thresholds T;
Covariances A+C+D+E   | H@A+C+Q@D _
            H@A+C+Q@D | A+C+D+E /
Options RSidual
Option Multiple
End

G25:  kenclinmzf.dat
#include kenclinmzf.dat
Begin Matrices= Group 2;
 T Full 2 1
End Matrices;
Thresholds T;
Matrix T 0.796055 0.796055 
Covariances A+C+D+E | A+C+D _
            A+C+D   | A+C+D+E /
Options RSidual
End 

G26:  kenclindzf.dat
#include kenclindzf.dat
Begin Matrices= Group 2;
 T Full 2 1
End Matrices;
Matrix T 0.796055 0.796055 
Thresholds T;
Covariances A+C+D+E   | H@A+C+Q@D _
            H@A+C+Q@D | A+C+D+E /
Options RSidual
Option Multiple
End

G27: kenclinmzm.dat
#include kenclinmzm.dat
Begin Matrices= Group 1;
 T Full 2 1
End Matrices;
Matrix T 1.140687 1.140687
Thresholds T;
Covariances A+C+D+E   | A+C+D _
            A+C+D | A+C+D+E /
Options RSidual
Option Multiple
End

G28:  kenclindzm.dat
#include kenclindzm.dat
Begin Matrices= Group 1;
 T Full 2 1
End Matrices;
Matrix T 1.140687 1.140687
Thresholds T;
Covariances A+C+D+E   | H@A+C+Q@D _
            H@A+C+Q@D | A+C+D+E /
Options RSidual
Option Multiple
End

G29: kenclindzof.dat
#include kenclindzof.dat
Begin Matrices= Group 2;
 T Full 2 1
 L Comp = A2
 M Comp = C2
 N Comp = E2
 O Comp = D2
 R Full 1 1 Free  ! parameter 201
 G Lower 1 1 = X2 
 I Lower 1 1 = Y2
 K Lower 1 1 = W2
End Matrices;
Bound 0 1 R 1 1
Matrix T 0.796055 1.140687 
Thresholds T;
Covariances A+C+D+E   |  R@H@G*X'+Y*I'+Q@W*K' _
            R@H@X*G'+I*Y'+Q@K*W' | L+M+N+O /
Options RSidual
Option Multiple
End

G30:  kenclindzom.dat
#include kenclindzom.dat
Begin Matrices= Group 2;
 T Full 2 1
 L Comp = A2
 M Comp = C2
 N Comp = E2
 O Comp = D2
 R Full 1 1 =R13  ! parameter 201
 G Lower 1 1 = X2 
 I Lower 1 1 = Y2
 K Lower 1 1 = W2
End Matrices;
Bound 0 1 R 1 1
Matrix T 0.796055 1.140687 
Thresholds T;
Covariances A+C+D+E   |  R@H@G*X'+Y*I'+Q@W*K' _
            R@H@X*G'+I*Y'+Q@K*W' | L+M+N+O /

Options RSidual
Option Multiple
End

G31: Summary
 Calculation
 Begin Matrices = Group 14;
  V Full 1 14
 End Matrices;
 Begin Algebra;
  B = A|C|E|D|L|M|N|O|R|V ;
 End Algebra;

 Option format=(15(F4.2,1x))
 Option append
 Option mxb=commun.mat
 Option issat
End