** Mx startup successful ** !@human; **MX-PC 1.63h** Job started on 03/05/07 at 09:35:18 ! 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 663 data vectors for analysis NOTE: Vectors contain a total of 8525 observations SELECT BMI1 BMI2 ; BEGIN MATRICES = GROUP 1; ! COPY ALL MATRICES FROM GROUP 1, BOTH DECLARED AND COMPUTED NOTE: Selection yields 661 data vectors for analysis NOTE: Vectors contain a total of 1298 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 400 data vectors for analysis NOTE: Vectors contain a total of 5135 observations SELECT BMI1 BMI2 ; BEGIN MATRICES = GROUP 1; NOTE: Selection yields 400 data vectors for analysis NOTE: Vectors contain a total of 780 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 647.0000 651.0000 Mean 21.9076 21.8684 Variance 0.9801 0.9528 Minimum 19.4150 18.5724 Maximum 26.0885 25.9592 Summary of VL file data for group 3 BMI1 BMI2 Code 1.0000 2.0000 Number 389.0000 391.0000 Mean 21.8207 21.8615 Variance 0.9155 1.0274 Minimum 19.4145 19.6650 Maximum 26.1492 25.6367 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.5613 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.1084 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.2995 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.3293 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 ADD_GEN 0.7492 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.5473 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.5613 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.1084 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.2995 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.9075 21.8676 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.3293 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 ADD_GEN 0.7492 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.5473 Vector of OBSERVED means BMI1 BMI2 Mean 21.9076 21.8684 Vector of EXPECTED means BMI1 BMI2 Mean 21.9075 21.8676 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX BMI1 BMI2 BMI1 0.9693 BMI2 0.6698 0.9693 Function value of this group: 3228.1306 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.5613 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.1084 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.2995 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.8219 21.8647 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.3293 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 ADD_GEN 0.7492 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.5473 Vector of OBSERVED means BMI1 BMI2 Mean 21.8207 21.8615 Vector of EXPECTED means BMI1 BMI2 Mean 21.8219 21.8647 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX BMI1 BMI2 BMI1 0.9693 BMI2 0.3078 0.9693 Function value of this group: 2150.2056 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.5613 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.1084 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.2995 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.5613 0.0000 0.2995 0.1084 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.5791 0.0000 0.3090 0.1119 MATRIX V This is a computed FULL matrix of order 1 by 1 [=A+C+E+D] 1 1 0.9693 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 DOM_GEN 0.3293 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 ADD_GEN 0.7492 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.5473 Your model has 7 estimated parameters and 2078 Observed statistics -2 times log-likelihood of data >>> 5378.336 Degrees of freedom >>>>>>>>>>>>>>>> 2071 Akaike's Information Criterion >>>> 1236.336 Bayesian Information Criterion >>>> -4525.126 Sample size Adjusted BIC >>>> -1236.202 Deviance Information Criterion >>>> -2622.005 Saturated model fit* >>>>>>>>>>> 4055.935 Saturated model df* >>>>>>>>>>> 1767 Difference Chi-squared >>>>>>>> 1322.401 Difference d.f. >>>>>>>>>>>>>>> 304 Probability >>>>>>>>>>>>>>>>>>>> 0.000 Akaike's Information Criterion > 714.401 * 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.5791 0.2131 0.7233 0 1 0 0 S 4 1 2 95.0 0.0000 0.0000 0.0000 0 0 0 1 S 4 1 3 95.0 0.3090 0.2737 0.3489 0 0 0 1 S 4 1 4 95.0 0.1119 0.0000 0.4809 0 0 0 1 This problem used 0.4% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 1.97 Execution 0: 0: 0:19.62 TOTAL 0: 0: 0:21.59 Total number of warnings issued: 0 ______________________________________________________________________________ ______________________________________________________________________________