**MX-PC 1.63h** Job started on 03/06/06 at 16:16:54 ! ESTIMATE GENETIC AND ENVIRONMENTAL COMPONENTS - ACED MODEL ! OZ BMI DATA - YOUNGER FEMALES The following MX script lines were read for group 1 #NGROUPS 4 Note: #NGroup set number of groups to 4 #DEFINE NVAR 1 #DEFINE NVAR2 2 TITLE G1: MODEL PARAMETERS CALCULATION BEGIN MATRICES; X LOWER NVAR NVAR FREE ! ADDITIVE GENETIC PATH, A Y LOWER NVAR NVAR ! COMMON ENVIRONMENTAL PATH, C Z LOWER NVAR NVAR FREE ! SPECIFIC ENVIRONMENTAL PATH, E W LOWER NVAR NVAR FREE ! DOMINANCE GENETIC PATH, D H FULL 1 1 ! SCALAR, 0.5 Q FULL 1 1 ! SCALAR, 0.25 END MATRICES; LABEL ROW X ADD_GEN LABEL ROW Y COMM_ENV LABEL ROW Z SPEC_ENV LABEL ROW W DOM_GEN MATRIX H .5 MATRIX Q .25 BEGIN ALGEBRA; A= X*X'; ! ADDITIVE GENETIC VARIANCE, A^2 C= Y*Y'; ! COMMON ENVIRONMENTAL VARIANCE, C^2 E= Z*Z'; ! SPECIFIC ENVIRONMENTAL VARIANCE, E^2 D= W*W'; ! DOMINANCE GENETIC VARIANCE, D^2 END ALGEBRA; END GROUP The following MX script lines were read for group 2 TITLE G2: MZ DATA #INCLUDE OZBMI2.DAT Note: Opening #include file 1 ozbmi2.dat DATA NINPUT=13 RECTANGULAR FILE=OZBMI2.REC Rectangular continuous data read initiated Note: Maximum ordinal/rectangular record length is: 1000 NOTE: Rectangular file contained 3808 records with data that contained a total of 48661 observations LABELS FAM AGECAT AGE ZYG PART WT1 WT2 HT1 HT2 HTWT1 HTWT2 BMI1 BMI2 Note: Closing #include file 1 SELECT IF ZYG=1 NOTE: Select if yields 1232 data vectors for analysis NOTE: Vectors contain a total of 15797 observations SELECT IF AGECAT=1 NOTE: Select if yields 569 data vectors for analysis NOTE: Vectors contain a total of 7272 observations SELECT BMI1 BMI2 ; BEGIN MATRICES = GROUP 1; ! COPY ALL MATRICES FROM GROUP 1, BOTH DECLARED AND COMPUTED NOTE: Selection yields 566 data vectors for analysis NOTE: Vectors contain a total of 1100 observations M FULL 1 NVAR2 FREE MEANS M; ! 1X2 VECTOR OF EXPECTED MEANS COVARIANCE ! 2X2 COVARIANCE MATRIX A+C+E+D | A+C+D _ A+C+D | A+C+E+D ; OPTION RSIDUALS; END The following MX script lines were read for group 3 TITLE G3: DZ DATA #INCLUDE OZBMI2.DAT Note: Opening #include file 1 ozbmi2.dat DATA NINPUT=13 RECTANGULAR FILE=OZBMI2.REC Rectangular continuous data read initiated Note: Maximum ordinal/rectangular record length is: 1000 NOTE: Rectangular file contained 3808 records with data that contained a total of 48661 observations LABELS FAM AGECAT AGE ZYG PART WT1 WT2 HT1 HT2 HTWT1 HTWT2 BMI1 BMI2 Note: Closing #include file 1 SELECT IF ZYG=3 NOTE: Select if yields 751 data vectors for analysis NOTE: Vectors contain a total of 9618 observations SELECT IF AGECAT=1 NOTE: Select if yields 351 data vectors for analysis NOTE: Vectors contain a total of 4483 observations SELECT BMI1 BMI2 ; BEGIN MATRICES = GROUP 1; NOTE: Selection yields 349 data vectors for analysis NOTE: Vectors contain a total of 677 observations M FULL 1 NVAR2 FREE END MATRICES; MEANS M; COVARIANCE A+C+E+D | H@A+C+Q@D _ H@A+C+Q@D | A+C+E+D; OPTION RSIDUALS END The following MX script lines were read for group 4 TITLE G4: STANDARDIZATION CALCULATION BEGIN MATRICES = GROUP 1; END MATRICES; START .6 ALL ! START ALL FREE PARAMETERS AT .6 START 20 M 2 1 1 - M 2 1 NVAR2 ! START FREE PARAMERS FOR MEANS AT 20 (OVERWRITE PREVIOUS STATEMENT) START 20 M 3 1 1 - M 3 1 NVAR2 BEGIN ALGEBRA; V=A+C+E+D; ! TOTAL VARIANCE P=A|C|E|D; ! CONCATENATE PARAMETER ESTIMATES S=P@V~; ! STANDARDIZED PARAMETER ESTIMATES END ALGEBRA; !ADE MODEL INTERVAL S 1 1 - S 1 4 ! 95% CONFIDENCE INTERVALS AROUND VARIANCE COMPONENTS OPTION NDECIMALS=4 OPTION SAT=4055.935,1767 ! LIKELIHOOD AND DF FROM SATURATED MODEL: OZBMIYFSAT.MXO END Summary of VL file data for group 2 BMI1 BMI2 Code 1.0000 2.0000 Number 550.0000 550.0000 Mean 21.3533 21.3389 Variance 0.7274 0.7905 Minimum 18.1125 19.3835 Maximum 25.5347 26.5016 Summary of VL file data for group 3 BMI1 BMI2 Code 1.0000 2.0000 Number 339.0000 338.0000 Mean 21.4481 21.4567 Variance 0.7702 0.8222 Minimum 18.9794 19.2262 Maximum 25.3704 26.8383 PARAMETER SPECIFICATIONS GROUP NUMBER: 1 Title G1: Model Parameters MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 1 by 1 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 Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 DOM_GEN 3 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 ADD_GEN 1 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 SPEC_ENV 2 GROUP NUMBER: 2 Title G2: MZ data MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 1 by 1 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 4 5 MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 DOM_GEN 3 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 ADD_GEN 1 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 SPEC_ENV 2 GROUP NUMBER: 3 Title G3: DZ data MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 1 by 1 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 6 7 MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 DOM_GEN 3 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 ADD_GEN 1 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 SPEC_ENV 2 GROUP NUMBER: 4 Title G4: Standardization MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 1 by 1 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 P This is a computed FULL matrix of order 1 by 4 It has no free parameters specified MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX S This is a computed FULL matrix of order 1 by 4 It has no free parameters specified MATRIX V This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 DOM_GEN 3 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 ADD_GEN 1 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 SPEC_ENV 2 Mx starting optimization; number of parameters = 7 MX PARAMETER ESTIMATES GROUP NUMBER: 1 Title G1: Model Parameters MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.3125 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.0000 MATRIX D This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 1 0.2952 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.1693 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 DOM_GEN 0.5433 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 ADD_GEN 0.5590 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 COMM_ENV 0.0000 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 SPEC_ENV -0.4115 GROUP NUMBER: 2 Title G2: MZ data MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.3125 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.0000 MATRIX D This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 1 0.2952 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.1693 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 21.3440 21.3486 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 DOM_GEN 0.5433 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 ADD_GEN 0.5590 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 COMM_ENV 0.0000 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 SPEC_ENV -0.4115 Vector of OBSERVED means BMI1 BMI2 Mean 21.3533 21.3389 Vector of EXPECTED means BMI1 BMI2 Mean 21.3440 21.3486 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX BMI1 BMI2 BMI1 0.7769 BMI2 0.6076 0.7769 Function value of this group: 2324.1103 Where the fit function is -2 * Log-likelihood of raw data GROUP NUMBER: 3 Title G3: DZ data MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.3125 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.0000 MATRIX D This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 1 0.2952 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.1693 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 21.4475 21.4578 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 DOM_GEN 0.5433 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 ADD_GEN 0.5590 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 COMM_ENV 0.0000 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 SPEC_ENV -0.4115 Vector of OBSERVED means BMI1 BMI2 Mean 21.4481 21.4567 Vector of EXPECTED means BMI1 BMI2 Mean 21.4475 21.4578 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX BMI1 BMI2 BMI1 0.7769 BMI2 0.2300 0.7769 Function value of this group: 1735.1002 Where the fit function is -2 * Log-likelihood of raw data GROUP NUMBER: 4 Title G4: Standardization MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.3125 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.0000 MATRIX D This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 1 0.2952 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.1693 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX P This is a computed FULL matrix of order 1 by 4 [=A|C|E|D] 1 2 3 4 1 0.3125 0.0000 0.1693 0.2952 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX S This is a computed FULL matrix of order 1 by 4 [=P@V~] 1 2 3 4 1 0.4022 0.0000 0.2179 0.3799 MATRIX V This is a computed FULL matrix of order 1 by 1 [=A+C+E+D] 1 1 0.7769 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 DOM_GEN 0.5433 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 ADD_GEN 0.5590 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 COMM_ENV 0.0000 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 SPEC_ENV -0.4115 Your model has 7 estimated parameters and 1777 Observed statistics -2 times log-likelihood of data >>> 4059.211 Degrees of freedom >>>>>>>>>>>>>>>> 1770 Akaike's Information Criterion >>>> 519.211 Bayesian Information Criterion >>>> -4005.143 Sample size Adjusted BIC >>>> -1194.497 Deviance Information Criterion >>>> -2378.621 Saturated model fit* >>>>>>>>>>> 4055.935 Saturated model df* >>>>>>>>>>> 1767 Difference Chi-squared >>>>>>>> 3.276 Difference d.f. >>>>>>>>>>>>>>> 3 Probability >>>>>>>>>>>>>>>>>>>> 0.351 Akaike's Information Criterion > -2.724 * Saturated model statistic supplied by user 4 Confidence intervals requested in group 4 Matrix Element Int. Estimate Lower Upper Lfail Ufail S 4 1 1 95.0 0.4022 0.0108 0.7590 0 1 0 1 S 4 1 2 95.0 0.0000 0.0000 0.0000 0 0 0 0 S 4 1 3 95.0 0.2179 0.1902 0.2500 0 0 0 0 S 4 1 4 95.0 0.3799 0.0228 0.7727 0 1 0 1 This problem used 11.0% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 1.73 Execution 0: 0: 0: 3.27 TOTAL 0: 0: 0: 5.00 Total number of warnings issued: 0 ______________________________________________________________________________