** Mx startup successful ** **MX-PC 1.65b** Job started on 03/05/08 at 23:15:12 ! FEMALE NEUROTICISM AT 12, 14 AND 16 YEARS ! MULTIVARIATE CHOLESKY The following MX script lines were read for group 1 #NGROUPS 5 Note: #NGroup set number of groups to 5 #DEFINE NVAR 3 ! VARIABLES #DEFINE NSIB 2 ! TWIN-1 & TWIN-2 #DEFINE NITEM 6 ! !_____________________________________________ ! G1: MODEL PARAMATERS CALCULATION BEGIN MATRICES; X LOWER NVAR NVAR FREE Y LOWER NVAR NVAR FREE Z LOWER NVAR NVAR FREE H FULL 1 1 END MATRICES; MATRIX H 0.5 START 1.26 X 1 1 1 START 1.26 X 1 2 2 START 1.26 X 1 3 3 START -2.72 Z 1 1 1 START -2.72 Z 1 2 2 START -2.72 Z 1 3 3 BEGIN ALGEBRA; A = X*X'; C = Y*Y'; E = Z*Z'; END ALGEBRA; END The following MX script lines were read for group 2 !____________________________________ ! G2: MZ TWINS DATA NI=33 RECTANGULAR FILE=JEPQ6.DAT Rectangular continuous data read initiated Note: Maximum ordinal/rectangular record length is: 1000 NOTE: Rectangular file contained 1215 records with data that contained a total of 25304 observations LABELS ZYG SEX SEX2 AGE1 AGE2 AGE3 AGE1B AGE2B AGE3B PSYCH12 PSYCH14 PSYCH16 EXTRA12 EXTRA14 EXTRA16 NEURO12 NEURO14 NEURO16 LIE12 LIE14 LIE16 PSYCH12B PSYCH14B PSYCH16B EXTRA12B EXTRA14B EXTRA16B NEURO12B NEURO14B NEURO16B LIE12B LIE14B LIE16B SELECT IF ZYG = 1 / NOTE: Select if yields 257 data vectors for analysis NOTE: Vectors contain a total of 5316 observations SELECT NEURO12 NEURO14 NEURO16 NEURO12B NEURO14B NEURO16B / MATRICES = GROUP 1 NOTE: Selection yields 192 data vectors for analysis NOTE: Vectors contain a total of 752 observations M FULL NSIB NVAR FREE !MEAN MATRIX END MATRICES START 19.09 M 1 1 START 19.09 M 1 2 START 19.09 M 1 3 SPECIFY M 80 81 82 80 81 82 / MEANS M / ! MODEL FOR MEANS COVARIANCES (A+C+E|A+C _ A+C|A+C+E ) / OP RS END The following MX script lines were read for group 3 !_______________________ ! G3: DZ TWINS DATA NI=33 RECTANGULAR FILE=JEPQ6.DAT Rectangular continuous data read initiated Note: Maximum ordinal/rectangular record length is: 1000 NOTE: Rectangular file contained 1215 records with data that contained a total of 25304 observations LABELS ZYG SEX SEX2 AGE1 AGE2 AGE3 AGE1B AGE2B AGE3B PSYCH12 PSYCH14 PSYCH16 EXTRA12 EXTRA14 EXTRA16 NEURO12 NEURO14 NEURO16 LIE12 LIE14 LIE16 PSYCH12B PSYCH14B PSYCH16B EXTRA12B EXTRA14B EXTRA16B NEURO12B NEURO14B NEURO16B LIE12B LIE14B LIE16B SELECT IF ZYG = 3 / NOTE: Select if yields 172 data vectors for analysis NOTE: Vectors contain a total of 3772 observations SELECT NEURO12 NEURO14 NEURO16 NEURO12B NEURO14B NEURO16B / MATRICES = GROUP 1 NOTE: Selection yields 131 data vectors for analysis NOTE: Vectors contain a total of 558 observations M FULL NSIB NVAR =M2 !MEAN MATRIX END MATRICES MEANS M / ! MODEL FOR MEANS COVARIANCES (A+C+E|H@A+C_ H@A+C|A+C+E ) / OP RS END The following MX script lines were read for group 4 !_____________________________________ ! G4: DZ TWINS OPP SEX DATA NI=33 RECTANGULAR FILE=JEPQ6.DAT Rectangular continuous data read initiated Note: Maximum ordinal/rectangular record length is: 1000 NOTE: Rectangular file contained 1215 records with data that contained a total of 25304 observations LABELS ZYG SEX SEX2 AGE1 AGE2 AGE3 AGE1B AGE2B AGE3B PSYCH12 PSYCH14 PSYCH16 EXTRA12 EXTRA14 EXTRA16 NEURO12 NEURO14 NEURO16 LIE12 LIE14 LIE16 PSYCH12B PSYCH14B PSYCH16B EXTRA12B EXTRA14B EXTRA16B NEURO12B NEURO14B NEURO16B LIE12B LIE14B LIE16B SELECT IF ZYG = 5 / NOTE: Select if yields 336 data vectors for analysis NOTE: Vectors contain a total of 7131 observations SELECT NEURO12 NEURO14 NEURO16 / MATRICES = GROUP 1 NOTE: Selection yields 248 data vectors for analysis NOTE: Vectors contain a total of 514 observations M FULL 1 NVAR !MEAN MATRIX END MATRICES SPECIFY M 80 81 82 / MEANS M / ! MODEL FOR MEANS COVARIANCES (A+C+E) / OP RS END The following MX script lines were read for group 5 !_____________________________________ ! G5: CALCULATE STANDARDISED SOLUTION AND GENETIC AND ENVIRONMENTAL CORRELATIONS CALC BEGIN MATRICES= G 1; I IDEN NVAR NVAR END MATRICES; BEGIN ALGEBRA; Q=\STND(A+C+E|A+C _ A+C|A+C+E ); ! ESTIMATES PHENOTYPIC CORRELATIONS I.E STANDARDIZES THE COVARIANCE MATRIX S=A%(A+C+E)| C%(A+C+E)|E%(A+C+E); !"PROPORTION OF OBSERVED PHENOTYPIC CORRELATION EXPLAINED BY A, C AND E" T=\STND(A) ; ! CORR OF A EFFECTS U=\STND(C) ; ! CORR OF C EFFECTS How am I supposed to take the square root of 0.000000000000000E+00 ? Diagonal elements are: 0.000000000000000E+00 0.000000000000000E+00 Matrix is 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 Matrix not fully standardized V=\STND(E) ; ! CORR OF E EFFECTS END ALGEBRA; OP RS OP MU OP ND=4 OP ISSAT END How am I supposed to take the square root of 0.000000000000000E+00 ? Diagonal elements are: 0.000000000000000E+00 0.000000000000000E+00 Matrix is 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 Matrix not fully standardized Summary of VL file data for group 2 NEURO12 NEURO14 NEURO16 NEURO12B NEURO14B NEURO16B Code 1.0000 2.0000 3.0000 4.0000 5.0000 6.0000 Number 128.0000 113.0000 134.0000 130.0000 113.0000 134.0000 Mean 10.7344 10.5221 10.5597 11.2385 10.7876 10.8731 Variance 23.6326 24.0371 22.3211 23.7662 23.1938 19.2003 Minimum 1.0000 1.0000 1.0000 1.0000 0.0000 0.0000 Maximum 20.0000 20.0000 19.0000 20.0000 20.0000 19.0000 Summary of VL file data for group 3 NEURO12 NEURO14 NEURO16 NEURO12B NEURO14B NEURO16B Code 1.0000 2.0000 3.0000 4.0000 5.0000 6.0000 Number 111.0000 93.0000 75.0000 111.0000 93.0000 75.0000 Mean 11.0541 10.1290 11.5600 11.1532 9.7849 11.1467 Variance 22.3034 24.6500 18.6731 20.8865 22.9860 17.7518 Minimum 0.0000 0.0000 3.0000 0.0000 0.0000 2.0000 Maximum 20.0000 20.0000 19.0000 20.0000 20.0000 20.0000 Summary of VL file data for group 4 NEURO12 NEURO14 NEURO16 Code 1.0000 2.0000 3.0000 Number 192.0000 170.0000 152.0000 Mean 10.5781 11.0824 11.9145 Variance 23.3064 22.5697 19.4598 Minimum 1.0000 0.0000 1.0000 Maximum 20.0000 20.0000 20.0000 PARAMETER SPECIFICATIONS GROUP NUMBER: 1 G1: Model paramaters MATRIX A This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 1 2 2 3 3 4 5 6 MATRIX Y This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 7 2 8 9 3 10 11 12 MATRIX Z This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 13 2 14 15 3 16 17 18 GROUP NUMBER: 2 G2: MZ twins MATRIX A This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX M This is a FULL matrix of order 2 by 3 1 2 3 1 80 81 82 2 80 81 82 MATRIX X This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 1 2 2 3 3 4 5 6 MATRIX Y This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 7 2 8 9 3 10 11 12 MATRIX Z This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 13 2 14 15 3 16 17 18 GROUP NUMBER: 3 G3: DZ twins MATRIX A This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX M This is a FULL matrix of order 2 by 3 1 2 3 1 80 81 82 2 80 81 82 MATRIX X This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 1 2 2 3 3 4 5 6 MATRIX Y This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 7 2 8 9 3 10 11 12 MATRIX Z This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 13 2 14 15 3 16 17 18 GROUP NUMBER: 4 G4: DZ twins opp sex MATRIX A This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX M This is a FULL matrix of order 1 by 3 1 2 3 1 80 81 82 MATRIX X This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 1 2 2 3 3 4 5 6 MATRIX Y This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 7 2 8 9 3 10 11 12 MATRIX Z This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 13 2 14 15 3 16 17 18 GROUP NUMBER: 5 G5: Calculate standardised solution and genetic and environmental correlations MATRIX A This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX I This is an IDENTITY matrix of order 3 by 3 MATRIX Q This is a computed FULL matrix of order 6 by 6 It has no free parameters specified MATRIX S This is a computed FULL matrix of order 3 by 9 It has no free parameters specified MATRIX T This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX U This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX V This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 1 2 2 3 3 4 5 6 MATRIX Y This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 7 2 8 9 3 10 11 12 MATRIX Z This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 13 2 14 15 3 16 17 18 Mx starting optimization; number of parameters = 21 *** WARNING! *** I am not sure I have found a solution that satisfies Kuhn-Tucker conditions for a minimum. NAG's IFAIL parameter is 1 We probably have a minimum here, but you might consider trying different starting values. You can randomize these with TH=n on the OU line, where n is the number of times you wish to do this. I STRONGLY recommend BOundaries to be set if you use TH MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1: Model paramaters MATRIX A This is a computed FULL matrix of order 3 by 3 [=X*X'] 1 2 3 1 6.5972 6.1950 3.3051 2 6.1950 10.3771 6.1876 3 3.3051 6.1876 3.7416 MATRIX C This is a computed FULL matrix of order 3 by 3 [=Y*Y'] 1 2 3 1 0.4917 0.4641 0.9260 2 0.4641 0.4381 0.8741 3 0.9260 0.8741 1.7439 MATRIX E This is a computed FULL matrix of order 3 by 3 [=Z*Z'] 1 2 3 1 15.7841 5.6562 5.5056 2 5.6562 12.6141 5.7016 3 5.5056 5.7016 14.7120 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX X This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 2.5685 2 2.4119 2.1354 3 1.2868 1.4442 0.0000 MATRIX Y This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.7012 2 0.6619 0.0000 3 1.3206 0.0000 0.0000 MATRIX Z This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 -3.9729 2 -1.4237 -3.2538 3 -1.3858 -1.1460 -3.3880 GROUP NUMBER: 2 G2: MZ twins MATRIX A This is a computed FULL matrix of order 3 by 3 [=X*X'] 1 2 3 1 6.5972 6.1950 3.3051 2 6.1950 10.3771 6.1876 3 3.3051 6.1876 3.7416 MATRIX C This is a computed FULL matrix of order 3 by 3 [=Y*Y'] 1 2 3 1 0.4917 0.4641 0.9260 2 0.4641 0.4381 0.8741 3 0.9260 0.8741 1.7439 MATRIX E This is a computed FULL matrix of order 3 by 3 [=Z*Z'] 1 2 3 1 15.7841 5.6562 5.5056 2 5.6562 12.6141 5.7016 3 5.5056 5.7016 14.7120 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX M This is a FULL matrix of order 2 by 3 1 2 3 1 10.8625 10.4872 11.2794 2 10.8625 10.4872 11.2794 MATRIX X This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 2.5685 2 2.4119 2.1354 3 1.2868 1.4442 0.0000 MATRIX Y This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.7012 2 0.6619 0.0000 3 1.3206 0.0000 0.0000 MATRIX Z This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 -3.9729 2 -1.4237 -3.2538 3 -1.3858 -1.1460 -3.3880 Vector of OBSERVED means NEURO12 NEURO14 NEURO16 NEURO12B NEURO14B NEURO16B Mean 10.7344 10.5221 10.5597 11.2385 10.7876 10.8731 Vector of EXPECTED means NEURO12 NEURO14 NEURO16 NEURO12B NEURO14B NEURO16B Mean 10.8625 10.4872 11.2794 10.8625 10.4872 11.2794 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX NEURO12 NEURO14 NEURO16 NEURO12B NEURO14B NEURO16B NEURO12 22.8729 NEURO14 12.3152 23.4292 NEURO16 9.7367 12.7632 20.1975 NEURO12B 7.0888 6.6591 4.2311 22.8729 NEURO14B 6.6591 10.8152 7.0616 12.3152 23.4292 NEURO16B 4.2311 7.0616 5.4855 9.7367 12.7632 20.1975 Function value of this group: 4258.8432 Where the fit function is -2 * Log-likelihood of raw data GROUP NUMBER: 3 G3: DZ twins MATRIX A This is a computed FULL matrix of order 3 by 3 [=X*X'] 1 2 3 1 6.5972 6.1950 3.3051 2 6.1950 10.3771 6.1876 3 3.3051 6.1876 3.7416 MATRIX C This is a computed FULL matrix of order 3 by 3 [=Y*Y'] 1 2 3 1 0.4917 0.4641 0.9260 2 0.4641 0.4381 0.8741 3 0.9260 0.8741 1.7439 MATRIX E This is a computed FULL matrix of order 3 by 3 [=Z*Z'] 1 2 3 1 15.7841 5.6562 5.5056 2 5.6562 12.6141 5.7016 3 5.5056 5.7016 14.7120 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX M This is a FULL matrix of order 2 by 3 1 2 3 1 10.8625 10.4872 11.2794 2 10.8625 10.4872 11.2794 MATRIX X This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 2.5685 2 2.4119 2.1354 3 1.2868 1.4442 0.0000 MATRIX Y This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.7012 2 0.6619 0.0000 3 1.3206 0.0000 0.0000 MATRIX Z This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 -3.9729 2 -1.4237 -3.2538 3 -1.3858 -1.1460 -3.3880 Vector of OBSERVED means NEURO12 NEURO14 NEURO16 NEURO12B NEURO14B NEURO16B Mean 11.0541 10.1290 11.5600 11.1532 9.7849 11.1467 Vector of EXPECTED means NEURO12 NEURO14 NEURO16 NEURO12B NEURO14B NEURO16B Mean 10.8625 10.4872 11.2794 10.8625 10.4872 11.2794 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX NEURO12 NEURO14 NEURO16 NEURO12B NEURO14B NEURO16B NEURO12 22.8729 NEURO14 12.3152 23.4292 NEURO16 9.7367 12.7632 20.1975 NEURO12B 3.7903 3.5616 2.5785 22.8729 NEURO14B 3.5616 5.6266 3.9679 12.3152 23.4292 NEURO16B 2.5785 3.9679 3.6147 9.7367 12.7632 20.1975 Function value of this group: 3187.7516 Where the fit function is -2 * Log-likelihood of raw data GROUP NUMBER: 4 G4: DZ twins opp sex MATRIX A This is a computed FULL matrix of order 3 by 3 [=X*X'] 1 2 3 1 6.5972 6.1950 3.3051 2 6.1950 10.3771 6.1876 3 3.3051 6.1876 3.7416 MATRIX C This is a computed FULL matrix of order 3 by 3 [=Y*Y'] 1 2 3 1 0.4917 0.4641 0.9260 2 0.4641 0.4381 0.8741 3 0.9260 0.8741 1.7439 MATRIX E This is a computed FULL matrix of order 3 by 3 [=Z*Z'] 1 2 3 1 15.7841 5.6562 5.5056 2 5.6562 12.6141 5.7016 3 5.5056 5.7016 14.7120 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX M This is a FULL matrix of order 1 by 3 1 2 3 1 10.8625 10.4872 11.2794 MATRIX X This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 2.5685 2 2.4119 2.1354 3 1.2868 1.4442 0.0000 MATRIX Y This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.7012 2 0.6619 0.0000 3 1.3206 0.0000 0.0000 MATRIX Z This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 -3.9729 2 -1.4237 -3.2538 3 -1.3858 -1.1460 -3.3880 Vector of OBSERVED means NEURO12 NEURO14 NEURO16 Mean 10.5781 11.0824 11.9145 Vector of EXPECTED means NEURO12 NEURO14 NEURO16 Mean 10.8625 10.4872 11.2794 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX NEURO12 NEURO14 NEURO16 NEURO12 22.8729 NEURO14 12.3152 23.4292 NEURO16 9.7367 12.7632 20.1975 Function value of this group: 2978.2886 Where the fit function is -2 * Log-likelihood of raw data GROUP NUMBER: 5 G5: Calculate standardised solution and genetic and environmental correlations MATRIX A This is a computed FULL matrix of order 3 by 3 [=X*X'] 1 2 3 1 6.5972 6.1950 3.3051 2 6.1950 10.3771 6.1876 3 3.3051 6.1876 3.7416 MATRIX C This is a computed FULL matrix of order 3 by 3 [=Y*Y'] 1 2 3 1 0.4917 0.4641 0.9260 2 0.4641 0.4381 0.8741 3 0.9260 0.8741 1.7439 MATRIX E This is a computed FULL matrix of order 3 by 3 [=Z*Z'] 1 2 3 1 15.7841 5.6562 5.5056 2 5.6562 12.6141 5.7016 3 5.5056 5.7016 14.7120 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX I This is an IDENTITY matrix of order 3 by 3 MATRIX Q This is a computed FULL matrix of order 6 by 6 [=\STND(A+C+E|A+C_A+C|A+C+E)] 1 2 3 4 5 6 1 1.0000 0.5320 0.4530 0.3099 0.2877 0.1969 2 0.5320 1.0000 0.5867 0.2877 0.4616 0.3246 3 0.4530 0.5867 1.0000 0.1969 0.3246 0.2716 4 0.3099 0.2877 0.1969 1.0000 0.5320 0.4530 5 0.2877 0.4616 0.3246 0.5320 1.0000 0.5867 6 0.1969 0.3246 0.2716 0.4530 0.5867 1.0000 MATRIX S This is a computed FULL matrix of order 3 by 9 [=A%(A+C+E)|C%(A+C+E)|E%(A+C+E)] 1 2 3 4 5 6 7 1 0.2884 0.5030 0.3395 0.0215 0.0377 0.0951 0.6901 2 0.5030 0.4429 0.4848 0.0377 0.0187 0.0685 0.4593 3 0.3395 0.4848 0.1853 0.0951 0.0685 0.0863 0.5654 8 9 1 0.4593 0.5654 2 0.5384 0.4467 3 0.4467 0.7284 MATRIX T This is a computed FULL matrix of order 3 by 3 [=\STND(A)] 1 2 3 1 1.0000 0.7487 0.6652 2 0.7487 1.0000 0.9930 3 0.6652 0.9930 1.0000 MATRIX U This is a computed FULL matrix of order 3 by 3 [=\STND(C)] 1 2 3 1 1.0000 1.0000 1.0000 2 1.0000 1.0000 1.0000 3 1.0000 1.0000 1.0000 MATRIX V This is a computed FULL matrix of order 3 by 3 [=\STND(E)] 1 2 3 1 1.0000 0.4009 0.3613 2 0.4009 1.0000 0.4185 3 0.3613 0.4185 1.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 2.5685 2 2.4119 2.1354 3 1.2868 1.4442 0.0000 MATRIX Y This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.7012 2 0.6619 0.0000 3 1.3206 0.0000 0.0000 MATRIX Z This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 -3.9729 2 -1.4237 -3.2538 3 -1.3858 -1.1460 -3.3880 EXPECTED MATRIX of this CALCULATION group *** WARNING! *** Minimization may not be successful. See above CODE GREEN - it probably was OK Your model has 21 estimated parameters and 1824 Observed statistics -2 times log-likelihood of data >>> 10424.883 Degrees of freedom >>>>>>>>>>>>>>>> 1803 Akaike's Information Criterion >>>> 6818.883 Bayesian Information Criterion >>>> -509.730 Sample size Adjusted BIC >>>> 2352.134 Deviance Information Criterion >>>> 1147.117 This problem used 3.3% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.67 Execution 0: 0: 0: 3.91 TOTAL 0: 0: 0: 4.58 Total number of warnings issued: 2 ______________________________________________________________________________ Multiple fit option in effect. The following MX script lines have been read: SAVE NF.MXS END Summary of VL file data for group 2 NEURO12 NEURO14 NEURO16 NEURO12B NEURO14B NEURO16B Code 1.0000 2.0000 3.0000 4.0000 5.0000 6.0000 Number 128.0000 113.0000 134.0000 130.0000 113.0000 134.0000 Mean 10.7344 10.5221 10.5597 11.2385 10.7876 10.8731 Variance 23.6326 24.0371 22.3211 23.7662 23.1938 19.2003 Minimum 1.0000 1.0000 1.0000 1.0000 0.0000 0.0000 Maximum 20.0000 20.0000 19.0000 20.0000 20.0000 19.0000 Summary of VL file data for group 3 NEURO12 NEURO14 NEURO16 NEURO12B NEURO14B NEURO16B Code 1.0000 2.0000 3.0000 4.0000 5.0000 6.0000 Number 111.0000 93.0000 75.0000 111.0000 93.0000 75.0000 Mean 11.0541 10.1290 11.5600 11.1532 9.7849 11.1467 Variance 22.3034 24.6500 18.6731 20.8865 22.9860 17.7518 Minimum 0.0000 0.0000 3.0000 0.0000 0.0000 2.0000 Maximum 20.0000 20.0000 19.0000 20.0000 20.0000 20.0000 Summary of VL file data for group 4 NEURO12 NEURO14 NEURO16 Code 1.0000 2.0000 3.0000 Number 192.0000 170.0000 152.0000 Mean 10.5781 11.0824 11.9145 Variance 23.3064 22.5697 19.4598 Minimum 1.0000 0.0000 1.0000 Maximum 20.0000 20.0000 20.0000 PARAMETER SPECIFICATIONS GROUP NUMBER: 1 G1: Model paramaters MATRIX A This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 1 2 2 3 3 4 5 6 MATRIX Y This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 7 2 8 9 3 10 11 12 MATRIX Z This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 13 2 14 15 3 16 17 18 GROUP NUMBER: 2 G2: MZ twins MATRIX A This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX M This is a FULL matrix of order 2 by 3 1 2 3 1 80 81 82 2 80 81 82 MATRIX X This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 1 2 2 3 3 4 5 6 MATRIX Y This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 7 2 8 9 3 10 11 12 MATRIX Z This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 13 2 14 15 3 16 17 18 GROUP NUMBER: 3 G3: DZ twins MATRIX A This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX M This is a FULL matrix of order 2 by 3 1 2 3 1 80 81 82 2 80 81 82 MATRIX X This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 1 2 2 3 3 4 5 6 MATRIX Y This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 7 2 8 9 3 10 11 12 MATRIX Z This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 13 2 14 15 3 16 17 18 GROUP NUMBER: 4 G4: DZ twins opp sex MATRIX A This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX M This is a FULL matrix of order 1 by 3 1 2 3 1 80 81 82 MATRIX X This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 1 2 2 3 3 4 5 6 MATRIX Y This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 7 2 8 9 3 10 11 12 MATRIX Z This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 13 2 14 15 3 16 17 18 GROUP NUMBER: 5 G5: Calculate standardised solution and genetic and environmental correlations MATRIX A This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX I This is an IDENTITY matrix of order 3 by 3 MATRIX Q This is a computed FULL matrix of order 6 by 6 It has no free parameters specified MATRIX S This is a computed FULL matrix of order 3 by 9 It has no free parameters specified MATRIX T This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX U This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX V This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 1 2 2 3 3 4 5 6 MATRIX Y This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 7 2 8 9 3 10 11 12 MATRIX Z This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 13 2 14 15 3 16 17 18 Mx starting optimization; number of parameters = 21 MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1: Model paramaters MATRIX A This is a computed FULL matrix of order 3 by 3 [=X*X'] 1 2 3 1 6.5971 6.1949 3.3051 2 6.1949 10.3770 6.1875 3 3.3051 6.1875 3.7416 MATRIX C This is a computed FULL matrix of order 3 by 3 [=Y*Y'] 1 2 3 1 0.4917 0.4642 0.9261 2 0.4642 0.4382 0.8742 3 0.9261 0.8742 1.7440 MATRIX E This is a computed FULL matrix of order 3 by 3 [=Z*Z'] 1 2 3 1 15.7841 5.6562 5.5056 2 5.6562 12.6141 5.7016 3 5.5056 5.7016 14.7120 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX X This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 2.5685 2 2.4119 2.1354 3 1.2868 1.4442 0.0000 MATRIX Y This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.7012 2 0.6619 0.0000 3 1.3206 0.0000 0.0000 MATRIX Z This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 -3.9729 2 -1.4237 -3.2538 3 -1.3858 -1.1459 -3.3880 GROUP NUMBER: 2 G2: MZ twins MATRIX A This is a computed FULL matrix of order 3 by 3 [=X*X'] 1 2 3 1 6.5971 6.1949 3.3051 2 6.1949 10.3770 6.1875 3 3.3051 6.1875 3.7416 MATRIX C This is a computed FULL matrix of order 3 by 3 [=Y*Y'] 1 2 3 1 0.4917 0.4642 0.9261 2 0.4642 0.4382 0.8742 3 0.9261 0.8742 1.7440 MATRIX E This is a computed FULL matrix of order 3 by 3 [=Z*Z'] 1 2 3 1 15.7841 5.6562 5.5056 2 5.6562 12.6141 5.7016 3 5.5056 5.7016 14.7120 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX M This is a FULL matrix of order 2 by 3 1 2 3 1 10.8625 10.4872 11.2794 2 10.8625 10.4872 11.2794 MATRIX X This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 2.5685 2 2.4119 2.1354 3 1.2868 1.4442 0.0000 MATRIX Y This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.7012 2 0.6619 0.0000 3 1.3206 0.0000 0.0000 MATRIX Z This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 -3.9729 2 -1.4237 -3.2538 3 -1.3858 -1.1459 -3.3880 Vector of OBSERVED means NEURO12 NEURO14 NEURO16 NEURO12B NEURO14B NEURO16B Mean 10.7344 10.5221 10.5597 11.2385 10.7876 10.8731 Vector of EXPECTED means NEURO12 NEURO14 NEURO16 NEURO12B NEURO14B NEURO16B Mean 10.8625 10.4872 11.2794 10.8625 10.4872 11.2794 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX NEURO12 NEURO14 NEURO16 NEURO12B NEURO14B NEURO16B NEURO12 22.8729 NEURO14 12.3152 23.4292 NEURO16 9.7367 12.7632 20.1975 NEURO12B 7.0888 6.6591 4.2311 22.8729 NEURO14B 6.6591 10.8152 7.0616 12.3152 23.4292 NEURO16B 4.2311 7.0616 5.4855 9.7367 12.7632 20.1975 Function value of this group: 4258.8432 Where the fit function is -2 * Log-likelihood of raw data GROUP NUMBER: 3 G3: DZ twins MATRIX A This is a computed FULL matrix of order 3 by 3 [=X*X'] 1 2 3 1 6.5971 6.1949 3.3051 2 6.1949 10.3770 6.1875 3 3.3051 6.1875 3.7416 MATRIX C This is a computed FULL matrix of order 3 by 3 [=Y*Y'] 1 2 3 1 0.4917 0.4642 0.9261 2 0.4642 0.4382 0.8742 3 0.9261 0.8742 1.7440 MATRIX E This is a computed FULL matrix of order 3 by 3 [=Z*Z'] 1 2 3 1 15.7841 5.6562 5.5056 2 5.6562 12.6141 5.7016 3 5.5056 5.7016 14.7120 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX M This is a FULL matrix of order 2 by 3 1 2 3 1 10.8625 10.4872 11.2794 2 10.8625 10.4872 11.2794 MATRIX X This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 2.5685 2 2.4119 2.1354 3 1.2868 1.4442 0.0000 MATRIX Y This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.7012 2 0.6619 0.0000 3 1.3206 0.0000 0.0000 MATRIX Z This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 -3.9729 2 -1.4237 -3.2538 3 -1.3858 -1.1459 -3.3880 Vector of OBSERVED means NEURO12 NEURO14 NEURO16 NEURO12B NEURO14B NEURO16B Mean 11.0541 10.1290 11.5600 11.1532 9.7849 11.1467 Vector of EXPECTED means NEURO12 NEURO14 NEURO16 NEURO12B NEURO14B NEURO16B Mean 10.8625 10.4872 11.2794 10.8625 10.4872 11.2794 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX NEURO12 NEURO14 NEURO16 NEURO12B NEURO14B NEURO16B NEURO12 22.8729 NEURO14 12.3152 23.4292 NEURO16 9.7367 12.7632 20.1975 NEURO12B 3.7903 3.5616 2.5786 22.8729 NEURO14B 3.5616 5.6267 3.9679 12.3152 23.4292 NEURO16B 2.5786 3.9679 3.6147 9.7367 12.7632 20.1975 Function value of this group: 3187.7516 Where the fit function is -2 * Log-likelihood of raw data GROUP NUMBER: 4 G4: DZ twins opp sex MATRIX A This is a computed FULL matrix of order 3 by 3 [=X*X'] 1 2 3 1 6.5971 6.1949 3.3051 2 6.1949 10.3770 6.1875 3 3.3051 6.1875 3.7416 MATRIX C This is a computed FULL matrix of order 3 by 3 [=Y*Y'] 1 2 3 1 0.4917 0.4642 0.9261 2 0.4642 0.4382 0.8742 3 0.9261 0.8742 1.7440 MATRIX E This is a computed FULL matrix of order 3 by 3 [=Z*Z'] 1 2 3 1 15.7841 5.6562 5.5056 2 5.6562 12.6141 5.7016 3 5.5056 5.7016 14.7120 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX M This is a FULL matrix of order 1 by 3 1 2 3 1 10.8625 10.4872 11.2794 MATRIX X This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 2.5685 2 2.4119 2.1354 3 1.2868 1.4442 0.0000 MATRIX Y This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.7012 2 0.6619 0.0000 3 1.3206 0.0000 0.0000 MATRIX Z This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 -3.9729 2 -1.4237 -3.2538 3 -1.3858 -1.1459 -3.3880 Vector of OBSERVED means NEURO12 NEURO14 NEURO16 Mean 10.5781 11.0824 11.9145 Vector of EXPECTED means NEURO12 NEURO14 NEURO16 Mean 10.8625 10.4872 11.2794 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX NEURO12 NEURO14 NEURO16 NEURO12 22.8729 NEURO14 12.3152 23.4292 NEURO16 9.7367 12.7632 20.1975 Function value of this group: 2978.2886 Where the fit function is -2 * Log-likelihood of raw data GROUP NUMBER: 5 G5: Calculate standardised solution and genetic and environmental correlations MATRIX A This is a computed FULL matrix of order 3 by 3 [=X*X'] 1 2 3 1 6.5971 6.1949 3.3051 2 6.1949 10.3770 6.1875 3 3.3051 6.1875 3.7416 MATRIX C This is a computed FULL matrix of order 3 by 3 [=Y*Y'] 1 2 3 1 0.4917 0.4642 0.9261 2 0.4642 0.4382 0.8742 3 0.9261 0.8742 1.7440 MATRIX E This is a computed FULL matrix of order 3 by 3 [=Z*Z'] 1 2 3 1 15.7841 5.6562 5.5056 2 5.6562 12.6141 5.7016 3 5.5056 5.7016 14.7120 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX I This is an IDENTITY matrix of order 3 by 3 MATRIX Q This is a computed FULL matrix of order 6 by 6 [=\STND(A+C+E|A+C_A+C|A+C+E)] 1 2 3 4 5 6 1 1.0000 0.5320 0.4530 0.3099 0.2877 0.1969 2 0.5320 1.0000 0.5867 0.2877 0.4616 0.3246 3 0.4530 0.5867 1.0000 0.1969 0.3246 0.2716 4 0.3099 0.2877 0.1969 1.0000 0.5320 0.4530 5 0.2877 0.4616 0.3246 0.5320 1.0000 0.5867 6 0.1969 0.3246 0.2716 0.4530 0.5867 1.0000 MATRIX S This is a computed FULL matrix of order 3 by 9 [=A%(A+C+E)|C%(A+C+E)|E%(A+C+E)] 1 2 3 4 5 6 7 1 0.2884 0.5030 0.3394 0.0215 0.0377 0.0951 0.6901 2 0.5030 0.4429 0.4848 0.0377 0.0187 0.0685 0.4593 3 0.3394 0.4848 0.1852 0.0951 0.0685 0.0863 0.5654 8 9 1 0.4593 0.5654 2 0.5384 0.4467 3 0.4467 0.7284 MATRIX T This is a computed FULL matrix of order 3 by 3 [=\STND(A)] 1 2 3 1 1.0000 0.7487 0.6652 2 0.7487 1.0000 0.9930 3 0.6652 0.9930 1.0000 MATRIX U This is a computed FULL matrix of order 3 by 3 [=\STND(C)] 1 2 3 1 1.0000 1.0000 1.0000 2 1.0000 1.0000 1.0000 3 1.0000 1.0000 1.0000 MATRIX V This is a computed FULL matrix of order 3 by 3 [=\STND(E)] 1 2 3 1 1.0000 0.4009 0.3613 2 0.4009 1.0000 0.4185 3 0.3613 0.4185 1.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 2.5685 2 2.4119 2.1354 3 1.2868 1.4442 0.0000 MATRIX Y This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.7012 2 0.6619 0.0000 3 1.3206 0.0000 0.0000 MATRIX Z This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 -3.9729 2 -1.4237 -3.2538 3 -1.3858 -1.1459 -3.3880 EXPECTED MATRIX of this CALCULATION group Your model has 21 estimated parameters and 1824 Observed statistics -2 times log-likelihood of data >>> 10424.883 Degrees of freedom >>>>>>>>>>>>>>>> 1803 Akaike's Information Criterion >>>> 6818.883 Bayesian Information Criterion >>>> -509.730 Sample size Adjusted BIC >>>> 2352.134 Deviance Information Criterion >>>> 1147.117 Saturated model fit* >>>>>>>>>>> 10424.883 Saturated model df* >>>>>>>>>>> 1803 Difference Chi-squared >>>>>>>> 0.000 Difference d.f. >>>>>>>>>>>>>>> 0 Probability incalculable Akaike's Information Criterion > 0.000 * Saturated model statistic computed earlier in this job This problem used 3.3% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.10 Execution 0: 0: 0: 0.95 TOTAL 0: 0: 0: 1.05 Total number of warnings issued: 2 ______________________________________________________________________________ !AE !GET NF.MXS !DROP 7 TO 12 !CE !GET NF.MXS !DROP 1 TO 6 !E !GET NF.MXS !DROP 1 TO 12