!raw data: choleski model iq age5 & age7 & age10 & age12 no sex diffs The following MX script lines were read for group 1 G1: CALCULATION GROUP DATA CALC NGROUPS=7 MATRICES X LOW 4 4 FREE !A Y LOW 4 4 FREE !C Z LOW 4 4 FREE !E M FULL 1 4 !MEANS H FULL 1 1 BEGIN ALGEBRA; A= X*X' ; C= Y*Y' ; E= Z*Z' ; END ALGEBRA; MATRIX H .5 SPECIFY M 104 105 106 107 START ALL 3 START 7 X 1 1 X 2 2 X 3 3 X 4 4 !Start 5 Y 1 1 Y 2 1 Y 3 1 Y 4 1 START 5 Y 1 1 Y 2 2 Y 3 3 Y 4 4 START 9 Z 1 1 Z 2 2 Z 3 3 Z 4 4 START 100 M 1 1 M 1 2 M 1 3 M 1 4 END The following MX script lines were read for group 2 G2: MALE MZ TWIN PAIRS DATA NINPUT_VARS=8 NOBSERVATIONS=0 MISSING=-1 NOTE: Missing value *string* set to '-1' REC FILE=IQ4MZM.REC NOTE: Rectangular file contained 42 records with data LABELS T1IQ5 T1IQ7 T1IQ10 T1IQ12 T2IQ5 T2IQ7 T2IQ10 T2IQ12 MATRICES = GROUP 1 COVARIANCES A+C+E | A+C _ A+C | A+C+E / MEANS M|M / OPTIONS NO_OUTPUT END The following MX script lines were read for group 3 G3: FEMALE MZ TWIN PAIRS DATA NINPUT_VARS=8 NOBSERVATIONS=0 MISSING=-1 NOTE: Missing value *string* set to '-1' REC FILE=IQ4MZF.REC NOTE: Rectangular file contained 47 records with data LABELS T1IQ5 T1IQ7 T1IQ10 T1IQ12 T2IQ5 T2IQ7 T2IQ10 T2IQ12 MATRICES = GROUP 1 COVARIANCES A+C+E | A+C _ A+C | A+C+E / MEANS M|M / OPTIONS NO_OUTPUT END The following MX script lines were read for group 4 G4: MALE DZ TWIN PAIRS DATA NINPUT_VARS=8 NOBSERVATIONS=0 MISSING=-1 NOTE: Missing value *string* set to '-1' REC FILE=IQ4DZM.REC NOTE: Rectangular file contained 44 records with data LABELS T1IQ5 T1IQ7 T1IQ10 T1IQ12 T2IQ5 T2IQ7 T2IQ10 T2IQ12 MATRICES = GROUP 1 COVARIANCES A+C+E | H@A+C_ H@A+C | A+C+E/ MEANS M|M / OPTIONS NO_OUTPUT END The following MX script lines were read for group 5 G5: FEMALE DZ TWIN PAIRS DATA NINPUT_VARS=8 NOBSERVATIONS=0 MISSING=-1 NOTE: Missing value *string* set to '-1' REC FILE=IQ4DZF.REC NOTE: Rectangular file contained 37 records with data LABELS T1IQ5 T1IQ7 T1IQ10 T1IQ12 T2IQ5 T2IQ7 T2IQ10 T2IQ12 MATRICES = GROUP 1 COVARIANCES A+C+E | H@A+C_ H@A+C | A+C+E/ MEANS M|M / OPTIONS NO_OUTPUT END The following MX script lines were read for group 6 G6: MALE-FEMALE DZ TWIN PAIRS DATA NINPUT_VARS=8 NOBSERVATIONS=0 MISSING=-1 NOTE: Missing value *string* set to '-1' REC FILE=IQ4DOS.REC NOTE: Rectangular file contained 39 records with data LABELS T1IQ5 T1IQ7 T1IQ10 T1IQ12 T2IQ5 T2IQ7 T2IQ10 T2IQ12 MATRICES = GROUP 1 COVARIANCES A+C+E | H@A+C_ H@A+C | A+C+E/ MEANS M|M / OPTIONS RSIDUALS OPTIONS NO_OUTPUT OPTIONS ND=3 IT=255 END The following MX script lines were read for group 7 G7: CALCULATE STANDARDIZED ESTIMATES DATA CALCULATION MATRICES = GROUP 1 BEGIN ALGEBRA; F=\STND(X*X'); G=\STND(Y*Y'); I=\STND(Z*Z'); END ALGEBRA; COMPUTE A%(A+C+E) | C%(A+C+E) | E%(A+C+E) / OPTIONS RS OPTION SAT=11078.7,1304 END Summary of VL file data for group 2 T1IQ5 T1IQ7 T1IQ10 T1IQ12 T2IQ5 T2IQ7 T2IQ10 Code 1.000 2.000 3.000 4.000 5.000 6.000 7.000 Number 42.000 37.000 38.000 32.000 42.000 37.000 38.000 Mean 103.000 104.432 105.737 98.656 100.905 102.892 107.184 Variance 189.619 193.975 228.089 133.038 177.181 158.961 206.519 Minimum 71.000 75.000 72.000 80.000 71.000 80.000 79.000 Maximum 137.000 139.000 145.000 128.000 129.000 130.000 142.000 T2IQ12 Code 8.000 Number 30.000 Mean 98.567 Variance 212.712 Minimum 67.000 Maximum 137.000 Summary of VL file data for group 3 T1IQ5 T1IQ7 T1IQ10 T1IQ12 T2IQ5 T2IQ7 T2IQ10 Code 1.000 2.000 3.000 4.000 5.000 6.000 7.000 Number 47.000 42.000 43.000 36.000 47.000 42.000 43.000 Mean 103.511 101.714 106.628 99.222 105.234 103.071 107.884 Variance 175.186 229.918 317.722 203.340 196.775 261.209 329.126 Minimum 75.000 72.000 73.000 70.000 75.000 74.000 63.000 Maximum 127.000 135.000 145.000 127.000 137.000 145.000 140.000 T2IQ12 Code 8.000 Number 36.000 Mean 100.333 Variance 187.500 Minimum 67.000 Maximum 126.000 Summary of VL file data for group 4 T1IQ5 T1IQ7 T1IQ10 T1IQ12 T2IQ5 T2IQ7 T2IQ10 Code 1.000 2.000 3.000 4.000 5.000 6.000 7.000 Number 44.000 41.000 41.000 33.000 43.000 41.000 41.000 Mean 98.818 100.463 106.244 100.303 102.186 101.268 107.098 Variance 152.194 242.883 186.623 171.969 137.175 212.928 131.454 Minimum 69.000 67.000 77.000 72.000 73.000 70.000 74.000 Maximum 126.000 138.000 143.000 138.000 139.000 125.000 131.000 T2IQ12 Code 8.000 Number 33.000 Mean 101.121 Variance 166.652 Minimum 80.000 Maximum 135.000 Summary of VL file data for group 5 T1IQ5 T1IQ7 T1IQ10 T1IQ12 T2IQ5 T2IQ7 T2IQ10 Code 1.000 2.000 3.000 4.000 5.000 6.000 7.000 Number 37.000 34.000 37.000 31.000 37.000 34.000 37.000 Mean 101.270 101.765 103.973 99.903 102.054 101.824 104.027 Variance 198.792 185.003 221.918 156.345 159.619 137.616 247.378 Minimum 74.000 62.000 69.000 61.000 80.000 79.000 70.000 Maximum 142.000 123.000 145.000 126.000 139.000 130.000 142.000 T2IQ12 Code 8.000 Number 31.000 Mean 99.935 Variance 149.996 Minimum 74.000 Maximum 124.000 Summary of VL file data for group 6 T1IQ5 T1IQ7 T1IQ10 T1IQ12 T2IQ5 T2IQ7 T2IQ10 Code 1.000 2.000 3.000 4.000 5.000 6.000 7.000 Number 39.000 38.000 37.000 34.000 39.000 38.000 37.000 Mean 107.205 106.237 112.811 102.294 102.359 104.526 107.973 Variance 173.702 229.286 255.613 181.502 132.897 251.881 220.134 Minimum 64.000 65.000 69.000 64.000 70.000 64.000 72.000 Maximum 128.000 137.000 140.000 125.000 121.000 133.000 137.000 T2IQ12 Code 8.000 Number 35.000 Mean 100.857 Variance 137.894 Minimum 66.000 Maximum 128.000 PARAMETER SPECIFICATIONS GROUP NUMBER: 1 G1: calculation group MATRIX A This is a computed FULL matrix of order 4 by 4 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 4 by 4 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 4 by 4 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 4 1 2 3 4 1 104 105 106 107 MATRIX X This is a LOWER TRIANGULAR matrix of order 4 by 4 1 2 3 4 1 1 2 2 3 3 4 5 6 4 7 8 9 10 MATRIX Y This is a LOWER TRIANGULAR matrix of order 4 by 4 1 2 3 4 1 11 2 12 13 3 14 15 16 4 17 18 19 20 MATRIX Z This is a LOWER TRIANGULAR matrix of order 4 by 4 1 2 3 4 1 21 2 22 23 3 24 25 26 4 27 28 29 30 GROUP NUMBER: 7 G7: calculate standardized estimates MATRIX A This is a computed FULL matrix of order 4 by 4 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 4 by 4 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 4 by 4 It has no free parameters specified MATRIX F This is a computed FULL matrix of order 4 by 4 It has no free parameters specified MATRIX G This is a computed FULL matrix of order 4 by 4 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 a computed FULL matrix of order 4 by 4 It has no free parameters specified MATRIX M This is a FULL matrix of order 1 by 4 1 2 3 4 1 104 105 106 107 MATRIX X This is a LOWER TRIANGULAR matrix of order 4 by 4 1 2 3 4 1 1 2 2 3 3 4 5 6 4 7 8 9 10 MATRIX Y This is a LOWER TRIANGULAR matrix of order 4 by 4 1 2 3 4 1 11 2 12 13 3 14 15 16 4 17 18 19 20 MATRIX Z This is a LOWER TRIANGULAR matrix of order 4 by 4 1 2 3 4 1 21 2 22 23 3 24 25 26 4 27 28 29 30 *** 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: calculation group MATRIX A This is a computed FULL matrix of order 4 by 4 [=X*X'] 1 2 3 4 1 43.926 61.688 67.943 51.834 2 61.688 89.179 97.648 79.814 3 67.943 97.648 107.047 86.326 4 51.834 79.814 86.326 80.510 MATRIX C This is a computed FULL matrix of order 4 by 4 [=Y*Y'] 1 2 3 4 1 87.131 59.123 59.720 58.098 2 59.123 68.386 52.302 51.035 3 59.720 52.302 76.815 67.816 4 58.098 51.035 67.816 60.821 MATRIX E This is a computed FULL matrix of order 4 by 4 [=Z*Z'] 1 2 3 4 1 42.364 6.703 4.013 -6.570 2 6.703 64.493 16.208 -1.466 3 4.013 16.208 54.632 -4.176 4 -6.570 -1.466 -4.176 34.949 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.500 MATRIX M This is a FULL matrix of order 1 by 4 1 2 3 4 1 102.656 102.338 106.607 99.333 MATRIX X This is a LOWER TRIANGULAR matrix of order 4 by 4 1 2 3 4 1 6.628 2 9.308 1.596 3 10.251 1.399 0.000 4 7.821 4.398 0.000 0.000 MATRIX Y This is a LOWER TRIANGULAR matrix of order 4 by 4 1 2 3 4 1 9.334 2 6.334 5.317 3 6.398 2.215 5.566 4 6.224 2.184 4.161 0.000 MATRIX Z This is a LOWER TRIANGULAR matrix of order 4 by 4 1 2 3 4 1 6.509 2 1.030 7.964 3 0.616 1.955 7.101 4 -1.009 -0.054 -0.486 5.804 GROUP NUMBER: 7 G7: calculate standardized estimates MATRIX A This is a computed FULL matrix of order 4 by 4 [=X*X'] 1 2 3 4 1 43.926 61.688 67.943 51.834 2 61.688 89.179 97.648 79.814 3 67.943 97.648 107.047 86.326 4 51.834 79.814 86.326 80.510 MATRIX C This is a computed FULL matrix of order 4 by 4 [=Y*Y'] 1 2 3 4 1 87.131 59.123 59.720 58.098 2 59.123 68.386 52.302 51.035 3 59.720 52.302 76.815 67.816 4 58.098 51.035 67.816 60.821 MATRIX E This is a computed FULL matrix of order 4 by 4 [=Z*Z'] 1 2 3 4 1 42.364 6.703 4.013 -6.570 2 6.703 64.493 16.208 -1.466 3 4.013 16.208 54.632 -4.176 4 -6.570 -1.466 -4.176 34.949 MATRIX F This is a computed FULL matrix of order 4 by 4 [=\STND(X*X')] 1 2 3 4 1 1.000 0.986 0.991 0.872 2 0.986 1.000 0.999 0.942 3 0.991 0.999 1.000 0.930 4 0.872 0.942 0.930 1.000 MATRIX G This is a computed FULL matrix of order 4 by 4 [=\STND(Y*Y')] 1 2 3 4 1 1.000 0.766 0.730 0.798 2 0.766 1.000 0.722 0.791 3 0.730 0.722 1.000 0.992 4 0.798 0.791 0.992 1.000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.500 MATRIX I This is a computed FULL matrix of order 4 by 4 [=\STND(Z*Z')] 1 2 3 4 1 1.000 0.128 0.083 -0.171 2 0.128 1.000 0.273 -0.031 3 0.083 0.273 1.000 -0.096 4 -0.171 -0.031 -0.096 1.000 MATRIX M This is a FULL matrix of order 1 by 4 1 2 3 4 1 102.656 102.338 106.607 99.333 MATRIX X This is a LOWER TRIANGULAR matrix of order 4 by 4 1 2 3 4 1 6.628 2 9.308 1.596 3 10.251 1.399 0.000 4 7.821 4.398 0.000 0.000 MATRIX Y This is a LOWER TRIANGULAR matrix of order 4 by 4 1 2 3 4 1 9.334 2 6.334 5.317 3 6.398 2.215 5.566 4 6.224 2.184 4.161 0.000 MATRIX Z This is a LOWER TRIANGULAR matrix of order 4 by 4 1 2 3 4 1 6.509 2 1.030 7.964 3 0.616 1.955 7.101 4 -1.009 -0.054 -0.486 5.804 EXPECTED MATRIX of this CALCULATION group 1 2 3 4 5 6 7 8 1 0.253 0.484 0.516 0.501 0.502 0.464 0.454 0.562 2 0.484 0.402 0.588 0.617 0.464 0.308 0.315 0.394 3 0.516 0.588 0.449 0.576 0.454 0.315 0.322 0.452 4 0.501 0.617 0.576 0.457 0.562 0.394 0.452 0.345 9 10 11 12 1 0.244 0.053 0.030 -0.064 2 0.053 0.290 0.098 -0.011 3 0.030 0.098 0.229 -0.028 4 -0.064 -0.011 -0.028 0.198 *** WARNING! *** Minimization may not be successful. See above CODE GREEN - it probably was OK Your model has 34 estimated parameters and 1524 Observed statistics -2 times log-likelihood of data >>> 11304.484 Degrees of freedom >>>>>>>>>>>>>>>> 1490 Saturated model fit* >>>>>>>>>>> 11078.700 Saturated model df* >>>>>>>>>>> 1304 Difference Chi-squared >>>>>>>> 225.784 Difference d.f. >>>>>>>>>>>>>>> 186 Probability >>>>>>>>>>>>>>>>>>>> 0.025 Akaike's Information Criterion > -146.216 * Saturated model statistic supplied by user This problem used 5.3% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.77 Execution 0: 0: 0:40.92 TOTAL 0: 0: 0:41.69 Total number of warnings issued: 2 ______________________________________________________________________________