! PD_META_ACE.MX	September 2000
! Meta-analysis of family and twin data for panic disorder
! Compare whether resemblance (S) in relatives predicted by FS data = TS data
! Note: Used raw numbers for FS data 
! NCS total lifetime prevalence (3.5%) used for thresholds (1.811) in FS


!***** GENERAL ACE MODEL FIRST ***********

#NGroups 18

G1: Model parameters
Calculation
Begin Matrices;
 X Full 1 1 Free        ! additive genetic factor
 Y Full 1 1 Free        ! shared environmental factor 
 Z Full 1 1 Free        ! specific environmental factor
 U Full 1 1 Free        ! family resemblance factor (1 for TS, 0 for FS)
 H Full 1 1
End Matrices;
 Specify X 1
 Specify Y 2
 Specify Z 3
 Specify U 5
 Matrix H .5
Begin Algebra;
 A= X*X';
 C= Y*Y';
 E= Z*Z';
 S= U*U';
End Algebra;
 Start .6 All
 Boundary 0.00001 .99999 all
 Boundary -5 5  20 30 50
Interval A 1 1
Interval C 1 1
Interval E 1 1
Interval S 1 1
End

G2: Constraint group to ensure a*a + c*c + e*e + s*s = 1
Data Constraint NInput=1
Begin Matrices = Group 1;
 I Iden 1 1
End Matrices;
Constraint I = A+C+E+S;
Option NO_Output
Option Multiple Issat
Option RSiduals
Option NDecimals=4 
Option Iterations=1000
End

G3: Scherrer Male MZ twin pairs
Data NInput_vars=2 
Lables STwin1 STwin2
Ctable 2 2	! Contingency Table (T1*T2)
 1673   83
   83   27
Begin Matrices= Group 1
 T Full 2 1
End Matrices
Specify T  20 20
Threshold T;
Covariances A+C+E+S | A+C+S _
            A+C+S   | A+C+E+S /
Option RSiduals
End 

G4: Scherrer Male DZ twin pairs
Data NInput_vars=2 
Lables STwin1 STwin2
Ctable 2 2	! Contingency Table (T1*T2)
 1325   82
   82    6
Begin Matrices= Group 1
 T Full 2 1
End Matrices
Specify T  20 20
Threshold T;
Covariances A+C+E+S   | H@A+C+S _
            H@A+C+S | A+C+E+S /
Option RSiduals
End 

G5: Kendler Female MZ twin pairs
Data NInput_vars=2 
Lables KTwin1 KTwin2
Ctable 2 2	! Contingency Table (T1*T2)
 511   35
  35    9
Begin Matrices= Group 1
 T Full 2 1
End Matrices
Specify T  30 30
Threshold T;
Covariances A+C+E+S | A+C+S _
            A+C+S   | A+C+E+S /
Option RSiduals
End 

G6: Kendler Female DZ twin pairs
Data NInput_vars=2 
Lables KTwin1 KTwin2
Ctable 2 2	! Contingency Table (T1*T2)
 376   29
  29    5
Begin Matrices= Group 1
 T Full 2 1
End Matrices
Specify T  30 30
Threshold T;
Covariances A+C+E+S   | H@A+C+S _
            H@A+C+S | A+C+E+S /
Option RSiduals
End 

G7: Perna MZ twin pairs
Data NInput_vars=2 
Lables PTwin1 PTwin2
Ctable 2 2	! Contingency Table (T1*T2)
 19   2
  1   4
Begin Matrices= Group 1
 T Full 2 1
End Matrices
Specify T  50 50
Threshold T;
Covariances A+C+E+S | A+C+S _
            A+C+S   | A+C+E+S /
Option RSiduals
End 

G8: Perna DZ twin pairs
Data NInput_vars=2 
Lables PTwin1 PTwin2
Ctable 2 2	! Contingency Table (T1*T2)
 28   3
  3   0
Begin Matrices= Group 1
 T Full 2 1
End Matrices
Specify T  50 50
Threshold T;
Covariances A+C+E+S   | H@A+C+S _
            H@A+C+S | A+C+E+S /
