** Mx startup successful ** **MX-PC 1.65b** Job started on 03/05/08 at 10:21:47 ! 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 Y LOWER NVAR NVAR FREE ! COMMON ENVIRONMENTAL PATH COEFFICIENT Z LOWER NVAR NVAR FREE ! UNIQUE ENVIRONMENTAL PATH COEFFICIENT 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 START 10 M 1 1 M 1 2 START -0.1 B 1 1 ! FOR THE VARIANCE COMPONENTS: DIAGONALS ! PHENOTYPIC VARIANCE ST 2.0 X 1 1 X 2 2 ST 2.0 Y 1 1 Y 2 2 ST 3.5 Z 1 1 Z 2 2 ! PHENOTYPIC COVARIANCE: OFF DIAGONALS ST 1.5 X 2 1 ST 1.5 Y 2 1 ST 1.5 Z 2 1 BEGIN ALGEBRA; ! **GO TO SLIDE** "TOTAL WITHIN-TWIN COVARIANCE" ! AIM IS TO ESTIMATE VARIANCE COMPONENTS ATTRIBUTABLE TO A, C AND E ! INSTRUCT MX HOW TO USE THE PATH COEFFICIENTS TO ESTIMATE VAR/COVARIANCE FOR A, C AND E 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 ! ***LOOK AT DATA FILE*** 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) MATRIX FOR TWIN 1 I.E. AGE FOR TWIN 1 AT TIMES 1 & 2 P FULL NDEF NVAR FIX ! DEFINITION VARIABLE(S) MATRIX 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 SELECT IF ZYG2 = 2 ; ! 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 ! **GO TO SLIDE** "PREDICTED MODEL" ! ***WRITE THE DZ COVARIANCE STATEMENT*** COVARIANCES 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 LATENT FACTOR 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 ; ! **GO TO SLIDE** "ESTIMATING CORRELATIONS BETWEEN LATENT GENETIC AND ENVIRONMENTAL FACTORS" 2 SLIDES ! T = \STND(A); ! ESTIMATE GENETIC CORRELATIONS (RG) OR T =\SQRT(I.A)~*A*\SQRT(I.A)~ ; ! U = ; ! ESTIMATE SHARED ENVIRONMENTAL CORRELATIONS (RC) ! V = ; ! ESTIMATE NON-SHARED ENVIRONMENTAL CORRELATIONS P = \STND(A+C+E | H@A+C _H@A+C | A+C+E); ! ESTIMATES PHENOTYPIC CORRELATIONS I.E STANDARDIZES THE COVARIANCE MATRIX ! **GO TO SLIDE** "PROPORTION OF OBSERVED PHENOTYPIC CORRELATION EXPLAINED BY A, C AND E" S = A%(A+C+E)| C%(A+C+E)|E%(A+C+E); ! ESTIMATE STANDARDIZED VARIANCE COMPONENTS, AND ! PROPORTION OF PHENOYPTIC CORRELATIONS EXPLAINED BY A, C AND E EFFECTS 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 latent factor 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 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.9675 0.1880 2 0.1880 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.8695 4.2807 2 4.2807 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.2236 MATRIX X This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.9836 2 0.1911 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.9675 0.1880 2 0.1880 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.8695 4.2807 2 4.2807 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.2236 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.9836 2 0.1911 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.4855 7.9818 -176.4855 (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.9675 0.1880 2 0.1880 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.8695 4.2807 2 4.2807 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.2236 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.9836 2 0.1911 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.4855 7.5781 -176.4855 (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 latent factor correlations MATRIX A This is a computed FULL matrix of order 2 by 2 [=X*X'] 1 2 1 0.9675 0.1880 2 0.1880 0.0365 MATRIX C This is a computed FULL matrix of order 2 by 2 [=Y*Y'] 1 2 1 5.8695 4.2807 2 4.2807 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.0461 0.0182 0.2799 0.4141 0.6739 0.5677 2 0.0182 0.0021 0.4141 0.2869 0.5677 0.7110 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.35 Execution 0: 0: 0:25.84 TOTAL 0: 0: 0:27.19 Total number of warnings issued: 0 ______________________________________________________________________________ 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 latent factor 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 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.9675 0.1880 2 0.1880 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.8695 4.2807 2 4.2807 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.2236 MATRIX X This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.9836 2 0.1911 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.9675 0.1880 2 0.1880 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.8695 4.2807 2 4.2807 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.2236 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.9836 2 0.1911 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.4853 7.9818 -176.4853 (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.9675 0.1880 2 0.1880 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.8695 4.2807 2 4.2807 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.2236 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.9836 2 0.1911 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.4853 7.5781 -176.4853 (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 latent factor correlations MATRIX A This is a computed FULL matrix of order 2 by 2 [=X*X'] 1 2 1 0.9675 0.1880 2 0.1880 0.0366 MATRIX C This is a computed FULL matrix of order 2 by 2 [=Y*Y'] 1 2 1 5.8695 4.2807 2 4.2807 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.0461 0.0182 0.2799 0.4141 0.6739 0.5677 2 0.0182 0.0021 0.4141 0.2869 0.5677 0.7110 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.15 Execution 0: 0: 0: 1.59 TOTAL 0: 0: 0: 1.74 Total number of warnings issued: 0 ______________________________________________________________________________ ! ***HOW MANY PARAMETERS DOES THIS MODEL ESTIMATE?*** ! ***RUN MODEL*** !AE DROP SHARED ENVIRONMENTAL EFFECTS Multiple fit option in effect. The following MX script lines have been read: GET BIVARIATE.MXS Note: 254684 binary records read DROP Y 1 1 1 - Y 1 NVAR NVAR 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 It has no free parameters specified 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 It has no free parameters specified 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 It has no free parameters specified 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 latent factor 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 Mx starting optimization; number of parameters = 9 *** 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 8.1089 5.3283 2 5.3283 6.1132 MATRIX B This is a FULL matrix of order 1 by 1 1 1 -0.1867 MATRIX C This is a computed FULL matrix of order 2 by 2 [=Y*Y'] 1 2 1 0.0000 0.0000 2 0.0000 0.0000 MATRIX E This is a computed FULL matrix of order 2 by 2 [=Z*Z'] 1 2 1 12.7064 4.9057 2 4.9057 11.4265 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.0182 10.2274 MATRIX X This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 2.8476 2 1.8711 1.6162 MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.0000 2 0.0000 0.0000 MATRIX Z This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 3.5646 2 1.3762 3.0875 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 8.1089 5.3283 2 5.3283 6.1132 MATRIX B This is a FULL matrix of order 1 by 1 1 1 -0.1867 MATRIX C This is a computed FULL matrix of order 2 by 2 [=Y*Y'] 1 2 1 0.0000 0.0000 2 0.0000 0.0000 MATRIX E This is a computed FULL matrix of order 2 by 2 [=Z*Z'] 1 2 1 12.7064 4.9057 2 4.9057 11.4265 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.0182 10.2274 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 2.8476 2 1.8711 1.6162 MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.0000 2 0.0000 0.0000 MATRIX Z This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 3.5646 2 1.3762 3.0875 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.9880 -176.2930 7.9880 -176.2930 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX NYAGG1.1 NYAGG3.1 NYAGG1.2 NYAGG3.2 NYAGG1.1 20.8153 NYAGG3.1 10.2340 17.5397 NYAGG1.2 8.1089 5.3283 20.8153 NYAGG3.2 5.3283 6.1132 10.2340 17.5397 Function value of this group: 11560.5356 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 8.1089 5.3283 2 5.3283 6.1132 MATRIX B This is a FULL matrix of order 1 by 1 1 1 -0.1867 MATRIX C This is a computed FULL matrix of order 2 by 2 [=Y*Y'] 1 2 1 0.0000 0.0000 2 0.0000 0.0000 MATRIX E This is a computed FULL matrix of order 2 by 2 [=Z*Z'] 1 2 1 12.7064 4.9057 2 4.9057 11.4265 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.0182 10.2274 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 2.8476 2 1.8711 1.6162 MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.0000 2 0.0000 0.0000 MATRIX Z This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 3.5646 2 1.3762 3.0875 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.5847 -176.2930 7.5847 -176.2930 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX NYAGG1.1 NYAGG3.1 NYAGG1.2 NYAGG3.2 NYAGG1.1 20.8153 NYAGG3.1 10.2340 17.5397 NYAGG1.2 4.0545 2.6641 20.8153 NYAGG3.2 2.6641 3.0566 10.2340 17.5397 Function value of this group: 25586.5449 Where the fit function is -2 * Log-likelihood of raw data GROUP NUMBER: 4 G4: Calculate A, C, and E latent factor correlations MATRIX A This is a computed FULL matrix of order 2 by 2 [=X*X'] 1 2 1 8.1089 5.3283 2 5.3283 6.1132 MATRIX C This is a computed FULL matrix of order 2 by 2 [=Y*Y'] 1 2 1 0.0000 0.0000 2 0.0000 0.0000 MATRIX E This is a computed FULL matrix of order 2 by 2 [=Z*Z'] 1 2 1 12.7064 4.9057 2 4.9057 11.4265 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.5356 0.1948 0.1394 2 0.5356 1.0000 0.1394 0.1743 3 0.1948 0.1394 1.0000 0.5356 4 0.1394 0.1743 0.5356 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.3896 0.5206 0.0000 0.0000 0.6104 0.4794 2 0.5206 0.3485 0.0000 0.0000 0.4794 0.6515 *** WARNING! *** Minimization may not be successful. See above CODE GREEN - it probably was OK Your model has 9 estimated parameters and 6561 Observed statistics -2 times log-likelihood of data >>> 37147.081 Degrees of freedom >>>>>>>>>>>>>>>> 6552 Akaike's Information Criterion >>>> 24043.081 Bayesian Information Criterion >>>> -6921.570 Sample size Adjusted BIC >>>> 3487.003 Deviance Information Criterion >>>> -900.685 Saturated model fit* >>>>>>>>>>> 37111.989 Saturated model df* >>>>>>>>>>> 6549 Difference Chi-squared >>>>>>>> 35.092 Difference d.f. >>>>>>>>>>>>>>> 3 Probability >>>>>>>>>>>>>>>>>>>> 0.000 Akaike's Information Criterion > 29.092 * 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.11 Execution 0: 0: 0:17.77 TOTAL 0: 0: 0:17.88 Total number of warnings issued: 2 ______________________________________________________________________________ !CE DROP ADDITIVE GENETIC EFFECTS ! ***WRITE SCRIPT*** Multiple fit option in effect. The following MX script lines have been read: GET BIVARIATE.MXS Note: 254684 binary records read DROP X 1 1 1 - X 1 NVAR NVAR 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 It has no free parameters specified 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 It has no free parameters specified 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 It has no free parameters specified 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 latent factor 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 Mx starting optimization; number of parameters = 9 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.0000 0.0000 2 0.0000 0.0000 MATRIX B This is a FULL matrix of order 1 by 1 1 1 -0.1870 MATRIX C This is a computed FULL matrix of order 2 by 2 [=Y*Y'] 1 2 1 6.5337 4.4041 2 4.4041 5.0813 MATRIX E This is a computed FULL matrix of order 2 by 2 [=Z*Z'] 1 2 1 14.4416 5.9368 2 5.9368 12.5509 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.0159 10.2244 MATRIX X This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.0000 2 0.0000 0.0000 MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 2.5561 2 1.7230 1.4535 MATRIX Z This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 3.8002 2 1.5622 3.1797 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.0000 0.0000 2 0.0000 0.0000 MATRIX B This is a FULL matrix of order 1 by 1 1 1 -0.1870 MATRIX C This is a computed FULL matrix of order 2 by 2 [=Y*Y'] 1 2 1 6.5337 4.4041 2 4.4041 5.0813 MATRIX E This is a computed FULL matrix of order 2 by 2 [=Z*Z'] 1 2 1 14.4416 5.9368 2 5.9368 12.5509 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.0159 10.2244 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.0000 2 0.0000 0.0000 MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 2.5561 2 1.7230 1.4535 MATRIX Z This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 3.8002 2 1.5622 3.1797 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.9815 -176.5508 7.9815 -176.5508 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX NYAGG1.1 NYAGG3.1 NYAGG1.2 NYAGG3.2 NYAGG1.1 20.9753 NYAGG3.1 10.3409 17.6321 NYAGG1.2 6.5337 4.4041 20.9753 NYAGG3.2 4.4041 5.0813 10.3409 17.6321 Function value of this group: 11551.7957 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.0000 0.0000 2 0.0000 0.0000 MATRIX B This is a FULL matrix of order 1 by 1 1 1 -0.1870 MATRIX C This is a computed FULL matrix of order 2 by 2 [=Y*Y'] 1 2 1 6.5337 4.4041 2 4.4041 5.0813 MATRIX E This is a computed FULL matrix of order 2 by 2 [=Z*Z'] 1 2 1 14.4416 5.9368 2 5.9368 12.5509 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.0159 10.2244 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.0000 2 0.0000 0.0000 MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 2.5561 2 1.7230 1.4535 MATRIX Z This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 3.8002 2 1.5622 3.1797 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.5777 -176.5508 7.5777 -176.5508 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX NYAGG1.1 NYAGG3.1 NYAGG1.2 NYAGG3.2 NYAGG1.1 20.9753 NYAGG3.1 10.3409 17.6321 NYAGG1.2 6.5337 4.4041 20.9753 NYAGG3.2 4.4041 5.0813 10.3409 17.6321 Function value of this group: 25560.4943 Where the fit function is -2 * Log-likelihood of raw data GROUP NUMBER: 4 G4: Calculate A, C, and E latent factor correlations MATRIX A This is a computed FULL matrix of order 2 by 2 [=X*X'] 1 2 1 0.0000 0.0000 2 0.0000 0.0000 MATRIX C This is a computed FULL matrix of order 2 by 2 [=Y*Y'] 1 2 1 6.5337 4.4041 2 4.4041 5.0813 MATRIX E This is a computed FULL matrix of order 2 by 2 [=Z*Z'] 1 2 1 14.4416 5.9368 2 5.9368 12.5509 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.3115 0.2290 2 0.5377 1.0000 0.2290 0.2882 3 0.3115 0.2290 1.0000 0.5377 4 0.2290 0.2882 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.0000 0.0000 0.3115 0.4259 0.6885 0.5741 2 0.0000 0.0000 0.4259 0.2882 0.5741 0.7118 Your model has 9 estimated parameters and 6561 Observed statistics -2 times log-likelihood of data >>> 37112.290 Degrees of freedom >>>>>>>>>>>>>>>> 6552 Akaike's Information Criterion >>>> 24008.290 Bayesian Information Criterion >>>> -6938.966 Sample size Adjusted BIC >>>> 3469.607 Deviance Information Criterion >>>> -918.081 Saturated model fit* >>>>>>>>>>> 37111.989 Saturated model df* >>>>>>>>>>> 6549 Difference Chi-squared >>>>>>>> 0.301 Difference d.f. >>>>>>>>>>>>>>> 3 Probability >>>>>>>>>>>>>>>>>>>> 0.960 Akaike's Information Criterion > -5.699 * 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.11 Execution 0: 0: 0:10.22 TOTAL 0: 0: 0:10.33 Total number of warnings issued: 2 ______________________________________________________________________________ !E DROP A AND C ! ***WRITE SCRIPT*** Multiple fit option in effect. The following MX script lines have been read: GET BIVARIATE.MXS Note: 254684 binary records read DROP Y 1 1 1 - Y 1 NVAR NVAR DROP X 1 1 1 - X 1 NVAR NVAR 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 It has no free parameters specified MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 It has no free parameters specified 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 It has no free parameters specified MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 It has no free parameters specified 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 It has no free parameters specified MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 It has no free parameters specified 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 latent factor 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 Mx starting optimization; number of parameters = 6 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.0000 0.0000 2 0.0000 0.0000 MATRIX B This is a FULL matrix of order 1 by 1 1 1 -0.1872 MATRIX C This is a computed FULL matrix of order 2 by 2 [=Y*Y'] 1 2 1 0.0000 0.0000 2 0.0000 0.0000 MATRIX E This is a computed FULL matrix of order 2 by 2 [=Z*Z'] 1 2 1 20.8778 10.2788 2 10.2788 17.6567 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.0138 10.2245 MATRIX X This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.0000 2 0.0000 0.0000 MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.0000 2 0.0000 0.0000 MATRIX Z This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 4.5692 2 2.2496 3.5491 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.0000 0.0000 2 0.0000 0.0000 MATRIX B This is a FULL matrix of order 1 by 1 1 1 -0.1872 MATRIX C This is a computed FULL matrix of order 2 by 2 [=Y*Y'] 1 2 1 0.0000 0.0000 2 0.0000 0.0000 MATRIX E This is a computed FULL matrix of order 2 by 2 [=Z*Z'] 1 2 1 20.8778 10.2788 2 10.2788 17.6567 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.0138 10.2245 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.0000 2 0.0000 0.0000 MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.0000 2 0.0000 0.0000 MATRIX Z This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 4.5692 2 2.2496 3.5491 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.9751 -176.8204 7.9751 -176.8204 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX NYAGG1.1 NYAGG3.1 NYAGG1.2 NYAGG3.2 NYAGG1.1 20.8778 NYAGG3.1 10.2788 17.6567 NYAGG1.2 0.0000 0.0000 20.8778 NYAGG3.2 0.0000 0.0000 10.2788 17.6567 Function value of this group: 11648.9171 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.0000 0.0000 2 0.0000 0.0000 MATRIX B This is a FULL matrix of order 1 by 1 1 1 -0.1872 MATRIX C This is a computed FULL matrix of order 2 by 2 [=Y*Y'] 1 2 1 0.0000 0.0000 2 0.0000 0.0000 MATRIX E This is a computed FULL matrix of order 2 by 2 [=Z*Z'] 1 2 1 20.8778 10.2788 2 10.2788 17.6567 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.0138 10.2245 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.0000 2 0.0000 0.0000 MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.0000 2 0.0000 0.0000 MATRIX Z This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 4.5692 2 2.2496 3.5491 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.5706 -176.8204 7.5706 -176.8204 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX NYAGG1.1 NYAGG3.1 NYAGG1.2 NYAGG3.2 NYAGG1.1 20.8778 NYAGG3.1 10.2788 17.6567 NYAGG1.2 0.0000 0.0000 20.8778 NYAGG3.2 0.0000 0.0000 10.2788 17.6567 Function value of this group: 25727.5852 Where the fit function is -2 * Log-likelihood of raw data GROUP NUMBER: 4 G4: Calculate A, C, and E latent factor correlations MATRIX A This is a computed FULL matrix of order 2 by 2 [=X*X'] 1 2 1 0.0000 0.0000 2 0.0000 0.0000 MATRIX C This is a computed FULL matrix of order 2 by 2 [=Y*Y'] 1 2 1 0.0000 0.0000 2 0.0000 0.0000 MATRIX E This is a computed FULL matrix of order 2 by 2 [=Z*Z'] 1 2 1 20.8778 10.2788 2 10.2788 17.6567 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.5354 0.0000 0.0000 2 0.5354 1.0000 0.0000 0.0000 3 0.0000 0.0000 1.0000 0.5354 4 0.0000 0.0000 0.5354 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.0000 0.0000 0.0000 0.0000 1.0000 1.0000 2 0.0000 0.0000 0.0000 0.0000 1.0000 1.0000 Your model has 6 estimated parameters and 6561 Observed statistics -2 times log-likelihood of data >>> 37376.502 Degrees of freedom >>>>>>>>>>>>>>>> 6555 Akaike's Information Criterion >>>> 24266.502 Bayesian Information Criterion >>>> -6818.533 Sample size Adjusted BIC >>>> 3594.806 Deviance Information Criterion >>>> -794.891 Saturated model fit* >>>>>>>>>>> 37111.989 Saturated model df* >>>>>>>>>>> 6549 Difference Chi-squared >>>>>>>> 264.513 Difference d.f. >>>>>>>>>>>>>>> 6 Probability >>>>>>>>>>>>>>>>>>>> 0.000 Akaike's Information Criterion > 252.513 * 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.11 Execution 0: 0: 0: 7.27 TOTAL 0: 0: 0: 7.38 Total number of warnings issued: 2 ______________________________________________________________________________ ! ***RUN MODEL***