**MX-PC 1.65b** Job started on 03/05/07 at 16:05:15 ! ESTIMATE CORRELATIONS FOR DZ TWINS BY IBD STATUS ! BRISBANE TWIN AND FAMILY STUDY OF COGNITION #DEFINE NVAR 1 #DEFINE NVARX2 2 The following MX script lines were read for group 1 #NGROUPS 3 Note: #NGroup set number of groups to 3 G1: DZ IBD2 TWINS DATA NINPUT=10 RECTANGULAR FILE=PIQDZ.REC Rectangular continuous data read initiated Note: Maximum ordinal/rectangular record length is: 1000 NOTE: Rectangular file contained 710 records with data that contained a total of 6536 observations LABELS FAM ID1 ID2 PIQ1 PIQ2 IBD0MNR IBD1MNR IBD2MNR PIHAT PICAT SELECT IF PICAT =2; NOTE: Select if yields 184 data vectors for analysis NOTE: Vectors contain a total of 1695 observations SELECT PIQ1 PIQ2 ; BEGIN MATRICES; NOTE: Selection yields 142 data vectors for analysis NOTE: Vectors contain a total of 223 observations M FULL NVAR NVARX2 FREE ! MEANS S DIAG NVARX2 NVARX2 FREE ! STANDARD DEVIATIONS R STND NVARX2 NVARX2 FREE ! CORRELATIONS END MATRICES; MATRIX M 110 110 ! STARTING VALUES MATRIX S 15 15 MATRIX R .5 MEANS M; COVARIANCE S*R*S'; OPTION RSIDUALS END The following MX script lines were read for group 2 G2: DZ IBD1 TWINS DATA NINPUT=10 RECTANGULAR FILE=PIQDZ.REC Rectangular continuous data read initiated Note: Maximum ordinal/rectangular record length is: 1000 NOTE: Rectangular file contained 710 records with data that contained a total of 6536 observations LABELS FAM ID1 ID2 PIQ1 PIQ2 IBD0MNR IBD1MNR IBD2MNR PIHAT PICAT SELECT IF PICAT =1; NOTE: Select if yields 366 data vectors for analysis NOTE: Vectors contain a total of 3350 observations SELECT PIQ1 PIQ2 ; BEGIN MATRICES; NOTE: Selection yields 278 data vectors for analysis NOTE: Vectors contain a total of 422 observations M FULL NVAR NVARX2 FREE ! MEANS S DIAG NVARX2 NVARX2 FREE ! STANDARD DEVIATIONS R STND NVARX2 NVARX2 FREE ! CORRELATIONS END MATRICES; MATRIX M 110 110 ! STARTING VALUES MATRIX S 15 15 MATRIX R .5 MEANS M; COVARIANCE S*R*S'; OPTION RSIDUALS END The following MX script lines were read for group 3 G3: DZ IBD0 TWINS DATA NINPUT=10 RECTANGULAR FILE=PIQDZ.REC Rectangular continuous data read initiated Note: Maximum ordinal/rectangular record length is: 1000 NOTE: Rectangular file contained 710 records with data that contained a total of 6536 observations LABELS FAM ID1 ID2 PIQ1 PIQ2 IBD0MNR IBD1MNR IBD2MNR PIHAT PICAT SELECT IF PICAT =0; NOTE: Select if yields 160 data vectors for analysis NOTE: Vectors contain a total of 1491 observations SELECT PIQ1 PIQ2 ; BEGIN MATRICES; NOTE: Selection yields 130 data vectors for analysis NOTE: Vectors contain a total of 211 observations M FULL NVAR NVARX2 FREE ! MEANS S DIAG NVARX2 NVARX2 FREE ! STANDARD DEVIATIONS R STND NVARX2 NVARX2 FREE ! CORRELATIONS END MATRICES; MATRIX M 110 110 ! STARTING VALUES MATRIX S 15 15 MATRIX R .5 MEANS M; COVARIANCE S*R*S'; OPTION RSIDUALS MX%P=PIQ.IL END Summary of VL file data for group 1 PIQ1 PIQ2 Code 1.0000 2.0000 Number 142.0000 81.0000 Mean 110.2817 112.3333 Variance 299.4136 282.3457 Minimum 68.0000 77.0000 Maximum 149.0000 150.0000 Summary of VL file data for group 2 PIQ1 PIQ2 Code 1.0000 2.0000 Number 276.0000 146.0000 Mean 112.3768 112.0753 Variance 259.3942 262.0423 Minimum 73.0000 68.0000 Maximum 149.0000 149.0000 Summary of VL file data for group 3 PIQ1 PIQ2 Code 1.0000 2.0000 Number 129.0000 82.0000 Mean 112.8527 109.7317 Variance 264.4202 309.4646 Minimum 77.0000 68.0000 Maximum 149.0000 143.0000 PARAMETER SPECIFICATIONS GROUP NUMBER: 1 G1: DZ IBD2 twins MATRIX M This is a FULL matrix of order 1 by 2 1 2 1 1 2 MATRIX R This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 5 0 MATRIX S This is a DIAGONAL matrix of order 2 by 2 1 2 1 3 2 0 4 GROUP NUMBER: 2 G2: DZ IBD1 twins MATRIX M This is a FULL matrix of order 1 by 2 1 2 1 6 7 MATRIX R This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 10 0 MATRIX S This is a DIAGONAL matrix of order 2 by 2 1 2 1 8 2 0 9 GROUP NUMBER: 3 G3: DZ IBD0 twins MATRIX M This is a FULL matrix of order 1 by 2 1 2 1 11 12 MATRIX R This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 15 0 MATRIX S This is a DIAGONAL matrix of order 2 by 2 1 2 1 13 2 0 14 Mx starting optimization; number of parameters = 15 MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1: DZ IBD2 twins MATRIX M This is a FULL matrix of order 1 by 2 1 2 1 110.2803 112.0023 MATRIX R This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.6004 1.0000 MATRIX S This is a DIAGONAL matrix of order 2 by 2 1 2 1 17.3032 2 0.0000 16.7683 Vector of OBSERVED means PIQ1 PIQ2 Mean 110.2817 112.3333 Vector of EXPECTED means PIQ1 PIQ2 Mean 110.2803 112.0023 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX PIQ1 PIQ2 PIQ1 299.4014 PIQ2 174.1952 281.1767 Function value of this group: 1863.0569 Where the fit function is -2 * Log-likelihood of raw data GROUP NUMBER: 2 G2: DZ IBD1 twins MATRIX M This is a FULL matrix of order 1 by 2 1 2 1 112.3779 111.8621 MATRIX R This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.2745 1.0000 MATRIX S This is a DIAGONAL matrix of order 2 by 2 1 2 1 16.1047 2 0.0000 16.2535 Vector of OBSERVED means PIQ1 PIQ2 Mean 112.3768 112.0753 Vector of EXPECTED means PIQ1 PIQ2 Mean 112.3779 111.8621 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX PIQ1 PIQ2 PIQ1 259.3604 PIQ2 71.8507 264.1756 Function value of this group: 3534.5464 Where the fit function is -2 * Log-likelihood of raw data GROUP NUMBER: 3 G3: DZ IBD0 twins Vector likelihood statistics written to piq.il Note: R code to read this file will be written to: piq.il.R MATRIX M This is a FULL matrix of order 1 by 2 1 2 1 112.8547 109.9328 MATRIX R This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.1559 1.0000 MATRIX S This is a DIAGONAL matrix of order 2 by 2 1 2 1 16.2575 2 0.0000 17.5953 Vector of OBSERVED means PIQ1 PIQ2 Mean 112.8527 109.7317 Vector of EXPECTED means PIQ1 PIQ2 Mean 112.8547 109.9328 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX PIQ1 PIQ2 PIQ1 264.3064 PIQ2 44.5876 309.5960 Function value of this group: 1786.5772 Where the fit function is -2 * Log-likelihood of raw data Your model has 15 estimated parameters and 856 Observed statistics -2 times log-likelihood of data >>> 7184.180 Degrees of freedom >>>>>>>>>>>>>>>> 841 Akaike's Information Criterion >>>> 5502.180 Bayesian Information Criterion >>>> 938.770 Sample size Adjusted BIC >>>> 2273.615 Deviance Information Criterion >>>> 1711.597 This problem used 2.3% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.40 Execution 0: 0: 0: 1.03 TOTAL 0: 0: 0: 1.43 Total number of warnings issued: 0 ______________________________________________________________________________