**MX-PC 1.65b** Job started on 03/05/08 at 10:29:34 ! MULTIVARIATE ANALYSIS OF DUTCH AGGRESSION DATA ! NYAGG1 = AGGRESSION IN '91 ! NYAGG3 = AGGRESSION IN '95 ! NYAGG4 = AGGRESSION IN '97 ! NYAGG5 = AGGRESSION IN 2000 ! ACE MODEL #DEFINE NVAR 2 ! NUMBER OF DEPENDENT VARIABLES PER INDIVIDUAL #DEFINE NDEF 1 ! NUMBER OF DEFINITION VARIABLES (AGE -VARIES PER YEAR-) The following MX script lines were read for group 1 #NGROUPS 4 Note: #NGroup set number of groups to 4 G1: CALCULATION GROUP DATA CALC BEGIN MATRICES; ! **GO TO SLIDE** "ADDITIVE GENETIC PATH COEFFICIENTS (X)" X LOWER NVAR NVAR FREE ! ADDITIVE GENETIC PATH COEFFICIENT (STANDARD DEVIATION) ! **GO TO SLIDE** "ADDITIVE GENETIC PATH COEFFICIENTS (Y)" Y LOWER NVAR NVAR FREE ! COMMON ENVIRONMENTAL PATH COEFFICIENT (STANDARD DEVIATION) Z LOWER NVAR NVAR FREE ! UNIQUE ENVIRONMENTAL PATH COEFFICIENT (STANDARD DEVIATION) H FULL 1 1 FIX ! .5 FOR DZ COV A **GO TO SLIDE** "ADDITIVE GENETIC CROSS-TWIN COVARIANCE (DZ)" M FULL 1 NVAR FREE ! GRAND MEAN PHENOTYPES B FULL 1 NDEF FREE ! AGE BETA END MATRICES; MATRIX H .5 ! MATRIX H IS FIXED AND CONTAINS .5 ! PROVIDE STARTING VALUES: ! FOR THE MEANS. PHENOTYPIC MEANS = 7.6 6.5 START 10 M 1 1 M 1 2 START -0.1 B 1 1 ! FOR THE VARIANCE COMPONENTS: DIAGONALS ! PHENOTYPIC VARIANCE ST 1.0 X 1 1 X 2 2 ST 2.5 Y 1 1 ST -1.5 Y 2 2 ST 3.5 Z 1 1 ST 3.0 Z 2 2 ! PHENOTYPIC COVARIANCE: OFF DIAGONALS ST 0.5 X 2 1 ST 1.7 Y 2 1 ST 1.6 Z 2 1 BEGIN ALGEBRA; ! **GO TO SLIDE** "TOTAL WITHIN-TWIN COVARIANCE" A= X*X' ; ! ADDITIVE GENETIC VARIANCE ***WRITE ALGEBRA*** C= Y*Y' ; ! COMMON ENVIRONMENTAL VARIANCE ***WRITE ALGEBRA*** E= Z*Z' ; ! UNIQUE ENVIRONMENTAL VARIANCE ***WRITE ALGEBRA*** END ALGEBRA; OPTIONS ND=2 END The following MX script lines were read for group 2 !__________________________________________________ ! G2: MZ TWINS, DATAGROUP DATA NINPUT_VARS= 19 ! NR OF INPUTVARS PER FAMILY MISSING=-1.00 NOTE: Missing value *string* set to '-1.00' RECTANGULAR FILE=AGGRESSION.DAT Rectangular continuous data read initiated Note: Maximum ordinal/rectangular record length is: 1000 NOTE: Rectangular file contained 3885 records with data that contained a total of 56380 observations LABELS ID TWZYG ZYG2 NYAGG1.1 NYAGG3.1 NYAGG4.1 NYAGG5.1 NYAGG1.2 NYAGG3.2 NYAGG4.2 NYAGG5.2 AGE1.1 AGE3.1 AGE4.1 AGE5.1 AGE1.2 AGE3.2 AGE4.2 AGE5.2 ! TELL MX WHAT IS FOUND IN THE DATAFILE SELECT IF ZYG2 = 1 ; ! SELECT MZ'S: ZYG2: 1=MZ 2=DZ NOTE: Select if yields 1077 data vectors for analysis NOTE: Vectors contain a total of 15599 observations SELECT NYAGG1.1 NYAGG3.1 ! NYAGG1.1 = AGGRESSION TWIN1 '91, NYAGG3.1 = AGGRESSION1 TWIN '95 NYAGG1.2 NYAGG3.2 ! NYAGG1.2 = AGGRESSION TWIN2 '91, NYAGG3.2 = AGGRESSION2 TWIN '95 AGE1.1 AGE3.1 ! AGE1.1 = AGE AT INTERVIEW TWIN1 '91, AGE3.1 = AGE AT INTERVIEW TWIN1 '95 AGE1.2 AGE3.2; ! AGE1.2 = AGE AT INTERVIEW TWIN2 '91, AGE3.2 = AGE AT INTERVIEW TWIN2 '95 DEFINITION NOTE: Selection yields 1077 data vectors for analysis NOTE: Vectors contain a total of 6309 observations AGE1.1 AGE3.1 AGE1.2 AGE3.2; MATRICES = GROUP 1 NOTE: Definition yields 743 data vectors for analysis NOTE: Vectors contain a total of 2001 observations O FULL NDEF NVAR FIX ! DEFINITION VARIABLE(S) FOR TWIN 1 I.E. AGE FOR TWIN 1 AT TIMES 1 & 2 P FULL NDEF NVAR FIX ! DEFINITION VARIABLE(S) FOR TWIN 2 I.E. AGE FOR TWIN 2 AT TIMES 1 & 2 END MATRICES; SPECIFY O AGE1.1 AGE3.1 ! AGE FOR TWIN 1 AT TIMES 1 & 2 SPECIFY P AGE1.2 AGE3.2 ! AGE FOR TWIN 2 AT TIMES 1 & 2 ! **GO TO SLIDE** "AGE EFFECTS ON MEAN " MEANS M + B*O| M + B*P; ! MODEL FOR MEANS = 1 BY 4 MATRIX ! TWIN1-VARIABLE-1 TWIN1-VARIABLE-2 | TWIN2-VARIABLE-1 TWIN2-VARIABLE-2 ! + + + + ! AGE EFFECT AGE EFFECT AGE EFFECT AGE EFFECT COVARIANCES ! MODEL FOR MZ VARIANCE/COVARIANCES A+C+E | A+C _ ! **GO TO SLIDE** "PREDICTED MODEL" A+C | A+C+E ; OPTIONS RSIDUALS END The following MX script lines were read for group 3 !__________________________________________________ ! G3: DZ TWINS, DATAGROUP DATA NINPUT_VARS= 19 ! NR OF INPUTVARS PER FAMILY MISSING=-1.00 NOTE: Missing value *string* set to '-1.00' RECTANGULAR FILE=AGGRESSION.DAT Rectangular continuous data read initiated Note: Maximum ordinal/rectangular record length is: 1000 NOTE: Rectangular file contained 3885 records with data that contained a total of 56380 observations LABELS ID TWZYG ZYG2 NYAGG1.1 NYAGG3.1 NYAGG4.1 NYAGG5.1 NYAGG1.2 NYAGG3.2 NYAGG4.2 NYAGG5.2 AGE1.1 AGE3.1 AGE4.1 AGE5.1 AGE1.2 AGE3.2 AGE4.2 AGE5.2 ! TELL MX WHAT IS FOUND IN THE DATAFILE ! WHAT'S WRONG HERE? SELECT IF ZYG2 = 2 ; ! SELECT MZ'S: ZYG2: 1=MZ 2=DZ NOTE: Select if yields 2808 data vectors for analysis NOTE: Vectors contain a total of 40781 observations SELECT NYAGG1.1 NYAGG3.1 NYAGG1.2 NYAGG3.2 AGE1.1 AGE3.1 AGE1.2 AGE3.2; DEFINITION NOTE: Selection yields 2808 data vectors for analysis NOTE: Vectors contain a total of 15792 observations AGE1.1 AGE3.1 AGE1.2 AGE3.2; Note: Global variable previously defined. Updating AGE1.1 Note: Global variable previously defined. Updating AGE3.1 Note: Global variable previously defined. Updating AGE1.2 Note: Global variable previously defined. Updating AGE3.2 MATRICES = GROUP 1 NOTE: Definition yields 1655 data vectors for analysis NOTE: Vectors contain a total of 4560 observations O FULL NDEF NVAR FIX ! DEFINITION VARIABLE(S) FOR TWIN 1 I.E. AGE FOR TWIN 1 AT TIMES 1 & 2 P FULL NDEF NVAR FIX ! DEFINITION VARIABLE(S) FOR TWIN 2 I.E. AGE FOR TWIN 2 AT TIMES 1 & 2 END MATRICES; SPECIFY O AGE1.1 AGE3.1 SPECIFY P AGE1.2 AGE3.2 MEANS M + B*O | M + B*P; ! MODEL FOR MEANS COVARIANCES ! ***WRITE THE DZ COVARIANCE STATEMENT*** A+C+E | H@A+C _ ! **GO TO SLIDE** "PREDICTED MODEL" H@A+C | A+C+E ; OPTIONS RSIDUALS END The following MX script lines were read for group 4 !__________________________________________________ ! G4: CALCULATE A, C, AND E CORRELATIONS DATA CALC BEGIN MATRICES A COMPUTED NVAR NVAR = A1 C COMPUTED NVAR NVAR = C1 E COMPUTED NVAR NVAR = E1 H FULL 1 1 I IDENT NVAR NVAR END MATRICES MA H 0.5 BEGIN ALGEBRA ; P = \STND(A+C+E | H@A+C _H@A+C | A+C+E); ! ESTIMATES PHENOTYPIC CORRELATIONS ! **GO TO SLIDE** T = \STND(A); ! ESTIMATE GENETIC CORRELATIONS U = \STND(C); ! ESTIMATE SHARED ENVIRONMENTAL CORRELATIONS V = \STND(E); ! ESTIMATE NON-SHARED ENVIRONMENTAL CORRELATIONS S = A%(A+C+E)| C%(A+C+E)|E%(A+C+E); ! !T =\SQRT(I.A)~*A*\SQRT(I.A)~ ; ! CALCULATES RG (OFF DIAGONAL) OR:T=\STND(A); !U = \SQRT(I.C)~*C*\SQRT(I.C)~ ; ! CALCULATES RC (OFF DIAGONAL) !V = \SQRT(I.E)~*E*\SQRT(I.E)~ ; ! CALCULATES RE (OFF DIAGONAL) END ALGEBRA ; OP ND=4 ! DECIMAL PLACES OP MULTIPLE OP ISSAT ! COMPARES ACE TO SUBMODELS BELOW END Summary of VL file data for group 2 AGE3.2 AGE1.2 AGE3.1 AGE1.1 NYAGG1.1 NYAGG3.1 NYAGG1.2 Code -4.0000 -3.0000 -2.0000 -1.0000 1.0000 2.0000 3.0000 Number 743.0000 743.0000 743.0000 743.0000 504.0000 501.0000 503.0000 Mean 336.0279 328.0273 329.4289 329.3497 7.6121 6.2330 7.4514 Variance **************************************** 25.6824 18.3844 20.3358 Minimum 14.2700 13.2500 14.3100 13.4000 0.0000 0.0000 0.0000 Maximum 999.0000 999.0000 999.0000 999.0000 27.4400 27.4400 25.3300 NYAGG3.2 Code 4.0000 Number 493.0000 Mean 6.4138 Variance 21.4903 Minimum 0.0000 Maximum 33.7800 Summary of VL file data for group 3 AGE3.2 AGE1.2 AGE3.1 AGE1.1 NYAGG1.1 NYAGG3.1 NYAGG1.2 Code -4.0000 -3.0000 -2.0000 -1.0000 1.0000 2.0000 3.0000 Number 1655.0000 1655.0000 1655.0000 1655.0000 1125.0000 1164.0000 1126.0000 Mean 300.2697 327.8228 296.7262 327.2353 7.6462 6.6732 7.7991 Variance **************************************** 19.4824 16.7076 19.7052 Minimum 14.2900 12.6000 14.2900 12.6000 0.0000 0.0000 0.0000 Maximum 999.0000 999.0000 999.0000 999.0000 32.7200 27.4400 29.5600 NYAGG3.2 Code 4.0000 Number 1145.0000 Mean 6.4779 Variance 17.5274 Minimum 0.0000 Maximum 26.3900 PARAMETER SPECIFICATIONS GROUP NUMBER: 1 G1: calculation group MATRIX A This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX B This is a FULL matrix of order 1 by 1 1 1 12 MATRIX C This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 2 by 2 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 2 1 2 1 10 11 MATRIX X This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 1 2 2 3 MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 4 2 5 6 MATRIX Z This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 7 2 8 9 GROUP NUMBER: 2 G2: MZ twins, datagroup MATRIX A This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX B This is a FULL matrix of order 1 by 1 1 1 12 MATRIX C This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 2 by 2 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 2 1 2 1 10 11 MATRIX O This is a FULL matrix of order 1 by 2 1 2 1 -1 -2 MATRIX P This is a FULL matrix of order 1 by 2 1 2 1 -3 -4 MATRIX X This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 1 2 2 3 MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 4 2 5 6 MATRIX Z This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 7 2 8 9 GROUP NUMBER: 3 G3: DZ twins, datagroup MATRIX A This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX B This is a FULL matrix of order 1 by 1 1 1 12 MATRIX C This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 2 by 2 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 2 1 2 1 10 11 MATRIX O This is a FULL matrix of order 1 by 2 1 2 1 -1 -2 MATRIX P This is a FULL matrix of order 1 by 2 1 2 1 -3 -4 MATRIX X This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 1 2 2 3 MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 4 2 5 6 MATRIX Z This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 7 2 8 9 GROUP NUMBER: 4 G4: Calculate A, C, and E correlations MATRIX A This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 2 by 2 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 2 by 2 MATRIX P This is a computed FULL matrix of order 4 by 4 It has no free parameters specified MATRIX S This is a computed FULL matrix of order 2 by 6 It has no free parameters specified MATRIX T This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX U This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX V This is a computed FULL matrix of order 2 by 2 It has no free parameters specified Mx starting optimization; number of parameters = 12 *** 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 2 by 2 [=X*X'] 1 2 1 0.9676 0.1881 2 0.1881 0.0365 MATRIX B This is a FULL matrix of order 1 by 1 1 1 -0.1869 MATRIX C This is a computed FULL matrix of order 2 by 2 [=Y*Y'] 1 2 1 5.8694 4.2806 2 4.2806 5.0582 MATRIX E This is a computed FULL matrix of order 2 by 2 [=Z*Z'] 1 2 1 14.1293 5.8690 2 5.8690 12.5367 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 2 1 2 1 11.0151 10.2235 MATRIX X This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.9837 2 0.1912 0.0000 MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 2.4227 2 1.7669 -1.3915 MATRIX Z This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 3.7589 2 1.5614 3.1779 GROUP NUMBER: 2 G2: MZ twins, datagroup MATRIX A This is a computed FULL matrix of order 2 by 2 [=X*X'] 1 2 1 0.9676 0.1881 2 0.1881 0.0365 MATRIX B This is a FULL matrix of order 1 by 1 1 1 -0.1869 MATRIX C This is a computed FULL matrix of order 2 by 2 [=Y*Y'] 1 2 1 5.8694 4.2806 2 4.2806 5.0582 MATRIX E This is a computed FULL matrix of order 2 by 2 [=Z*Z'] 1 2 1 14.1293 5.8690 2 5.8690 12.5367 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 2 1 2 1 11.0151 10.2235 MATRIX O This is a FULL matrix of order 1 by 2 1 2 1 16.2300 999.0000 MATRIX P This is a FULL matrix of order 1 by 2 1 2 1 16.2300 999.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.9837 2 0.1912 0.0000 MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 2.4227 2 1.7669 -1.3915 MATRIX Z This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 3.7589 2 1.5614 3.1779 Vector of OBSERVED means NYAGG1.1 NYAGG3.1 NYAGG1.2 NYAGG3.2 Mean 7.6121 6.2330 7.4514 6.4138 Vector of EXPECTED means NYAGG1.1 NYAGG3.1 NYAGG1.2 NYAGG3.2 Mean 7.9818 -176.4838 7.9818 -176.4838 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX NYAGG1.1 NYAGG3.1 NYAGG1.2 NYAGG3.2 NYAGG1.1 20.9663 NYAGG3.1 10.3377 17.6314 NYAGG1.2 6.8370 4.4687 20.9663 NYAGG3.2 4.4687 5.0948 10.3377 17.6314 Function value of this group: 11551.6650 Where the fit function is -2 * Log-likelihood of raw data GROUP NUMBER: 3 G3: DZ twins, datagroup MATRIX A This is a computed FULL matrix of order 2 by 2 [=X*X'] 1 2 1 0.9676 0.1881 2 0.1881 0.0365 MATRIX B This is a FULL matrix of order 1 by 1 1 1 -0.1869 MATRIX C This is a computed FULL matrix of order 2 by 2 [=Y*Y'] 1 2 1 5.8694 4.2806 2 4.2806 5.0582 MATRIX E This is a computed FULL matrix of order 2 by 2 [=Z*Z'] 1 2 1 14.1293 5.8690 2 5.8690 12.5367 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 2 1 2 1 11.0151 10.2235 MATRIX O This is a FULL matrix of order 1 by 2 1 2 1 18.3900 999.0000 MATRIX P This is a FULL matrix of order 1 by 2 1 2 1 18.3900 999.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.9837 2 0.1912 0.0000 MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 2.4227 2 1.7669 -1.3915 MATRIX Z This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 3.7589 2 1.5614 3.1779 Vector of OBSERVED means NYAGG1.1 NYAGG3.1 NYAGG1.2 NYAGG3.2 Mean 7.6462 6.6732 7.7991 6.4779 Vector of EXPECTED means NYAGG1.1 NYAGG3.1 NYAGG1.2 NYAGG3.2 Mean 7.5781 -176.4838 7.5781 -176.4838 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX NYAGG1.1 NYAGG3.1 NYAGG1.2 NYAGG3.2 NYAGG1.1 20.9663 NYAGG3.1 10.3377 17.6314 NYAGG1.2 6.3532 4.3747 20.9663 NYAGG3.2 4.3747 5.0765 10.3377 17.6314 Function value of this group: 25560.3239 Where the fit function is -2 * Log-likelihood of raw data GROUP NUMBER: 4 G4: Calculate A, C, and E correlations MATRIX A This is a computed FULL matrix of order 2 by 2 [=X*X'] 1 2 1 0.9676 0.1881 2 0.1881 0.0365 MATRIX C This is a computed FULL matrix of order 2 by 2 [=Y*Y'] 1 2 1 5.8694 4.2806 2 4.2806 5.0582 MATRIX E This is a computed FULL matrix of order 2 by 2 [=Z*Z'] 1 2 1 14.1293 5.8690 2 5.8690 12.5367 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 2 by 2 MATRIX P This is a computed FULL matrix of order 4 by 4 [=\STND(A+C+E|H@A+C_H@A+C|A+C+E)] 1 2 3 4 1 1.0000 0.5377 0.3030 0.2275 2 0.5377 1.0000 0.2275 0.2879 3 0.3030 0.2275 1.0000 0.5377 4 0.2275 0.2879 0.5377 1.0000 MATRIX S This is a computed FULL matrix of order 2 by 6 [=A%(A+C+E)|C%(A+C+E)|E%(A+C+E)] 1 2 3 4 5 6 1 0.0462 0.0182 0.2799 0.4141 0.6739 0.5677 2 0.0182 0.0021 0.4141 0.2869 0.5677 0.7110 MATRIX T This is a computed FULL matrix of order 2 by 2 [=\STND(A)] 1 2 1 1.0000 1.0000 2 1.0000 1.0000 MATRIX U This is a computed FULL matrix of order 2 by 2 [=\STND(C)] 1 2 1 1.0000 0.7856 2 0.7856 1.0000 MATRIX V This is a computed FULL matrix of order 2 by 2 [=\STND(E)] 1 2 1 1.0000 0.4410 2 0.4410 1.0000 *** WARNING! *** Minimization may not be successful. See above CODE GREEN - it probably was OK Your model has 12 estimated parameters and 6561 Observed statistics -2 times log-likelihood of data >>> 37111.989 Degrees of freedom >>>>>>>>>>>>>>>> 6549 Akaike's Information Criterion >>>> 24013.989 Bayesian Information Criterion >>>> -6927.443 Sample size Adjusted BIC >>>> 3476.365 Deviance Information Criterion >>>> -909.314 This problem used 5.0% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 1.33 Execution 0: 0: 0:22.23 TOTAL 0: 0: 0:23.56 Total number of warnings issued: 2 ______________________________________________________________________________ Multiple fit option in effect. The following MX script lines have been read: SAVE BIVARIATE.MXS END Summary of VL file data for group 2 AGE3.2 AGE1.2 AGE3.1 AGE1.1 NYAGG1.1 NYAGG3.1 NYAGG1.2 Code -4.0000 -3.0000 -2.0000 -1.0000 1.0000 2.0000 3.0000 Number 743.0000 743.0000 743.0000 743.0000 504.0000 501.0000 503.0000 Mean 336.0279 328.0273 329.4289 329.3497 7.6121 6.2330 7.4514 Variance **************************************** 25.6824 18.3844 20.3358 Minimum 14.2700 13.2500 14.3100 13.4000 0.0000 0.0000 0.0000 Maximum 999.0000 999.0000 999.0000 999.0000 27.4400 27.4400 25.3300 NYAGG3.2 Code 4.0000 Number 493.0000 Mean 6.4138 Variance 21.4903 Minimum 0.0000 Maximum 33.7800 Summary of VL file data for group 3 AGE3.2 AGE1.2 AGE3.1 AGE1.1 NYAGG1.1 NYAGG3.1 NYAGG1.2 Code -4.0000 -3.0000 -2.0000 -1.0000 1.0000 2.0000 3.0000 Number 1655.0000 1655.0000 1655.0000 1655.0000 1125.0000 1164.0000 1126.0000 Mean 300.2697 327.8228 296.7262 327.2353 7.6462 6.6732 7.7991 Variance **************************************** 19.4824 16.7076 19.7052 Minimum 14.2900 12.6000 14.2900 12.6000 0.0000 0.0000 0.0000 Maximum 999.0000 999.0000 999.0000 999.0000 32.7200 27.4400 29.5600 NYAGG3.2 Code 4.0000 Number 1145.0000 Mean 6.4779 Variance 17.5274 Minimum 0.0000 Maximum 26.3900 PARAMETER SPECIFICATIONS GROUP NUMBER: 1 G1: calculation group MATRIX A This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX B This is a FULL matrix of order 1 by 1 1 1 12 MATRIX C This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 2 by 2 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 2 1 2 1 10 11 MATRIX X This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 1 2 2 3 MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 4 2 5 6 MATRIX Z This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 7 2 8 9 GROUP NUMBER: 2 G2: MZ twins, datagroup MATRIX A This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX B This is a FULL matrix of order 1 by 1 1 1 12 MATRIX C This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 2 by 2 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 2 1 2 1 10 11 MATRIX O This is a FULL matrix of order 1 by 2 1 2 1 -1 -2 MATRIX P This is a FULL matrix of order 1 by 2 1 2 1 -3 -4 MATRIX X This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 1 2 2 3 MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 4 2 5 6 MATRIX Z This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 7 2 8 9 GROUP NUMBER: 3 G3: DZ twins, datagroup MATRIX A This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX B This is a FULL matrix of order 1 by 1 1 1 12 MATRIX C This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 2 by 2 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 2 1 2 1 10 11 MATRIX O This is a FULL matrix of order 1 by 2 1 2 1 -1 -2 MATRIX P This is a FULL matrix of order 1 by 2 1 2 1 -3 -4 MATRIX X This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 1 2 2 3 MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 4 2 5 6 MATRIX Z This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 7 2 8 9 GROUP NUMBER: 4 G4: Calculate A, C, and E correlations MATRIX A This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 2 by 2 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 2 by 2 MATRIX P This is a computed FULL matrix of order 4 by 4 It has no free parameters specified MATRIX S This is a computed FULL matrix of order 2 by 6 It has no free parameters specified MATRIX T This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX U This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX V This is a computed FULL matrix of order 2 by 2 It has no free parameters specified Mx starting optimization; number of parameters = 12 MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1: calculation group MATRIX A This is a computed FULL matrix of order 2 by 2 [=X*X'] 1 2 1 0.9676 0.1881 2 0.1881 0.0366 MATRIX B This is a FULL matrix of order 1 by 1 1 1 -0.1869 MATRIX C This is a computed FULL matrix of order 2 by 2 [=Y*Y'] 1 2 1 5.8694 4.2806 2 4.2806 5.0582 MATRIX E This is a computed FULL matrix of order 2 by 2 [=Z*Z'] 1 2 1 14.1293 5.8690 2 5.8690 12.5367 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 2 1 2 1 11.0151 10.2235 MATRIX X This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.9837 2 0.1912 0.0100 MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 2.4227 2 1.7669 -1.3915 MATRIX Z This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 3.7589 2 1.5614 3.1779 GROUP NUMBER: 2 G2: MZ twins, datagroup MATRIX A This is a computed FULL matrix of order 2 by 2 [=X*X'] 1 2 1 0.9676 0.1881 2 0.1881 0.0366 MATRIX B This is a FULL matrix of order 1 by 1 1 1 -0.1869 MATRIX C This is a computed FULL matrix of order 2 by 2 [=Y*Y'] 1 2 1 5.8694 4.2806 2 4.2806 5.0582 MATRIX E This is a computed FULL matrix of order 2 by 2 [=Z*Z'] 1 2 1 14.1293 5.8690 2 5.8690 12.5367 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 2 1 2 1 11.0151 10.2235 MATRIX O This is a FULL matrix of order 1 by 2 1 2 1 16.2300 999.0000 MATRIX P This is a FULL matrix of order 1 by 2 1 2 1 16.2300 999.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.9837 2 0.1912 0.0100 MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 2.4227 2 1.7669 -1.3915 MATRIX Z This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 3.7589 2 1.5614 3.1779 Vector of OBSERVED means NYAGG1.1 NYAGG3.1 NYAGG1.2 NYAGG3.2 Mean 7.6121 6.2330 7.4514 6.4138 Vector of EXPECTED means NYAGG1.1 NYAGG3.1 NYAGG1.2 NYAGG3.2 Mean 7.9818 -176.4837 7.9818 -176.4837 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX NYAGG1.1 NYAGG3.1 NYAGG1.2 NYAGG3.2 NYAGG1.1 20.9663 NYAGG3.1 10.3377 17.6315 NYAGG1.2 6.8370 4.4687 20.9663 NYAGG3.2 4.4687 5.0949 10.3377 17.6315 Function value of this group: 11551.6643 Where the fit function is -2 * Log-likelihood of raw data GROUP NUMBER: 3 G3: DZ twins, datagroup MATRIX A This is a computed FULL matrix of order 2 by 2 [=X*X'] 1 2 1 0.9676 0.1881 2 0.1881 0.0366 MATRIX B This is a FULL matrix of order 1 by 1 1 1 -0.1869 MATRIX C This is a computed FULL matrix of order 2 by 2 [=Y*Y'] 1 2 1 5.8694 4.2806 2 4.2806 5.0582 MATRIX E This is a computed FULL matrix of order 2 by 2 [=Z*Z'] 1 2 1 14.1293 5.8690 2 5.8690 12.5367 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 2 1 2 1 11.0151 10.2235 MATRIX O This is a FULL matrix of order 1 by 2 1 2 1 18.3900 999.0000 MATRIX P This is a FULL matrix of order 1 by 2 1 2 1 18.3900 999.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.9837 2 0.1912 0.0100 MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 2.4227 2 1.7669 -1.3915 MATRIX Z This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 3.7589 2 1.5614 3.1779 Vector of OBSERVED means NYAGG1.1 NYAGG3.1 NYAGG1.2 NYAGG3.2 Mean 7.6462 6.6732 7.7991 6.4779 Vector of EXPECTED means NYAGG1.1 NYAGG3.1 NYAGG1.2 NYAGG3.2 Mean 7.5781 -176.4837 7.5781 -176.4837 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX NYAGG1.1 NYAGG3.1 NYAGG1.2 NYAGG3.2 NYAGG1.1 20.9663 NYAGG3.1 10.3377 17.6315 NYAGG1.2 6.3532 4.3747 20.9663 NYAGG3.2 4.3747 5.0765 10.3377 17.6315 Function value of this group: 25560.3246 Where the fit function is -2 * Log-likelihood of raw data GROUP NUMBER: 4 G4: Calculate A, C, and E correlations MATRIX A This is a computed FULL matrix of order 2 by 2 [=X*X'] 1 2 1 0.9676 0.1881 2 0.1881 0.0366 MATRIX C This is a computed FULL matrix of order 2 by 2 [=Y*Y'] 1 2 1 5.8694 4.2806 2 4.2806 5.0582 MATRIX E This is a computed FULL matrix of order 2 by 2 [=Z*Z'] 1 2 1 14.1293 5.8690 2 5.8690 12.5367 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 2 by 2 MATRIX P This is a computed FULL matrix of order 4 by 4 [=\STND(A+C+E|H@A+C_H@A+C|A+C+E)] 1 2 3 4 1 1.0000 0.5377 0.3030 0.2275 2 0.5377 1.0000 0.2275 0.2879 3 0.3030 0.2275 1.0000 0.5377 4 0.2275 0.2879 0.5377 1.0000 MATRIX S This is a computed FULL matrix of order 2 by 6 [=A%(A+C+E)|C%(A+C+E)|E%(A+C+E)] 1 2 3 4 5 6 1 0.0462 0.0182 0.2799 0.4141 0.6739 0.5677 2 0.0182 0.0021 0.4141 0.2869 0.5677 0.7110 MATRIX T This is a computed FULL matrix of order 2 by 2 [=\STND(A)] 1 2 1 1.0000 0.9986 2 0.9986 1.0000 MATRIX U This is a computed FULL matrix of order 2 by 2 [=\STND(C)] 1 2 1 1.0000 0.7856 2 0.7856 1.0000 MATRIX V This is a computed FULL matrix of order 2 by 2 [=\STND(E)] 1 2 1 1.0000 0.4410 2 0.4410 1.0000 Your model has 12 estimated parameters and 6561 Observed statistics -2 times log-likelihood of data >>> 37111.989 Degrees of freedom >>>>>>>>>>>>>>>> 6549 Akaike's Information Criterion >>>> 24013.989 Bayesian Information Criterion >>>> -6927.443 Sample size Adjusted BIC >>>> 3476.365 Deviance Information Criterion >>>> -909.314 Saturated model fit* >>>>>>>>>>> 37111.989 Saturated model df* >>>>>>>>>>> 6549 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 5.0% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.13 Execution 0: 0: 0: 1.60 TOTAL 0: 0: 0: 1.73 Total number of warnings issued: 2 ______________________________________________________________________________ !AE DROP SHARED ENVIRONMENTAL EFFECTS !GET BIVARIATE.MXS !DROP Y 1 1 1 - Y 1 NVAR NVAR !CE DROP ADDITIVE GENETIC EFFECTS ! ***WRITE SCRIPT*** !GET BIVARIATE.MXS !DROP X 1 1 1 X 1 2 2 X 1 2 1 !E DROP A AND C ! ***WRITE SCRIPT*** !GET BIVARIATE.MXS !DROP Y 1 1 1 - Y 1 NVAR NVAR !DROP X 1 1 1 - X 1 NVAR NVAR