! Correlations for Total mole count (ln(x+1)) #define nvar 1 #NGroups 4 G1: DZ IBD2 twins Data NInput=8 Rectangular File=moledz.dat Labels mole1 mole2 zyg sex1 sex2 p0 p1 p2 Select if p2 =1; Select mole1 mole2 ; Matrices = Group 1 S diag 2 2 free R stnd 2 2 free M full 2 1 free End Matrices; Matrix S .5 .5 Matrix R .5 Matrix M 4 4 Means M / Covariance S*R*S' / End G2: DZ IBD1 twins Data NInput=8 Rectangular File=moledz.dat Labels mole1 mole2 zyg sex1 sex2 p0 p1 p2 Select if p1 =1; Select mole1 mole2; Matrices S diag 2 2 free R stnd 2 2 free M full 2 1 free End Matrices; Matrix S .5 .5 Matrix R .5 Matrix M 4 4 Means M / Covariance S*R*S' / End G3: DZ IBD0 twins Data NInput=8 Rectangular File=moledz.dat Labels mole1 mole2 zyg sex1 sex2 p0 p1 p2 Select if p0 =1; Select mole1 mole2 ; Matrices S diag 2 2 free R stnd 2 2 free M full 2 1 free End Matrices; Matrix S .5 .5 Matrix R .5 Matrix M 4 4 Means M / Covariance S*R*S' / End G4: MZ twins Data NInput=5 Rectangular File=molemz.dat Labels mole1 mole2 zyg sex1 sex2 Select mole1 mole2 ; Matrices S diag 2 2 free R stnd 2 2 free M full 2 1 free End Matrices; Matrix S .5 .5 Matrix R .5 Matrix M 4 4 Means M / Covariance S*R*S' / Option Multiple Issat End Save molecor.mxs ! Test for linkage ! Set 3 DZ IBD correlations equal Equate R 1 2 1 R 2 2 1 R 3 2 1 End Get molecor.mxs ! Test for residual polygenic variance ! Set DZ IBD2 correlation equal to MZ correlation Equate R 1 2 1 R 4 2 1 End