**MX-PC 1.65b** Job started on 03/05/08 at 10:38:00 ! 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 3 ! 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 M 1 3 START -0.1 B 1 1 ! FOR THE VARIANCE COMPONENTS: DIAGONALS ! PHENOTYPIC VARIANCE ST 2.0 X 1 1 X 2 2 X 3 3 ST 2.0 Y 1 1 Y 2 2 Y 2 2 ST 3.5 Z 1 1 Z 2 2 Z 3 3 ! 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 NYAGG4.1 ! NYAGG1.1 = AGGRESSION TWIN1 '91, NYAGG3.1 = AGGRESSION1 TWIN '95 NYAGG1.2 NYAGG3.2 NYAGG4.2 ! NYAGG1.2 = AGGRESSION TWIN2 '91, NYAGG3.2 = AGGRESSION2 TWIN '95 AGE1.1 AGE3.1 AGE4.1 ! AGE1.1 = AGE AT INTERVIEW TWIN1 '91, AGE3.1 = AGE AT INTERVIEW TWIN1 '95 AGE1.2 AGE3.2 AGE4.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 9262 observations AGE1.1 AGE3.1 AGE4.1 AGE1.2 AGE3.2 AGE4.2 ; NOTE: Definition yields 877 data vectors for analysis NOTE: Vectors contain a total of 2800 observations MATRICES = GROUP 1 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 AGE4.1 ! AGE FOR TWIN 1 AT TIMES 1 & 2 SPECIFY P AGE1.2 AGE3.2 AGE4.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 NYAGG4.1 ! NYAGG1.1 = AGGRESSION TWIN1 '91, NYAGG3.1 = AGGRESSION1 TWIN '95 NYAGG1.2 NYAGG3.2 NYAGG4.2 ! NYAGG1.2 = AGGRESSION TWIN2 '91, NYAGG3.2 = AGGRESSION2 TWIN '95 AGE1.1 AGE3.1 AGE4.1 ! AGE1.1 = AGE AT INTERVIEW TWIN1 '91, AGE3.1 = AGE AT INTERVIEW TWIN1 '95 AGE1.2 AGE3.2 AGE4.2 ; ! AGE1.2 = AGE AT INTERVIEW TWIN2 '91, AGE3.2 = AGE AT INTERVIEW TWIN2 '95 DEFINITION NOTE: Selection yields 2808 data vectors for analysis NOTE: Vectors contain a total of 23593 observations AGE1.1 AGE3.1 AGE4.1 AGE1.2 AGE3.2 AGE4.2 ; Note: Global variable previously defined. Updating AGE1.1 Note: Global variable previously defined. Updating AGE3.1 Note: Global variable previously defined. Updating AGE4.1 Note: Global variable previously defined. Updating AGE1.2 Note: Global variable previously defined. Updating AGE3.2 Note: Global variable previously defined. Updating AGE4.2 NOTE: Definition yields 2061 data vectors for analysis NOTE: Vectors contain a total of 6745 observations MATRICES = GROUP 1 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 AGE4.1 ! AGE FOR TWIN 1 AT TIMES 1 & 2 SPECIFY P AGE1.2 AGE3.2 AGE4.2 ! AGE FOR TWIN 2 AT TIMES 1 & 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 = \STND(C) ; ! ESTIMATE SHARED ENVIRONMENTAL CORRELATIONS (RC) How am I supposed to take the square root of 0.000000000000000E+00 ? Diagonal elements are: 0.000000000000000E+00 4.000000000000000 Matrix is 0.4000D+01 0.3000D+01 0.0000D+00 0.3000D+01 0.6250D+01 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 Matrix not fully standardized V = \STND(E) ; ! 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 How am I supposed to take the square root of 0.000000000000000E+00 ? Diagonal elements are: 0.000000000000000E+00 4.000000000000000 Matrix is 0.4000D+01 0.3000D+01 0.0000D+00 0.3000D+01 0.6250D+01 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 Matrix not fully standardized Summary of VL file data for group 2 AGE4.2 AGE3.2 AGE1.2 AGE4.1 AGE3.1 AGE1.1 NYAGG1.1 Code -6.0000 -5.0000 -4.0000 -3.0000 -2.0000 -1.0000 1.0000 Number 877.0000 877.0000 877.0000 877.0000 877.0000 877.0000 504.0000 Mean 567.2520 435.0872 430.5477 537.2522 430.6173 431.6680 7.6121 Variance ************************************************************ 25.6824 Minimum 12.5700 14.2700 13.2500 12.5700 14.3100 13.4000 0.0000 Maximum 999.0000 999.0000 999.0000 999.0000 999.0000 999.0000 27.4400 NYAGG3.1 NYAGG4.1 NYAGG1.2 NYAGG3.2 NYAGG4.2 Code 2.0000 3.0000 4.0000 5.0000 6.0000 Number 501.0000 413.0000 503.0000 493.0000 386.0000 Mean 6.2330 5.4632 7.4514 6.4138 5.3900 Variance 18.3844 16.4659 20.3358 21.4903 16.7190 Minimum 0.0000 0.0000 0.0000 0.0000 0.0000 Maximum 27.4400 22.1700 25.3300 33.7800 20.0600 Summary of VL file data for group 3 AGE4.2 AGE3.2 AGE1.2 AGE4.1 AGE3.1 AGE1.1 NYAGG1.1 Code -6.0000 -5.0000 -4.0000 -3.0000 -2.0000 -1.0000 1.0000 Number 2061.0000 2061.0000 2061.0000 2061.0000 2061.0000 2061.0000 1125.0000 Mean 490.2664 437.9138 460.0392 466.2491 435.0684 459.5674 7.6462 Variance ************************************************************ 19.4824 Minimum 11.5000 14.2900 12.6000 11.6000 14.2900 12.6000 0.0000 Maximum 999.0000 999.0000 999.0000 999.0000 999.0000 999.0000 32.7200 NYAGG3.1 NYAGG4.1 NYAGG1.2 NYAGG3.2 NYAGG4.2 Code 2.0000 3.0000 4.0000 5.0000 6.0000 Number 1164.0000 1115.0000 1126.0000 1145.0000 1070.0000 Mean 6.6732 5.5484 7.7991 6.4779 5.3727 Variance 16.7076 14.0051 19.7052 17.5274 14.9030 Minimum 0.0000 0.0000 0.0000 0.0000 0.0000 Maximum 27.4400 21.1100 29.5600 26.3900 22.1700 PARAMETER SPECIFICATIONS GROUP NUMBER: 1 G1: calculation group MATRIX A This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX B This is a FULL matrix of order 1 by 1 1 1 22 MATRIX C This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX M This is a FULL matrix of order 1 by 3 1 2 3 1 19 20 21 MATRIX X This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 1 2 2 3 3 4 5 6 MATRIX Y This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 7 2 8 9 3 10 11 12 MATRIX Z This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 13 2 14 15 3 16 17 18 GROUP NUMBER: 2 G2: MZ twins, datagroup MATRIX A This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX B This is a FULL matrix of order 1 by 1 1 1 22 MATRIX C This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX M This is a FULL matrix of order 1 by 3 1 2 3 1 19 20 21 MATRIX O This is a FULL matrix of order 1 by 3 1 2 3 1 -1 -2 -3 MATRIX P This is a FULL matrix of order 1 by 3 1 2 3 1 -4 -5 -6 MATRIX X This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 1 2 2 3 3 4 5 6 MATRIX Y This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 7 2 8 9 3 10 11 12 MATRIX Z This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 13 2 14 15 3 16 17 18 GROUP NUMBER: 3 G3: DZ twins, datagroup MATRIX A This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX B This is a FULL matrix of order 1 by 1 1 1 22 MATRIX C This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX M This is a FULL matrix of order 1 by 3 1 2 3 1 19 20 21 MATRIX O This is a FULL matrix of order 1 by 3 1 2 3 1 -1 -2 -3 MATRIX P This is a FULL matrix of order 1 by 3 1 2 3 1 -4 -5 -6 MATRIX X This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 1 2 2 3 3 4 5 6 MATRIX Y This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 7 2 8 9 3 10 11 12 MATRIX Z This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 13 2 14 15 3 16 17 18 GROUP NUMBER: 4 G4: Calculate A, C, and E latent factor correlations MATRIX A This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX I This is an IDENTITY matrix of order 3 by 3 MATRIX P This is a computed FULL matrix of order 6 by 6 It has no free parameters specified MATRIX S This is a computed FULL matrix of order 3 by 9 It has no free parameters specified MATRIX T This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX U This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX V This is a computed FULL matrix of order 3 by 3 It has no free parameters specified Mx starting optimization; number of parameters = 22 *** 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 3 by 3 [=X*X'] 1 2 3 1 1.3735 0.2609 1.7843 2 0.2609 0.0963 0.4990 3 1.7843 0.4990 2.8658 MATRIX B This is a FULL matrix of order 1 by 1 1 1 -0.1132 MATRIX C This is a computed FULL matrix of order 3 by 3 [=Y*Y'] 1 2 3 1 5.7383 4.2931 2.5900 2 4.2931 5.0204 3.4508 3 2.5900 3.4508 2.4349 MATRIX E This is a computed FULL matrix of order 3 by 3 [=Z*Z'] 1 2 3 1 13.9799 5.7241 3.9215 2 5.7241 12.5660 5.9709 3 3.9215 5.9709 9.4127 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX M This is a FULL matrix of order 1 by 3 1 2 3 1 9.8322 8.8438 8.1585 MATRIX X This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 1.1720 2 0.2227 0.2162 3 1.5225 0.7402 -0.0003 MATRIX Y This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 2.3955 2 1.7922 1.3448 3 1.0812 1.1251 -0.0002 MATRIX Z This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 3.7390 2 1.5309 3.1972 3 1.0488 1.3653 2.5394 GROUP NUMBER: 2 G2: MZ twins, datagroup MATRIX A This is a computed FULL matrix of order 3 by 3 [=X*X'] 1 2 3 1 1.3735 0.2609 1.7843 2 0.2609 0.0963 0.4990 3 1.7843 0.4990 2.8658 MATRIX B This is a FULL matrix of order 1 by 1 1 1 -0.1132 MATRIX C This is a computed FULL matrix of order 3 by 3 [=Y*Y'] 1 2 3 1 5.7383 4.2931 2.5900 2 4.2931 5.0204 3.4508 3 2.5900 3.4508 2.4349 MATRIX E This is a computed FULL matrix of order 3 by 3 [=Z*Z'] 1 2 3 1 13.9799 5.7241 3.9215 2 5.7241 12.5660 5.9709 3 3.9215 5.9709 9.4127 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX M This is a FULL matrix of order 1 by 3 1 2 3 1 9.8322 8.8438 8.1585 MATRIX O This is a FULL matrix of order 1 by 3 1 2 3 1 16.2300 999.0000 999.0000 MATRIX P This is a FULL matrix of order 1 by 3 1 2 3 1 16.2300 999.0000 999.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 1.1720 2 0.2227 0.2162 3 1.5225 0.7402 -0.0003 MATRIX Y This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 2.3955 2 1.7922 1.3448 3 1.0812 1.1251 -0.0002 MATRIX Z This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 3.7390 2 1.5309 3.1972 3 1.0488 1.3653 2.5394 Vector of OBSERVED means NYAGG1.1 NYAGG3.1 NYAGG4.1 NYAGG1.2 NYAGG3.2 NYAGG4.2 Mean 7.6121 6.2330 5.4632 7.4514 6.4138 5.3900 Vector of EXPECTED means NYAGG1.1 NYAGG3.1 NYAGG4.1 NYAGG1.2 NYAGG3.2 NYAGG4.2 Mean 7.9944 -104.2780 -104.9633 7.9944 -104.2780 -104.9633 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX NYAGG1.1 NYAGG3.1 NYAGG4.1 NYAGG1.2 NYAGG3.2 NYAGG4.2 NYAGG1.1 21.0917 NYAGG3.1 10.2782 17.6828 NYAGG4.1 8.2958 9.9207 14.7133 NYAGG1.2 7.1118 4.5540 4.3743 21.0917 NYAGG3.2 4.5540 5.1167 3.9498 10.2782 17.6828 NYAGG4.2 4.3743 3.9498 5.3006 8.2958 9.9207 14.7133 Function value of this group: 15783.4696 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 3 by 3 [=X*X'] 1 2 3 1 1.3735 0.2609 1.7843 2 0.2609 0.0963 0.4990 3 1.7843 0.4990 2.8658 MATRIX B This is a FULL matrix of order 1 by 1 1 1 -0.1132 MATRIX C This is a computed FULL matrix of order 3 by 3 [=Y*Y'] 1 2 3 1 5.7383 4.2931 2.5900 2 4.2931 5.0204 3.4508 3 2.5900 3.4508 2.4349 MATRIX E This is a computed FULL matrix of order 3 by 3 [=Z*Z'] 1 2 3 1 13.9799 5.7241 3.9215 2 5.7241 12.5660 5.9709 3 3.9215 5.9709 9.4127 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX M This is a FULL matrix of order 1 by 3 1 2 3 1 9.8322 8.8438 8.1585 MATRIX O This is a FULL matrix of order 1 by 3 1 2 3 1 18.3900 999.0000 999.0000 MATRIX P This is a FULL matrix of order 1 by 3 1 2 3 1 18.3900 999.0000 999.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 1.1720 2 0.2227 0.2162 3 1.5225 0.7402 -0.0003 MATRIX Y This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 2.3955 2 1.7922 1.3448 3 1.0812 1.1251 -0.0002 MATRIX Z This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 3.7390 2 1.5309 3.1972 3 1.0488 1.3653 2.5394 Vector of OBSERVED means NYAGG1.1 NYAGG3.1 NYAGG4.1 NYAGG1.2 NYAGG3.2 NYAGG4.2 Mean 7.6462 6.6732 5.5484 7.7991 6.4779 5.3727 Vector of EXPECTED means NYAGG1.1 NYAGG3.1 NYAGG4.1 NYAGG1.2 NYAGG3.2 NYAGG4.2 Mean 7.7498 -104.2780 -104.9633 7.7498 -104.2780 -104.9633 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX NYAGG1.1 NYAGG3.1 NYAGG4.1 NYAGG1.2 NYAGG3.2 NYAGG4.2 NYAGG1.1 21.0917 NYAGG3.1 10.2782 17.6828 NYAGG4.1 8.2958 9.9207 14.7133 NYAGG1.2 6.4250 4.4236 3.4822 21.0917 NYAGG3.2 4.4236 5.0686 3.7003 10.2782 17.6828 NYAGG4.2 3.4822 3.7003 3.8677 8.2958 9.9207 14.7133 Function value of this group: 36822.4247 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 3 by 3 [=X*X'] 1 2 3 1 1.3735 0.2609 1.7843 2 0.2609 0.0963 0.4990 3 1.7843 0.4990 2.8658 MATRIX C This is a computed FULL matrix of order 3 by 3 [=Y*Y'] 1 2 3 1 5.7383 4.2931 2.5900 2 4.2931 5.0204 3.4508 3 2.5900 3.4508 2.4349 MATRIX E This is a computed FULL matrix of order 3 by 3 [=Z*Z'] 1 2 3 1 13.9799 5.7241 3.9215 2 5.7241 12.5660 5.9709 3 3.9215 5.9709 9.4127 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX I This is an IDENTITY matrix of order 3 by 3 MATRIX P This is a computed FULL matrix of order 6 by 6 [=\STND(A+C+E|H@A+C_H@A+C|A+C+E)] 1 2 3 4 5 6 1 1.0000 0.5322 0.4709 0.3046 0.2291 0.1977 2 0.5322 1.0000 0.6151 0.2291 0.2866 0.2294 3 0.4709 0.6151 1.0000 0.1977 0.2294 0.2629 4 0.3046 0.2291 0.1977 1.0000 0.5322 0.4709 5 0.2291 0.2866 0.2294 0.5322 1.0000 0.6151 6 0.1977 0.2294 0.2629 0.4709 0.6151 1.0000 MATRIX S This is a computed FULL matrix of order 3 by 9 [=A%(A+C+E)|C%(A+C+E)|E%(A+C+E)] 1 2 3 4 5 6 7 1 0.0651 0.0254 0.2151 0.2721 0.4177 0.3122 0.6628 2 0.0254 0.0054 0.0503 0.4177 0.2839 0.3478 0.5569 3 0.2151 0.0503 0.1948 0.3122 0.3478 0.1655 0.4727 8 9 1 0.5569 0.4727 2 0.7106 0.6019 3 0.6019 0.6397 MATRIX T This is a computed FULL matrix of order 3 by 3 [=\STND(A)] 1 2 3 1 1.0000 0.7174 0.8994 2 0.7174 1.0000 0.9498 3 0.8994 0.9498 1.0000 MATRIX U This is a computed FULL matrix of order 3 by 3 [=\STND(C)] 1 2 3 1 1.0000 0.7999 0.6929 2 0.7999 1.0000 0.9870 3 0.6929 0.9870 1.0000 MATRIX V This is a computed FULL matrix of order 3 by 3 [=\STND(E)] 1 2 3 1 1.0000 0.4319 0.3419 2 0.4319 1.0000 0.5490 3 0.3419 0.5490 1.0000 *** WARNING! *** Minimization may not be successful. See above CODE GREEN - it probably was OK Your model has 22 estimated parameters and 9545 Observed statistics -2 times log-likelihood of data >>> 52605.894 Degrees of freedom >>>>>>>>>>>>>>>> 9523 Akaike's Information Criterion >>>> 33559.894 Bayesian Information Criterion >>>>-11719.937 Sample size Adjusted BIC >>>> 3409.126 Deviance Information Criterion >>>> -2968.885 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: 2:10.63 TOTAL 0: 0: 2:11.98 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 AGE4.2 AGE3.2 AGE1.2 AGE4.1 AGE3.1 AGE1.1 NYAGG1.1 Code -6.0000 -5.0000 -4.0000 -3.0000 -2.0000 -1.0000 1.0000 Number 877.0000 877.0000 877.0000 877.0000 877.0000 877.0000 504.0000 Mean 567.2520 435.0872 430.5477 537.2522 430.6173 431.6680 7.6121 Variance ************************************************************ 25.6824 Minimum 12.5700 14.2700 13.2500 12.5700 14.3100 13.4000 0.0000 Maximum 999.0000 999.0000 999.0000 999.0000 999.0000 999.0000 27.4400 NYAGG3.1 NYAGG4.1 NYAGG1.2 NYAGG3.2 NYAGG4.2 Code 2.0000 3.0000 4.0000 5.0000 6.0000 Number 501.0000 413.0000 503.0000 493.0000 386.0000 Mean 6.2330 5.4632 7.4514 6.4138 5.3900 Variance 18.3844 16.4659 20.3358 21.4903 16.7190 Minimum 0.0000 0.0000 0.0000 0.0000 0.0000 Maximum 27.4400 22.1700 25.3300 33.7800 20.0600 Summary of VL file data for group 3 AGE4.2 AGE3.2 AGE1.2 AGE4.1 AGE3.1 AGE1.1 NYAGG1.1 Code -6.0000 -5.0000 -4.0000 -3.0000 -2.0000 -1.0000 1.0000 Number 2061.0000 2061.0000 2061.0000 2061.0000 2061.0000 2061.0000 1125.0000 Mean 490.2664 437.9138 460.0392 466.2491 435.0684 459.5674 7.6462 Variance ************************************************************ 19.4824 Minimum 11.5000 14.2900 12.6000 11.6000 14.2900 12.6000 0.0000 Maximum 999.0000 999.0000 999.0000 999.0000 999.0000 999.0000 32.7200 NYAGG3.1 NYAGG4.1 NYAGG1.2 NYAGG3.2 NYAGG4.2 Code 2.0000 3.0000 4.0000 5.0000 6.0000 Number 1164.0000 1115.0000 1126.0000 1145.0000 1070.0000 Mean 6.6732 5.5484 7.7991 6.4779 5.3727 Variance 16.7076 14.0051 19.7052 17.5274 14.9030 Minimum 0.0000 0.0000 0.0000 0.0000 0.0000 Maximum 27.4400 21.1100 29.5600 26.3900 22.1700 PARAMETER SPECIFICATIONS GROUP NUMBER: 1 G1: calculation group MATRIX A This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX B This is a FULL matrix of order 1 by 1 1 1 22 MATRIX C This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX M This is a FULL matrix of order 1 by 3 1 2 3 1 19 20 21 MATRIX X This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 1 2 2 3 3 4 5 6 MATRIX Y This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 7 2 8 9 3 10 11 12 MATRIX Z This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 13 2 14 15 3 16 17 18 GROUP NUMBER: 2 G2: MZ twins, datagroup MATRIX A This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX B This is a FULL matrix of order 1 by 1 1 1 22 MATRIX C This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX M This is a FULL matrix of order 1 by 3 1 2 3 1 19 20 21 MATRIX O This is a FULL matrix of order 1 by 3 1 2 3 1 -1 -2 -3 MATRIX P This is a FULL matrix of order 1 by 3 1 2 3 1 -4 -5 -6 MATRIX X This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 1 2 2 3 3 4 5 6 MATRIX Y This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 7 2 8 9 3 10 11 12 MATRIX Z This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 13 2 14 15 3 16 17 18 GROUP NUMBER: 3 G3: DZ twins, datagroup MATRIX A This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX B This is a FULL matrix of order 1 by 1 1 1 22 MATRIX C This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX M This is a FULL matrix of order 1 by 3 1 2 3 1 19 20 21 MATRIX O This is a FULL matrix of order 1 by 3 1 2 3 1 -1 -2 -3 MATRIX P This is a FULL matrix of order 1 by 3 1 2 3 1 -4 -5 -6 MATRIX X This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 1 2 2 3 3 4 5 6 MATRIX Y This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 7 2 8 9 3 10 11 12 MATRIX Z This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 13 2 14 15 3 16 17 18 GROUP NUMBER: 4 G4: Calculate A, C, and E latent factor correlations MATRIX A This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX I This is an IDENTITY matrix of order 3 by 3 MATRIX P This is a computed FULL matrix of order 6 by 6 It has no free parameters specified MATRIX S This is a computed FULL matrix of order 3 by 9 It has no free parameters specified MATRIX T This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX U This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX V This is a computed FULL matrix of order 3 by 3 It has no free parameters specified Mx starting optimization; number of parameters = 22 MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1: calculation group MATRIX A This is a computed FULL matrix of order 3 by 3 [=X*X'] 1 2 3 1 1.3735 0.2609 1.7843 2 0.2609 0.0963 0.4990 3 1.7843 0.4990 2.8659 MATRIX B This is a FULL matrix of order 1 by 1 1 1 -0.1132 MATRIX C This is a computed FULL matrix of order 3 by 3 [=Y*Y'] 1 2 3 1 5.7383 4.2931 2.5900 2 4.2931 5.0204 3.4508 3 2.5900 3.4508 2.4350 MATRIX E This is a computed FULL matrix of order 3 by 3 [=Z*Z'] 1 2 3 1 13.9799 5.7241 3.9215 2 5.7241 12.5660 5.9709 3 3.9215 5.9709 9.4127 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX M This is a FULL matrix of order 1 by 3 1 2 3 1 9.8322 8.8438 8.1585 MATRIX X This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 1.1720 2 0.2227 0.2162 3 1.5225 0.7402 0.0100 MATRIX Y This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 2.3955 2 1.7922 1.3448 3 1.0812 1.1251 0.0100 MATRIX Z This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 3.7390 2 1.5309 3.1972 3 1.0488 1.3653 2.5394 GROUP NUMBER: 2 G2: MZ twins, datagroup MATRIX A This is a computed FULL matrix of order 3 by 3 [=X*X'] 1 2 3 1 1.3735 0.2609 1.7843 2 0.2609 0.0963 0.4990 3 1.7843 0.4990 2.8659 MATRIX B This is a FULL matrix of order 1 by 1 1 1 -0.1132 MATRIX C This is a computed FULL matrix of order 3 by 3 [=Y*Y'] 1 2 3 1 5.7383 4.2931 2.5900 2 4.2931 5.0204 3.4508 3 2.5900 3.4508 2.4350 MATRIX E This is a computed FULL matrix of order 3 by 3 [=Z*Z'] 1 2 3 1 13.9799 5.7241 3.9215 2 5.7241 12.5660 5.9709 3 3.9215 5.9709 9.4127 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX M This is a FULL matrix of order 1 by 3 1 2 3 1 9.8322 8.8438 8.1585 MATRIX O This is a FULL matrix of order 1 by 3 1 2 3 1 16.2300 999.0000 999.0000 MATRIX P This is a FULL matrix of order 1 by 3 1 2 3 1 16.2300 999.0000 999.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 1.1720 2 0.2227 0.2162 3 1.5225 0.7402 0.0100 MATRIX Y This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 2.3955 2 1.7922 1.3448 3 1.0812 1.1251 0.0100 MATRIX Z This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 3.7390 2 1.5309 3.1972 3 1.0488 1.3653 2.5394 Vector of OBSERVED means NYAGG1.1 NYAGG3.1 NYAGG4.1 NYAGG1.2 NYAGG3.2 NYAGG4.2 Mean 7.6121 6.2330 5.4632 7.4514 6.4138 5.3900 Vector of EXPECTED means NYAGG1.1 NYAGG3.1 NYAGG4.1 NYAGG1.2 NYAGG3.2 NYAGG4.2 Mean 7.9944 -104.2779 -104.9631 7.9944 -104.2779 -104.9631 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX NYAGG1.1 NYAGG3.1 NYAGG4.1 NYAGG1.2 NYAGG3.2 NYAGG4.2 NYAGG1.1 21.0917 NYAGG3.1 10.2782 17.6828 NYAGG4.1 8.2958 9.9207 14.7135 NYAGG1.2 7.1118 4.5540 4.3743 21.0917 NYAGG3.2 4.5540 5.1167 3.9498 10.2782 17.6828 NYAGG4.2 4.3743 3.9498 5.3008 8.2958 9.9207 14.7135 Function value of this group: 15783.4679 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 3 by 3 [=X*X'] 1 2 3 1 1.3735 0.2609 1.7843 2 0.2609 0.0963 0.4990 3 1.7843 0.4990 2.8659 MATRIX B This is a FULL matrix of order 1 by 1 1 1 -0.1132 MATRIX C This is a computed FULL matrix of order 3 by 3 [=Y*Y'] 1 2 3 1 5.7383 4.2931 2.5900 2 4.2931 5.0204 3.4508 3 2.5900 3.4508 2.4350 MATRIX E This is a computed FULL matrix of order 3 by 3 [=Z*Z'] 1 2 3 1 13.9799 5.7241 3.9215 2 5.7241 12.5660 5.9709 3 3.9215 5.9709 9.4127 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX M This is a FULL matrix of order 1 by 3 1 2 3 1 9.8322 8.8438 8.1585 MATRIX O This is a FULL matrix of order 1 by 3 1 2 3 1 18.3900 999.0000 999.0000 MATRIX P This is a FULL matrix of order 1 by 3 1 2 3 1 18.3900 999.0000 999.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 1.1720 2 0.2227 0.2162 3 1.5225 0.7402 0.0100 MATRIX Y This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 2.3955 2 1.7922 1.3448 3 1.0812 1.1251 0.0100 MATRIX Z This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 3.7390 2 1.5309 3.1972 3 1.0488 1.3653 2.5394 Vector of OBSERVED means NYAGG1.1 NYAGG3.1 NYAGG4.1 NYAGG1.2 NYAGG3.2 NYAGG4.2 Mean 7.6462 6.6732 5.5484 7.7991 6.4779 5.3727 Vector of EXPECTED means NYAGG1.1 NYAGG3.1 NYAGG4.1 NYAGG1.2 NYAGG3.2 NYAGG4.2 Mean 7.7498 -104.2779 -104.9631 7.7498 -104.2779 -104.9631 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX NYAGG1.1 NYAGG3.1 NYAGG4.1 NYAGG1.2 NYAGG3.2 NYAGG4.2 NYAGG1.1 21.0917 NYAGG3.1 10.2782 17.6828 NYAGG4.1 8.2958 9.9207 14.7135 NYAGG1.2 6.4250 4.4236 3.4822 21.0917 NYAGG3.2 4.4236 5.0686 3.7003 10.2782 17.6828 NYAGG4.2 3.4822 3.7003 3.8679 8.2958 9.9207 14.7135 Function value of this group: 36822.4264 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 3 by 3 [=X*X'] 1 2 3 1 1.3735 0.2609 1.7843 2 0.2609 0.0963 0.4990 3 1.7843 0.4990 2.8659 MATRIX C This is a computed FULL matrix of order 3 by 3 [=Y*Y'] 1 2 3 1 5.7383 4.2931 2.5900 2 4.2931 5.0204 3.4508 3 2.5900 3.4508 2.4350 MATRIX E This is a computed FULL matrix of order 3 by 3 [=Z*Z'] 1 2 3 1 13.9799 5.7241 3.9215 2 5.7241 12.5660 5.9709 3 3.9215 5.9709 9.4127 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX I This is an IDENTITY matrix of order 3 by 3 MATRIX P This is a computed FULL matrix of order 6 by 6 [=\STND(A+C+E|H@A+C_H@A+C|A+C+E)] 1 2 3 4 5 6 1 1.0000 0.5322 0.4709 0.3046 0.2291 0.1977 2 0.5322 1.0000 0.6150 0.2291 0.2866 0.2294 3 0.4709 0.6150 1.0000 0.1977 0.2294 0.2629 4 0.3046 0.2291 0.1977 1.0000 0.5322 0.4709 5 0.2291 0.2866 0.2294 0.5322 1.0000 0.6150 6 0.1977 0.2294 0.2629 0.4709 0.6150 1.0000 MATRIX S This is a computed FULL matrix of order 3 by 9 [=A%(A+C+E)|C%(A+C+E)|E%(A+C+E)] 1 2 3 4 5 6 7 1 0.0651 0.0254 0.2151 0.2721 0.4177 0.3122 0.6628 2 0.0254 0.0054 0.0503 0.4177 0.2839 0.3478 0.5569 3 0.2151 0.0503 0.1948 0.3122 0.3478 0.1655 0.4727 8 9 1 0.5569 0.4727 2 0.7106 0.6019 3 0.6019 0.6397 MATRIX T This is a computed FULL matrix of order 3 by 3 [=\STND(A)] 1 2 3 1 1.0000 0.7174 0.8993 2 0.7174 1.0000 0.9498 3 0.8993 0.9498 1.0000 MATRIX U This is a computed FULL matrix of order 3 by 3 [=\STND(C)] 1 2 3 1 1.0000 0.7999 0.6929 2 0.7999 1.0000 0.9870 3 0.6929 0.9870 1.0000 MATRIX V This is a computed FULL matrix of order 3 by 3 [=\STND(E)] 1 2 3 1 1.0000 0.4319 0.3419 2 0.4319 1.0000 0.5490 3 0.3419 0.5490 1.0000 Your model has 22 estimated parameters and 9545 Observed statistics -2 times log-likelihood of data >>> 52605.894 Degrees of freedom >>>>>>>>>>>>>>>> 9523 Akaike's Information Criterion >>>> 33559.894 Bayesian Information Criterion >>>>-11719.937 Sample size Adjusted BIC >>>> 3409.126 Deviance Information Criterion >>>> -2968.885 Saturated model fit* >>>>>>>>>>> 52605.894 Saturated model df* >>>>>>>>>>> 9523 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: 4.17 TOTAL 0: 0: 0: 4.30 Total number of warnings issued: 2 ______________________________________________________________________________ ! ***HOW MANY PARAMETERS DOES THIS MODEL ESTIMATE?*** ! ***RUN MODEL*** !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 NVAR NVAR !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 ! ***RUN MODEL***