Option RSiduals
End 

G9: FS1 (Noyes) - probands
Data Ninput=2
Labels Rel1 Rel2
Ctable  2 2
 -1  -1
205  36
Begin Matrices= Group 1
 T Full 2 1
End Matrices
Matrix T  1.811 1.811
Threshold T;
Covariances A+C+E   | H@A+C _
            H@A+C | A+C+E /
Option RSiduals
End 


G10: FS1 (Noyes) - controls
Data Ninput=2
Labels Rel1 Rel2
Ctable  2 2
109   4
 -1  -1
Begin Matrices= Group 1
 T Full 2 1
End Matrices
Matrix T  1.811 1.811
Threshold T;
Covariances A+C+E   | H@A+C _
            H@A+C | A+C+E /
Option RSiduals
End 


G11: FS2 (Mendlewicz) - probands
Data Ninput=2
Labels Rel1 Rel2
Ctable  2 2
 -1  -1
109  13
Begin Matrices= Group 1
 T Full 2 1
End Matrices
Matrix T  1.811 1.811
Threshold T;
Covariances A+C+E   | H@A+C _
            H@A+C | A+C+E /
Option RSiduals
End 

G12: FS2 (Mendlewicz) - controls
Data Ninput=2
Labels Rel1 Rel2
Ctable  2 2
129   1
 -1  -1
Begin Matrices= Group 1
 T Full 2 1
End Matrices
Matrix T  1.811 1.811
Threshold T;
Covariances A+C+E   | H@A+C _
            H@A+C | A+C+E /
Option RSiduals
End 
 

G13: FS3 (Maier) - probands
Data Ninput=2
Labels Rel1 Rel2
Ctable  2 2
 -1  -1
164  10
Begin Matrices= Group 1
 T Full 2 1
End Matrices
Matrix T  1.811 1.811
Threshold T;
Covariances A+C+E   | H@A+C _
            H@A+C | A+C+E /
Option RSiduals
End 

G14: FS3 (Maier) - controls
Data Ninput=2
Labels Rel1 Rel2
Ctable  2 2
303   6
 -1  -1
Begin Matrices= Group 1
 T Full 2 1
End Matrices
Matrix T  1.811 1.811
Threshold T;
Covariances A+C+E   | H@A+C _
            H@A+C | A+C+E /
Option RSiduals
End 
 

G15: FS4 (Horwath) - probands
Data Ninput=2
Labels Rel1 Rel2
Ctable  2 2
 -1  -1
530  53
Begin Matrices= Group 1
 T Full 2 1
End Matrices
Matrix T  1.811 1.811
Threshold T;
Covariances A+C+E   | H@A+C _
            H@A+C | A+C+E /
Option RSiduals
End 

G16: FS4 (Horwath) - controls
Data Ninput=2
Labels Rel1 Rel2
Ctable  2 2
252   3
 -1  -1
Begin Matrices= Group 1
 T Full 2 1
End Matrices
Matrix T  1.811 1.811
Threshold T;
Covariances A+C+E   | H@A+C _
            H@A+C | A+C+E /
Option RSiduals
End 
 

G17: FS5 (Fyer) - probands
Data Ninput=2
Labels Rel1 Rel2
Ctable  2 2
 -1  -1
212  24
Begin Matrices= Group 1
 T Full 2 1
End Matrices
Matrix T  1.811 1.811
Threshold T;
Covariances A+C+E   | H@A+C _
            H@A+C | A+C+E /
Option RSiduals
End 

G18: FS5 (Fyer) - controls
Data Ninput=2
Labels Rel1 Rel2
Ctable  2 2
369  11
 -1  -1
Begin Matrices= Group 1
 T Full 2 1
End Matrices
Matrix T  1.811 1.811
Threshold T;
Covariances A+C+E   | H@A+C _
            H@A+C | A+C+E /
Option RSiduals
End 
 
! ** Fit model equating family resemblance between TS and FS (S=0) **
Drop 5
End

!***** NOW FIT AE SUB-MODEL ********
Drop 2
End