**MX-PC 1.66b** Job started on 08/11/08 at 09:43:33 ! SCRIPT NAME : RAWVC4A.MX (DP) ! GOAL : TO EVALUATE BEST MODEL FOR VARIANCE COMPONENTS ! DATA : CONTINUOUS ! INPUT : RAW DATA ! UNI/BI/MULTI : TRIVARIATE ! DATA-GROUPS : MZ DZ ! MEANS MODEL : GRAND MEAN ! VARIANCE COVARIANCE MODEL(S) : 1.ACE INCLUDING RG, RC, RE ! ! DOWNLOADING MX SOFTWARE: HTTP://VIEWS.VCU.EDU/MX ! MX SCRIPT'S LIBRARY: HTTP://WWW.PSY.VU.NL/MXBIB !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! #DEFINE NVAR 3 ! NUMBER OF DEPENDENT VARIABLES PER INDIVIDUAL The following MX script lines were read for group 1 G1: CALCULATION GROUP DATA CALC NGROUPS=4 BEGIN MATRICES; X LOWER NVAR NVAR FREE ! ADDITIVE GENETIC PATH COEFFICIENT (STANDARD DEVIATION) 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 G FULL 1 NVAR FREE ! GRAND MEAN PHENOTYPES END MATRICES; MATRIX H .5 ! MATRIX H IS FIXED AND CONTAINS .5 ST 10 X 1 1 1 ! PROVIDE STARTING VALUES ST 45 X 1 2 2 ST 10 X 1 3 3 ST 1.5 Y 1 1 1 ST 0 Y 1 2 2 ST 0 Y 1 3 3 ST 7 Z 1 1 1 ST 20 Z 1 2 2 ST 18 Z 1 3 3 ST 60 G 1 1 1 ST 700 G 1 1 2 ST 380 G 1 1 3 BEGIN ALGEBRA; A= X*X' ; ! ADDITIVE GENETIC VARIANCE C= Y*Y' ; ! COMMON ENVIRONMENTAL VARIANCE E= Z*Z' ; ! UNIQUE ENVIRONMENTAL VARIANCE V = A%(A+C+E) | C%(A+C+E) | E%(A+C+E) ; !STANDARDIZED VARIANCE COMPONENTS END ALGEBRA; LABELS ROW V1 WM BBGM BBWM LABELS COLUMN V1 A1 A2 A3 C1 C2 C3 E1 E2 E3 END The following MX script lines were read for group 2 G2: MZ TWINS, DATAGROUP DA NINPUT_VARS=8 ! NR OF INPUTVARS PER FAMILY MISSING=-200.00 ! MISSING VALUES = -200.00 NOTE: Missing value *string* set to '-200.00' RECTANGULAR FILE=MRI.DAT ! READ IN RAW DATAFILE, MUST BE IN SAME DIRECTORY Rectangular continuous data read initiated Note: Maximum ordinal/rectangular record length is: 1000 Note: It be increased by maxrec= parameter on the data line. NOTE: Rectangular file contained 322 records with data that contained a total of 1588 observations LABELS FAMNR ZYG WM1 BBGMT1 BBWMT1 WM2 BBGMT2 BBWMT2 SELECT IF ZYG =1 ; NOTE: Select if yields 68 data vectors for analysis NOTE: Vectors contain a total of 366 observations SELECT WM1 BBGMT1 BBWMT1 WM2 BBGMT2 BBWMT2 ; ! SELECT ONE PHENOTYPE PER TWIN MEMBER, AND 2 COVARIATES PER TWIN MEMBER MATRICES = GROUP 1 NOTE: Selection yields 67 data vectors for analysis NOTE: Vectors contain a total of 230 observations END MATRICES; MEANS G| G; ! MODEL FOR MEANS, ASSUMING GRAND MEAN TWIN 1 = GRAND MEAN TWIN 2 COVARIANCES ! MODEL FOR MZ VARIANCE/COVARIANCES A+C+E | A+C _ A+C | A+C+E ; OPTIONS RSIDUALS END The following MX script lines were read for group 3 G3: DZ TWINS, DATAGROUP DA NINPUT_VARS=8 ! NR OF INPUTVARS PER FAMILY MISSING=-200.00 ! MISSING VALUES = -200.00 NOTE: Missing value *string* set to '-200.00' RECTANGULAR FILE=MRI.DAT ! READ IN RAW DATAFILE, MUST BE IN SAME DIRECTORY Rectangular continuous data read initiated Note: Maximum ordinal/rectangular record length is: 1000 Note: It be increased by maxrec= parameter on the data line. NOTE: Rectangular file contained 322 records with data that contained a total of 1588 observations LABELS FAMNR ZYG WM1 BBGMT1 BBWMT1 WM2 BBGMT2 BBWMT2 SELECT IF ZYG =3 ; NOTE: Select if yields 43 data vectors for analysis NOTE: Vectors contain a total of 218 observations SELECT WM1 BBGMT1 BBWMT1 WM2 BBGMT2 BBWMT2 ; ! SELECT ONE PHENOTYPE PER TWIN MEMBER, AND 2 COVARIATES PER TWIN MEMBER MATRICES = GROUP 1 NOTE: Selection yields 43 data vectors for analysis NOTE: Vectors contain a total of 132 observations END MATRICES; MEANS G| G; ! MODEL FOR MEANS, ASSUMING GRAND MEAN TWIN 1 = GRAND MEAN TWIN 2 COVARIANCES ! MODEL FOR DZ VARIANCE/COVARIANCES A+C+E | H@A+C _ H@A+C | A+C+E ; OPTIONS RSIDUALS END The following MX script lines were read for group 4 G4: SUMMARIZE RELEVANT STATISTICS DATA CALC BEGIN MATRICES ; A COMPUTED NVAR NVAR = A1 C COMPUTED NVAR NVAR = C1 E COMPUTED NVAR NVAR = E1 V COMPUTED NVAR 9 = V1 END MATRICES ; BEGIN ALGEBRA ; K = \STND(A) ; ! CALCULATES RG (OFF DIAGONAL) L = \STND(C) ; ! CALCULATES RC (OFF DIAGONAL) How am I supposed to take the square root of 0.000000000000000E+00 ? Diagonal elements are: 0.000000000000000E+00 2.250000000000000 Matrix is 0.2250D+01 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 Matrix not fully standardized M = \STND(E) ; ! CALCULATES RE (OFF DIAGONAL) END ALGEBRA ; ! INTERVAL @95 V 1 1 1 V 1 1 3 V 1 1 5 ! REQUEST CI'S ON %A, %C, %E FOR FIRST PHENOTYPE ! INTERVAL @95 V 1 2 2 V 1 2 4 V 1 2 6 ! REQUEST CI'S ON %A, %C, %E FOR SECOND PHENOTYPE ! INTERVAL @95 K 4 2 1 L 4 2 1 M 4 2 1 ! REQUEST CI'S ON RG, RC, RE END How am I supposed to take the square root of 0.000000000000000E+00 ? Diagonal elements are: 0.000000000000000E+00 2.250000000000000 Matrix is 0.2250D+01 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 Matrix not fully standardized Summary of VL file data for group 2 WM1 BBGMT1 BBWMT1 WM2 BBGMT2 BBWMT2 Code 1.0000 2.0000 3.0000 4.0000 5.0000 6.0000 Number 48.0000 33.0000 33.0000 50.0000 33.0000 33.0000 Mean 64.0554 705.2912 393.2188 63.9522 699.5758 389.2467 Variance 153.0001 3076.3244 2689.1675 236.4818 3508.8562 2789.0215 Minimum 34.8900 570.6100 221.6100 27.7400 523.6600 207.2600 Maximum 87.1900 860.5200 479.7300 91.7700 836.6400 516.8400 Summary of VL file data for group 3 WM1 BBGMT1 BBWMT1 WM2 BBGMT2 BBWMT2 Code 1.0000 2.0000 3.0000 4.0000 5.0000 6.0000 Number 33.0000 17.0000 17.0000 31.0000 17.0000 17.0000 Mean 60.0403 690.4318 402.3047 62.1271 655.9288 373.9806 Variance 157.8822 1440.8459 1176.7411 96.4641 1726.0283 1863.3563 Minimum 24.9400 603.9300 347.5700 39.9700 603.7300 260.1600 Maximum 82.5000 750.2500 489.4900 85.4800 771.7100 426.3600 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 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 G This is a FULL matrix of order 1 by 3 1 2 3 1 19 20 21 MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX V This is a computed FULL matrix of order 3 by 9 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 1 2 2 3 3 4 5 6 MATRIX Y This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 7 2 8 9 3 10 11 12 MATRIX Z This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 13 2 14 15 3 16 17 18 GROUP NUMBER: 2 G2: MZ twins, datagroup 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 G This is a FULL matrix of order 1 by 3 1 2 3 1 19 20 21 MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX V This is a computed FULL matrix of order 3 by 9 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 1 2 2 3 3 4 5 6 MATRIX Y This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 7 2 8 9 3 10 11 12 MATRIX Z This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 13 2 14 15 3 16 17 18 GROUP NUMBER: 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 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 G This is a FULL matrix of order 1 by 3 1 2 3 1 19 20 21 MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX V This is a computed FULL matrix of order 3 by 9 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 1 2 2 3 3 4 5 6 MATRIX Y This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 7 2 8 9 3 10 11 12 MATRIX Z This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 13 2 14 15 3 16 17 18 GROUP NUMBER: 4 G4: summarize relevant statistics 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 K This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX L This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX M 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 9 It has no free parameters specified Mx starting optimization; number of parameters = 21 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 109.5421 148.9947 92.7224 2 148.9947 2322.0396 1912.3678 3 92.7224 1912.3678 1737.4462 MATRIX C This is a computed FULL matrix of order 3 by 3 [=Y*Y'] 1 2 3 1 2.3468 -24.1240 27.4001 2 -24.1240 248.2309 -281.8857 3 27.4001 -281.8857 320.1257 MATRIX E This is a computed FULL matrix of order 3 by 3 [=Z*Z'] 1 2 3 1 57.7516 -24.9571 15.5015 2 -24.9571 453.1544 -101.3189 3 15.5015 -101.3189 353.0115 MATRIX G This is a FULL matrix of order 1 by 3 1 2 3 1 62.7153 691.4319 389.4945 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX V This is a computed FULL matrix of order 3 by 9 [=A%(A+C+E)|C%(A+C+E)|E%(A+C+E)] A1 A2 A3 C1 C2 C3 WM 0.6457 1.4912 0.6837 0.0138 -0.2414 0.2020 BBGM 1.4912 0.7680 1.2506 -0.2414 0.0821 -0.1843 BBWM 0.6837 1.2506 0.7208 0.2020 -0.1843 0.1328 E1 E2 E3 WM 0.3404 -0.2498 0.1143 BBGM -0.2498 0.1499 -0.0663 BBWM 0.1143 -0.0663 0.1464 MATRIX X This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 10.4662 2 14.2358 46.0368 3 8.8592 38.8005 12.3887 MATRIX Y This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 1.5319 2 -15.7475 -0.4976 3 17.8860 0.4536 0.0994 MATRIX Z This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 7.5994 2 -3.2841 21.0326 3 2.0398 -4.4987 18.1277 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 109.5421 148.9947 92.7224 2 148.9947 2322.0396 1912.3678 3 92.7224 1912.3678 1737.4462 MATRIX C This is a computed FULL matrix of order 3 by 3 [=Y*Y'] 1 2 3 1 2.3468 -24.1240 27.4001 2 -24.1240 248.2309 -281.8857 3 27.4001 -281.8857 320.1257 MATRIX E This is a computed FULL matrix of order 3 by 3 [=Z*Z'] 1 2 3 1 57.7516 -24.9571 15.5015 2 -24.9571 453.1544 -101.3189 3 15.5015 -101.3189 353.0115 MATRIX G This is a FULL matrix of order 1 by 3 1 2 3 1 62.7153 691.4319 389.4945 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX V This is a computed FULL matrix of order 3 by 9 [=A%(A+C+E)|C%(A+C+E)|E%(A+C+E)] A1 A2 A3 C1 C2 C3 WM 0.6457 1.4912 0.6837 0.0138 -0.2414 0.2020 BBGM 1.4912 0.7680 1.2506 -0.2414 0.0821 -0.1843 BBWM 0.6837 1.2506 0.7208 0.2020 -0.1843 0.1328 E1 E2 E3 WM 0.3404 -0.2498 0.1143 BBGM -0.2498 0.1499 -0.0663 BBWM 0.1143 -0.0663 0.1464 MATRIX X This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 10.4662 2 14.2358 46.0368 3 8.8592 38.8005 12.3887 MATRIX Y This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 1.5319 2 -15.7475 -0.4976 3 17.8860 0.4536 0.0994 MATRIX Z This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 7.5994 2 -3.2841 21.0326 3 2.0398 -4.4987 18.1277 Vector of OBSERVED means WM1 BBGMT1 BBWMT1 WM2 BBGMT2 BBWMT2 Mean 64.0554 705.2912 393.2188 63.9522 699.5758 389.2467 Vector of EXPECTED means WM1 BBGMT1 BBWMT1 WM2 BBGMT2 BBWMT2 Mean 62.7153 691.4319 389.4945 62.7153 691.4319 389.4945 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX WM1 BBGMT1 BBWMT1 WM2 BBGMT2 BBWMT2 WM1 169.6404 BBGMT1 99.9136 3023.4250 BBWMT1 135.6240 1529.1632 2410.5834 WM2 111.8889 124.8708 120.1225 169.6404 BBGMT2 124.8708 2570.2706 1630.4821 99.9136 3023.4250 BBWMT2 120.1225 1630.4821 2057.5719 135.6240 1529.1632 2410.5834 Function value of this group: 2081.2014 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 109.5421 148.9947 92.7224 2 148.9947 2322.0396 1912.3678 3 92.7224 1912.3678 1737.4462 MATRIX C This is a computed FULL matrix of order 3 by 3 [=Y*Y'] 1 2 3 1 2.3468 -24.1240 27.4001 2 -24.1240 248.2309 -281.8857 3 27.4001 -281.8857 320.1257 MATRIX E This is a computed FULL matrix of order 3 by 3 [=Z*Z'] 1 2 3 1 57.7516 -24.9571 15.5015 2 -24.9571 453.1544 -101.3189 3 15.5015 -101.3189 353.0115 MATRIX G This is a FULL matrix of order 1 by 3 1 2 3 1 62.7153 691.4319 389.4945 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX V This is a computed FULL matrix of order 3 by 9 [=A%(A+C+E)|C%(A+C+E)|E%(A+C+E)] A1 A2 A3 C1 C2 C3 WM 0.6457 1.4912 0.6837 0.0138 -0.2414 0.2020 BBGM 1.4912 0.7680 1.2506 -0.2414 0.0821 -0.1843 BBWM 0.6837 1.2506 0.7208 0.2020 -0.1843 0.1328 E1 E2 E3 WM 0.3404 -0.2498 0.1143 BBGM -0.2498 0.1499 -0.0663 BBWM 0.1143 -0.0663 0.1464 MATRIX X This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 10.4662 2 14.2358 46.0368 3 8.8592 38.8005 12.3887 MATRIX Y This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 1.5319 2 -15.7475 -0.4976 3 17.8860 0.4536 0.0994 MATRIX Z This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 7.5994 2 -3.2841 21.0326 3 2.0398 -4.4987 18.1277 Vector of OBSERVED means WM1 BBGMT1 BBWMT1 WM2 BBGMT2 BBWMT2 Mean 60.0403 690.4318 402.3047 62.1271 655.9288 373.9806 Vector of EXPECTED means WM1 BBGMT1 BBWMT1 WM2 BBGMT2 BBWMT2 Mean 62.7153 691.4319 389.4945 62.7153 691.4319 389.4945 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX WM1 BBGMT1 BBWMT1 WM2 BBGMT2 BBWMT2 WM1 169.6404 BBGMT1 99.9136 3023.4250 BBWMT1 135.6240 1529.1632 2410.5834 WM2 57.1178 50.3734 73.7613 169.6404 BBGMT2 50.3734 1409.2508 674.2982 99.9136 3023.4250 BBWMT2 73.7613 674.2982 1188.8488 135.6240 1529.1632 2410.5834 Function value of this group: 1198.1693 Where the fit function is -2 * Log-likelihood of raw data GROUP NUMBER: 4 G4: summarize relevant statistics MATRIX A This is a computed FULL matrix of order 3 by 3 [=X*X'] 1 2 3 1 109.5421 148.9947 92.7224 2 148.9947 2322.0396 1912.3678 3 92.7224 1912.3678 1737.4462 MATRIX C This is a computed FULL matrix of order 3 by 3 [=Y*Y'] 1 2 3 1 2.3468 -24.1240 27.4001 2 -24.1240 248.2309 -281.8857 3 27.4001 -281.8857 320.1257 MATRIX E This is a computed FULL matrix of order 3 by 3 [=Z*Z'] 1 2 3 1 57.7516 -24.9571 15.5015 2 -24.9571 453.1544 -101.3189 3 15.5015 -101.3189 353.0115 MATRIX K This is a computed FULL matrix of order 3 by 3 [=\STND(A)] 1 2 3 1 1.0000 0.2954 0.2125 2 0.2954 1.0000 0.9521 3 0.2125 0.9521 1.0000 MATRIX L This is a computed FULL matrix of order 3 by 3 [=\STND(C)] 1 2 3 1 1.0000 -0.9995 0.9997 2 -0.9995 1.0000 -1.0000 3 0.9997 -1.0000 1.0000 MATRIX M This is a computed FULL matrix of order 3 by 3 [=\STND(E)] 1 2 3 1 1.0000 -0.1543 0.1086 2 -0.1543 1.0000 -0.2533 3 0.1086 -0.2533 1.0000 MATRIX V This is a computed FULL matrix of order 3 by 9 [=A%(A+C+E)|C%(A+C+E)|E%(A+C+E)] A1 A2 A3 C1 C2 C3 WM 0.6457 1.4912 0.6837 0.0138 -0.2414 0.2020 BBGM 1.4912 0.7680 1.2506 -0.2414 0.0821 -0.1843 BBWM 0.6837 1.2506 0.7208 0.2020 -0.1843 0.1328 E1 E2 E3 WM 0.3404 -0.2498 0.1143 BBGM -0.2498 0.1499 -0.0663 BBWM 0.1143 -0.0663 0.1464 Your model has 21 estimated parameters and 362 Observed statistics -2 times log-likelihood of data >>> 3279.371 Degrees of freedom >>>>>>>>>>>>>>>> 341 Akaike's Information Criterion >>>> 2597.371 Bayesian Information Criterion >>>> 838.253 Sample size Adjusted BIC >>>> 1377.039 Deviance Information Criterion >>>> 1151.611 This problem used 0.5% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.22 Execution 0: 0: 0: 0.41 TOTAL 0: 0: 0: 0.63 Total number of warnings issued: 0 ______________________________________________________________________________