** Mx startup successful ** **MX-PC 1.50** Job started on 03/06/01 at 13:53:19 !Correlations for 3 DZ IBD and MZ groups !'Partitioned Twin Analysis' The following MX script lines were read for group 1 #DEFINE NVAR 1 #NGROUPS 4 Note: #NGroup set number of groups to 4 G1: DZ IBD2 TWINS DATA NINPUT=13 RECTANGULAR FILE=MARKER.DAT NOTE: Rectangular file contained 519 records with data LABELS FAMID INDID1 INDID2 PIBD0 PIBD1 PIBD2 ZYGOSITY PLT11 PLT12 PLT13 PLT21 PLT22 PLT23 SELECT IF ZYGOSITY >2; NOTE: Select if yields 232 data vectors for analysis NOTE: Vectors contain a total of 2662 observations SELECT IF PIBD2 >.80; NOTE: Select if yields 43 data vectors for analysis NOTE: Vectors contain a total of 500 observations SELECT PLT11 PLT21; MATRICES = GROUP 1 NOTE: Selection yields 38 data vectors for analysis NOTE: Vectors contain a total of 75 observations 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 The following MX script lines were read for group 2 G2: DZ IBD1 TWINS DATA NINPUT=13 RECTANGULAR FILE=MARKER.DAT NOTE: Rectangular file contained 519 records with data LABELS FAMID INDID1 INDID2 PIBD0 PIBD1 PIBD2 ZYGOSITY PLT11 PLT12 PLT13 PLT21 PLT22 PLT23 SELECT IF ZYGOSITY >2; NOTE: Select if yields 232 data vectors for analysis NOTE: Vectors contain a total of 2662 observations SELECT IF PIBD1 >.80; NOTE: Select if yields 67 data vectors for analysis NOTE: Vectors contain a total of 771 observations SELECT PLT11 PLT21; MATRICES NOTE: Selection yields 58 data vectors for analysis NOTE: Vectors contain a total of 114 observations 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 The following MX script lines were read for group 3 G3: DZ IBD0 TWINS DATA NINPUT=13 RECTANGULAR FILE=MARKER.DAT NOTE: Rectangular file contained 519 records with data LABELS FAMID INDID1 INDID2 PIBD0 PIBD1 PIBD2 ZYGOSITY PLT11 PLT12 PLT13 PLT21 PLT22 PLT23 SELECT IF ZYGOSITY >2; NOTE: Select if yields 232 data vectors for analysis NOTE: Vectors contain a total of 2662 observations SELECT IF PIBD0 >.80; NOTE: Select if yields 37 data vectors for analysis NOTE: Vectors contain a total of 419 observations SELECT PLT11 PLT21; MATRICES NOTE: Selection yields 32 data vectors for analysis NOTE: Vectors contain a total of 64 observations 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 The following MX script lines were read for group 4 G4: MZ TWINS DATA NINPUT=13 RECTANGULAR FILE=MARKER.DAT NOTE: Rectangular file contained 519 records with data LABELS FAMID INDID1 INDID2 PIBD0 PIBD1 PIBD2 ZYGOSITY PLT11 PLT12 PLT13 PLT21 PLT22 PLT23 SELECT IF ZYGOSITY <3; NOTE: Select if yields 287 data vectors for analysis NOTE: Vectors contain a total of 2962 observations SELECT PLT11 PLT21; MATRICES NOTE: Selection yields 166 data vectors for analysis NOTE: Vectors contain a total of 328 observations 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 Summary of VL file data for group 1 PLT11 PLT21 Code 1.0000 2.0000 Number 38.0000 37.0000 Mean 2.7950 2.7276 Variance 0.2711 0.2431 Minimum 1.9200 1.7200 Maximum 4.2600 4.0200 Summary of VL file data for group 2 PLT11 PLT21 Code 1.0000 2.0000 Number 58.0000 56.0000 Mean 2.9200 2.8020 Variance 0.2255 0.2195 Minimum 1.9800 1.7700 Maximum 3.9400 3.9100 Summary of VL file data for group 3 PLT11 PLT21 Code 1.0000 2.0000 Number 32.0000 32.0000 Mean 2.8994 2.7906 Variance 0.2188 0.3680 Minimum 2.0000 1.6600 Maximum 4.2900 4.5900 Summary of VL file data for group 4 PLT11 PLT21 Code 1.0000 2.0000 Number 165.0000 163.0000 Mean 2.8489 2.8618 Variance 0.2925 0.3018 Minimum 1.8000 1.8000 Maximum 4.4800 4.5800 PARAMETER SPECIFICATIONS GROUP NUMBER: 1 G1: DZ IBD2 twins MATRIX M This is a FULL matrix of order 2 by 1 1 1 4 2 5 MATRIX R This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 3 0 MATRIX S This is a DIAGONAL matrix of order 2 by 2 1 2 1 1 2 0 2 GROUP NUMBER: 2 G2: DZ IBD1 twins MATRIX M This is a FULL matrix of order 2 by 1 1 1 9 2 10 MATRIX R This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 8 0 MATRIX S This is a DIAGONAL matrix of order 2 by 2 1 2 1 6 2 0 7 GROUP NUMBER: 3 G3: DZ IBD0 twins MATRIX M This is a FULL matrix of order 2 by 1 1 1 14 2 15 MATRIX R This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 13 0 MATRIX S This is a DIAGONAL matrix of order 2 by 2 1 2 1 11 2 0 12 GROUP NUMBER: 4 G4: MZ twins MATRIX M This is a FULL matrix of order 2 by 1 1 1 19 2 20 MATRIX R This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 18 0 MATRIX S This is a DIAGONAL matrix of order 2 by 2 1 2 1 16 2 0 17 MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1: DZ IBD2 twins MATRIX M This is a FULL matrix of order 2 by 1 1 1 2.7950 2 2.7285 MATRIX R This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.6811 1.0000 MATRIX S This is a DIAGONAL matrix of order 2 by 2 1 2 1 0.5206 2 0.0000 0.4900 GROUP NUMBER: 2 G2: DZ IBD1 twins MATRIX M This is a FULL matrix of order 2 by 1 1 1 2.9200 2 2.8113 MATRIX R This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.3424 1.0000 MATRIX S This is a DIAGONAL matrix of order 2 by 2 1 2 1 0.4748 2 0.0000 0.4705 GROUP NUMBER: 3 G3: DZ IBD0 twins MATRIX M This is a FULL matrix of order 2 by 1 1 1 2.8994 2 2.7906 MATRIX R This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.1797 1.0000 MATRIX S This is a DIAGONAL matrix of order 2 by 2 1 2 1 0.4677 2 0.0000 0.6066 GROUP NUMBER: 4 G4: MZ twins MATRIX M This is a FULL matrix of order 2 by 1 1 1 2.8495 2 2.8677 MATRIX R This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.8508 1.0000 MATRIX S This is a DIAGONAL matrix of order 2 by 2 1 2 1 0.5397 2 0.0000 0.5504 Your model has 20 estimated parameters and 581 Observed statistics -2 times log-likelihood of data >>> 657.202 Degrees of freedom >>>>>>>>>>>>>>>> 561 This problem used 26.8% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 2.58 Execution 0: 0: 0: 1.97 TOTAL 0: 0: 0: 4.55 Total number of warnings issued: 0 ______________________________________________________________________________ Multiple fit option in effect. The following MX script lines have been read: 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 Summary of VL file data for group 1 PLT11 PLT21 Code 1.0000 2.0000 Number 38.0000 37.0000 Mean 2.7950 2.7276 Variance 0.2711 0.2431 Minimum 1.9200 1.7200 Maximum 4.2600 4.0200 Summary of VL file data for group 2 PLT11 PLT21 Code 1.0000 2.0000 Number 58.0000 56.0000 Mean 2.9200 2.8020 Variance 0.2255 0.2195 Minimum 1.9800 1.7700 Maximum 3.9400 3.9100 Summary of VL file data for group 3 PLT11 PLT21 Code 1.0000 2.0000 Number 32.0000 32.0000 Mean 2.8994 2.7906 Variance 0.2188 0.3680 Minimum 2.0000 1.6600 Maximum 4.2900 4.5900 Summary of VL file data for group 4 PLT11 PLT21 Code 1.0000 2.0000 Number 165.0000 163.0000 Mean 2.8489 2.8618 Variance 0.2925 0.3018 Minimum 1.8000 1.8000 Maximum 4.4800 4.5800 PARAMETER SPECIFICATIONS GROUP NUMBER: 1 G1: DZ IBD2 twins MATRIX M This is a FULL matrix of order 2 by 1 1 1 4 2 5 MATRIX R This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 3 0 MATRIX S This is a DIAGONAL matrix of order 2 by 2 1 2 1 1 2 0 2 GROUP NUMBER: 2 G2: DZ IBD1 twins MATRIX M This is a FULL matrix of order 2 by 1 1 1 9 2 10 MATRIX R This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 3 0 MATRIX S This is a DIAGONAL matrix of order 2 by 2 1 2 1 6 2 0 7 GROUP NUMBER: 3 G3: DZ IBD0 twins MATRIX M This is a FULL matrix of order 2 by 1 1 1 14 2 15 MATRIX R This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 3 0 MATRIX S This is a DIAGONAL matrix of order 2 by 2 1 2 1 11 2 0 12 GROUP NUMBER: 4 G4: MZ twins MATRIX M This is a FULL matrix of order 2 by 1 1 1 19 2 20 MATRIX R This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 18 0 MATRIX S This is a DIAGONAL matrix of order 2 by 2 1 2 1 16 2 0 17 MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1: DZ IBD2 twins MATRIX M This is a FULL matrix of order 2 by 1 1 1 2.7950 2 2.7282 MATRIX R This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.4249 1.0000 MATRIX S This is a DIAGONAL matrix of order 2 by 2 1 2 1 0.4851 2 0.0000 0.4576 GROUP NUMBER: 2 G2: DZ IBD1 twins MATRIX M This is a FULL matrix of order 2 by 1 1 1 2.9200 2 2.8136 MATRIX R This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.4249 1.0000 MATRIX S This is a DIAGONAL matrix of order 2 by 2 1 2 1 0.4840 2 0.0000 0.4809 GROUP NUMBER: 3 G3: DZ IBD0 twins MATRIX M This is a FULL matrix of order 2 by 1 1 1 2.8994 2 2.7906 MATRIX R This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.4249 1.0000 MATRIX S This is a DIAGONAL matrix of order 2 by 2 1 2 1 0.4966 2 0.0000 0.6441 GROUP NUMBER: 4 G4: MZ twins MATRIX M This is a FULL matrix of order 2 by 1 1 1 2.8495 2 2.8677 MATRIX R This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.8508 1.0000 MATRIX S This is a DIAGONAL matrix of order 2 by 2 1 2 1 0.5397 2 0.0000 0.5504 Your model has 18 estimated parameters and 581 Observed statistics -2 times log-likelihood of data >>> 665.288 Degrees of freedom >>>>>>>>>>>>>>>> 563 Saturated model fit* >>>>>>>>>>> 657.202 Saturated model df* >>>>>>>>>>> 561 Difference Chi-squared >>>>>>>> 8.086 Difference d.f. >>>>>>>>>>>>>>> 2 Probability >>>>>>>>>>>>>>>>>>>> 0.018 Akaike's Information Criterion > 4.086 * Saturated model statistic computed earlier in this job This problem used 26.8% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.33 Execution 0: 0: 0: 0.99 TOTAL 0: 0: 0: 1.32 Total number of warnings issued: 0 ______________________________________________________________________________ Multiple fit option in effect. The following MX script lines have been read: GET COR.MXS Note: 53224 binary records read ! Test for residual polygenic variance ! Set DZ IBD2 correlation equal to MZ correlation EQUATE R 1 2 1 R 4 2 1 END Summary of VL file data for group 1 PLT11 PLT21 Code 1.0000 2.0000 Number 38.0000 37.0000 Mean 2.7950 2.7276 Variance 0.2711 0.2431 Minimum 1.9200 1.7200 Maximum 4.2600 4.0200 Summary of VL file data for group 2 PLT11 PLT21 Code 1.0000 2.0000 Number 58.0000 56.0000 Mean 2.9200 2.8020 Variance 0.2255 0.2195 Minimum 1.9800 1.7700 Maximum 3.9400 3.9100 Summary of VL file data for group 3 PLT11 PLT21 Code 1.0000 2.0000 Number 32.0000 32.0000 Mean 2.8994 2.7906 Variance 0.2188 0.3680 Minimum 2.0000 1.6600 Maximum 4.2900 4.5900 Summary of VL file data for group 4 PLT11 PLT21 Code 1.0000 2.0000 Number 165.0000 163.0000 Mean 2.8489 2.8618 Variance 0.2925 0.3018 Minimum 1.8000 1.8000 Maximum 4.4800 4.5800 PARAMETER SPECIFICATIONS GROUP NUMBER: 1 G1: DZ IBD2 twins MATRIX M This is a FULL matrix of order 2 by 1 1 1 4 2 5 MATRIX R This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 3 0 MATRIX S This is a DIAGONAL matrix of order 2 by 2 1 2 1 1 2 0 2 GROUP NUMBER: 2 G2: DZ IBD1 twins MATRIX M This is a FULL matrix of order 2 by 1 1 1 9 2 10 MATRIX R This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 8 0 MATRIX S This is a DIAGONAL matrix of order 2 by 2 1 2 1 6 2 0 7 GROUP NUMBER: 3 G3: DZ IBD0 twins MATRIX M This is a FULL matrix of order 2 by 1 1 1 14 2 15 MATRIX R This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 13 0 MATRIX S This is a DIAGONAL matrix of order 2 by 2 1 2 1 11 2 0 12 GROUP NUMBER: 4 G4: MZ twins MATRIX M This is a FULL matrix of order 2 by 1 1 1 19 2 20 MATRIX R This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 3 0 MATRIX S This is a DIAGONAL matrix of order 2 by 2 1 2 1 16 2 0 17 MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1: DZ IBD2 twins MATRIX M This is a FULL matrix of order 2 by 1 1 1 2.7950 2 2.7287 MATRIX R This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.8277 1.0000 MATRIX S This is a DIAGONAL matrix of order 2 by 2 1 2 1 0.6123 2 0.0000 0.5756 GROUP NUMBER: 2 G2: DZ IBD1 twins MATRIX M This is a FULL matrix of order 2 by 1 1 1 2.9200 2 2.8113 MATRIX R This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.3424 1.0000 MATRIX S This is a DIAGONAL matrix of order 2 by 2 1 2 1 0.4748 2 0.0000 0.4705 GROUP NUMBER: 3 G3: DZ IBD0 twins MATRIX M This is a FULL matrix of order 2 by 1 1 1 2.8994 2 2.7906 MATRIX R This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.1797 1.0000 MATRIX S This is a DIAGONAL matrix of order 2 by 2 1 2 1 0.4677 2 0.0000 0.6066 GROUP NUMBER: 4 G4: MZ twins MATRIX M This is a FULL matrix of order 2 by 1 1 1 2.8495 2 2.8676 MATRIX R This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.8277 1.0000 MATRIX S This is a DIAGONAL matrix of order 2 by 2 1 2 1 0.5234 2 0.0000 0.5335 Your model has 19 estimated parameters and 581 Observed statistics -2 times log-likelihood of data >>> 662.730 Degrees of freedom >>>>>>>>>>>>>>>> 562 Saturated model fit* >>>>>>>>>>> 657.202 Saturated model df* >>>>>>>>>>> 561 Difference Chi-squared >>>>>>>> 5.528 Difference d.f. >>>>>>>>>>>>>>> 1 Probability >>>>>>>>>>>>>>>>>>>> 0.019 Akaike's Information Criterion > 3.528 * Saturated model statistic computed earlier in this job This problem used 26.8% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.06 Execution 0: 0: 0: 0.99 TOTAL 0: 0: 0: 1.05 Total number of warnings issued: 0 ______________________________________________________________________________