!Correlations for 3 DZ IBD and MZ groups !'Partitioned Twin Analysis' #define nvar 1 #NGroups 4 G1: DZ IBD2 twins Data NInput=13 Rectangular File=marker.dat Labels famid indid1 indid2 pibd0 pibd1 pibd2 zygosity plt11 plt12 plt13 plt21 plt22 plt23 Select if zygosity >2; Select if pibd2 >.80; Select plt11 plt21; Matrices = Group 1 S Diag 2 2 free R Stnd 2 2 free M Full 2 1 free End Matrices; Matrix S 1 1 Matrix R .5 Matrix M 3 3 Means M; Covariance S*R*S'; End G2: DZ IBD1 twins Data NInput=13 Rectangular File=marker.dat Labels famid indid1 indid2 pibd0 pibd1 pibd2 zygosity plt11 plt12 plt13 plt21 plt22 plt23 Select if zygosity >2; Select if pibd1 >.80; Select plt11 plt21; Matrices S Diag 2 2 free R Stnd 2 2 free M Full 2 1 free End Matrices; Matrix S 1 1 Matrix R .5 Matrix M 3 3 Means M; Covariance S*R*S'; End G3: DZ IBD0 twins Data NInput=13 Rectangular File=marker.dat Labels famid indid1 indid2 pibd0 pibd1 pibd2 zygosity plt11 plt12 plt13 plt21 plt22 plt23 Select if zygosity >2; Select if pibd0 >.80; Select plt11 plt21; Matrices S Diag 2 2 free R Stnd 2 2 free M Full 2 1 free End Matrices; Matrix S 1 1 Matrix R .5 Matrix M 3 3 Means M; Covariance S*R*S'; End G4: MZ twins Data NInput=13 Rectangular File=marker.dat Labels famid indid1 indid2 pibd0 pibd1 pibd2 zygosity plt11 plt12 plt13 plt21 plt22 plt23 Select if zygosity <3; Select plt11 plt21; Matrices S Diag 2 2 free R Stnd 2 2 free M Full 2 1 free End Matrices; Matrix S 1 1 Matrix R .5 Matrix M 3 3 Means M; Covariance S*R*S'; Option Multiple Issat End Save cor.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 cor.mxs ! Test for residual polygenic variance ! Set DZ IBD2 correlation equal to MZ correlation Equate R 1 2 1 R 4 2 1 End