** Mx startup successful ** **MX-PC 1.66b** Job started on 08/11/08 at 09:44:03 ! SCRIPT NAME : RAWVC4A.MX (DP) ! GOAL : TO EVALUATE BEST MODEL FOR VARIANCE COMPONENTS ! DATA : CONTINUOUS ! INPUT : RAW DATA ! UNI/BI/MULTI : BIVARIATE ! 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 2 ! 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 5.5 X 1 1 1 ! PROVIDE STARTING VALUES ST 7.5 X 1 2 2 ST 5.5 Y 1 1 1 ST 7.5 Y 1 2 2 ST 5.5 Z 1 1 1 ST 7.5 Z 1 2 2 ST 60 G 1 1 1 ST 700 G 1 1 2 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 LABELS COLUMN V1 A1 A2 C1 C2 E1 E2 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 WM2 BBGMT2 ; ! 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 164 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 WM2 BBGMT2 ; ! 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 98 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 6 = V1 END MATRICES ; BEGIN ALGEBRA ; K = \STND(A) ; ! CALCULATES RG (OFF DIAGONAL) L = \STND(C) ; ! CALCULATES RC (OFF DIAGONAL) 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 Summary of VL file data for group 2 WM1 BBGMT1 WM2 BBGMT2 Code 1.0000 2.0000 3.0000 4.0000 Number 48.0000 33.0000 50.0000 33.0000 Mean 64.0554 705.2912 63.9522 699.5758 Variance 153.0001 3076.3244 236.4818 3508.8562 Minimum 34.8900 570.6100 27.7400 523.6600 Maximum 87.1900 860.5200 91.7700 836.6400 Summary of VL file data for group 3 WM1 BBGMT1 WM2 BBGMT2 Code 1.0000 2.0000 3.0000 4.0000 Number 33.0000 17.0000 31.0000 17.0000 Mean 60.0403 690.4318 62.1271 655.9288 Variance 157.8822 1440.8459 96.4641 1726.0283 Minimum 24.9400 603.9300 39.9700 603.7300 Maximum 82.5000 750.2500 85.4800 771.7100 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 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 G This is a FULL matrix of order 1 by 2 1 2 1 10 11 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 2 by 6 It has no free parameters specified 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 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 G This is a FULL matrix of order 1 by 2 1 2 1 10 11 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 2 by 6 It has no free parameters specified 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 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 G This is a FULL matrix of order 1 by 2 1 2 1 10 11 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 2 by 6 It has no free parameters specified 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: summarize relevant statistics 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 K This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX L This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX V This is a computed FULL matrix of order 2 by 6 It has no free parameters specified Mx starting optimization; number of parameters = 11 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 111.5853 130.4175 2 130.4175 2487.1875 MATRIX C This is a computed FULL matrix of order 2 by 2 [=Y*Y'] 1 2 1 5.6312E-10 7.3848E-08 2 7.3848E-08 1.0280E-05 MATRIX E This is a computed FULL matrix of order 2 by 2 [=Z*Z'] 1 2 1 57.6137 -31.6614 2 -31.6614 470.8741 MATRIX G This is a FULL matrix of order 1 by 2 1 2 1 62.6772 690.9162 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 2 by 6 [=A%(A+C+E)|C%(A+C+E)|E%(A+C+E)] A1 A2 C1 C2 E1 WM 6.5949E-01 1.3206E+00 3.3282E-12 7.4778E-10 3.4051E-01 BBGM 1.3206E+00 8.4082E-01 7.4778E-10 3.4751E-09 -3.2060E-01 E2 WM -3.2060E-01 BBGM 1.5918E-01 MATRIX X This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 10.5634 2 12.3462 48.3193 MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 2.3730E-05 2 3.1120E-03 7.7148E-04 MATRIX Z This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 7.5904 2 -4.1713 21.2949 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 111.5853 130.4175 2 130.4175 2487.1875 MATRIX C This is a computed FULL matrix of order 2 by 2 [=Y*Y'] 1 2 1 5.6312E-10 7.3848E-08 2 7.3848E-08 1.0280E-05 MATRIX E This is a computed FULL matrix of order 2 by 2 [=Z*Z'] 1 2 1 57.6137 -31.6614 2 -31.6614 470.8741 MATRIX G This is a FULL matrix of order 1 by 2 1 2 1 62.6772 690.9162 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 2 by 6 [=A%(A+C+E)|C%(A+C+E)|E%(A+C+E)] A1 A2 C1 C2 E1 WM 6.5949E-01 1.3206E+00 3.3282E-12 7.4778E-10 3.4051E-01 BBGM 1.3206E+00 8.4082E-01 7.4778E-10 3.4751E-09 -3.2060E-01 E2 WM -3.2060E-01 BBGM 1.5918E-01 MATRIX X This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 10.5634 2 12.3462 48.3193 MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 2.3730E-05 2 3.1120E-03 7.7148E-04 MATRIX Z This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 7.5904 2 -4.1713 21.2949 Vector of OBSERVED means WM1 BBGMT1 WM2 BBGMT2 Mean 64.0554 705.2912 63.9522 699.5758 Vector of EXPECTED means WM1 BBGMT1 WM2 BBGMT2 Mean 62.6772 690.9162 62.6772 690.9162 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX WM1 BBGMT1 WM2 BBGMT2 WM1 169.1990 BBGMT1 98.7562 2958.0616 WM2 111.5853 130.4175 169.1990 BBGMT2 130.4175 2487.1875 98.7562 2958.0616 Function value of this group: 1439.2184 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 111.5853 130.4175 2 130.4175 2487.1875 MATRIX C This is a computed FULL matrix of order 2 by 2 [=Y*Y'] 1 2 1 5.6312E-10 7.3848E-08 2 7.3848E-08 1.0280E-05 MATRIX E This is a computed FULL matrix of order 2 by 2 [=Z*Z'] 1 2 1 57.6137 -31.6614 2 -31.6614 470.8741 MATRIX G This is a FULL matrix of order 1 by 2 1 2 1 62.6772 690.9162 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 2 by 6 [=A%(A+C+E)|C%(A+C+E)|E%(A+C+E)] A1 A2 C1 C2 E1 WM 6.5949E-01 1.3206E+00 3.3282E-12 7.4778E-10 3.4051E-01 BBGM 1.3206E+00 8.4082E-01 7.4778E-10 3.4751E-09 -3.2060E-01 E2 WM -3.2060E-01 BBGM 1.5918E-01 MATRIX X This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 10.5634 2 12.3462 48.3193 MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 2.3730E-05 2 3.1120E-03 7.7148E-04 MATRIX Z This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 7.5904 2 -4.1713 21.2949 Vector of OBSERVED means WM1 BBGMT1 WM2 BBGMT2 Mean 60.0403 690.4318 62.1271 655.9288 Vector of EXPECTED means WM1 BBGMT1 WM2 BBGMT2 Mean 62.6772 690.9162 62.6772 690.9162 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX WM1 BBGMT1 WM2 BBGMT2 WM1 169.1990 BBGMT1 98.7562 2958.0616 WM2 55.7927 65.2088 169.1990 BBGMT2 65.2088 1243.5937 98.7562 2958.0616 Function value of this group: 859.3373 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 2 by 2 [=X*X'] 1 2 1 111.5853 130.4175 2 130.4175 2487.1875 MATRIX C This is a computed FULL matrix of order 2 by 2 [=Y*Y'] 1 2 1 5.6312E-10 7.3848E-08 2 7.3848E-08 1.0280E-05 MATRIX E This is a computed FULL matrix of order 2 by 2 [=Z*Z'] 1 2 1 57.6137 -31.6614 2 -31.6614 470.8741 MATRIX K This is a computed FULL matrix of order 2 by 2 [=\STND(A)] 1 2 1 1.0000 0.2476 2 0.2476 1.0000 MATRIX L This is a computed FULL matrix of order 2 by 2 [=\STND(C)] 1 2 1 1.0000 0.9706 2 0.9706 1.0000 MATRIX M This is a computed FULL matrix of order 2 by 2 [=\STND(E)] 1 2 1 1.0000 -0.1922 2 -0.1922 1.0000 MATRIX V This is a computed FULL matrix of order 2 by 6 [=A%(A+C+E)|C%(A+C+E)|E%(A+C+E)] A1 A2 C1 C2 E1 WM 6.5949E-01 1.3206E+00 3.3282E-12 7.4778E-10 3.4051E-01 BBGM 1.3206E+00 8.4082E-01 7.4778E-10 3.4751E-09 -3.2060E-01 E2 WM -3.2060E-01 BBGM 1.5918E-01 Your model has 11 estimated parameters and 262 Observed statistics -2 times log-likelihood of data >>> 2298.556 Degrees of freedom >>>>>>>>>>>>>>>> 251 Akaike's Information Criterion >>>> 1796.556 Bayesian Information Criterion >>>> 559.368 Sample size Adjusted BIC >>>> 955.952 Deviance Information Criterion >>>> 790.021 9 Confidence intervals requested in group 4 Matrix Element Int. Estimate Lower Upper Lfail Ufail V 1 1 1 95.0 0.6595 0.2325 0.7812 0 1 0 0 V 1 1 3 95.0 0.0000 0.0000 0.3839 0 0 0 0 V 1 1 5 95.0 0.3405 0.2188 0.5221 0 0 0 0 V 1 2 2 95.0 0.8408 0.3935 0.9104 0 0 0 0 V 1 2 4 95.0 0.0000 0.0000 0.4387 1 0 0 0 V 1 2 6 95.0 0.1592 0.0896 0.2882 0 0 0 0 K 4 2 1 95.0 0.2476 -0.3492 0.7678 0 1 0 1 L 4 2 1 95.0 0.9706 -1.0000 1.0000 0 0 0 0 M 4 2 1 95.0 -0.1922 -0.5299 0.2185 0 0 0 0 This problem used 0.6% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.20 Execution 0: 0: 0: 1.91 TOTAL 0: 0: 0: 2.11 Total number of warnings issued: 0 ______________________________________________________________________________