** Mx startup successful ** **MX-PC 1.54** Job started on 03/03/03 at 16:08:44 ! Correlations for lipid levels: 3 DZ and 1 MZ group The following MX script lines were read for group 1 #DEFINE NVAR 3 #NGROUPS 4 Note: #NGroup set number of groups to 4 G1: DZ IBD2 TWINS DATA NINPUT=23 RECTANGULAR FILE=PARTIONEDADULTDUTCH3.DAT Rectangular continuous data read initiated NOTE: Rectangular file contained 117 records with data that contained a total of 2681 observations MISSING=. NOTE: Missing value *string* set to '.' 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; NOTE: Select if yields 19 data vectors for analysis NOTE: Vectors contain a total of 433 observations SELECT LDL3 APOB3 LNAPOE3 LDL4 APOB4 LNAPOE4 ; MATRICES NOTE: Selection yields 19 data vectors for analysis NOTE: Vectors contain a total of 112 observations 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 The following MX script lines were read for group 2 G2: DZ IBD1 TWINS DATA NINPUT=23 RECTANGULAR FILE= PARTIONEDADULTDUTCH3.DAT Rectangular continuous data read initiated NOTE: Rectangular file contained 117 records with data that contained a total of 2681 observations MISSING=. NOTE: Missing value *string* set to '.' 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; NOTE: Select if yields 67 data vectors for analysis NOTE: Vectors contain a total of 1535 observations SELECT LDL3 APOB3 LNAPOE3 LDL4 APOB4 LNAPOE4 ; MATRICES NOTE: Selection yields 67 data vectors for analysis NOTE: Vectors contain a total of 398 observations 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 The following MX script lines were read for group 3 G3: DZ IBD0 TWINS DATA NINPUT=23 RECTANGULAR FILE= PARTIONEDADULTDUTCH3.DAT Rectangular continuous data read initiated NOTE: Rectangular file contained 117 records with data that contained a total of 2681 observations MISSING=. NOTE: Missing value *string* set to '.' 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; NOTE: Select if yields 31 data vectors for analysis NOTE: Vectors contain a total of 713 observations SELECT LDL3 APOB3 LNAPOE3 LDL4 APOB4 LNAPOE4 ; MATRICES NOTE: Selection yields 31 data vectors for analysis NOTE: Vectors contain a total of 186 observations 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 The following MX script lines were read for group 4 G4: MZ TWINS DATA NINPUT=14 RECTANGULAR FILE=ADULTDUTCHMZ3.DAT Rectangular continuous data read initiated NOTE: Rectangular file contained 92 records with data that contained a total of 1284 observations MISSING=. NOTE: Missing value *string* set to '.' 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 NOTE: Selection yields 92 data vectors for analysis NOTE: Vectors contain a total of 548 observations 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 Summary of VL file data for group 1 LDL3 APOB3 LNAPOE3 LDL4 APOB4 LNAPOE4 Code 1.0000 2.0000 3.0000 4.0000 5.0000 6.0000 Number 19.0000 19.0000 18.0000 19.0000 19.0000 18.0000 Mean 3.5079 1.1695 10.4689 3.4826 1.1889 9.6817 Variance 1.1139 0.1154 6.4771 1.3891 0.1916 5.9205 Minimum 0.5800 0.4100 5.6700 0.6000 0.4100 5.0000 Maximum 5.6300 1.7200 14.3500 6.2700 2.3000 13.9700 Summary of VL file data for group 2 LDL3 APOB3 LNAPOE3 LDL4 APOB4 LNAPOE4 Code 1.0000 2.0000 3.0000 4.0000 5.0000 6.0000 Number 66.0000 67.0000 66.0000 66.0000 67.0000 66.0000 Mean 3.5670 1.2099 11.2700 3.6852 1.2081 11.3709 Variance 0.9857 0.1106 8.3283 1.2128 0.1660 5.7552 Minimum 2.0300 0.6800 3.4400 1.4100 0.4800 5.6700 Maximum 5.6500 1.9600 20.9000 6.3900 2.2500 16.5400 Summary of VL file data for group 3 LDL3 APOB3 LNAPOE3 LDL4 APOB4 LNAPOE4 Code 1.0000 2.0000 3.0000 4.0000 5.0000 6.0000 Number 31.0000 31.0000 31.0000 31.0000 31.0000 31.0000 Mean 3.6829 1.2539 11.0058 3.7277 1.2313 11.0726 Variance 0.6123 0.0602 5.8027 0.5042 0.0850 4.7699 Minimum 2.2500 0.6900 6.2800 1.8900 0.6200 5.0000 Maximum 5.3900 1.7800 16.5400 5.2100 1.9200 15.0500 Summary of VL file data for group 4 LDL3 APOB3 LNAPOE3 LDL4 APOB4 LNAPOE4 Code 1.0000 2.0000 3.0000 4.0000 5.0000 6.0000 Number 91.0000 91.0000 90.0000 92.0000 92.0000 92.0000 Mean 3.6877 1.2524 11.0157 3.7190 1.2485 10.9265 Variance 0.9640 0.1097 7.3560 0.7488 0.0998 8.3479 Minimum 1.6600 0.4800 4.2600 1.7300 0.6200 0.1500 Maximum 5.9800 2.1400 17.4200 5.9600 2.0400 17.9900 PARAMETER SPECIFICATIONS GROUP NUMBER: 1 G1: DZ IBD2 twins MATRIX M This is a FULL matrix of order 6 by 1 1 1 22 2 23 3 24 4 25 5 26 6 27 MATRIX R This is a STANDARDISED matrix of order 6 by 6 1 2 3 4 5 6 1 0 2 7 0 3 8 9 0 4 10 11 12 0 5 13 14 15 16 0 6 17 18 19 20 21 0 MATRIX S This is a DIAGONAL matrix of order 6 by 6 1 2 3 4 5 6 1 1 2 0 2 3 0 0 3 4 0 0 0 4 5 0 0 0 0 5 6 0 0 0 0 0 6 GROUP NUMBER: 2 G2: DZ IBD1 twins MATRIX M This is a FULL matrix of order 6 by 1 1 1 49 2 50 3 51 4 52 5 53 6 54 MATRIX R This is a STANDARDISED matrix of order 6 by 6 1 2 3 4 5 6 1 0 2 34 0 3 35 36 0 4 37 38 39 0 5 40 41 42 43 0 6 44 45 46 47 48 0 MATRIX S This is a DIAGONAL matrix of order 6 by 6 1 2 3 4 5 6 1 28 2 0 29 3 0 0 30 4 0 0 0 31 5 0 0 0 0 32 6 0 0 0 0 0 33 GROUP NUMBER: 3 G3: DZ IBD0 twins MATRIX M This is a FULL matrix of order 6 by 1 1 1 76 2 77 3 78 4 79 5 80 6 81 MATRIX R This is a STANDARDISED matrix of order 6 by 6 1 2 3 4 5 6 1 0 2 61 0 3 62 63 0 4 64 65 66 0 5 67 68 69 70 0 6 71 72 73 74 75 0 MATRIX S This is a DIAGONAL matrix of order 6 by 6 1 2 3 4 5 6 1 55 2 0 56 3 0 0 57 4 0 0 0 58 5 0 0 0 0 59 6 0 0 0 0 0 60 GROUP NUMBER: 4 G4: MZ twins MATRIX M This is a FULL matrix of order 6 by 1 1 1 103 2 104 3 105 4 106 5 107 6 108 MATRIX R This is a STANDARDISED matrix of order 6 by 6 1 2 3 4 5 6 1 0 2 88 0 3 89 90 0 4 91 92 93 0 5 94 95 96 97 0 6 98 99 100 101 102 0 MATRIX S This is a DIAGONAL matrix of order 6 by 6 1 2 3 4 5 6 1 82 2 0 83 3 0 0 84 4 0 0 0 85 5 0 0 0 0 86 6 0 0 0 0 0 87 Mx starting optimization; number of parameters = 108 MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1: DZ IBD2 twins MATRIX M This is a FULL matrix of order 6 by 1 1 1 3.5080 2 1.1695 3 10.5837 4 3.4827 5 1.1890 6 9.7275 MATRIX R This is a STANDARDISED matrix of order 6 by 6 1 2 3 4 5 6 1 1.0000 2 0.9004 1.0000 3 0.3961 0.1954 1.0000 4 0.8125 0.7722 0.0826 1.0000 5 0.6504 0.6387 -0.0391 0.9384 1.0000 6 0.4615 0.2625 0.8343 0.3225 0.3182 1.0000 MATRIX S This is a DIAGONAL matrix of order 6 by 6 1 2 3 4 5 6 1 1.0554 2 0.0000 0.3398 3 0.0000 0.0000 2.5637 4 0.0000 0.0000 0.0000 1.1787 5 0.0000 0.0000 0.0000 0.0000 0.4377 6 0.0000 0.0000 0.0000 0.0000 0.0000 2.4128 GROUP NUMBER: 2 G2: DZ IBD1 twins MATRIX M This is a FULL matrix of order 6 by 1 1 1 3.5702 2 1.2099 3 11.2623 4 3.6848 5 1.2081 6 11.3604 MATRIX R This is a STANDARDISED matrix of order 6 by 6 1 2 3 4 5 6 1 1.0000 2 0.8857 1.0000 3 0.0497 0.0222 1.0000 4 0.4967 0.5139 -0.0526 1.0000 5 0.4394 0.5065 0.0052 0.9423 1.0000 6 0.0867 -0.0162 0.5555 0.1487 0.1214 1.0000 MATRIX S This is a DIAGONAL matrix of order 6 by 6 1 2 3 4 5 6 1 0.9871 2 0.0000 0.3325 3 0.0000 0.0000 2.8856 4 0.0000 0.0000 0.0000 1.0938 5 0.0000 0.0000 0.0000 0.0000 0.4074 6 0.0000 0.0000 0.0000 0.0000 0.0000 2.3992 GROUP NUMBER: 3 G3: DZ IBD0 twins MATRIX M This is a FULL matrix of order 6 by 1 1 1 3.6829 2 1.2539 3 11.0054 4 3.7277 5 1.2313 6 11.0723 MATRIX R This is a STANDARDISED matrix of order 6 by 6 1 2 3 4 5 6 1 1.0000 2 0.6901 1.0000 3 0.1387 -0.0626 1.0000 4 -0.2163 -0.0951 -0.0148 1.0000 5 -0.2057 0.0243 -0.0804 0.8514 1.0000 6 -0.3176 -0.1925 0.1436 0.1725 0.3250 1.0000 MATRIX S This is a DIAGONAL matrix of order 6 by 6 1 2 3 4 5 6 1 0.7826 2 0.0000 0.2454 3 0.0000 0.0000 2.4088 4 0.0000 0.0000 0.0000 0.7101 5 0.0000 0.0000 0.0000 0.0000 0.2915 6 0.0000 0.0000 0.0000 0.0000 0.0000 2.1840 GROUP NUMBER: 4 G4: MZ twins MATRIX M This is a FULL matrix of order 6 by 1 1 1 3.6823 2 1.2499 3 10.8799 4 3.7190 5 1.2485 6 10.9262 MATRIX R This is a STANDARDISED matrix of order 6 by 6 1 2 3 4 5 6 1 1.0000 2 0.8917 1.0000 3 0.4667 0.4738 1.0000 4 0.7841 0.7753 0.4917 1.0000 5 0.7028 0.7888 0.4146 0.8850 1.0000 6 0.3477 0.3565 0.8859 0.5113 0.4375 1.0000 MATRIX S This is a DIAGONAL matrix of order 6 by 6 1 2 3 4 5 6 1 0.9798 2 0.0000 0.3308 3 0.0000 0.0000 2.8816 4 0.0000 0.0000 0.0000 0.8654 5 0.0000 0.0000 0.0000 0.0000 0.3160 6 0.0000 0.0000 0.0000 0.0000 0.0000 2.8894 Your model has 108 estimated parameters and 1244 Observed statistics -2 times log-likelihood of data >>> 2243.422 Degrees of freedom >>>>>>>>>>>>>>>> 1136 This problem used 1.2% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 3.27 Execution 0: 0: 1:44.38 TOTAL 0: 0: 1:47.65 Total number of warnings issued: 0 ______________________________________________________________________________ !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