! Correlations for lipid levels: 3 DZ and 1 MZ group #define nvar 3 #NGroups 4 G1: DZ IBD2 twins Data NInput=23 Rectangular File=partionedAdultDutch3.dat Missing=. Labels fam zyg1_5 nih_id3 lipidme3 sex3 age3 ldl3 apob3 apoe3 lnapoe3 nih_id4 lipidme4 sex4 age4 ldl4 apob4 apoe4 lnapoe4 z0_65 z1_65 z2_65 pihat65 pi65cat Select if pi65cat = 2; Select ldl3 apob3 lnapoe3 ldl4 apob4 lnapoe4 ; Matrices S diag 6 6 free !standard deviations R stnd 6 6 free !correlations M full 6 1 free !means End Matrices; Matrix M 3 1 9 3 1 9 !starting values Matrix S .5 .5 .9 .5 .5 .9 Matrix R .5 .4 .5 .4 .4 .5 .4 .4 .4 .5 .4 .4 .4 .4 .5 Means M / Covariance S*R*S' / End G2: DZ IBD1 twins Data NInput=23 Rectangular File= partionedAdultDutch3.dat Missing=. Labels fam zyg1_5 nih_id3 lipidme3 sex3 age3 ldl3 apob3 apoe3 lnapoe3 nih_id4 lipidme4 sex4 age4 ldl4 apob4 apoe4 lnapoe4 z0_65 z1_65 z2_65 pihat65 pi65cat Select if pi65cat = 1; Select ldl3 apob3 lnapoe3 ldl4 apob4 lnapoe4 ; Matrices S diag 6 6 free !standard deviations R stnd 6 6 free !correlations M full 6 1 free !means End Matrices; Matrix M 3 1 9 3 1 9 !starting values Matrix S .5 .5 .9 .5 .5 .9 Matrix R .5 .4 .5 .4 .4 .5 .4 .4 .4 .5 .4 .4 .4 .4 .5 Means M / Covariance S*R*S' / End G3: DZ IBD0 twins Data NInput=23 Rectangular File= partionedAdultDutch3.dat Missing=. Labels fam zyg1_5 nih_id3 lipidme3 sex3 age3 ldl3 apob3 apoe3 lnapoe3 nih_id4 lipidme4 sex4 age4 ldl4 apob4 apoe4 lnapoe4 z0_65 z1_65 z2_65 pihat65 pi65cat Select if pi65cat = 0; Select ldl3 apob3 lnapoe3 ldl4 apob4 lnapoe4 ; Matrices S diag 6 6 free !standard deviations R stnd 6 6 free !correlations M full 6 1 free !means End Matrices; Matrix M 3 1 9 3 1 9 !starting values Matrix S .5 .5 .9 .5 .5 .9 Matrix R .5 .4 .5 .4 .4 .5 .4 .4 .4 .5 .4 .4 .4 .4 .5 Means M / Covariance S*R*S' / End G4: MZ twins Data NInput=14 Rectangular File=AdultDutchMZ3.dat missing=. Labels fam zyg1_5 nih_id3 sex3 age3 ldl3 apob3 lnapoe3 nih_id4 sex4 age4 ldl4 apob4 lnapoe4 Select ldl3 apob3 lnapoe3 ldl4 apob4 lnapoe4 ; Matrices S diag 6 6 free !standard deviations R stnd 6 6 free !correlations M full 6 1 free !means End Matrices; Matrix M 3 1 9 3 1 9 !starting values Matrix S .5 .5 .5 .5 .5 .5 Matrix R .5 .4 .5 .4 .4 .5 .4 .4 .4 .5 .4 .4 .4 .4 .5 Means M / Covariance S*R*S' / ! Option Multiple Issat End !Save lipidcor.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 lipidcor.mxs ! Test for residual polygenic variance ! Set DZ IBD2 correlation equal to MZ correlation !Equate R 1 2 1 R 4 2 1 !End