** Mx startup successful ** !@human; **MX-PC 1.63h** Job started on 03/10/06 at 13:55:35 ! ESTIMATE GENETIC AND ENVIRONMENTAL COMPONENTS - ACED MODEL ! OZ FINGERPRINT DATA - MALE AND FEMALE CORRECTED FOR SEX !#LOOP $VARIABLE START STOP INCREMENT !EG Note: #loop section begins here 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 L FULL NVAR 1 FREE ! QTL GENETIC PATH, Q H FULL 1 1 ! SCALAR, 0.5 I FULL 1 1 ! SCALAR, 0.25 M FULL 1 1 FREE END MATRICES; LABEL ROW X ADD_GEN LABEL ROW Y COMM_ENV LABEL ROW Z SPEC_ENV LABEL ROW W DOM_GEN LABEL ROW L QTL MATRIX H .5 MATRIX I .25 START .6 ALL ! START ALL FREE PARAMETERS AT .6 START 0 M 1 1 1 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 Q= L*L'; ! QTL VARIANCE, Q^2 END ALGEBRA; OPTION NO_OUTPUT END The following MX script lines were read for group 2 TITLE G2: MZ DATA DATA NINPUT=2 RECTANGULAR FILE=MZDAT.REC Rectangular continuous data read initiated Note: Maximum ordinal/rectangular record length is: 1000 NOTE: Rectangular file contained 125 records with data that contained a total of 250 observations LABELS ARC1 ARC2 SELECT ARC1 ARC2 ; BEGIN MATRICES = GROUP 1; ! COPY ALL MATRICES FROM GROUP 1, BOTH DECLARED AND COMPUTED NOTE: Selection yields 125 data vectors for analysis NOTE: Vectors contain a total of 250 observations MEANS M|M ; ! 1X2 VECTOR OF EXPECTED MEANS COVARIANCE ! 2X2 COVARIANCE MATRIX A+C+E+D+Q | A+C+D+Q _ A+C+D+Q | A+C+E+D+Q ; OPTION RSIDUALS; OPTION NO_OUTPUT END The following MX script lines were read for group 3 TITLE G3: DZ DATA DATA NINPUT=15 RECTANGULAR FILE=GRID39.REC Rectangular continuous data read initiated Note: Maximum ordinal/rectangular record length is: 1000 NOTE: Rectangular file contained 392 records with data that contained a total of 4868 observations LABELS PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ARC1 ARC2 ARC3 ARC4 ARC5 SELECT PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ARC1 ARC2 ARC3 ARC4 ARC5 ; DEFINITON PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ; NOTE: Selection yields 392 data vectors for analysis NOTE: Vectors contain a total of 4868 observations NOTE: Definition yields 392 data vectors for analysis NOTE: Vectors contain a total of 948 observations BEGIN MATRICES ; A COMPUTED NVAR NVAR =A1 C COMPUTED NVAR NVAR =C1 D COMPUTED NVAR NVAR =D1 E COMPUTED NVAR NVAR =E1 Q COMPUTED NVAR NVAR =Q1 M FULL 1 1 =M1 H FULL 1 1 =H1 I FULL 1 1 =I1 P FULL 1 1 R FULL 1 1 S FULL 1 1 T FULL 1 1 U FULL 1 1 V FULL 1 1 W FULL 1 1 X FULL 1 1 Y FULL 1 1 Z FULL 1 1 END MATRICES; SPECIFY P PAB SPECIFY R PAC SPECIFY S PAD SPECIFY T PAE SPECIFY U PBC SPECIFY V PBD SPECIFY W PBE SPECIFY X PCD SPECIFY Y PCE SPECIFY Z PDE MEANS M|M|M|M|M ; COVARIANCE A+C+E+D+Q | H@A+C+I@D+P@Q | H@A+C+I@D+R@Q | H@A+C+I@D+S@Q | H@A+C+I@D+T@Q _ H@A+C+I@D+P@Q | A+C+E+D+Q | H@A+C+I@D+U@Q | H@A+C+I@D+V@Q | H@A+C+I@D+W@Q _ H@A+C+I@D+R@Q | H@A+C+I@D+U@Q | A+C+E+D+Q | H@A+C+I@D+X@Q | H@A+C+I@D+Y@Q _ H@A+C+I@D+S@Q | H@A+C+I@D+V@Q | H@A+C+I@D+X@Q | A+C+E+D+Q | H@A+C+I@D+Z@Q _ H@A+C+I@D+T@Q | H@A+C+I@D+W@Q | H@A+C+I@D+Y@Q | H@A+C+I@D+Z@Q | A+C+E+D+Q ; OPTION NO_OUTPUT END The following MX script lines were read for group 4 TITLE G4: STANDARDIZATION CALCULATION BEGIN MATRICES = GROUP 1; END MATRICES; BEGIN ALGEBRA; V=A+C+E+D+Q; ! TOTAL VARIANCE P=A|C|E|D|Q; ! CONCATENATE PARAMETER ESTIMATES S=P@V~; ! STANDARDIZED PARAMETER ESTIMATES END ALGEBRA; OPTION NDECIMALS=4 OPTION SUB=3176.157,1194 ! LIKELIHOOD AND DF FROM SATURATED MODEL: OZBMIYFSAT.MXO !DROP L 1 1 1 END Summary of VL file data for group 2 ARC1 ARC2 Code 1.0000 2.0000 Number 125.0000 125.0000 Mean -0.0090 -0.0190 Variance 0.8861 0.8757 Minimum -2.3713 -2.9867 Maximum 2.6268 2.2563 Summary of VL file data for group 3 PDE PCE PCD PBE PBD PBC PAE Code -10.0000 -9.0000 -8.0000 -7.0000 -6.0000 -5.0000 -4.0000 Number 392.0000 392.0000 392.0000 392.0000 392.0000 392.0000 392.0000 Mean 0.4988 0.5012 0.5066 0.4998 0.4947 0.5160 0.5011 Variance 0.0005 0.0006 0.0081 0.0000 0.0092 0.0313 0.0005 Minimum 0.0400 0.4650 0.0020 0.4180 0.0010 0.0020 0.4920 Maximum 0.5030 0.9970 0.9980 0.5000 0.9980 0.9980 0.9590 PAD PAC PAB ARC1 ARC2 ARC3 ARC4 Code -3.0000 -2.0000 -1.0000 1.0000 2.0000 3.0000 4.0000 Number 392.0000 392.0000 392.0000 392.0000 392.0000 132.0000 29.0000 Mean 0.5028 0.5069 0.4811 0.0844 0.0470 0.2260 0.2899 Variance 0.0074 0.0379 0.1050 0.9724 1.0117 1.1402 1.0366 Minimum 0.0010 0.0010 0.0010 -2.0600 -2.2000 -1.9900 -1.8400 Maximum 0.9980 0.9980 0.9980 3.1560 3.4310 3.5340 2.7000 ARC5 Code 5.0000 Number 3.0000 Mean -0.0407 Variance 0.5965 Minimum -0.9500 Maximum 0.9380 PARAMETER SPECIFICATIONS 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 I This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a FULL matrix of order 1 by 1 1 QTL 4 MATRIX M This is a FULL matrix of order 1 by 1 1 1 5 MATRIX P This is a computed FULL matrix of order 1 by 5 It has no free parameters specified MATRIX Q This is a computed 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 5 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 = 5 MX PARAMETER ESTIMATES GROUP NUMBER: 4 TITLE G4: STANDARDIZATION MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.7485 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.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.2384 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX I This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX L This is a FULL matrix of order 1 by 1 1 QTL 0.0000 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0656 MATRIX P This is a computed FULL matrix of order 1 by 5 [=A|C|E|D|Q] 1 2 3 4 5 1 0.7485 0.0000 0.2384 0.0000 0.0000 MATRIX Q This is a computed FULL matrix of order 1 by 1 [=L*L'] 1 1 0.0000 MATRIX S This is a computed FULL matrix of order 1 by 5 [=P@V~] 1 2 3 4 5 1 0.7584 0.0000 0.2416 0.0000 0.0000 MATRIX V This is a computed FULL matrix of order 1 by 1 [=A+C+E+D+Q] 1 1 0.9869 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 DOM_GEN 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 ADD_GEN 0.8651 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.4883 Your model has 5 estimated parameters and 1198 Observed statistics -2 times log-likelihood of data >>> 3176.157 Degrees of freedom >>>>>>>>>>>>>>>> 1193 Akaike's Information Criterion >>>> 790.157 Bayesian Information Criterion >>>> -2138.879 Sample size Adjusted BIC >>>> -245.473 Deviance Information Criterion >>>> -1042.586 Submodel fit* >>>>>>>>>>> 3176.157 Submodel df* >>>>>>>>>>> 1194 Difference Chi-squared >>>>>>>> 0.000 Difference d.f. >>>>>>>>>>>>>>> 1 Probability >>>>>>>>>>>>>>>>>>>> 0.986 Akaike's Information Criterion > -2.000 * Submodel statistic supplied by user This problem used 0.0% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.23 Execution 0: 0: 0:15.72 TOTAL 0: 0: 0:15.95 Total number of warnings issued: 0 ______________________________________________________________________________ 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 L FULL NVAR 1 FREE ! QTL GENETIC PATH, Q H FULL 1 1 ! SCALAR, 0.5 I FULL 1 1 ! SCALAR, 0.25 M FULL 1 1 FREE END MATRICES; LABEL ROW X ADD_GEN LABEL ROW Y COMM_ENV LABEL ROW Z SPEC_ENV LABEL ROW W DOM_GEN LABEL ROW L QTL MATRIX H .5 MATRIX I .25 START .6 ALL ! START ALL FREE PARAMETERS AT .6 START 0 M 1 1 1 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 Q= L*L'; ! QTL VARIANCE, Q^2 END ALGEBRA; OPTION NO_OUTPUT END The following MX script lines were read for group 2 TITLE G2: MZ DATA DATA NINPUT=2 RECTANGULAR FILE=MZDAT.REC Rectangular continuous data read initiated Note: Maximum ordinal/rectangular record length is: 1000 NOTE: Rectangular file contained 125 records with data that contained a total of 250 observations LABELS ARC1 ARC2 SELECT ARC1 ARC2 ; BEGIN MATRICES = GROUP 1; ! COPY ALL MATRICES FROM GROUP 1, BOTH DECLARED AND COMPUTED NOTE: Selection yields 125 data vectors for analysis NOTE: Vectors contain a total of 250 observations MEANS M|M ; ! 1X2 VECTOR OF EXPECTED MEANS COVARIANCE ! 2X2 COVARIANCE MATRIX A+C+E+D+Q | A+C+D+Q _ A+C+D+Q | A+C+E+D+Q ; OPTION RSIDUALS; OPTION NO_OUTPUT END The following MX script lines were read for group 3 TITLE G3: DZ DATA DATA NINPUT=15 RECTANGULAR FILE=GRID40.REC Rectangular continuous data read initiated Note: Maximum ordinal/rectangular record length is: 1000 NOTE: Rectangular file contained 392 records with data that contained a total of 4868 observations LABELS PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ARC1 ARC2 ARC3 ARC4 ARC5 SELECT PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ARC1 ARC2 ARC3 ARC4 ARC5 ; DEFINITON PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ; NOTE: Selection yields 392 data vectors for analysis NOTE: Vectors contain a total of 4868 observations Note: Global variable previously defined. Updating PAB Note: Global variable previously defined. Updating PAC Note: Global variable previously defined. Updating PAD Note: Global variable previously defined. Updating PAE Note: Global variable previously defined. Updating PBC Note: Global variable previously defined. Updating PBD Note: Global variable previously defined. Updating PBE Note: Global variable previously defined. Updating PCD Note: Global variable previously defined. Updating PCE Note: Global variable previously defined. Updating PDE NOTE: Definition yields 392 data vectors for analysis NOTE: Vectors contain a total of 948 observations BEGIN MATRICES ; A COMPUTED NVAR NVAR =A1 C COMPUTED NVAR NVAR =C1 D COMPUTED NVAR NVAR =D1 E COMPUTED NVAR NVAR =E1 Q COMPUTED NVAR NVAR =Q1 M FULL 1 1 =M1 H FULL 1 1 =H1 I FULL 1 1 =I1 P FULL 1 1 R FULL 1 1 S FULL 1 1 T FULL 1 1 U FULL 1 1 V FULL 1 1 W FULL 1 1 X FULL 1 1 Y FULL 1 1 Z FULL 1 1 END MATRICES; SPECIFY P PAB SPECIFY R PAC SPECIFY S PAD SPECIFY T PAE SPECIFY U PBC SPECIFY V PBD SPECIFY W PBE SPECIFY X PCD SPECIFY Y PCE SPECIFY Z PDE MEANS M|M|M|M|M ; COVARIANCE A+C+E+D+Q | H@A+C+I@D+P@Q | H@A+C+I@D+R@Q | H@A+C+I@D+S@Q | H@A+C+I@D+T@Q _ H@A+C+I@D+P@Q | A+C+E+D+Q | H@A+C+I@D+U@Q | H@A+C+I@D+V@Q | H@A+C+I@D+W@Q _ H@A+C+I@D+R@Q | H@A+C+I@D+U@Q | A+C+E+D+Q | H@A+C+I@D+X@Q | H@A+C+I@D+Y@Q _ H@A+C+I@D+S@Q | H@A+C+I@D+V@Q | H@A+C+I@D+X@Q | A+C+E+D+Q | H@A+C+I@D+Z@Q _ H@A+C+I@D+T@Q | H@A+C+I@D+W@Q | H@A+C+I@D+Y@Q | H@A+C+I@D+Z@Q | A+C+E+D+Q ; OPTION NO_OUTPUT END The following MX script lines were read for group 4 TITLE G4: STANDARDIZATION CALCULATION BEGIN MATRICES = GROUP 1; END MATRICES; BEGIN ALGEBRA; V=A+C+E+D+Q; ! TOTAL VARIANCE P=A|C|E|D|Q; ! CONCATENATE PARAMETER ESTIMATES S=P@V~; ! STANDARDIZED PARAMETER ESTIMATES END ALGEBRA; OPTION NDECIMALS=4 OPTION SUB=3176.157,1194 ! LIKELIHOOD AND DF FROM SATURATED MODEL: OZBMIYFSAT.MXO !DROP L 1 1 1 END Summary of VL file data for group 2 ARC1 ARC2 Code 1.0000 2.0000 Number 125.0000 125.0000 Mean -0.0090 -0.0190 Variance 0.8861 0.8757 Minimum -2.3713 -2.9867 Maximum 2.6268 2.2563 Summary of VL file data for group 3 PDE PCE PCD PBE PBD PBC PAE Code -10.0000 -9.0000 -8.0000 -7.0000 -6.0000 -5.0000 -4.0000 Number 392.0000 392.0000 392.0000 392.0000 392.0000 392.0000 392.0000 Mean 0.4994 0.5006 0.5032 0.4993 0.4941 0.5153 0.5004 Variance 0.0001 0.0008 0.0066 0.0002 0.0072 0.0247 0.0001 Minimum 0.2830 0.2370 0.0080 0.2090 0.0050 0.0060 0.4460 Maximum 0.5000 0.9930 0.9940 0.5000 0.9940 0.9940 0.7230 PAD PAC PAB ARC1 ARC2 ARC3 ARC4 Code -3.0000 -2.0000 -1.0000 1.0000 2.0000 3.0000 4.0000 Number 392.0000 392.0000 392.0000 392.0000 392.0000 132.0000 29.0000 Mean 0.5007 0.5081 0.4901 0.0844 0.0470 0.2260 0.2899 Variance 0.0062 0.0284 0.0878 0.9724 1.0117 1.1402 1.0366 Minimum 0.0050 0.0050 0.0050 -2.0600 -2.2000 -1.9900 -1.8400 Maximum 0.9930 0.9930 0.9940 3.1560 3.4310 3.5340 2.7000 ARC5 Code 5.0000 Number 3.0000 Mean -0.0407 Variance 0.5965 Minimum -0.9500 Maximum 0.9380 PARAMETER SPECIFICATIONS 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 I This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a FULL matrix of order 1 by 1 1 QTL 4 MATRIX M This is a FULL matrix of order 1 by 1 1 1 5 MATRIX P This is a computed FULL matrix of order 1 by 5 It has no free parameters specified MATRIX Q This is a computed 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 5 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 = 5 MX PARAMETER ESTIMATES GROUP NUMBER: 4 TITLE G4: STANDARDIZATION MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.7443 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.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.2384 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX I This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX L This is a FULL matrix of order 1 by 1 1 QTL -0.0645 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0656 MATRIX P This is a computed FULL matrix of order 1 by 5 [=A|C|E|D|Q] 1 2 3 4 5 1 0.7443 0.0000 0.2384 0.0000 0.0042 MATRIX Q This is a computed FULL matrix of order 1 by 1 [=L*L'] 1 1 0.0042 MATRIX S This is a computed FULL matrix of order 1 by 5 [=P@V~] 1 2 3 4 5 1 0.7542 0.0000 0.2416 0.0000 0.0042 MATRIX V This is a computed FULL matrix of order 1 by 1 [=A+C+E+D+Q] 1 1 0.9869 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 DOM_GEN 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 ADD_GEN 0.8627 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.4883 Your model has 5 estimated parameters and 1198 Observed statistics -2 times log-likelihood of data >>> 3176.155 Degrees of freedom >>>>>>>>>>>>>>>> 1193 Akaike's Information Criterion >>>> 790.155 Bayesian Information Criterion >>>> -2138.880 Sample size Adjusted BIC >>>> -245.474 Deviance Information Criterion >>>> -1042.587 Submodel fit* >>>>>>>>>>> 3176.157 Submodel df* >>>>>>>>>>> 1194 Difference Chi-squared >>>>>>>> 0.002 Difference d.f. >>>>>>>>>>>>>>> 1 Probability >>>>>>>>>>>>>>>>>>>> 0.961 Akaike's Information Criterion > -1.998 * Submodel statistic supplied by user This problem used 0.0% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.14 Execution 0: 0: 0:12.25 TOTAL 0: 0: 0:12.39 Total number of warnings issued: 0 ______________________________________________________________________________ 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 L FULL NVAR 1 FREE ! QTL GENETIC PATH, Q H FULL 1 1 ! SCALAR, 0.5 I FULL 1 1 ! SCALAR, 0.25 M FULL 1 1 FREE END MATRICES; LABEL ROW X ADD_GEN LABEL ROW Y COMM_ENV LABEL ROW Z SPEC_ENV LABEL ROW W DOM_GEN LABEL ROW L QTL MATRIX H .5 MATRIX I .25 START .6 ALL ! START ALL FREE PARAMETERS AT .6 START 0 M 1 1 1 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 Q= L*L'; ! QTL VARIANCE, Q^2 END ALGEBRA; OPTION NO_OUTPUT END The following MX script lines were read for group 2 TITLE G2: MZ DATA DATA NINPUT=2 RECTANGULAR FILE=MZDAT.REC Rectangular continuous data read initiated Note: Maximum ordinal/rectangular record length is: 1000 NOTE: Rectangular file contained 125 records with data that contained a total of 250 observations LABELS ARC1 ARC2 SELECT ARC1 ARC2 ; BEGIN MATRICES = GROUP 1; ! COPY ALL MATRICES FROM GROUP 1, BOTH DECLARED AND COMPUTED NOTE: Selection yields 125 data vectors for analysis NOTE: Vectors contain a total of 250 observations MEANS M|M ; ! 1X2 VECTOR OF EXPECTED MEANS COVARIANCE ! 2X2 COVARIANCE MATRIX A+C+E+D+Q | A+C+D+Q _ A+C+D+Q | A+C+E+D+Q ; OPTION RSIDUALS; OPTION NO_OUTPUT END The following MX script lines were read for group 3 TITLE G3: DZ DATA DATA NINPUT=15 RECTANGULAR FILE=GRID41.REC Rectangular continuous data read initiated Note: Maximum ordinal/rectangular record length is: 1000 NOTE: Rectangular file contained 392 records with data that contained a total of 4868 observations LABELS PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ARC1 ARC2 ARC3 ARC4 ARC5 SELECT PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ARC1 ARC2 ARC3 ARC4 ARC5 ; DEFINITON PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ; NOTE: Selection yields 392 data vectors for analysis NOTE: Vectors contain a total of 4868 observations Note: Global variable previously defined. Updating PAB Note: Global variable previously defined. Updating PAC Note: Global variable previously defined. Updating PAD Note: Global variable previously defined. Updating PAE Note: Global variable previously defined. Updating PBC Note: Global variable previously defined. Updating PBD Note: Global variable previously defined. Updating PBE Note: Global variable previously defined. Updating PCD Note: Global variable previously defined. Updating PCE Note: Global variable previously defined. Updating PDE NOTE: Definition yields 392 data vectors for analysis NOTE: Vectors contain a total of 948 observations BEGIN MATRICES ; A COMPUTED NVAR NVAR =A1 C COMPUTED NVAR NVAR =C1 D COMPUTED NVAR NVAR =D1 E COMPUTED NVAR NVAR =E1 Q COMPUTED NVAR NVAR =Q1 M FULL 1 1 =M1 H FULL 1 1 =H1 I FULL 1 1 =I1 P FULL 1 1 R FULL 1 1 S FULL 1 1 T FULL 1 1 U FULL 1 1 V FULL 1 1 W FULL 1 1 X FULL 1 1 Y FULL 1 1 Z FULL 1 1 END MATRICES; SPECIFY P PAB SPECIFY R PAC SPECIFY S PAD SPECIFY T PAE SPECIFY U PBC SPECIFY V PBD SPECIFY W PBE SPECIFY X PCD SPECIFY Y PCE SPECIFY Z PDE MEANS M|M|M|M|M ; COVARIANCE A+C+E+D+Q | H@A+C+I@D+P@Q | H@A+C+I@D+R@Q | H@A+C+I@D+S@Q | H@A+C+I@D+T@Q _ H@A+C+I@D+P@Q | A+C+E+D+Q | H@A+C+I@D+U@Q | H@A+C+I@D+V@Q | H@A+C+I@D+W@Q _ H@A+C+I@D+R@Q | H@A+C+I@D+U@Q | A+C+E+D+Q | H@A+C+I@D+X@Q | H@A+C+I@D+Y@Q _ H@A+C+I@D+S@Q | H@A+C+I@D+V@Q | H@A+C+I@D+X@Q | A+C+E+D+Q | H@A+C+I@D+Z@Q _ H@A+C+I@D+T@Q | H@A+C+I@D+W@Q | H@A+C+I@D+Y@Q | H@A+C+I@D+Z@Q | A+C+E+D+Q ; OPTION NO_OUTPUT END The following MX script lines were read for group 4 TITLE G4: STANDARDIZATION CALCULATION BEGIN MATRICES = GROUP 1; END MATRICES; BEGIN ALGEBRA; V=A+C+E+D+Q; ! TOTAL VARIANCE P=A|C|E|D|Q; ! CONCATENATE PARAMETER ESTIMATES S=P@V~; ! STANDARDIZED PARAMETER ESTIMATES END ALGEBRA; OPTION NDECIMALS=4 OPTION SUB=3176.157,1194 ! LIKELIHOOD AND DF FROM SATURATED MODEL: OZBMIYFSAT.MXO !DROP L 1 1 1 END Summary of VL file data for group 2 ARC1 ARC2 Code 1.0000 2.0000 Number 125.0000 125.0000 Mean -0.0090 -0.0190 Variance 0.8861 0.8757 Minimum -2.3713 -2.9867 Maximum 2.6268 2.2563 Summary of VL file data for group 3 PDE PCE PCD PBE PBD PBC PAE Code -10.0000 -9.0000 -8.0000 -7.0000 -6.0000 -5.0000 -4.0000 Number 392.0000 392.0000 392.0000 392.0000 392.0000 392.0000 392.0000 Mean 0.5000 0.5000 0.5005 0.4987 0.4935 0.5166 0.4997 Variance 0.0000 0.0013 0.0080 0.0006 0.0104 0.0315 0.0000 Minimum 0.4950 0.0000 0.0000 0.0000 0.0000 0.0000 0.4020 Maximum 0.5000 0.9990 0.9990 0.5000 0.9990 0.9990 0.5000 PAD PAC PAB ARC1 ARC2 ARC3 ARC4 Code -3.0000 -2.0000 -1.0000 1.0000 2.0000 3.0000 4.0000 Number 392.0000 392.0000 392.0000 392.0000 392.0000 132.0000 29.0000 Mean 0.4984 0.5124 0.5447 0.0844 0.0470 0.2260 0.2899 Variance 0.0087 0.0361 0.1290 0.9724 1.0117 1.1402 1.0366 Minimum 0.0000 0.0000 0.0000 -2.0600 -2.2000 -1.9900 -1.8400 Maximum 0.9990 0.9990 1.5000 3.1560 3.4310 3.5340 2.7000 ARC5 Code 5.0000 Number 3.0000 Mean -0.0407 Variance 0.5965 Minimum -0.9500 Maximum 0.9380 PARAMETER SPECIFICATIONS 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 I This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a FULL matrix of order 1 by 1 1 QTL 4 MATRIX M This is a FULL matrix of order 1 by 1 1 1 5 MATRIX P This is a computed FULL matrix of order 1 by 5 It has no free parameters specified MATRIX Q This is a computed 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 5 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 = 5 MX PARAMETER ESTIMATES GROUP NUMBER: 4 TITLE G4: STANDARDIZATION MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.6659 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.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.2398 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX I This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX L This is a FULL matrix of order 1 by 1 1 QTL 0.2861 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0651 MATRIX P This is a computed FULL matrix of order 1 by 5 [=A|C|E|D|Q] 1 2 3 4 5 1 0.6659 0.0000 0.2398 0.0000 0.0818 MATRIX Q This is a computed FULL matrix of order 1 by 1 [=L*L'] 1 1 0.0818 MATRIX S This is a computed FULL matrix of order 1 by 5 [=P@V~] 1 2 3 4 5 1 0.6743 0.0000 0.2429 0.0000 0.0829 MATRIX V This is a computed FULL matrix of order 1 by 1 [=A+C+E+D+Q] 1 1 0.9875 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 DOM_GEN 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 ADD_GEN 0.8160 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.4897 Your model has 5 estimated parameters and 1198 Observed statistics -2 times log-likelihood of data >>> 3175.095 Degrees of freedom >>>>>>>>>>>>>>>> 1193 Akaike's Information Criterion >>>> 789.095 Bayesian Information Criterion >>>> -2139.410 Sample size Adjusted BIC >>>> -246.004 Deviance Information Criterion >>>> -1043.116 Submodel fit* >>>>>>>>>>> 3176.157 Submodel df* >>>>>>>>>>> 1194 Difference Chi-squared >>>>>>>> 1.062 Difference d.f. >>>>>>>>>>>>>>> 1 Probability >>>>>>>>>>>>>>>>>>>> 0.303 Akaike's Information Criterion > -0.938 * Submodel statistic supplied by user This problem used 0.0% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.14 Execution 0: 0: 0:10.41 TOTAL 0: 0: 0:10.55 Total number of warnings issued: 0 ______________________________________________________________________________ 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 L FULL NVAR 1 FREE ! QTL GENETIC PATH, Q H FULL 1 1 ! SCALAR, 0.5 I FULL 1 1 ! SCALAR, 0.25 M FULL 1 1 FREE END MATRICES; LABEL ROW X ADD_GEN LABEL ROW Y COMM_ENV LABEL ROW Z SPEC_ENV LABEL ROW W DOM_GEN LABEL ROW L QTL MATRIX H .5 MATRIX I .25 START .6 ALL ! START ALL FREE PARAMETERS AT .6 START 0 M 1 1 1 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 Q= L*L'; ! QTL VARIANCE, Q^2 END ALGEBRA; OPTION NO_OUTPUT END The following MX script lines were read for group 2 TITLE G2: MZ DATA DATA NINPUT=2 RECTANGULAR FILE=MZDAT.REC Rectangular continuous data read initiated Note: Maximum ordinal/rectangular record length is: 1000 NOTE: Rectangular file contained 125 records with data that contained a total of 250 observations LABELS ARC1 ARC2 SELECT ARC1 ARC2 ; BEGIN MATRICES = GROUP 1; ! COPY ALL MATRICES FROM GROUP 1, BOTH DECLARED AND COMPUTED NOTE: Selection yields 125 data vectors for analysis NOTE: Vectors contain a total of 250 observations MEANS M|M ; ! 1X2 VECTOR OF EXPECTED MEANS COVARIANCE ! 2X2 COVARIANCE MATRIX A+C+E+D+Q | A+C+D+Q _ A+C+D+Q | A+C+E+D+Q ; OPTION RSIDUALS; OPTION NO_OUTPUT END The following MX script lines were read for group 3 TITLE G3: DZ DATA DATA NINPUT=15 RECTANGULAR FILE=GRID42.REC Rectangular continuous data read initiated Note: Maximum ordinal/rectangular record length is: 1000 NOTE: Rectangular file contained 392 records with data that contained a total of 4868 observations LABELS PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ARC1 ARC2 ARC3 ARC4 ARC5 SELECT PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ARC1 ARC2 ARC3 ARC4 ARC5 ; DEFINITON PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ; NOTE: Selection yields 392 data vectors for analysis NOTE: Vectors contain a total of 4868 observations Note: Global variable previously defined. Updating PAB Note: Global variable previously defined. Updating PAC Note: Global variable previously defined. Updating PAD Note: Global variable previously defined. Updating PAE Note: Global variable previously defined. Updating PBC Note: Global variable previously defined. Updating PBD Note: Global variable previously defined. Updating PBE Note: Global variable previously defined. Updating PCD Note: Global variable previously defined. Updating PCE Note: Global variable previously defined. Updating PDE NOTE: Definition yields 392 data vectors for analysis NOTE: Vectors contain a total of 948 observations BEGIN MATRICES ; A COMPUTED NVAR NVAR =A1 C COMPUTED NVAR NVAR =C1 D COMPUTED NVAR NVAR =D1 E COMPUTED NVAR NVAR =E1 Q COMPUTED NVAR NVAR =Q1 M FULL 1 1 =M1 H FULL 1 1 =H1 I FULL 1 1 =I1 P FULL 1 1 R FULL 1 1 S FULL 1 1 T FULL 1 1 U FULL 1 1 V FULL 1 1 W FULL 1 1 X FULL 1 1 Y FULL 1 1 Z FULL 1 1 END MATRICES; SPECIFY P PAB SPECIFY R PAC SPECIFY S PAD SPECIFY T PAE SPECIFY U PBC SPECIFY V PBD SPECIFY W PBE SPECIFY X PCD SPECIFY Y PCE SPECIFY Z PDE MEANS M|M|M|M|M ; COVARIANCE A+C+E+D+Q | H@A+C+I@D+P@Q | H@A+C+I@D+R@Q | H@A+C+I@D+S@Q | H@A+C+I@D+T@Q _ H@A+C+I@D+P@Q | A+C+E+D+Q | H@A+C+I@D+U@Q | H@A+C+I@D+V@Q | H@A+C+I@D+W@Q _ H@A+C+I@D+R@Q | H@A+C+I@D+U@Q | A+C+E+D+Q | H@A+C+I@D+X@Q | H@A+C+I@D+Y@Q _ H@A+C+I@D+S@Q | H@A+C+I@D+V@Q | H@A+C+I@D+X@Q | A+C+E+D+Q | H@A+C+I@D+Z@Q _ H@A+C+I@D+T@Q | H@A+C+I@D+W@Q | H@A+C+I@D+Y@Q | H@A+C+I@D+Z@Q | A+C+E+D+Q ; OPTION NO_OUTPUT END The following MX script lines were read for group 4 TITLE G4: STANDARDIZATION CALCULATION BEGIN MATRICES = GROUP 1; END MATRICES; BEGIN ALGEBRA; V=A+C+E+D+Q; ! TOTAL VARIANCE P=A|C|E|D|Q; ! CONCATENATE PARAMETER ESTIMATES S=P@V~; ! STANDARDIZED PARAMETER ESTIMATES END ALGEBRA; OPTION NDECIMALS=4 OPTION SUB=3176.157,1194 ! LIKELIHOOD AND DF FROM SATURATED MODEL: OZBMIYFSAT.MXO !DROP L 1 1 1 END Summary of VL file data for group 2 ARC1 ARC2 Code 1.0000 2.0000 Number 125.0000 125.0000 Mean -0.0090 -0.0190 Variance 0.8861 0.8757 Minimum -2.3713 -2.9867 Maximum 2.6268 2.2563 Summary of VL file data for group 3 PDE PCE PCD PBE PBD PBC PAE Code -10.0000 -9.0000 -8.0000 -7.0000 -6.0000 -5.0000 -4.0000 Number 392.0000 392.0000 392.0000 392.0000 392.0000 392.0000 392.0000 Mean 0.4983 0.5000 0.5006 0.5004 0.4953 0.5221 0.5013 Variance 0.0006 0.0010 0.0084 0.0004 0.0101 0.0274 0.0004 Minimum 0.1630 0.1630 0.0020 0.3380 0.0020 0.0020 0.5000 Maximum 0.5000 0.9970 0.9970 0.8360 0.9970 0.9970 0.8360 PAD PAC PAB ARC1 ARC2 ARC3 ARC4 Code -3.0000 -2.0000 -1.0000 1.0000 2.0000 3.0000 4.0000 Number 392.0000 392.0000 392.0000 392.0000 392.0000 132.0000 29.0000 Mean 0.5005 0.5120 0.4954 0.0844 0.0470 0.2260 0.2899 Variance 0.0066 0.0337 0.0949 0.9724 1.0117 1.1402 1.0366 Minimum 0.0020 0.0020 0.0020 -2.0600 -2.2000 -1.9900 -1.8400 Maximum 0.9970 0.9970 0.9970 3.1560 3.4310 3.5340 2.7000 ARC5 Code 5.0000 Number 3.0000 Mean -0.0407 Variance 0.5965 Minimum -0.9500 Maximum 0.9380 PARAMETER SPECIFICATIONS 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 I This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a FULL matrix of order 1 by 1 1 QTL 4 MATRIX M This is a FULL matrix of order 1 by 1 1 1 5 MATRIX P This is a computed FULL matrix of order 1 by 5 It has no free parameters specified MATRIX Q This is a computed 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 5 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 = 5 MX PARAMETER ESTIMATES GROUP NUMBER: 4 TITLE G4: STANDARDIZATION MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.6265 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.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.2395 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX I This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX L This is a FULL matrix of order 1 by 1 1 QTL 0.3478 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0650 MATRIX P This is a computed FULL matrix of order 1 by 5 [=A|C|E|D|Q] 1 2 3 4 5 1 0.6265 0.0000 0.2395 0.0000 0.1210 MATRIX Q This is a computed FULL matrix of order 1 by 1 [=L*L'] 1 1 0.1210 MATRIX S This is a computed FULL matrix of order 1 by 5 [=P@V~] 1 2 3 4 5 1 0.6348 0.0000 0.2427 0.0000 0.1226 MATRIX V This is a computed FULL matrix of order 1 by 1 [=A+C+E+D+Q] 1 1 0.9870 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 DOM_GEN 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 ADD_GEN 0.7915 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.4894 Your model has 5 estimated parameters and 1198 Observed statistics -2 times log-likelihood of data >>> 3174.427 Degrees of freedom >>>>>>>>>>>>>>>> 1193 Akaike's Information Criterion >>>> 788.427 Bayesian Information Criterion >>>> -2139.744 Sample size Adjusted BIC >>>> -246.338 Deviance Information Criterion >>>> -1043.450 Submodel fit* >>>>>>>>>>> 3176.157 Submodel df* >>>>>>>>>>> 1194 Difference Chi-squared >>>>>>>> 1.730 Difference d.f. >>>>>>>>>>>>>>> 1 Probability >>>>>>>>>>>>>>>>>>>> 0.188 Akaike's Information Criterion > -0.270 * Submodel statistic supplied by user This problem used 0.0% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.14 Execution 0: 0: 0: 9.94 TOTAL 0: 0: 0:10.08 Total number of warnings issued: 0 ______________________________________________________________________________ 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 L FULL NVAR 1 FREE ! QTL GENETIC PATH, Q H FULL 1 1 ! SCALAR, 0.5 I FULL 1 1 ! SCALAR, 0.25 M FULL 1 1 FREE END MATRICES; LABEL ROW X ADD_GEN LABEL ROW Y COMM_ENV LABEL ROW Z SPEC_ENV LABEL ROW W DOM_GEN LABEL ROW L QTL MATRIX H .5 MATRIX I .25 START .6 ALL ! START ALL FREE PARAMETERS AT .6 START 0 M 1 1 1 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 Q= L*L'; ! QTL VARIANCE, Q^2 END ALGEBRA; OPTION NO_OUTPUT END The following MX script lines were read for group 2 TITLE G2: MZ DATA DATA NINPUT=2 RECTANGULAR FILE=MZDAT.REC Rectangular continuous data read initiated Note: Maximum ordinal/rectangular record length is: 1000 NOTE: Rectangular file contained 125 records with data that contained a total of 250 observations LABELS ARC1 ARC2 SELECT ARC1 ARC2 ; BEGIN MATRICES = GROUP 1; ! COPY ALL MATRICES FROM GROUP 1, BOTH DECLARED AND COMPUTED NOTE: Selection yields 125 data vectors for analysis NOTE: Vectors contain a total of 250 observations MEANS M|M ; ! 1X2 VECTOR OF EXPECTED MEANS COVARIANCE ! 2X2 COVARIANCE MATRIX A+C+E+D+Q | A+C+D+Q _ A+C+D+Q | A+C+E+D+Q ; OPTION RSIDUALS; OPTION NO_OUTPUT END The following MX script lines were read for group 3 TITLE G3: DZ DATA DATA NINPUT=15 RECTANGULAR FILE=GRID43.REC Rectangular continuous data read initiated Note: Maximum ordinal/rectangular record length is: 1000 NOTE: Rectangular file contained 392 records with data that contained a total of 4868 observations LABELS PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ARC1 ARC2 ARC3 ARC4 ARC5 SELECT PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ARC1 ARC2 ARC3 ARC4 ARC5 ; DEFINITON PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ; NOTE: Selection yields 392 data vectors for analysis NOTE: Vectors contain a total of 4868 observations Note: Global variable previously defined. Updating PAB Note: Global variable previously defined. Updating PAC Note: Global variable previously defined. Updating PAD Note: Global variable previously defined. Updating PAE Note: Global variable previously defined. Updating PBC Note: Global variable previously defined. Updating PBD Note: Global variable previously defined. Updating PBE Note: Global variable previously defined. Updating PCD Note: Global variable previously defined. Updating PCE Note: Global variable previously defined. Updating PDE NOTE: Definition yields 392 data vectors for analysis NOTE: Vectors contain a total of 948 observations BEGIN MATRICES ; A COMPUTED NVAR NVAR =A1 C COMPUTED NVAR NVAR =C1 D COMPUTED NVAR NVAR =D1 E COMPUTED NVAR NVAR =E1 Q COMPUTED NVAR NVAR =Q1 M FULL 1 1 =M1 H FULL 1 1 =H1 I FULL 1 1 =I1 P FULL 1 1 R FULL 1 1 S FULL 1 1 T FULL 1 1 U FULL 1 1 V FULL 1 1 W FULL 1 1 X FULL 1 1 Y FULL 1 1 Z FULL 1 1 END MATRICES; SPECIFY P PAB SPECIFY R PAC SPECIFY S PAD SPECIFY T PAE SPECIFY U PBC SPECIFY V PBD SPECIFY W PBE SPECIFY X PCD SPECIFY Y PCE SPECIFY Z PDE MEANS M|M|M|M|M ; COVARIANCE A+C+E+D+Q | H@A+C+I@D+P@Q | H@A+C+I@D+R@Q | H@A+C+I@D+S@Q | H@A+C+I@D+T@Q _ H@A+C+I@D+P@Q | A+C+E+D+Q | H@A+C+I@D+U@Q | H@A+C+I@D+V@Q | H@A+C+I@D+W@Q _ H@A+C+I@D+R@Q | H@A+C+I@D+U@Q | A+C+E+D+Q | H@A+C+I@D+X@Q | H@A+C+I@D+Y@Q _ H@A+C+I@D+S@Q | H@A+C+I@D+V@Q | H@A+C+I@D+X@Q | A+C+E+D+Q | H@A+C+I@D+Z@Q _ H@A+C+I@D+T@Q | H@A+C+I@D+W@Q | H@A+C+I@D+Y@Q | H@A+C+I@D+Z@Q | A+C+E+D+Q ; OPTION NO_OUTPUT END The following MX script lines were read for group 4 TITLE G4: STANDARDIZATION CALCULATION BEGIN MATRICES = GROUP 1; END MATRICES; BEGIN ALGEBRA; V=A+C+E+D+Q; ! TOTAL VARIANCE P=A|C|E|D|Q; ! CONCATENATE PARAMETER ESTIMATES S=P@V~; ! STANDARDIZED PARAMETER ESTIMATES END ALGEBRA; OPTION NDECIMALS=4 OPTION SUB=3176.157,1194 ! LIKELIHOOD AND DF FROM SATURATED MODEL: OZBMIYFSAT.MXO !DROP L 1 1 1 END Summary of VL file data for group 2 ARC1 ARC2 Code 1.0000 2.0000 Number 125.0000 125.0000 Mean -0.0090 -0.0190 Variance 0.8861 0.8757 Minimum -2.3713 -2.9867 Maximum 2.6268 2.2563 Summary of VL file data for group 3 PDE PCE PCD PBE PBD PBC PAE Code -10.0000 -9.0000 -8.0000 -7.0000 -6.0000 -5.0000 -4.0000 Number 392.0000 392.0000 392.0000 392.0000 392.0000 392.0000 392.0000 Mean 0.4975 0.5000 0.5031 0.5013 0.4996 0.5167 0.5013 Variance 0.0013 0.0013 0.0080 0.0006 0.0086 0.0286 0.0006 Minimum 0.0010 0.0010 0.0070 0.4970 0.0010 0.0010 0.4990 Maximum 0.5020 0.9960 0.9980 0.9960 0.9980 0.9980 0.9980 PAD PAC PAB ARC1 ARC2 ARC3 ARC4 Code -3.0000 -2.0000 -1.0000 1.0000 2.0000 3.0000 4.0000 Number 392.0000 392.0000 392.0000 392.0000 392.0000 132.0000 29.0000 Mean 0.5010 0.5072 0.4948 0.0844 0.0470 0.2260 0.2899 Variance 0.0057 0.0325 0.0938 0.9724 1.0117 1.1402 1.0366 Minimum 0.0010 0.0010 0.0010 -2.0600 -2.2000 -1.9900 -1.8400 Maximum 0.9980 0.9980 0.9980 3.1560 3.4310 3.5340 2.7000 ARC5 Code 5.0000 Number 3.0000 Mean -0.0407 Variance 0.5965 Minimum -0.9500 Maximum 0.9380 PARAMETER SPECIFICATIONS 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 I This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a FULL matrix of order 1 by 1 1 QTL 4 MATRIX M This is a FULL matrix of order 1 by 1 1 1 5 MATRIX P This is a computed FULL matrix of order 1 by 5 It has no free parameters specified MATRIX Q This is a computed 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 5 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 = 5 MX PARAMETER ESTIMATES GROUP NUMBER: 4 TITLE G4: STANDARDIZATION MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.6029 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.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.2395 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX I This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX L This is a FULL matrix of order 1 by 1 1 QTL 0.3800 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0651 MATRIX P This is a computed FULL matrix of order 1 by 5 [=A|C|E|D|Q] 1 2 3 4 5 1 0.6029 0.0000 0.2395 0.0000 0.1444 MATRIX Q This is a computed FULL matrix of order 1 by 1 [=L*L'] 1 1 0.1444 MATRIX S This is a computed FULL matrix of order 1 by 5 [=P@V~] 1 2 3 4 5 1 0.6110 0.0000 0.2427 0.0000 0.1464 MATRIX V This is a computed FULL matrix of order 1 by 1 [=A+C+E+D+Q] 1 1 0.9868 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 DOM_GEN 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 ADD_GEN 0.7765 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.4893 Your model has 5 estimated parameters and 1198 Observed statistics -2 times log-likelihood of data >>> 3173.610 Degrees of freedom >>>>>>>>>>>>>>>> 1193 Akaike's Information Criterion >>>> 787.610 Bayesian Information Criterion >>>> -2140.152 Sample size Adjusted BIC >>>> -246.746 Deviance Information Criterion >>>> -1043.859 Submodel fit* >>>>>>>>>>> 3176.157 Submodel df* >>>>>>>>>>> 1194 Difference Chi-squared >>>>>>>> 2.547 Difference d.f. >>>>>>>>>>>>>>> 1 Probability >>>>>>>>>>>>>>>>>>>> 0.111 Akaike's Information Criterion > 0.547 * Submodel statistic supplied by user This problem used 0.0% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.14 Execution 0: 0: 0: 9.44 TOTAL 0: 0: 0: 9.58 Total number of warnings issued: 0 ______________________________________________________________________________ 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 L FULL NVAR 1 FREE ! QTL GENETIC PATH, Q H FULL 1 1 ! SCALAR, 0.5 I FULL 1 1 ! SCALAR, 0.25 M FULL 1 1 FREE END MATRICES; LABEL ROW X ADD_GEN LABEL ROW Y COMM_ENV LABEL ROW Z SPEC_ENV LABEL ROW W DOM_GEN LABEL ROW L QTL MATRIX H .5 MATRIX I .25 START .6 ALL ! START ALL FREE PARAMETERS AT .6 START 0 M 1 1 1 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 Q= L*L'; ! QTL VARIANCE, Q^2 END ALGEBRA; OPTION NO_OUTPUT END The following MX script lines were read for group 2 TITLE G2: MZ DATA DATA NINPUT=2 RECTANGULAR FILE=MZDAT.REC Rectangular continuous data read initiated Note: Maximum ordinal/rectangular record length is: 1000 NOTE: Rectangular file contained 125 records with data that contained a total of 250 observations LABELS ARC1 ARC2 SELECT ARC1 ARC2 ; BEGIN MATRICES = GROUP 1; ! COPY ALL MATRICES FROM GROUP 1, BOTH DECLARED AND COMPUTED NOTE: Selection yields 125 data vectors for analysis NOTE: Vectors contain a total of 250 observations MEANS M|M ; ! 1X2 VECTOR OF EXPECTED MEANS COVARIANCE ! 2X2 COVARIANCE MATRIX A+C+E+D+Q | A+C+D+Q _ A+C+D+Q | A+C+E+D+Q ; OPTION RSIDUALS; OPTION NO_OUTPUT END The following MX script lines were read for group 3 TITLE G3: DZ DATA DATA NINPUT=15 RECTANGULAR FILE=GRID44.REC Rectangular continuous data read initiated Note: Maximum ordinal/rectangular record length is: 1000 NOTE: Rectangular file contained 392 records with data that contained a total of 4868 observations LABELS PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ARC1 ARC2 ARC3 ARC4 ARC5 SELECT PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ARC1 ARC2 ARC3 ARC4 ARC5 ; DEFINITON PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ; NOTE: Selection yields 392 data vectors for analysis NOTE: Vectors contain a total of 4868 observations Note: Global variable previously defined. Updating PAB Note: Global variable previously defined. Updating PAC Note: Global variable previously defined. Updating PAD Note: Global variable previously defined. Updating PAE Note: Global variable previously defined. Updating PBC Note: Global variable previously defined. Updating PBD Note: Global variable previously defined. Updating PBE Note: Global variable previously defined. Updating PCD Note: Global variable previously defined. Updating PCE Note: Global variable previously defined. Updating PDE NOTE: Definition yields 392 data vectors for analysis NOTE: Vectors contain a total of 948 observations BEGIN MATRICES ; A COMPUTED NVAR NVAR =A1 C COMPUTED NVAR NVAR =C1 D COMPUTED NVAR NVAR =D1 E COMPUTED NVAR NVAR =E1 Q COMPUTED NVAR NVAR =Q1 M FULL 1 1 =M1 H FULL 1 1 =H1 I FULL 1 1 =I1 P FULL 1 1 R FULL 1 1 S FULL 1 1 T FULL 1 1 U FULL 1 1 V FULL 1 1 W FULL 1 1 X FULL 1 1 Y FULL 1 1 Z FULL 1 1 END MATRICES; SPECIFY P PAB SPECIFY R PAC SPECIFY S PAD SPECIFY T PAE SPECIFY U PBC SPECIFY V PBD SPECIFY W PBE SPECIFY X PCD SPECIFY Y PCE SPECIFY Z PDE MEANS M|M|M|M|M ; COVARIANCE A+C+E+D+Q | H@A+C+I@D+P@Q | H@A+C+I@D+R@Q | H@A+C+I@D+S@Q | H@A+C+I@D+T@Q _ H@A+C+I@D+P@Q | A+C+E+D+Q | H@A+C+I@D+U@Q | H@A+C+I@D+V@Q | H@A+C+I@D+W@Q _ H@A+C+I@D+R@Q | H@A+C+I@D+U@Q | A+C+E+D+Q | H@A+C+I@D+X@Q | H@A+C+I@D+Y@Q _ H@A+C+I@D+S@Q | H@A+C+I@D+V@Q | H@A+C+I@D+X@Q | A+C+E+D+Q | H@A+C+I@D+Z@Q _ H@A+C+I@D+T@Q | H@A+C+I@D+W@Q | H@A+C+I@D+Y@Q | H@A+C+I@D+Z@Q | A+C+E+D+Q ; OPTION NO_OUTPUT END The following MX script lines were read for group 4 TITLE G4: STANDARDIZATION CALCULATION BEGIN MATRICES = GROUP 1; END MATRICES; BEGIN ALGEBRA; V=A+C+E+D+Q; ! TOTAL VARIANCE P=A|C|E|D|Q; ! CONCATENATE PARAMETER ESTIMATES S=P@V~; ! STANDARDIZED PARAMETER ESTIMATES END ALGEBRA; OPTION NDECIMALS=4 OPTION SUB=3176.157,1194 ! LIKELIHOOD AND DF FROM SATURATED MODEL: OZBMIYFSAT.MXO !DROP L 1 1 1 END Summary of VL file data for group 2 ARC1 ARC2 Code 1.0000 2.0000 Number 125.0000 125.0000 Mean -0.0090 -0.0190 Variance 0.8861 0.8757 Minimum -2.3713 -2.9867 Maximum 2.6268 2.2563 Summary of VL file data for group 3 PDE PCE PCD PBE PBD PBC PAE Code -10.0000 -9.0000 -8.0000 -7.0000 -6.0000 -5.0000 -4.0000 Number 392.0000 392.0000 392.0000 392.0000 392.0000 392.0000 392.0000 Mean 0.4986 0.5008 0.5034 0.5021 0.5033 0.5113 0.5004 Variance 0.0006 0.0007 0.0074 0.0009 0.0085 0.0301 0.0001 Minimum 0.0040 0.3170 0.0030 0.5000 0.0000 0.0000 0.4960 Maximum 0.5000 0.9880 0.9980 0.9950 0.9990 0.9980 0.6790 PAD PAC PAB ARC1 ARC2 ARC3 ARC4 Code -3.0000 -2.0000 -1.0000 1.0000 2.0000 3.0000 4.0000 Number 392.0000 392.0000 392.0000 392.0000 392.0000 132.0000 29.0000 Mean 0.4990 0.5012 0.4937 0.0844 0.0470 0.2260 0.2899 Variance 0.0071 0.0311 0.0921 0.9724 1.0117 1.1402 1.0366 Minimum 0.0010 0.0000 0.0000 -2.0600 -2.2000 -1.9900 -1.8400 Maximum 0.9990 0.9990 0.9990 3.1560 3.4310 3.5340 2.7000 ARC5 Code 5.0000 Number 3.0000 Mean -0.0407 Variance 0.5965 Minimum -0.9500 Maximum 0.9380 PARAMETER SPECIFICATIONS 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 I This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a FULL matrix of order 1 by 1 1 QTL 4 MATRIX M This is a FULL matrix of order 1 by 1 1 1 5 MATRIX P This is a computed FULL matrix of order 1 by 5 It has no free parameters specified MATRIX Q This is a computed 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 5 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 = 5 *** 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: 4 TITLE G4: STANDARDIZATION MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.5838 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.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.2391 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX I This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX L This is a FULL matrix of order 1 by 1 1 QTL 0.4049 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0658 MATRIX P This is a computed FULL matrix of order 1 by 5 [=A|C|E|D|Q] 1 2 3 4 5 1 0.5838 0.0000 0.2391 0.0000 0.1639 MATRIX Q This is a computed FULL matrix of order 1 by 1 [=L*L'] 1 1 0.1639 MATRIX S This is a computed FULL matrix of order 1 by 5 [=P@V~] 1 2 3 4 5 1 0.5916 0.0000 0.2423 0.0000 0.1661 MATRIX V This is a computed FULL matrix of order 1 by 1 [=A+C+E+D+Q] 1 1 0.9867 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 DOM_GEN 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 ADD_GEN 0.7640 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.4889 *** WARNING! *** Minimization may not be successful. See above CODE GREEN - it probably was OK Your model has 5 estimated parameters and 1198 Observed statistics -2 times log-likelihood of data >>> 3173.026 Degrees of freedom >>>>>>>>>>>>>>>> 1193 Akaike's Information Criterion >>>> 787.026 Bayesian Information Criterion >>>> -2140.445 Sample size Adjusted BIC >>>> -247.039 Deviance Information Criterion >>>> -1044.151 Submodel fit* >>>>>>>>>>> 3176.157 Submodel df* >>>>>>>>>>> 1194 Difference Chi-squared >>>>>>>> 3.131 Difference d.f. >>>>>>>>>>>>>>> 1 Probability >>>>>>>>>>>>>>>>>>>> 0.077 Akaike's Information Criterion > 1.131 * Submodel statistic supplied by user This problem used 0.0% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.14 Execution 0: 0: 0: 8.95 TOTAL 0: 0: 0: 9.09 Total number of warnings issued: 2 ______________________________________________________________________________ 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 L FULL NVAR 1 FREE ! QTL GENETIC PATH, Q H FULL 1 1 ! SCALAR, 0.5 I FULL 1 1 ! SCALAR, 0.25 M FULL 1 1 FREE END MATRICES; LABEL ROW X ADD_GEN LABEL ROW Y COMM_ENV LABEL ROW Z SPEC_ENV LABEL ROW W DOM_GEN LABEL ROW L QTL MATRIX H .5 MATRIX I .25 START .6 ALL ! START ALL FREE PARAMETERS AT .6 START 0 M 1 1 1 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 Q= L*L'; ! QTL VARIANCE, Q^2 END ALGEBRA; OPTION NO_OUTPUT END The following MX script lines were read for group 2 TITLE G2: MZ DATA DATA NINPUT=2 RECTANGULAR FILE=MZDAT.REC Rectangular continuous data read initiated Note: Maximum ordinal/rectangular record length is: 1000 NOTE: Rectangular file contained 125 records with data that contained a total of 250 observations LABELS ARC1 ARC2 SELECT ARC1 ARC2 ; BEGIN MATRICES = GROUP 1; ! COPY ALL MATRICES FROM GROUP 1, BOTH DECLARED AND COMPUTED NOTE: Selection yields 125 data vectors for analysis NOTE: Vectors contain a total of 250 observations MEANS M|M ; ! 1X2 VECTOR OF EXPECTED MEANS COVARIANCE ! 2X2 COVARIANCE MATRIX A+C+E+D+Q | A+C+D+Q _ A+C+D+Q | A+C+E+D+Q ; OPTION RSIDUALS; OPTION NO_OUTPUT END The following MX script lines were read for group 3 TITLE G3: DZ DATA DATA NINPUT=15 RECTANGULAR FILE=GRID45.REC Rectangular continuous data read initiated Note: Maximum ordinal/rectangular record length is: 1000 NOTE: Rectangular file contained 392 records with data that contained a total of 4868 observations LABELS PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ARC1 ARC2 ARC3 ARC4 ARC5 SELECT PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ARC1 ARC2 ARC3 ARC4 ARC5 ; DEFINITON PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ; NOTE: Selection yields 392 data vectors for analysis NOTE: Vectors contain a total of 4868 observations Note: Global variable previously defined. Updating PAB Note: Global variable previously defined. Updating PAC Note: Global variable previously defined. Updating PAD Note: Global variable previously defined. Updating PAE Note: Global variable previously defined. Updating PBC Note: Global variable previously defined. Updating PBD Note: Global variable previously defined. Updating PBE Note: Global variable previously defined. Updating PCD Note: Global variable previously defined. Updating PCE Note: Global variable previously defined. Updating PDE NOTE: Definition yields 392 data vectors for analysis NOTE: Vectors contain a total of 948 observations BEGIN MATRICES ; A COMPUTED NVAR NVAR =A1 C COMPUTED NVAR NVAR =C1 D COMPUTED NVAR NVAR =D1 E COMPUTED NVAR NVAR =E1 Q COMPUTED NVAR NVAR =Q1 M FULL 1 1 =M1 H FULL 1 1 =H1 I FULL 1 1 =I1 P FULL 1 1 R FULL 1 1 S FULL 1 1 T FULL 1 1 U FULL 1 1 V FULL 1 1 W FULL 1 1 X FULL 1 1 Y FULL 1 1 Z FULL 1 1 END MATRICES; SPECIFY P PAB SPECIFY R PAC SPECIFY S PAD SPECIFY T PAE SPECIFY U PBC SPECIFY V PBD SPECIFY W PBE SPECIFY X PCD SPECIFY Y PCE SPECIFY Z PDE MEANS M|M|M|M|M ; COVARIANCE A+C+E+D+Q | H@A+C+I@D+P@Q | H@A+C+I@D+R@Q | H@A+C+I@D+S@Q | H@A+C+I@D+T@Q _ H@A+C+I@D+P@Q | A+C+E+D+Q | H@A+C+I@D+U@Q | H@A+C+I@D+V@Q | H@A+C+I@D+W@Q _ H@A+C+I@D+R@Q | H@A+C+I@D+U@Q | A+C+E+D+Q | H@A+C+I@D+X@Q | H@A+C+I@D+Y@Q _ H@A+C+I@D+S@Q | H@A+C+I@D+V@Q | H@A+C+I@D+X@Q | A+C+E+D+Q | H@A+C+I@D+Z@Q _ H@A+C+I@D+T@Q | H@A+C+I@D+W@Q | H@A+C+I@D+Y@Q | H@A+C+I@D+Z@Q | A+C+E+D+Q ; OPTION NO_OUTPUT END The following MX script lines were read for group 4 TITLE G4: STANDARDIZATION CALCULATION BEGIN MATRICES = GROUP 1; END MATRICES; BEGIN ALGEBRA; V=A+C+E+D+Q; ! TOTAL VARIANCE P=A|C|E|D|Q; ! CONCATENATE PARAMETER ESTIMATES S=P@V~; ! STANDARDIZED PARAMETER ESTIMATES END ALGEBRA; OPTION NDECIMALS=4 OPTION SUB=3176.157,1194 ! LIKELIHOOD AND DF FROM SATURATED MODEL: OZBMIYFSAT.MXO !DROP L 1 1 1 END Summary of VL file data for group 2 ARC1 ARC2 Code 1.0000 2.0000 Number 125.0000 125.0000 Mean -0.0090 -0.0190 Variance 0.8861 0.8757 Minimum -2.3713 -2.9867 Maximum 2.6268 2.2563 Summary of VL file data for group 3 PDE PCE PCD PBE PBD PBC PAE Code -10.0000 -9.0000 -8.0000 -7.0000 -6.0000 -5.0000 -4.0000 Number 392.0000 392.0000 392.0000 392.0000 392.0000 392.0000 392.0000 Mean 0.4997 0.5012 0.5021 0.5026 0.5061 0.5059 0.5013 Variance 0.0011 0.0006 0.0065 0.0013 0.0085 0.0308 0.0004 Minimum 0.0020 0.4950 0.0010 0.5000 0.0010 0.0010 0.4860 Maximum 0.9080 0.9860 0.9980 0.9970 0.9980 0.9980 0.8320 PAD PAC PAB ARC1 ARC2 ARC3 ARC4 Code -3.0000 -2.0000 -1.0000 1.0000 2.0000 3.0000 4.0000 Number 392.0000 392.0000 392.0000 392.0000 392.0000 132.0000 29.0000 Mean 0.5030 0.5001 0.4870 0.0844 0.0470 0.2260 0.2899 Variance 0.0075 0.0309 0.0952 0.9724 1.0117 1.1402 1.0366 Minimum 0.0040 0.0010 0.0010 -2.0600 -2.2000 -1.9900 -1.8400 Maximum 0.9980 0.9980 0.9980 3.1560 3.4310 3.5340 2.7000 ARC5 Code 5.0000 Number 3.0000 Mean -0.0407 Variance 0.5965 Minimum -0.9500 Maximum 0.9380 PARAMETER SPECIFICATIONS 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 I This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a FULL matrix of order 1 by 1 1 QTL 4 MATRIX M This is a FULL matrix of order 1 by 1 1 1 5 MATRIX P This is a computed FULL matrix of order 1 by 5 It has no free parameters specified MATRIX Q This is a computed 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 5 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 = 5 MX PARAMETER ESTIMATES GROUP NUMBER: 4 TITLE G4: STANDARDIZATION MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.5879 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.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.2383 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX I This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX L This is a FULL matrix of order 1 by 1 1 QTL 0.4006 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0655 MATRIX P This is a computed FULL matrix of order 1 by 5 [=A|C|E|D|Q] 1 2 3 4 5 1 0.5879 0.0000 0.2383 0.0000 0.1605 MATRIX Q This is a computed FULL matrix of order 1 by 1 [=L*L'] 1 1 0.1605 MATRIX S This is a computed FULL matrix of order 1 by 5 [=P@V~] 1 2 3 4 5 1 0.5958 0.0000 0.2415 0.0000 0.1627 MATRIX V This is a computed FULL matrix of order 1 by 1 [=A+C+E+D+Q] 1 1 0.9867 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 DOM_GEN 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 ADD_GEN 0.7668 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.4882 Your model has 5 estimated parameters and 1198 Observed statistics -2 times log-likelihood of data >>> 3172.976 Degrees of freedom >>>>>>>>>>>>>>>> 1193 Akaike's Information Criterion >>>> 786.976 Bayesian Information Criterion >>>> -2140.470 Sample size Adjusted BIC >>>> -247.064 Deviance Information Criterion >>>> -1044.176 Submodel fit* >>>>>>>>>>> 3176.157 Submodel df* >>>>>>>>>>> 1194 Difference Chi-squared >>>>>>>> 3.181 Difference d.f. >>>>>>>>>>>>>>> 1 Probability >>>>>>>>>>>>>>>>>>>> 0.074 Akaike's Information Criterion > 1.181 * Submodel statistic supplied by user This problem used 0.0% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.16 Execution 0: 0: 0: 8.89 TOTAL 0: 0: 0: 9.05 Total number of warnings issued: 2 ______________________________________________________________________________ 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 L FULL NVAR 1 FREE ! QTL GENETIC PATH, Q H FULL 1 1 ! SCALAR, 0.5 I FULL 1 1 ! SCALAR, 0.25 M FULL 1 1 FREE END MATRICES; LABEL ROW X ADD_GEN LABEL ROW Y COMM_ENV LABEL ROW Z SPEC_ENV LABEL ROW W DOM_GEN LABEL ROW L QTL MATRIX H .5 MATRIX I .25 START .6 ALL ! START ALL FREE PARAMETERS AT .6 START 0 M 1 1 1 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 Q= L*L'; ! QTL VARIANCE, Q^2 END ALGEBRA; OPTION NO_OUTPUT END The following MX script lines were read for group 2 TITLE G2: MZ DATA DATA NINPUT=2 RECTANGULAR FILE=MZDAT.REC Rectangular continuous data read initiated Note: Maximum ordinal/rectangular record length is: 1000 NOTE: Rectangular file contained 125 records with data that contained a total of 250 observations LABELS ARC1 ARC2 SELECT ARC1 ARC2 ; BEGIN MATRICES = GROUP 1; ! COPY ALL MATRICES FROM GROUP 1, BOTH DECLARED AND COMPUTED NOTE: Selection yields 125 data vectors for analysis NOTE: Vectors contain a total of 250 observations MEANS M|M ; ! 1X2 VECTOR OF EXPECTED MEANS COVARIANCE ! 2X2 COVARIANCE MATRIX A+C+E+D+Q | A+C+D+Q _ A+C+D+Q | A+C+E+D+Q ; OPTION RSIDUALS; OPTION NO_OUTPUT END The following MX script lines were read for group 3 TITLE G3: DZ DATA DATA NINPUT=15 RECTANGULAR FILE=GRID46.REC Rectangular continuous data read initiated Note: Maximum ordinal/rectangular record length is: 1000 NOTE: Rectangular file contained 392 records with data that contained a total of 4868 observations LABELS PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ARC1 ARC2 ARC3 ARC4 ARC5 SELECT PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ARC1 ARC2 ARC3 ARC4 ARC5 ; DEFINITON PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ; NOTE: Selection yields 392 data vectors for analysis NOTE: Vectors contain a total of 4868 observations Note: Global variable previously defined. Updating PAB Note: Global variable previously defined. Updating PAC Note: Global variable previously defined. Updating PAD Note: Global variable previously defined. Updating PAE Note: Global variable previously defined. Updating PBC Note: Global variable previously defined. Updating PBD Note: Global variable previously defined. Updating PBE Note: Global variable previously defined. Updating PCD Note: Global variable previously defined. Updating PCE Note: Global variable previously defined. Updating PDE NOTE: Definition yields 392 data vectors for analysis NOTE: Vectors contain a total of 948 observations BEGIN MATRICES ; A COMPUTED NVAR NVAR =A1 C COMPUTED NVAR NVAR =C1 D COMPUTED NVAR NVAR =D1 E COMPUTED NVAR NVAR =E1 Q COMPUTED NVAR NVAR =Q1 M FULL 1 1 =M1 H FULL 1 1 =H1 I FULL 1 1 =I1 P FULL 1 1 R FULL 1 1 S FULL 1 1 T FULL 1 1 U FULL 1 1 V FULL 1 1 W FULL 1 1 X FULL 1 1 Y FULL 1 1 Z FULL 1 1 END MATRICES; SPECIFY P PAB SPECIFY R PAC SPECIFY S PAD SPECIFY T PAE SPECIFY U PBC SPECIFY V PBD SPECIFY W PBE SPECIFY X PCD SPECIFY Y PCE SPECIFY Z PDE MEANS M|M|M|M|M ; COVARIANCE A+C+E+D+Q | H@A+C+I@D+P@Q | H@A+C+I@D+R@Q | H@A+C+I@D+S@Q | H@A+C+I@D+T@Q _ H@A+C+I@D+P@Q | A+C+E+D+Q | H@A+C+I@D+U@Q | H@A+C+I@D+V@Q | H@A+C+I@D+W@Q _ H@A+C+I@D+R@Q | H@A+C+I@D+U@Q | A+C+E+D+Q | H@A+C+I@D+X@Q | H@A+C+I@D+Y@Q _ H@A+C+I@D+S@Q | H@A+C+I@D+V@Q | H@A+C+I@D+X@Q | A+C+E+D+Q | H@A+C+I@D+Z@Q _ H@A+C+I@D+T@Q | H@A+C+I@D+W@Q | H@A+C+I@D+Y@Q | H@A+C+I@D+Z@Q | A+C+E+D+Q ; OPTION NO_OUTPUT END The following MX script lines were read for group 4 TITLE G4: STANDARDIZATION CALCULATION BEGIN MATRICES = GROUP 1; END MATRICES; BEGIN ALGEBRA; V=A+C+E+D+Q; ! TOTAL VARIANCE P=A|C|E|D|Q; ! CONCATENATE PARAMETER ESTIMATES S=P@V~; ! STANDARDIZED PARAMETER ESTIMATES END ALGEBRA; OPTION NDECIMALS=4 OPTION SUB=3176.157,1194 ! LIKELIHOOD AND DF FROM SATURATED MODEL: OZBMIYFSAT.MXO !DROP L 1 1 1 END Summary of VL file data for group 2 ARC1 ARC2 Code 1.0000 2.0000 Number 125.0000 125.0000 Mean -0.0090 -0.0190 Variance 0.8861 0.8757 Minimum -2.3713 -2.9867 Maximum 2.6268 2.2563 Summary of VL file data for group 3 PDE PCE PCD PBE PBD PBC PAE Code -10.0000 -9.0000 -8.0000 -7.0000 -6.0000 -5.0000 -4.0000 Number 392.0000 392.0000 392.0000 392.0000 392.0000 392.0000 392.0000 Mean 0.5000 0.5012 0.4985 0.5026 0.5072 0.5030 0.5023 Variance 0.0013 0.0006 0.0066 0.0013 0.0098 0.0314 0.0011 Minimum 0.0000 0.4960 0.0000 0.5000 0.0000 0.0000 0.4760 Maximum 0.9950 0.9860 0.9980 0.9990 0.9990 0.9990 0.9970 PAD PAC PAB ARC1 ARC2 ARC3 ARC4 Code -3.0000 -2.0000 -1.0000 1.0000 2.0000 3.0000 4.0000 Number 392.0000 392.0000 392.0000 392.0000 392.0000 132.0000 29.0000 Mean 0.5031 0.5017 0.4734 0.0844 0.0470 0.2260 0.2899 Variance 0.0084 0.0293 0.0977 0.9724 1.0117 1.1402 1.0366 Minimum 0.0020 0.0000 0.0000 -2.0600 -2.2000 -1.9900 -1.8400 Maximum 0.9980 0.9990 0.9990 3.1560 3.4310 3.5340 2.7000 ARC5 Code 5.0000 Number 3.0000 Mean -0.0407 Variance 0.5965 Minimum -0.9500 Maximum 0.9380 PARAMETER SPECIFICATIONS 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 I This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a FULL matrix of order 1 by 1 1 QTL 4 MATRIX M This is a FULL matrix of order 1 by 1 1 1 5 MATRIX P This is a computed FULL matrix of order 1 by 5 It has no free parameters specified MATRIX Q This is a computed 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 5 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 = 5 MX PARAMETER ESTIMATES GROUP NUMBER: 4 TITLE G4: STANDARDIZATION MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.6055 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.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.2380 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX I This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX L This is a FULL matrix of order 1 by 1 1 QTL 0.3779 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0654 MATRIX P This is a computed FULL matrix of order 1 by 5 [=A|C|E|D|Q] 1 2 3 4 5 1 0.6055 0.0000 0.2380 0.0000 0.1428 MATRIX Q This is a computed FULL matrix of order 1 by 1 [=L*L'] 1 1 0.1428 MATRIX S This is a computed FULL matrix of order 1 by 5 [=P@V~] 1 2 3 4 5 1 0.6139 0.0000 0.2413 0.0000 0.1448 MATRIX V This is a computed FULL matrix of order 1 by 1 [=A+C+E+D+Q] 1 1 0.9863 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 DOM_GEN 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 ADD_GEN 0.7781 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.4878 Your model has 5 estimated parameters and 1198 Observed statistics -2 times log-likelihood of data >>> 3173.500 Degrees of freedom >>>>>>>>>>>>>>>> 1193 Akaike's Information Criterion >>>> 787.500 Bayesian Information Criterion >>>> -2140.208 Sample size Adjusted BIC >>>> -246.802 Deviance Information Criterion >>>> -1043.914 Submodel fit* >>>>>>>>>>> 3176.157 Submodel df* >>>>>>>>>>> 1194 Difference Chi-squared >>>>>>>> 2.657 Difference d.f. >>>>>>>>>>>>>>> 1 Probability >>>>>>>>>>>>>>>>>>>> 0.103 Akaike's Information Criterion > 0.657 * Submodel statistic supplied by user This problem used 0.0% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.14 Execution 0: 0: 0: 9.79 TOTAL 0: 0: 0: 9.93 Total number of warnings issued: 2 ______________________________________________________________________________ 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 L FULL NVAR 1 FREE ! QTL GENETIC PATH, Q H FULL 1 1 ! SCALAR, 0.5 I FULL 1 1 ! SCALAR, 0.25 M FULL 1 1 FREE END MATRICES; LABEL ROW X ADD_GEN LABEL ROW Y COMM_ENV LABEL ROW Z SPEC_ENV LABEL ROW W DOM_GEN LABEL ROW L QTL MATRIX H .5 MATRIX I .25 START .6 ALL ! START ALL FREE PARAMETERS AT .6 START 0 M 1 1 1 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 Q= L*L'; ! QTL VARIANCE, Q^2 END ALGEBRA; OPTION NO_OUTPUT END The following MX script lines were read for group 2 TITLE G2: MZ DATA DATA NINPUT=2 RECTANGULAR FILE=MZDAT.REC Rectangular continuous data read initiated Note: Maximum ordinal/rectangular record length is: 1000 NOTE: Rectangular file contained 125 records with data that contained a total of 250 observations LABELS ARC1 ARC2 SELECT ARC1 ARC2 ; BEGIN MATRICES = GROUP 1; ! COPY ALL MATRICES FROM GROUP 1, BOTH DECLARED AND COMPUTED NOTE: Selection yields 125 data vectors for analysis NOTE: Vectors contain a total of 250 observations MEANS M|M ; ! 1X2 VECTOR OF EXPECTED MEANS COVARIANCE ! 2X2 COVARIANCE MATRIX A+C+E+D+Q | A+C+D+Q _ A+C+D+Q | A+C+E+D+Q ; OPTION RSIDUALS; OPTION NO_OUTPUT END The following MX script lines were read for group 3 TITLE G3: DZ DATA DATA NINPUT=15 RECTANGULAR FILE=GRID47.REC Rectangular continuous data read initiated Note: Maximum ordinal/rectangular record length is: 1000 NOTE: Rectangular file contained 392 records with data that contained a total of 4868 observations LABELS PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ARC1 ARC2 ARC3 ARC4 ARC5 SELECT PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ARC1 ARC2 ARC3 ARC4 ARC5 ; DEFINITON PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ; NOTE: Selection yields 392 data vectors for analysis NOTE: Vectors contain a total of 4868 observations Note: Global variable previously defined. Updating PAB Note: Global variable previously defined. Updating PAC Note: Global variable previously defined. Updating PAD Note: Global variable previously defined. Updating PAE Note: Global variable previously defined. Updating PBC Note: Global variable previously defined. Updating PBD Note: Global variable previously defined. Updating PBE Note: Global variable previously defined. Updating PCD Note: Global variable previously defined. Updating PCE Note: Global variable previously defined. Updating PDE NOTE: Definition yields 392 data vectors for analysis NOTE: Vectors contain a total of 948 observations BEGIN MATRICES ; A COMPUTED NVAR NVAR =A1 C COMPUTED NVAR NVAR =C1 D COMPUTED NVAR NVAR =D1 E COMPUTED NVAR NVAR =E1 Q COMPUTED NVAR NVAR =Q1 M FULL 1 1 =M1 H FULL 1 1 =H1 I FULL 1 1 =I1 P FULL 1 1 R FULL 1 1 S FULL 1 1 T FULL 1 1 U FULL 1 1 V FULL 1 1 W FULL 1 1 X FULL 1 1 Y FULL 1 1 Z FULL 1 1 END MATRICES; SPECIFY P PAB SPECIFY R PAC SPECIFY S PAD SPECIFY T PAE SPECIFY U PBC SPECIFY V PBD SPECIFY W PBE SPECIFY X PCD SPECIFY Y PCE SPECIFY Z PDE MEANS M|M|M|M|M ; COVARIANCE A+C+E+D+Q | H@A+C+I@D+P@Q | H@A+C+I@D+R@Q | H@A+C+I@D+S@Q | H@A+C+I@D+T@Q _ H@A+C+I@D+P@Q | A+C+E+D+Q | H@A+C+I@D+U@Q | H@A+C+I@D+V@Q | H@A+C+I@D+W@Q _ H@A+C+I@D+R@Q | H@A+C+I@D+U@Q | A+C+E+D+Q | H@A+C+I@D+X@Q | H@A+C+I@D+Y@Q _ H@A+C+I@D+S@Q | H@A+C+I@D+V@Q | H@A+C+I@D+X@Q | A+C+E+D+Q | H@A+C+I@D+Z@Q _ H@A+C+I@D+T@Q | H@A+C+I@D+W@Q | H@A+C+I@D+Y@Q | H@A+C+I@D+Z@Q | A+C+E+D+Q ; OPTION NO_OUTPUT END The following MX script lines were read for group 4 TITLE G4: STANDARDIZATION CALCULATION BEGIN MATRICES = GROUP 1; END MATRICES; BEGIN ALGEBRA; V=A+C+E+D+Q; ! TOTAL VARIANCE P=A|C|E|D|Q; ! CONCATENATE PARAMETER ESTIMATES S=P@V~; ! STANDARDIZED PARAMETER ESTIMATES END ALGEBRA; OPTION NDECIMALS=4 OPTION SUB=3176.157,1194 ! LIKELIHOOD AND DF FROM SATURATED MODEL: OZBMIYFSAT.MXO !DROP L 1 1 1 END Summary of VL file data for group 2 ARC1 ARC2 Code 1.0000 2.0000 Number 125.0000 125.0000 Mean -0.0090 -0.0190 Variance 0.8861 0.8757 Minimum -2.3713 -2.9867 Maximum 2.6268 2.2563 Summary of VL file data for group 3 PDE PCE PCD PBE PBD PBC PAE Code -10.0000 -9.0000 -8.0000 -7.0000 -6.0000 -5.0000 -4.0000 Number 392.0000 392.0000 392.0000 392.0000 392.0000 392.0000 392.0000 Mean 0.5000 0.5013 0.4991 0.5024 0.5077 0.4994 0.5015 Variance 0.0013 0.0006 0.0075 0.0013 0.0108 0.0354 0.0007 Minimum 0.0010 0.5000 0.0000 0.4300 0.0000 0.0000 0.4890 Maximum 0.9990 0.9930 0.9990 0.9990 0.9990 0.9990 0.9990 PAD PAC PAB ARC1 ARC2 ARC3 ARC4 Code -3.0000 -2.0000 -1.0000 1.0000 2.0000 3.0000 4.0000 Number 392.0000 392.0000 392.0000 392.0000 392.0000 132.0000 29.0000 Mean 0.5011 0.5035 0.4745 0.0844 0.0470 0.2260 0.2899 Variance 0.0095 0.0339 0.1080 0.9724 1.0117 1.1402 1.0366 Minimum 0.0000 0.0000 0.0000 -2.0600 -2.2000 -1.9900 -1.8400 Maximum 0.9990 0.9990 0.9990 3.1560 3.4310 3.5340 2.7000 ARC5 Code 5.0000 Number 3.0000 Mean -0.0407 Variance 0.5965 Minimum -0.9500 Maximum 0.9380 PARAMETER SPECIFICATIONS 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 I This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a FULL matrix of order 1 by 1 1 QTL 4 MATRIX M This is a FULL matrix of order 1 by 1 1 1 5 MATRIX P This is a computed FULL matrix of order 1 by 5 It has no free parameters specified MATRIX Q This is a computed 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 5 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 = 5 MX PARAMETER ESTIMATES GROUP NUMBER: 4 TITLE G4: STANDARDIZATION MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.6239 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.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.2377 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX I This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX L This is a FULL matrix of order 1 by 1 1 QTL 0.3529 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0657 MATRIX P This is a computed FULL matrix of order 1 by 5 [=A|C|E|D|Q] 1 2 3 4 5 1 0.6239 0.0000 0.2377 0.0000 0.1246 MATRIX Q This is a computed FULL matrix of order 1 by 1 [=L*L'] 1 1 0.1246 MATRIX S This is a computed FULL matrix of order 1 by 5 [=P@V~] 1 2 3 4 5 1 0.6327 0.0000 0.2410 0.0000 0.1263 MATRIX V This is a computed FULL matrix of order 1 by 1 [=A+C+E+D+Q] 1 1 0.9862 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 DOM_GEN 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 ADD_GEN 0.7899 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.4876 Your model has 5 estimated parameters and 1198 Observed statistics -2 times log-likelihood of data >>> 3173.628 Degrees of freedom >>>>>>>>>>>>>>>> 1193 Akaike's Information Criterion >>>> 787.628 Bayesian Information Criterion >>>> -2140.144 Sample size Adjusted BIC >>>> -246.738 Deviance Information Criterion >>>> -1043.850 Submodel fit* >>>>>>>>>>> 3176.157 Submodel df* >>>>>>>>>>> 1194 Difference Chi-squared >>>>>>>> 2.529 Difference d.f. >>>>>>>>>>>>>>> 1 Probability >>>>>>>>>>>>>>>>>>>> 0.112 Akaike's Information Criterion > 0.529 * Submodel statistic supplied by user This problem used 0.0% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.11 Execution 0: 0: 0:10.27 TOTAL 0: 0: 0:10.38 Total number of warnings issued: 2 ______________________________________________________________________________ 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 L FULL NVAR 1 FREE ! QTL GENETIC PATH, Q H FULL 1 1 ! SCALAR, 0.5 I FULL 1 1 ! SCALAR, 0.25 M FULL 1 1 FREE END MATRICES; LABEL ROW X ADD_GEN LABEL ROW Y COMM_ENV LABEL ROW Z SPEC_ENV LABEL ROW W DOM_GEN LABEL ROW L QTL MATRIX H .5 MATRIX I .25 START .6 ALL ! START ALL FREE PARAMETERS AT .6 START 0 M 1 1 1 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 Q= L*L'; ! QTL VARIANCE, Q^2 END ALGEBRA; OPTION NO_OUTPUT END The following MX script lines were read for group 2 TITLE G2: MZ DATA DATA NINPUT=2 RECTANGULAR FILE=MZDAT.REC Rectangular continuous data read initiated Note: Maximum ordinal/rectangular record length is: 1000 NOTE: Rectangular file contained 125 records with data that contained a total of 250 observations LABELS ARC1 ARC2 SELECT ARC1 ARC2 ; BEGIN MATRICES = GROUP 1; ! COPY ALL MATRICES FROM GROUP 1, BOTH DECLARED AND COMPUTED NOTE: Selection yields 125 data vectors for analysis NOTE: Vectors contain a total of 250 observations MEANS M|M ; ! 1X2 VECTOR OF EXPECTED MEANS COVARIANCE ! 2X2 COVARIANCE MATRIX A+C+E+D+Q | A+C+D+Q _ A+C+D+Q | A+C+E+D+Q ; OPTION RSIDUALS; OPTION NO_OUTPUT END The following MX script lines were read for group 3 TITLE G3: DZ DATA DATA NINPUT=15 RECTANGULAR FILE=GRID48.REC Rectangular continuous data read initiated Note: Maximum ordinal/rectangular record length is: 1000 NOTE: Rectangular file contained 392 records with data that contained a total of 4868 observations LABELS PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ARC1 ARC2 ARC3 ARC4 ARC5 SELECT PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ARC1 ARC2 ARC3 ARC4 ARC5 ; DEFINITON PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ; NOTE: Selection yields 392 data vectors for analysis NOTE: Vectors contain a total of 4868 observations Note: Global variable previously defined. Updating PAB Note: Global variable previously defined. Updating PAC Note: Global variable previously defined. Updating PAD Note: Global variable previously defined. Updating PAE Note: Global variable previously defined. Updating PBC Note: Global variable previously defined. Updating PBD Note: Global variable previously defined. Updating PBE Note: Global variable previously defined. Updating PCD Note: Global variable previously defined. Updating PCE Note: Global variable previously defined. Updating PDE NOTE: Definition yields 392 data vectors for analysis NOTE: Vectors contain a total of 948 observations BEGIN MATRICES ; A COMPUTED NVAR NVAR =A1 C COMPUTED NVAR NVAR =C1 D COMPUTED NVAR NVAR =D1 E COMPUTED NVAR NVAR =E1 Q COMPUTED NVAR NVAR =Q1 M FULL 1 1 =M1 H FULL 1 1 =H1 I FULL 1 1 =I1 P FULL 1 1 R FULL 1 1 S FULL 1 1 T FULL 1 1 U FULL 1 1 V FULL 1 1 W FULL 1 1 X FULL 1 1 Y FULL 1 1 Z FULL 1 1 END MATRICES; SPECIFY P PAB SPECIFY R PAC SPECIFY S PAD SPECIFY T PAE SPECIFY U PBC SPECIFY V PBD SPECIFY W PBE SPECIFY X PCD SPECIFY Y PCE SPECIFY Z PDE MEANS M|M|M|M|M ; COVARIANCE A+C+E+D+Q | H@A+C+I@D+P@Q | H@A+C+I@D+R@Q | H@A+C+I@D+S@Q | H@A+C+I@D+T@Q _ H@A+C+I@D+P@Q | A+C+E+D+Q | H@A+C+I@D+U@Q | H@A+C+I@D+V@Q | H@A+C+I@D+W@Q _ H@A+C+I@D+R@Q | H@A+C+I@D+U@Q | A+C+E+D+Q | H@A+C+I@D+X@Q | H@A+C+I@D+Y@Q _ H@A+C+I@D+S@Q | H@A+C+I@D+V@Q | H@A+C+I@D+X@Q | A+C+E+D+Q | H@A+C+I@D+Z@Q _ H@A+C+I@D+T@Q | H@A+C+I@D+W@Q | H@A+C+I@D+Y@Q | H@A+C+I@D+Z@Q | A+C+E+D+Q ; OPTION NO_OUTPUT END The following MX script lines were read for group 4 TITLE G4: STANDARDIZATION CALCULATION BEGIN MATRICES = GROUP 1; END MATRICES; BEGIN ALGEBRA; V=A+C+E+D+Q; ! TOTAL VARIANCE P=A|C|E|D|Q; ! CONCATENATE PARAMETER ESTIMATES S=P@V~; ! STANDARDIZED PARAMETER ESTIMATES END ALGEBRA; OPTION NDECIMALS=4 OPTION SUB=3176.157,1194 ! LIKELIHOOD AND DF FROM SATURATED MODEL: OZBMIYFSAT.MXO !DROP L 1 1 1 END Summary of VL file data for group 2 ARC1 ARC2 Code 1.0000 2.0000 Number 125.0000 125.0000 Mean -0.0090 -0.0190 Variance 0.8861 0.8757 Minimum -2.3713 -2.9867 Maximum 2.6268 2.2563 Summary of VL file data for group 3 PDE PCE PCD PBE PBD PBC PAE Code -10.0000 -9.0000 -8.0000 -7.0000 -6.0000 -5.0000 -4.0000 Number 392.0000 392.0000 392.0000 392.0000 392.0000 392.0000 392.0000 Mean 0.4999 0.5025 0.4978 0.5025 0.5024 0.5034 0.5000 Variance 0.0012 0.0012 0.0099 0.0013 0.0089 0.0391 0.0013 Minimum 0.0000 0.5000 0.0000 0.4840 0.0000 0.0000 0.0000 Maximum 0.9770 0.9990 0.9990 0.9990 0.9990 0.9990 0.9940 PAD PAC PAB ARC1 ARC2 ARC3 ARC4 Code -3.0000 -2.0000 -1.0000 1.0000 2.0000 3.0000 4.0000 Number 392.0000 392.0000 392.0000 392.0000 392.0000 132.0000 29.0000 Mean 0.5005 0.5030 0.4894 0.0844 0.0470 0.2260 0.2899 Variance 0.0102 0.0358 0.1041 0.9724 1.0117 1.1402 1.0366 Minimum 0.0000 0.0000 0.0000 -2.0600 -2.2000 -1.9900 -1.8400 Maximum 0.9990 0.9990 0.9990 3.1560 3.4310 3.5340 2.7000 ARC5 Code 5.0000 Number 3.0000 Mean -0.0407 Variance 0.5965 Minimum -0.9500 Maximum 0.9380 PARAMETER SPECIFICATIONS 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 I This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a FULL matrix of order 1 by 1 1 QTL 4 MATRIX M This is a FULL matrix of order 1 by 1 1 1 5 MATRIX P This is a computed FULL matrix of order 1 by 5 It has no free parameters specified MATRIX Q This is a computed 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 5 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 = 5 MX PARAMETER ESTIMATES GROUP NUMBER: 4 TITLE G4: STANDARDIZATION MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.5608 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.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.2379 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX I This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX L This is a FULL matrix of order 1 by 1 1 QTL 0.4318 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0644 MATRIX P This is a computed FULL matrix of order 1 by 5 [=A|C|E|D|Q] 1 2 3 4 5 1 0.5608 0.0000 0.2379 0.0000 0.1864 MATRIX Q This is a computed FULL matrix of order 1 by 1 [=L*L'] 1 1 0.1864 MATRIX S This is a computed FULL matrix of order 1 by 5 [=P@V~] 1 2 3 4 5 1 0.5693 0.0000 0.2415 0.0000 0.1892 MATRIX V This is a computed FULL matrix of order 1 by 1 [=A+C+E+D+Q] 1 1 0.9852 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 DOM_GEN 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 ADD_GEN 0.7489 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.4878 Your model has 5 estimated parameters and 1198 Observed statistics -2 times log-likelihood of data >>> 3169.467 Degrees of freedom >>>>>>>>>>>>>>>> 1193 Akaike's Information Criterion >>>> 783.467 Bayesian Information Criterion >>>> -2142.224 Sample size Adjusted BIC >>>> -248.818 Deviance Information Criterion >>>> -1045.930 Submodel fit* >>>>>>>>>>> 3176.157 Submodel df* >>>>>>>>>>> 1194 Difference Chi-squared >>>>>>>> 6.690 Difference d.f. >>>>>>>>>>>>>>> 1 Probability >>>>>>>>>>>>>>>>>>>> 0.010 Akaike's Information Criterion > 4.690 * Submodel statistic supplied by user This problem used 0.0% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.14 Execution 0: 0: 0: 8.55 TOTAL 0: 0: 0: 8.69 Total number of warnings issued: 2 ______________________________________________________________________________ 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 L FULL NVAR 1 FREE ! QTL GENETIC PATH, Q H FULL 1 1 ! SCALAR, 0.5 I FULL 1 1 ! SCALAR, 0.25 M FULL 1 1 FREE END MATRICES; LABEL ROW X ADD_GEN LABEL ROW Y COMM_ENV LABEL ROW Z SPEC_ENV LABEL ROW W DOM_GEN LABEL ROW L QTL MATRIX H .5 MATRIX I .25 START .6 ALL ! START ALL FREE PARAMETERS AT .6 START 0 M 1 1 1 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 Q= L*L'; ! QTL VARIANCE, Q^2 END ALGEBRA; OPTION NO_OUTPUT END The following MX script lines were read for group 2 TITLE G2: MZ DATA DATA NINPUT=2 RECTANGULAR FILE=MZDAT.REC Rectangular continuous data read initiated Note: Maximum ordinal/rectangular record length is: 1000 NOTE: Rectangular file contained 125 records with data that contained a total of 250 observations LABELS ARC1 ARC2 SELECT ARC1 ARC2 ; BEGIN MATRICES = GROUP 1; ! COPY ALL MATRICES FROM GROUP 1, BOTH DECLARED AND COMPUTED NOTE: Selection yields 125 data vectors for analysis NOTE: Vectors contain a total of 250 observations MEANS M|M ; ! 1X2 VECTOR OF EXPECTED MEANS COVARIANCE ! 2X2 COVARIANCE MATRIX A+C+E+D+Q | A+C+D+Q _ A+C+D+Q | A+C+E+D+Q ; OPTION RSIDUALS; OPTION NO_OUTPUT END The following MX script lines were read for group 3 TITLE G3: DZ DATA DATA NINPUT=15 RECTANGULAR FILE=GRID49.REC Rectangular continuous data read initiated Note: Maximum ordinal/rectangular record length is: 1000 NOTE: Rectangular file contained 392 records with data that contained a total of 4868 observations LABELS PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ARC1 ARC2 ARC3 ARC4 ARC5 SELECT PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ARC1 ARC2 ARC3 ARC4 ARC5 ; DEFINITON PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ; NOTE: Selection yields 392 data vectors for analysis NOTE: Vectors contain a total of 4868 observations Note: Global variable previously defined. Updating PAB Note: Global variable previously defined. Updating PAC Note: Global variable previously defined. Updating PAD Note: Global variable previously defined. Updating PAE Note: Global variable previously defined. Updating PBC Note: Global variable previously defined. Updating PBD Note: Global variable previously defined. Updating PBE Note: Global variable previously defined. Updating PCD Note: Global variable previously defined. Updating PCE Note: Global variable previously defined. Updating PDE NOTE: Definition yields 392 data vectors for analysis NOTE: Vectors contain a total of 948 observations BEGIN MATRICES ; A COMPUTED NVAR NVAR =A1 C COMPUTED NVAR NVAR =C1 D COMPUTED NVAR NVAR =D1 E COMPUTED NVAR NVAR =E1 Q COMPUTED NVAR NVAR =Q1 M FULL 1 1 =M1 H FULL 1 1 =H1 I FULL 1 1 =I1 P FULL 1 1 R FULL 1 1 S FULL 1 1 T FULL 1 1 U FULL 1 1 V FULL 1 1 W FULL 1 1 X FULL 1 1 Y FULL 1 1 Z FULL 1 1 END MATRICES; SPECIFY P PAB SPECIFY R PAC SPECIFY S PAD SPECIFY T PAE SPECIFY U PBC SPECIFY V PBD SPECIFY W PBE SPECIFY X PCD SPECIFY Y PCE SPECIFY Z PDE MEANS M|M|M|M|M ; COVARIANCE A+C+E+D+Q | H@A+C+I@D+P@Q | H@A+C+I@D+R@Q | H@A+C+I@D+S@Q | H@A+C+I@D+T@Q _ H@A+C+I@D+P@Q | A+C+E+D+Q | H@A+C+I@D+U@Q | H@A+C+I@D+V@Q | H@A+C+I@D+W@Q _ H@A+C+I@D+R@Q | H@A+C+I@D+U@Q | A+C+E+D+Q | H@A+C+I@D+X@Q | H@A+C+I@D+Y@Q _ H@A+C+I@D+S@Q | H@A+C+I@D+V@Q | H@A+C+I@D+X@Q | A+C+E+D+Q | H@A+C+I@D+Z@Q _ H@A+C+I@D+T@Q | H@A+C+I@D+W@Q | H@A+C+I@D+Y@Q | H@A+C+I@D+Z@Q | A+C+E+D+Q ; OPTION NO_OUTPUT END The following MX script lines were read for group 4 TITLE G4: STANDARDIZATION CALCULATION BEGIN MATRICES = GROUP 1; END MATRICES; BEGIN ALGEBRA; V=A+C+E+D+Q; ! TOTAL VARIANCE P=A|C|E|D|Q; ! CONCATENATE PARAMETER ESTIMATES S=P@V~; ! STANDARDIZED PARAMETER ESTIMATES END ALGEBRA; OPTION NDECIMALS=4 OPTION SUB=3176.157,1194 ! LIKELIHOOD AND DF FROM SATURATED MODEL: OZBMIYFSAT.MXO !DROP L 1 1 1 END Summary of VL file data for group 2 ARC1 ARC2 Code 1.0000 2.0000 Number 125.0000 125.0000 Mean -0.0090 -0.0190 Variance 0.8861 0.8757 Minimum -2.3713 -2.9867 Maximum 2.6268 2.2563 Summary of VL file data for group 3 PDE PCE PCD PBE PBD PBC PAE Code -10.0000 -9.0000 -8.0000 -7.0000 -6.0000 -5.0000 -4.0000 Number 392.0000 392.0000 392.0000 392.0000 392.0000 392.0000 392.0000 Mean 0.4994 0.5024 0.4938 0.5024 0.5017 0.5030 0.4995 Variance 0.0008 0.0012 0.0094 0.0012 0.0081 0.0339 0.0009 Minimum 0.0040 0.4820 0.0040 0.4880 0.0040 0.0040 0.0040 Maximum 0.7800 0.9950 0.9950 0.9950 0.9950 0.9950 0.8240 PAD PAC PAB ARC1 ARC2 ARC3 ARC4 Code -3.0000 -2.0000 -1.0000 1.0000 2.0000 3.0000 4.0000 Number 392.0000 392.0000 392.0000 392.0000 392.0000 132.0000 29.0000 Mean 0.5000 0.5021 0.4900 0.0844 0.0470 0.2260 0.2899 Variance 0.0094 0.0316 0.0901 0.9724 1.0117 1.1402 1.0366 Minimum 0.0040 0.0040 0.0040 -2.0600 -2.2000 -1.9900 -1.8400 Maximum 0.9950 0.9950 0.9930 3.1560 3.4310 3.5340 2.7000 ARC5 Code 5.0000 Number 3.0000 Mean -0.0407 Variance 0.5965 Minimum -0.9500 Maximum 0.9380 PARAMETER SPECIFICATIONS 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 I This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a FULL matrix of order 1 by 1 1 QTL 4 MATRIX M This is a FULL matrix of order 1 by 1 1 1 5 MATRIX P This is a computed FULL matrix of order 1 by 5 It has no free parameters specified MATRIX Q This is a computed 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 5 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 = 5 MX PARAMETER ESTIMATES GROUP NUMBER: 4 TITLE G4: STANDARDIZATION MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.5087 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.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.2376 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX I This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX L This is a FULL matrix of order 1 by 1 1 QTL 0.4881 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0633 MATRIX P This is a computed FULL matrix of order 1 by 5 [=A|C|E|D|Q] 1 2 3 4 5 1 0.5087 0.0000 0.2376 0.0000 0.2383 MATRIX Q This is a computed FULL matrix of order 1 by 1 [=L*L'] 1 1 0.2383 MATRIX S This is a computed FULL matrix of order 1 by 5 [=P@V~] 1 2 3 4 5 1 0.5167 0.0000 0.2413 0.0000 0.2420 MATRIX V This is a computed FULL matrix of order 1 by 1 [=A+C+E+D+Q] 1 1 0.9846 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 DOM_GEN 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 ADD_GEN 0.7132 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.4874 Your model has 5 estimated parameters and 1198 Observed statistics -2 times log-likelihood of data >>> 3166.774 Degrees of freedom >>>>>>>>>>>>>>>> 1193 Akaike's Information Criterion >>>> 780.774 Bayesian Information Criterion >>>> -2143.571 Sample size Adjusted BIC >>>> -250.165 Deviance Information Criterion >>>> -1047.277 Submodel fit* >>>>>>>>>>> 3176.157 Submodel df* >>>>>>>>>>> 1194 Difference Chi-squared >>>>>>>> 9.383 Difference d.f. >>>>>>>>>>>>>>> 1 Probability >>>>>>>>>>>>>>>>>>>> 0.002 Akaike's Information Criterion > 7.383 * Submodel statistic supplied by user This problem used 0.0% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.12 Execution 0: 0: 0: 9.32 TOTAL 0: 0: 0: 9.44 Total number of warnings issued: 2 ______________________________________________________________________________ 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 L FULL NVAR 1 FREE ! QTL GENETIC PATH, Q H FULL 1 1 ! SCALAR, 0.5 I FULL 1 1 ! SCALAR, 0.25 M FULL 1 1 FREE END MATRICES; LABEL ROW X ADD_GEN LABEL ROW Y COMM_ENV LABEL ROW Z SPEC_ENV LABEL ROW W DOM_GEN LABEL ROW L QTL MATRIX H .5 MATRIX I .25 START .6 ALL ! START ALL FREE PARAMETERS AT .6 START 0 M 1 1 1 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 Q= L*L'; ! QTL VARIANCE, Q^2 END ALGEBRA; OPTION NO_OUTPUT END The following MX script lines were read for group 2 TITLE G2: MZ DATA DATA NINPUT=2 RECTANGULAR FILE=MZDAT.REC Rectangular continuous data read initiated Note: Maximum ordinal/rectangular record length is: 1000 NOTE: Rectangular file contained 125 records with data that contained a total of 250 observations LABELS ARC1 ARC2 SELECT ARC1 ARC2 ; BEGIN MATRICES = GROUP 1; ! COPY ALL MATRICES FROM GROUP 1, BOTH DECLARED AND COMPUTED NOTE: Selection yields 125 data vectors for analysis NOTE: Vectors contain a total of 250 observations MEANS M|M ; ! 1X2 VECTOR OF EXPECTED MEANS COVARIANCE ! 2X2 COVARIANCE MATRIX A+C+E+D+Q | A+C+D+Q _ A+C+D+Q | A+C+E+D+Q ; OPTION RSIDUALS; OPTION NO_OUTPUT END The following MX script lines were read for group 3 TITLE G3: DZ DATA DATA NINPUT=15 RECTANGULAR FILE=GRID50.REC Rectangular continuous data read initiated Note: Maximum ordinal/rectangular record length is: 1000 NOTE: Rectangular file contained 392 records with data that contained a total of 4868 observations LABELS PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ARC1 ARC2 ARC3 ARC4 ARC5 SELECT PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ARC1 ARC2 ARC3 ARC4 ARC5 ; DEFINITON PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ; NOTE: Selection yields 392 data vectors for analysis NOTE: Vectors contain a total of 4868 observations Note: Global variable previously defined. Updating PAB Note: Global variable previously defined. Updating PAC Note: Global variable previously defined. Updating PAD Note: Global variable previously defined. Updating PAE Note: Global variable previously defined. Updating PBC Note: Global variable previously defined. Updating PBD Note: Global variable previously defined. Updating PBE Note: Global variable previously defined. Updating PCD Note: Global variable previously defined. Updating PCE Note: Global variable previously defined. Updating PDE NOTE: Definition yields 392 data vectors for analysis NOTE: Vectors contain a total of 948 observations BEGIN MATRICES ; A COMPUTED NVAR NVAR =A1 C COMPUTED NVAR NVAR =C1 D COMPUTED NVAR NVAR =D1 E COMPUTED NVAR NVAR =E1 Q COMPUTED NVAR NVAR =Q1 M FULL 1 1 =M1 H FULL 1 1 =H1 I FULL 1 1 =I1 P FULL 1 1 R FULL 1 1 S FULL 1 1 T FULL 1 1 U FULL 1 1 V FULL 1 1 W FULL 1 1 X FULL 1 1 Y FULL 1 1 Z FULL 1 1 END MATRICES; SPECIFY P PAB SPECIFY R PAC SPECIFY S PAD SPECIFY T PAE SPECIFY U PBC SPECIFY V PBD SPECIFY W PBE SPECIFY X PCD SPECIFY Y PCE SPECIFY Z PDE MEANS M|M|M|M|M ; COVARIANCE A+C+E+D+Q | H@A+C+I@D+P@Q | H@A+C+I@D+R@Q | H@A+C+I@D+S@Q | H@A+C+I@D+T@Q _ H@A+C+I@D+P@Q | A+C+E+D+Q | H@A+C+I@D+U@Q | H@A+C+I@D+V@Q | H@A+C+I@D+W@Q _ H@A+C+I@D+R@Q | H@A+C+I@D+U@Q | A+C+E+D+Q | H@A+C+I@D+X@Q | H@A+C+I@D+Y@Q _ H@A+C+I@D+S@Q | H@A+C+I@D+V@Q | H@A+C+I@D+X@Q | A+C+E+D+Q | H@A+C+I@D+Z@Q _ H@A+C+I@D+T@Q | H@A+C+I@D+W@Q | H@A+C+I@D+Y@Q | H@A+C+I@D+Z@Q | A+C+E+D+Q ; OPTION NO_OUTPUT END The following MX script lines were read for group 4 TITLE G4: STANDARDIZATION CALCULATION BEGIN MATRICES = GROUP 1; END MATRICES; BEGIN ALGEBRA; V=A+C+E+D+Q; ! TOTAL VARIANCE P=A|C|E|D|Q; ! CONCATENATE PARAMETER ESTIMATES S=P@V~; ! STANDARDIZED PARAMETER ESTIMATES END ALGEBRA; OPTION NDECIMALS=4 OPTION SUB=3176.157,1194 ! LIKELIHOOD AND DF FROM SATURATED MODEL: OZBMIYFSAT.MXO !DROP L 1 1 1 END Summary of VL file data for group 2 ARC1 ARC2 Code 1.0000 2.0000 Number 125.0000 125.0000 Mean -0.0090 -0.0190 Variance 0.8861 0.8757 Minimum -2.3713 -2.9867 Maximum 2.6268 2.2563 Summary of VL file data for group 3 PDE PCE PCD PBE PBD PBC PAE Code -10.0000 -9.0000 -8.0000 -7.0000 -6.0000 -5.0000 -4.0000 Number 392.0000 392.0000 392.0000 392.0000 392.0000 392.0000 392.0000 Mean 0.4989 0.5025 0.4917 0.5025 0.5010 0.4989 0.4989 Variance 0.0006 0.0012 0.0096 0.0012 0.0082 0.0331 0.0007 Minimum 0.0030 0.4870 0.0030 0.4900 0.0030 0.0030 0.0030 Maximum 0.5830 0.9960 0.9960 0.9960 0.9960 0.9960 0.5960 PAD PAC PAB ARC1 ARC2 ARC3 ARC4 Code -3.0000 -2.0000 -1.0000 1.0000 2.0000 3.0000 4.0000 Number 392.0000 392.0000 392.0000 392.0000 392.0000 132.0000 29.0000 Mean 0.4983 0.5008 0.4914 0.0844 0.0470 0.2260 0.2899 Variance 0.0096 0.0312 0.0925 0.9724 1.0117 1.1402 1.0366 Minimum 0.0030 0.0030 0.0030 -2.0600 -2.2000 -1.9900 -1.8400 Maximum 0.9960 0.9960 0.9950 3.1560 3.4310 3.5340 2.7000 ARC5 Code 5.0000 Number 3.0000 Mean -0.0407 Variance 0.5965 Minimum -0.9500 Maximum 0.9380 PARAMETER SPECIFICATIONS 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 I This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a FULL matrix of order 1 by 1 1 QTL 4 MATRIX M This is a FULL matrix of order 1 by 1 1 1 5 MATRIX P This is a computed FULL matrix of order 1 by 5 It has no free parameters specified MATRIX Q This is a computed 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 5 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 = 5 MX PARAMETER ESTIMATES GROUP NUMBER: 4 TITLE G4: STANDARDIZATION MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.4844 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.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.2377 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX I This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX L This is a FULL matrix of order 1 by 1 1 QTL 0.5113 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0615 MATRIX P This is a computed FULL matrix of order 1 by 5 [=A|C|E|D|Q] 1 2 3 4 5 1 0.4844 0.0000 0.2377 0.0000 0.2615 MATRIX Q This is a computed FULL matrix of order 1 by 1 [=L*L'] 1 1 0.2615 MATRIX S This is a computed FULL matrix of order 1 by 5 [=P@V~] 1 2 3 4 5 1 0.4925 0.0000 0.2417 0.0000 0.2658 MATRIX V This is a computed FULL matrix of order 1 by 1 [=A+C+E+D+Q] 1 1 0.9835 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 DOM_GEN 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 ADD_GEN 0.6960 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.4875 Your model has 5 estimated parameters and 1198 Observed statistics -2 times log-likelihood of data >>> 3164.700 Degrees of freedom >>>>>>>>>>>>>>>> 1193 Akaike's Information Criterion >>>> 778.700 Bayesian Information Criterion >>>> -2144.608 Sample size Adjusted BIC >>>> -251.202 Deviance Information Criterion >>>> -1048.314 Submodel fit* >>>>>>>>>>> 3176.157 Submodel df* >>>>>>>>>>> 1194 Difference Chi-squared >>>>>>>> 11.457 Difference d.f. >>>>>>>>>>>>>>> 1 Probability >>>>>>>>>>>>>>>>>>>> 0.001 Akaike's Information Criterion > 9.457 * Submodel statistic supplied by user This problem used 0.0% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.12 Execution 0: 0: 0: 8.95 TOTAL 0: 0: 0: 9.07 Total number of warnings issued: 2 ______________________________________________________________________________ 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 L FULL NVAR 1 FREE ! QTL GENETIC PATH, Q H FULL 1 1 ! SCALAR, 0.5 I FULL 1 1 ! SCALAR, 0.25 M FULL 1 1 FREE END MATRICES; LABEL ROW X ADD_GEN LABEL ROW Y COMM_ENV LABEL ROW Z SPEC_ENV LABEL ROW W DOM_GEN LABEL ROW L QTL MATRIX H .5 MATRIX I .25 START .6 ALL ! START ALL FREE PARAMETERS AT .6 START 0 M 1 1 1 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 Q= L*L'; ! QTL VARIANCE, Q^2 END ALGEBRA; OPTION NO_OUTPUT END The following MX script lines were read for group 2 TITLE G2: MZ DATA DATA NINPUT=2 RECTANGULAR FILE=MZDAT.REC Rectangular continuous data read initiated Note: Maximum ordinal/rectangular record length is: 1000 NOTE: Rectangular file contained 125 records with data that contained a total of 250 observations LABELS ARC1 ARC2 SELECT ARC1 ARC2 ; BEGIN MATRICES = GROUP 1; ! COPY ALL MATRICES FROM GROUP 1, BOTH DECLARED AND COMPUTED NOTE: Selection yields 125 data vectors for analysis NOTE: Vectors contain a total of 250 observations MEANS M|M ; ! 1X2 VECTOR OF EXPECTED MEANS COVARIANCE ! 2X2 COVARIANCE MATRIX A+C+E+D+Q | A+C+D+Q _ A+C+D+Q | A+C+E+D+Q ; OPTION RSIDUALS; OPTION NO_OUTPUT END The following MX script lines were read for group 3 TITLE G3: DZ DATA DATA NINPUT=15 RECTANGULAR FILE=GRID51.REC Rectangular continuous data read initiated Note: Maximum ordinal/rectangular record length is: 1000 NOTE: Rectangular file contained 392 records with data that contained a total of 4868 observations LABELS PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ARC1 ARC2 ARC3 ARC4 ARC5 SELECT PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ARC1 ARC2 ARC3 ARC4 ARC5 ; DEFINITON PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ; NOTE: Selection yields 392 data vectors for analysis NOTE: Vectors contain a total of 4868 observations Note: Global variable previously defined. Updating PAB Note: Global variable previously defined. Updating PAC Note: Global variable previously defined. Updating PAD Note: Global variable previously defined. Updating PAE Note: Global variable previously defined. Updating PBC Note: Global variable previously defined. Updating PBD Note: Global variable previously defined. Updating PBE Note: Global variable previously defined. Updating PCD Note: Global variable previously defined. Updating PCE Note: Global variable previously defined. Updating PDE NOTE: Definition yields 392 data vectors for analysis NOTE: Vectors contain a total of 948 observations BEGIN MATRICES ; A COMPUTED NVAR NVAR =A1 C COMPUTED NVAR NVAR =C1 D COMPUTED NVAR NVAR =D1 E COMPUTED NVAR NVAR =E1 Q COMPUTED NVAR NVAR =Q1 M FULL 1 1 =M1 H FULL 1 1 =H1 I FULL 1 1 =I1 P FULL 1 1 R FULL 1 1 S FULL 1 1 T FULL 1 1 U FULL 1 1 V FULL 1 1 W FULL 1 1 X FULL 1 1 Y FULL 1 1 Z FULL 1 1 END MATRICES; SPECIFY P PAB SPECIFY R PAC SPECIFY S PAD SPECIFY T PAE SPECIFY U PBC SPECIFY V PBD SPECIFY W PBE SPECIFY X PCD SPECIFY Y PCE SPECIFY Z PDE MEANS M|M|M|M|M ; COVARIANCE A+C+E+D+Q | H@A+C+I@D+P@Q | H@A+C+I@D+R@Q | H@A+C+I@D+S@Q | H@A+C+I@D+T@Q _ H@A+C+I@D+P@Q | A+C+E+D+Q | H@A+C+I@D+U@Q | H@A+C+I@D+V@Q | H@A+C+I@D+W@Q _ H@A+C+I@D+R@Q | H@A+C+I@D+U@Q | A+C+E+D+Q | H@A+C+I@D+X@Q | H@A+C+I@D+Y@Q _ H@A+C+I@D+S@Q | H@A+C+I@D+V@Q | H@A+C+I@D+X@Q | A+C+E+D+Q | H@A+C+I@D+Z@Q _ H@A+C+I@D+T@Q | H@A+C+I@D+W@Q | H@A+C+I@D+Y@Q | H@A+C+I@D+Z@Q | A+C+E+D+Q ; OPTION NO_OUTPUT END The following MX script lines were read for group 4 TITLE G4: STANDARDIZATION CALCULATION BEGIN MATRICES = GROUP 1; END MATRICES; BEGIN ALGEBRA; V=A+C+E+D+Q; ! TOTAL VARIANCE P=A|C|E|D|Q; ! CONCATENATE PARAMETER ESTIMATES S=P@V~; ! STANDARDIZED PARAMETER ESTIMATES END ALGEBRA; OPTION NDECIMALS=4 OPTION SUB=3176.157,1194 ! LIKELIHOOD AND DF FROM SATURATED MODEL: OZBMIYFSAT.MXO !DROP L 1 1 1 END Summary of VL file data for group 2 ARC1 ARC2 Code 1.0000 2.0000 Number 125.0000 125.0000 Mean -0.0090 -0.0190 Variance 0.8861 0.8757 Minimum -2.3713 -2.9867 Maximum 2.6268 2.2563 Summary of VL file data for group 3 PDE PCE PCD PBE PBD PBC PAE Code -10.0000 -9.0000 -8.0000 -7.0000 -6.0000 -5.0000 -4.0000 Number 392.0000 392.0000 392.0000 392.0000 392.0000 392.0000 392.0000 Mean 0.4990 0.5025 0.4929 0.5025 0.4963 0.4969 0.4985 Variance 0.0006 0.0012 0.0100 0.0012 0.0079 0.0310 0.0006 Minimum 0.0270 0.5000 0.0000 0.4950 0.0010 0.0010 0.0270 Maximum 0.5660 0.9900 0.9980 0.9970 0.9980 0.9980 0.5000 PAD PAC PAB ARC1 ARC2 ARC3 ARC4 Code -3.0000 -2.0000 -1.0000 1.0000 2.0000 3.0000 4.0000 Number 392.0000 392.0000 392.0000 392.0000 392.0000 132.0000 29.0000 Mean 0.4968 0.4968 0.5003 0.0844 0.0470 0.2260 0.2899 Variance 0.0108 0.0334 0.0948 0.9724 1.0117 1.1402 1.0366 Minimum 0.0010 0.0010 0.0010 -2.0600 -2.2000 -1.9900 -1.8400 Maximum 0.9980 0.9980 0.9980 3.1560 3.4310 3.5340 2.7000 ARC5 Code 5.0000 Number 3.0000 Mean -0.0407 Variance 0.5965 Minimum -0.9500 Maximum 0.9380 PARAMETER SPECIFICATIONS 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 I This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a FULL matrix of order 1 by 1 1 QTL 4 MATRIX M This is a FULL matrix of order 1 by 1 1 1 5 MATRIX P This is a computed FULL matrix of order 1 by 5 It has no free parameters specified MATRIX Q This is a computed 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 5 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 = 5 MX PARAMETER ESTIMATES GROUP NUMBER: 4 TITLE G4: STANDARDIZATION MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.4588 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.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.2377 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX I This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX L This is a FULL matrix of order 1 by 1 1 QTL 0.5361 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0612 MATRIX P This is a computed FULL matrix of order 1 by 5 [=A|C|E|D|Q] 1 2 3 4 5 1 0.4588 0.0000 0.2377 0.0000 0.2874 MATRIX Q This is a computed FULL matrix of order 1 by 1 [=L*L'] 1 1 0.2874 MATRIX S This is a computed FULL matrix of order 1 by 5 [=P@V~] 1 2 3 4 5 1 0.4663 0.0000 0.2416 0.0000 0.2921 MATRIX V This is a computed FULL matrix of order 1 by 1 [=A+C+E+D+Q] 1 1 0.9840 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 DOM_GEN 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 ADD_GEN 0.6774 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.4876 Your model has 5 estimated parameters and 1198 Observed statistics -2 times log-likelihood of data >>> 3163.384 Degrees of freedom >>>>>>>>>>>>>>>> 1193 Akaike's Information Criterion >>>> 777.384 Bayesian Information Criterion >>>> -2145.266 Sample size Adjusted BIC >>>> -251.860 Deviance Information Criterion >>>> -1048.972 Submodel fit* >>>>>>>>>>> 3176.157 Submodel df* >>>>>>>>>>> 1194 Difference Chi-squared >>>>>>>> 12.773 Difference d.f. >>>>>>>>>>>>>>> 1 Probability >>>>>>>>>>>>>>>>>>>> 0.000 Akaike's Information Criterion > 10.773 * Submodel statistic supplied by user This problem used 0.0% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.12 Execution 0: 0: 0: 8.51 TOTAL 0: 0: 0: 8.63 Total number of warnings issued: 2 ______________________________________________________________________________ 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 L FULL NVAR 1 FREE ! QTL GENETIC PATH, Q H FULL 1 1 ! SCALAR, 0.5 I FULL 1 1 ! SCALAR, 0.25 M FULL 1 1 FREE END MATRICES; LABEL ROW X ADD_GEN LABEL ROW Y COMM_ENV LABEL ROW Z SPEC_ENV LABEL ROW W DOM_GEN LABEL ROW L QTL MATRIX H .5 MATRIX I .25 START .6 ALL ! START ALL FREE PARAMETERS AT .6 START 0 M 1 1 1 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 Q= L*L'; ! QTL VARIANCE, Q^2 END ALGEBRA; OPTION NO_OUTPUT END The following MX script lines were read for group 2 TITLE G2: MZ DATA DATA NINPUT=2 RECTANGULAR FILE=MZDAT.REC Rectangular continuous data read initiated Note: Maximum ordinal/rectangular record length is: 1000 NOTE: Rectangular file contained 125 records with data that contained a total of 250 observations LABELS ARC1 ARC2 SELECT ARC1 ARC2 ; BEGIN MATRICES = GROUP 1; ! COPY ALL MATRICES FROM GROUP 1, BOTH DECLARED AND COMPUTED NOTE: Selection yields 125 data vectors for analysis NOTE: Vectors contain a total of 250 observations MEANS M|M ; ! 1X2 VECTOR OF EXPECTED MEANS COVARIANCE ! 2X2 COVARIANCE MATRIX A+C+E+D+Q | A+C+D+Q _ A+C+D+Q | A+C+E+D+Q ; OPTION RSIDUALS; OPTION NO_OUTPUT END The following MX script lines were read for group 3 TITLE G3: DZ DATA DATA NINPUT=15 RECTANGULAR FILE=GRID52.REC Rectangular continuous data read initiated Note: Maximum ordinal/rectangular record length is: 1000 NOTE: Rectangular file contained 392 records with data that contained a total of 4868 observations LABELS PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ARC1 ARC2 ARC3 ARC4 ARC5 SELECT PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ARC1 ARC2 ARC3 ARC4 ARC5 ; DEFINITON PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ; NOTE: Selection yields 392 data vectors for analysis NOTE: Vectors contain a total of 4868 observations Note: Global variable previously defined. Updating PAB Note: Global variable previously defined. Updating PAC Note: Global variable previously defined. Updating PAD Note: Global variable previously defined. Updating PAE Note: Global variable previously defined. Updating PBC Note: Global variable previously defined. Updating PBD Note: Global variable previously defined. Updating PBE Note: Global variable previously defined. Updating PCD Note: Global variable previously defined. Updating PCE Note: Global variable previously defined. Updating PDE NOTE: Definition yields 392 data vectors for analysis NOTE: Vectors contain a total of 948 observations BEGIN MATRICES ; A COMPUTED NVAR NVAR =A1 C COMPUTED NVAR NVAR =C1 D COMPUTED NVAR NVAR =D1 E COMPUTED NVAR NVAR =E1 Q COMPUTED NVAR NVAR =Q1 M FULL 1 1 =M1 H FULL 1 1 =H1 I FULL 1 1 =I1 P FULL 1 1 R FULL 1 1 S FULL 1 1 T FULL 1 1 U FULL 1 1 V FULL 1 1 W FULL 1 1 X FULL 1 1 Y FULL 1 1 Z FULL 1 1 END MATRICES; SPECIFY P PAB SPECIFY R PAC SPECIFY S PAD SPECIFY T PAE SPECIFY U PBC SPECIFY V PBD SPECIFY W PBE SPECIFY X PCD SPECIFY Y PCE SPECIFY Z PDE MEANS M|M|M|M|M ; COVARIANCE A+C+E+D+Q | H@A+C+I@D+P@Q | H@A+C+I@D+R@Q | H@A+C+I@D+S@Q | H@A+C+I@D+T@Q _ H@A+C+I@D+P@Q | A+C+E+D+Q | H@A+C+I@D+U@Q | H@A+C+I@D+V@Q | H@A+C+I@D+W@Q _ H@A+C+I@D+R@Q | H@A+C+I@D+U@Q | A+C+E+D+Q | H@A+C+I@D+X@Q | H@A+C+I@D+Y@Q _ H@A+C+I@D+S@Q | H@A+C+I@D+V@Q | H@A+C+I@D+X@Q | A+C+E+D+Q | H@A+C+I@D+Z@Q _ H@A+C+I@D+T@Q | H@A+C+I@D+W@Q | H@A+C+I@D+Y@Q | H@A+C+I@D+Z@Q | A+C+E+D+Q ; OPTION NO_OUTPUT END The following MX script lines were read for group 4 TITLE G4: STANDARDIZATION CALCULATION BEGIN MATRICES = GROUP 1; END MATRICES; BEGIN ALGEBRA; V=A+C+E+D+Q; ! TOTAL VARIANCE P=A|C|E|D|Q; ! CONCATENATE PARAMETER ESTIMATES S=P@V~; ! STANDARDIZED PARAMETER ESTIMATES END ALGEBRA; OPTION NDECIMALS=4 OPTION SUB=3176.157,1194 ! LIKELIHOOD AND DF FROM SATURATED MODEL: OZBMIYFSAT.MXO !DROP L 1 1 1 END Summary of VL file data for group 2 ARC1 ARC2 Code 1.0000 2.0000 Number 125.0000 125.0000 Mean -0.0090 -0.0190 Variance 0.8861 0.8757 Minimum -2.3713 -2.9867 Maximum 2.6268 2.2563 Summary of VL file data for group 3 PDE PCE PCD PBE PBD PBC PAE Code -10.0000 -9.0000 -8.0000 -7.0000 -6.0000 -5.0000 -4.0000 Number 392.0000 392.0000 392.0000 392.0000 392.0000 392.0000 392.0000 Mean 0.4996 0.5021 0.4928 0.5021 0.4935 0.5016 0.4979 Variance 0.0003 0.0009 0.0096 0.0009 0.0084 0.0311 0.0009 Minimum 0.1810 0.5000 0.0010 0.5000 0.0010 0.0010 0.0100 Maximum 0.6850 0.9870 0.9870 0.9940 0.9980 0.9980 0.5000 PAD PAC PAB ARC1 ARC2 ARC3 ARC4 Code -3.0000 -2.0000 -1.0000 1.0000 2.0000 3.0000 4.0000 Number 392.0000 392.0000 392.0000 392.0000 392.0000 132.0000 29.0000 Mean 0.4954 0.4919 0.5092 0.0844 0.0470 0.2260 0.2899 Variance 0.0111 0.0342 0.0990 0.9724 1.0117 1.1402 1.0366 Minimum 0.0040 0.0010 0.0010 -2.0600 -2.2000 -1.9900 -1.8400 Maximum 0.9980 0.9980 0.9980 3.1560 3.4310 3.5340 2.7000 ARC5 Code 5.0000 Number 3.0000 Mean -0.0407 Variance 0.5965 Minimum -0.9500 Maximum 0.9380 PARAMETER SPECIFICATIONS 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 I This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a FULL matrix of order 1 by 1 1 QTL 4 MATRIX M This is a FULL matrix of order 1 by 1 1 1 5 MATRIX P This is a computed FULL matrix of order 1 by 5 It has no free parameters specified MATRIX Q This is a computed 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 5 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 = 5 MX PARAMETER ESTIMATES GROUP NUMBER: 4 TITLE G4: STANDARDIZATION MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.5014 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.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.2382 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX I This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX L This is a FULL matrix of order 1 by 1 1 QTL 0.4954 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0618 MATRIX P This is a computed FULL matrix of order 1 by 5 [=A|C|E|D|Q] 1 2 3 4 5 1 0.5014 0.0000 0.2382 0.0000 0.2454 MATRIX Q This is a computed FULL matrix of order 1 by 1 [=L*L'] 1 1 0.2454 MATRIX S This is a computed FULL matrix of order 1 by 5 [=P@V~] 1 2 3 4 5 1 0.5090 0.0000 0.2418 0.0000 0.2492 MATRIX V This is a computed FULL matrix of order 1 by 1 [=A+C+E+D+Q] 1 1 0.9851 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 DOM_GEN 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 ADD_GEN 0.7081 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.4881 Your model has 5 estimated parameters and 1198 Observed statistics -2 times log-likelihood of data >>> 3166.970 Degrees of freedom >>>>>>>>>>>>>>>> 1193 Akaike's Information Criterion >>>> 780.970 Bayesian Information Criterion >>>> -2143.473 Sample size Adjusted BIC >>>> -250.067 Deviance Information Criterion >>>> -1047.179 Submodel fit* >>>>>>>>>>> 3176.157 Submodel df* >>>>>>>>>>> 1194 Difference Chi-squared >>>>>>>> 9.187 Difference d.f. >>>>>>>>>>>>>>> 1 Probability >>>>>>>>>>>>>>>>>>>> 0.002 Akaike's Information Criterion > 7.187 * Submodel statistic supplied by user This problem used 0.0% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.12 Execution 0: 0: 0: 8.52 TOTAL 0: 0: 0: 8.64 Total number of warnings issued: 2 ______________________________________________________________________________ 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 L FULL NVAR 1 FREE ! QTL GENETIC PATH, Q H FULL 1 1 ! SCALAR, 0.5 I FULL 1 1 ! SCALAR, 0.25 M FULL 1 1 FREE END MATRICES; LABEL ROW X ADD_GEN LABEL ROW Y COMM_ENV LABEL ROW Z SPEC_ENV LABEL ROW W DOM_GEN LABEL ROW L QTL MATRIX H .5 MATRIX I .25 START .6 ALL ! START ALL FREE PARAMETERS AT .6 START 0 M 1 1 1 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 Q= L*L'; ! QTL VARIANCE, Q^2 END ALGEBRA; OPTION NO_OUTPUT END The following MX script lines were read for group 2 TITLE G2: MZ DATA DATA NINPUT=2 RECTANGULAR FILE=MZDAT.REC Rectangular continuous data read initiated Note: Maximum ordinal/rectangular record length is: 1000 NOTE: Rectangular file contained 125 records with data that contained a total of 250 observations LABELS ARC1 ARC2 SELECT ARC1 ARC2 ; BEGIN MATRICES = GROUP 1; ! COPY ALL MATRICES FROM GROUP 1, BOTH DECLARED AND COMPUTED NOTE: Selection yields 125 data vectors for analysis NOTE: Vectors contain a total of 250 observations MEANS M|M ; ! 1X2 VECTOR OF EXPECTED MEANS COVARIANCE ! 2X2 COVARIANCE MATRIX A+C+E+D+Q | A+C+D+Q _ A+C+D+Q | A+C+E+D+Q ; OPTION RSIDUALS; OPTION NO_OUTPUT END The following MX script lines were read for group 3 TITLE G3: DZ DATA DATA NINPUT=15 RECTANGULAR FILE=GRID53.REC Rectangular continuous data read initiated Note: Maximum ordinal/rectangular record length is: 1000 NOTE: Rectangular file contained 392 records with data that contained a total of 4868 observations LABELS PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ARC1 ARC2 ARC3 ARC4 ARC5 SELECT PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ARC1 ARC2 ARC3 ARC4 ARC5 ; DEFINITON PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ; NOTE: Selection yields 392 data vectors for analysis NOTE: Vectors contain a total of 4868 observations Note: Global variable previously defined. Updating PAB Note: Global variable previously defined. Updating PAC Note: Global variable previously defined. Updating PAD Note: Global variable previously defined. Updating PAE Note: Global variable previously defined. Updating PBC Note: Global variable previously defined. Updating PBD Note: Global variable previously defined. Updating PBE Note: Global variable previously defined. Updating PCD Note: Global variable previously defined. Updating PCE Note: Global variable previously defined. Updating PDE NOTE: Definition yields 392 data vectors for analysis NOTE: Vectors contain a total of 948 observations BEGIN MATRICES ; A COMPUTED NVAR NVAR =A1 C COMPUTED NVAR NVAR =C1 D COMPUTED NVAR NVAR =D1 E COMPUTED NVAR NVAR =E1 Q COMPUTED NVAR NVAR =Q1 M FULL 1 1 =M1 H FULL 1 1 =H1 I FULL 1 1 =I1 P FULL 1 1 R FULL 1 1 S FULL 1 1 T FULL 1 1 U FULL 1 1 V FULL 1 1 W FULL 1 1 X FULL 1 1 Y FULL 1 1 Z FULL 1 1 END MATRICES; SPECIFY P PAB SPECIFY R PAC SPECIFY S PAD SPECIFY T PAE SPECIFY U PBC SPECIFY V PBD SPECIFY W PBE SPECIFY X PCD SPECIFY Y PCE SPECIFY Z PDE MEANS M|M|M|M|M ; COVARIANCE A+C+E+D+Q | H@A+C+I@D+P@Q | H@A+C+I@D+R@Q | H@A+C+I@D+S@Q | H@A+C+I@D+T@Q _ H@A+C+I@D+P@Q | A+C+E+D+Q | H@A+C+I@D+U@Q | H@A+C+I@D+V@Q | H@A+C+I@D+W@Q _ H@A+C+I@D+R@Q | H@A+C+I@D+U@Q | A+C+E+D+Q | H@A+C+I@D+X@Q | H@A+C+I@D+Y@Q _ H@A+C+I@D+S@Q | H@A+C+I@D+V@Q | H@A+C+I@D+X@Q | A+C+E+D+Q | H@A+C+I@D+Z@Q _ H@A+C+I@D+T@Q | H@A+C+I@D+W@Q | H@A+C+I@D+Y@Q | H@A+C+I@D+Z@Q | A+C+E+D+Q ; OPTION NO_OUTPUT END The following MX script lines were read for group 4 TITLE G4: STANDARDIZATION CALCULATION BEGIN MATRICES = GROUP 1; END MATRICES; BEGIN ALGEBRA; V=A+C+E+D+Q; ! TOTAL VARIANCE P=A|C|E|D|Q; ! CONCATENATE PARAMETER ESTIMATES S=P@V~; ! STANDARDIZED PARAMETER ESTIMATES END ALGEBRA; OPTION NDECIMALS=4 OPTION SUB=3176.157,1194 ! LIKELIHOOD AND DF FROM SATURATED MODEL: OZBMIYFSAT.MXO !DROP L 1 1 1 END Summary of VL file data for group 2 ARC1 ARC2 Code 1.0000 2.0000 Number 125.0000 125.0000 Mean -0.0090 -0.0190 Variance 0.8861 0.8757 Minimum -2.3713 -2.9867 Maximum 2.6268 2.2563 Summary of VL file data for group 3 PDE PCE PCD PBE PBD PBC PAE Code -10.0000 -9.0000 -8.0000 -7.0000 -6.0000 -5.0000 -4.0000 Number 392.0000 392.0000 392.0000 392.0000 392.0000 392.0000 392.0000 Mean 0.5008 0.5011 0.4921 0.5017 0.4942 0.5033 0.4983 Variance 0.0002 0.0006 0.0085 0.0007 0.0070 0.0256 0.0007 Minimum 0.4950 0.4550 0.0090 0.5000 0.0060 0.0060 0.0150 Maximum 0.7940 0.9830 0.9810 0.9900 0.9910 0.9920 0.5000 PAD PAC PAB ARC1 ARC2 ARC3 ARC4 Code -3.0000 -2.0000 -1.0000 1.0000 2.0000 3.0000 4.0000 Number 392.0000 392.0000 392.0000 392.0000 392.0000 132.0000 29.0000 Mean 0.4966 0.4926 0.5140 0.0844 0.0470 0.2260 0.2899 Variance 0.0085 0.0263 0.0838 0.9724 1.0117 1.1402 1.0366 Minimum 0.0080 0.0070 0.0060 -2.0600 -2.2000 -1.9900 -1.8400 Maximum 0.9930 0.9910 0.9920 3.1560 3.4310 3.5340 2.7000 ARC5 Code 5.0000 Number 3.0000 Mean -0.0407 Variance 0.5965 Minimum -0.9500 Maximum 0.9380 PARAMETER SPECIFICATIONS 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 I This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a FULL matrix of order 1 by 1 1 QTL 4 MATRIX M This is a FULL matrix of order 1 by 1 1 1 5 MATRIX P This is a computed FULL matrix of order 1 by 5 It has no free parameters specified MATRIX Q This is a computed 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 5 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 = 5 MX PARAMETER ESTIMATES GROUP NUMBER: 4 TITLE G4: STANDARDIZATION MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.5290 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.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.2389 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX I This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX L This is a FULL matrix of order 1 by 1 1 QTL 0.4663 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0621 MATRIX P This is a computed FULL matrix of order 1 by 5 [=A|C|E|D|Q] 1 2 3 4 5 1 0.5290 0.0000 0.2389 0.0000 0.2174 MATRIX Q This is a computed FULL matrix of order 1 by 1 [=L*L'] 1 1 0.2174 MATRIX S This is a computed FULL matrix of order 1 by 5 [=P@V~] 1 2 3 4 5 1 0.5369 0.0000 0.2424 0.0000 0.2207 MATRIX V This is a computed FULL matrix of order 1 by 1 [=A+C+E+D+Q] 1 1 0.9853 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 DOM_GEN 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 ADD_GEN 0.7274 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.4887 Your model has 5 estimated parameters and 1198 Observed statistics -2 times log-likelihood of data >>> 3169.886 Degrees of freedom >>>>>>>>>>>>>>>> 1193 Akaike's Information Criterion >>>> 783.886 Bayesian Information Criterion >>>> -2142.015 Sample size Adjusted BIC >>>> -248.609 Deviance Information Criterion >>>> -1045.721 Submodel fit* >>>>>>>>>>> 3176.157 Submodel df* >>>>>>>>>>> 1194 Difference Chi-squared >>>>>>>> 6.271 Difference d.f. >>>>>>>>>>>>>>> 1 Probability >>>>>>>>>>>>>>>>>>>> 0.012 Akaike's Information Criterion > 4.271 * Submodel statistic supplied by user This problem used 0.0% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.14 Execution 0: 0: 0: 9.25 TOTAL 0: 0: 0: 9.39 Total number of warnings issued: 2 ______________________________________________________________________________ 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 L FULL NVAR 1 FREE ! QTL GENETIC PATH, Q H FULL 1 1 ! SCALAR, 0.5 I FULL 1 1 ! SCALAR, 0.25 M FULL 1 1 FREE END MATRICES; LABEL ROW X ADD_GEN LABEL ROW Y COMM_ENV LABEL ROW Z SPEC_ENV LABEL ROW W DOM_GEN LABEL ROW L QTL MATRIX H .5 MATRIX I .25 START .6 ALL ! START ALL FREE PARAMETERS AT .6 START 0 M 1 1 1 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 Q= L*L'; ! QTL VARIANCE, Q^2 END ALGEBRA; OPTION NO_OUTPUT END The following MX script lines were read for group 2 TITLE G2: MZ DATA DATA NINPUT=2 RECTANGULAR FILE=MZDAT.REC Rectangular continuous data read initiated Note: Maximum ordinal/rectangular record length is: 1000 NOTE: Rectangular file contained 125 records with data that contained a total of 250 observations LABELS ARC1 ARC2 SELECT ARC1 ARC2 ; BEGIN MATRICES = GROUP 1; ! COPY ALL MATRICES FROM GROUP 1, BOTH DECLARED AND COMPUTED NOTE: Selection yields 125 data vectors for analysis NOTE: Vectors contain a total of 250 observations MEANS M|M ; ! 1X2 VECTOR OF EXPECTED MEANS COVARIANCE ! 2X2 COVARIANCE MATRIX A+C+E+D+Q | A+C+D+Q _ A+C+D+Q | A+C+E+D+Q ; OPTION RSIDUALS; OPTION NO_OUTPUT END The following MX script lines were read for group 3 TITLE G3: DZ DATA DATA NINPUT=15 RECTANGULAR FILE=GRID54.REC Rectangular continuous data read initiated Note: Maximum ordinal/rectangular record length is: 1000 NOTE: Rectangular file contained 392 records with data that contained a total of 4868 observations LABELS PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ARC1 ARC2 ARC3 ARC4 ARC5 SELECT PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ARC1 ARC2 ARC3 ARC4 ARC5 ; DEFINITON PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ; NOTE: Selection yields 392 data vectors for analysis NOTE: Vectors contain a total of 4868 observations Note: Global variable previously defined. Updating PAB Note: Global variable previously defined. Updating PAC Note: Global variable previously defined. Updating PAD Note: Global variable previously defined. Updating PAE Note: Global variable previously defined. Updating PBC Note: Global variable previously defined. Updating PBD Note: Global variable previously defined. Updating PBE Note: Global variable previously defined. Updating PCD Note: Global variable previously defined. Updating PCE Note: Global variable previously defined. Updating PDE NOTE: Definition yields 392 data vectors for analysis NOTE: Vectors contain a total of 948 observations BEGIN MATRICES ; A COMPUTED NVAR NVAR =A1 C COMPUTED NVAR NVAR =C1 D COMPUTED NVAR NVAR =D1 E COMPUTED NVAR NVAR =E1 Q COMPUTED NVAR NVAR =Q1 M FULL 1 1 =M1 H FULL 1 1 =H1 I FULL 1 1 =I1 P FULL 1 1 R FULL 1 1 S FULL 1 1 T FULL 1 1 U FULL 1 1 V FULL 1 1 W FULL 1 1 X FULL 1 1 Y FULL 1 1 Z FULL 1 1 END MATRICES; SPECIFY P PAB SPECIFY R PAC SPECIFY S PAD SPECIFY T PAE SPECIFY U PBC SPECIFY V PBD SPECIFY W PBE SPECIFY X PCD SPECIFY Y PCE SPECIFY Z PDE MEANS M|M|M|M|M ; COVARIANCE A+C+E+D+Q | H@A+C+I@D+P@Q | H@A+C+I@D+R@Q | H@A+C+I@D+S@Q | H@A+C+I@D+T@Q _ H@A+C+I@D+P@Q | A+C+E+D+Q | H@A+C+I@D+U@Q | H@A+C+I@D+V@Q | H@A+C+I@D+W@Q _ H@A+C+I@D+R@Q | H@A+C+I@D+U@Q | A+C+E+D+Q | H@A+C+I@D+X@Q | H@A+C+I@D+Y@Q _ H@A+C+I@D+S@Q | H@A+C+I@D+V@Q | H@A+C+I@D+X@Q | A+C+E+D+Q | H@A+C+I@D+Z@Q _ H@A+C+I@D+T@Q | H@A+C+I@D+W@Q | H@A+C+I@D+Y@Q | H@A+C+I@D+Z@Q | A+C+E+D+Q ; OPTION NO_OUTPUT END The following MX script lines were read for group 4 TITLE G4: STANDARDIZATION CALCULATION BEGIN MATRICES = GROUP 1; END MATRICES; BEGIN ALGEBRA; V=A+C+E+D+Q; ! TOTAL VARIANCE P=A|C|E|D|Q; ! CONCATENATE PARAMETER ESTIMATES S=P@V~; ! STANDARDIZED PARAMETER ESTIMATES END ALGEBRA; OPTION NDECIMALS=4 OPTION SUB=3176.157,1194 ! LIKELIHOOD AND DF FROM SATURATED MODEL: OZBMIYFSAT.MXO !DROP L 1 1 1 END Summary of VL file data for group 2 ARC1 ARC2 Code 1.0000 2.0000 Number 125.0000 125.0000 Mean -0.0090 -0.0190 Variance 0.8861 0.8757 Minimum -2.3713 -2.9867 Maximum 2.6268 2.2563 Summary of VL file data for group 3 PDE PCE PCD PBE PBD PBC PAE Code -10.0000 -9.0000 -8.0000 -7.0000 -6.0000 -5.0000 -4.0000 Number 392.0000 392.0000 392.0000 392.0000 392.0000 392.0000 392.0000 Mean 0.5021 0.5002 0.4901 0.5013 0.4945 0.5047 0.4987 Variance 0.0009 0.0010 0.0106 0.0006 0.0077 0.0310 0.0006 Minimum 0.4980 0.0890 0.0030 0.5000 0.0020 0.0020 0.0090 Maximum 0.9110 0.9890 0.9940 0.9960 0.9970 0.9970 0.5000 PAD PAC PAB ARC1 ARC2 ARC3 ARC4 Code -3.0000 -2.0000 -1.0000 1.0000 2.0000 3.0000 4.0000 Number 392.0000 392.0000 392.0000 392.0000 392.0000 132.0000 29.0000 Mean 0.4973 0.4931 0.5189 0.0844 0.0470 0.2260 0.2899 Variance 0.0089 0.0302 0.0976 0.9724 1.0117 1.1402 1.0366 Minimum 0.0020 0.0020 0.0020 -2.0600 -2.2000 -1.9900 -1.8400 Maximum 0.9970 0.9950 0.9970 3.1560 3.4310 3.5340 2.7000 ARC5 Code 5.0000 Number 3.0000 Mean -0.0407 Variance 0.5965 Minimum -0.9500 Maximum 0.9380 PARAMETER SPECIFICATIONS 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 I This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a FULL matrix of order 1 by 1 1 QTL 4 MATRIX M This is a FULL matrix of order 1 by 1 1 1 5 MATRIX P This is a computed FULL matrix of order 1 by 5 It has no free parameters specified MATRIX Q This is a computed 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 5 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 = 5 MX PARAMETER ESTIMATES GROUP NUMBER: 4 TITLE G4: STANDARDIZATION MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.6091 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.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.2390 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX I This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX L This is a FULL matrix of order 1 by 1 1 QTL 0.3715 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0631 MATRIX P This is a computed FULL matrix of order 1 by 5 [=A|C|E|D|Q] 1 2 3 4 5 1 0.6091 0.0000 0.2390 0.0000 0.1380 MATRIX Q This is a computed FULL matrix of order 1 by 1 [=L*L'] 1 1 0.1380 MATRIX S This is a computed FULL matrix of order 1 by 5 [=P@V~] 1 2 3 4 5 1 0.6177 0.0000 0.2424 0.0000 0.1399 MATRIX V This is a computed FULL matrix of order 1 by 1 [=A+C+E+D+Q] 1 1 0.9862 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 DOM_GEN 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 ADD_GEN 0.7805 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.4889 Your model has 5 estimated parameters and 1198 Observed statistics -2 times log-likelihood of data >>> 3173.299 Degrees of freedom >>>>>>>>>>>>>>>> 1193 Akaike's Information Criterion >>>> 787.299 Bayesian Information Criterion >>>> -2140.308 Sample size Adjusted BIC >>>> -246.902 Deviance Information Criterion >>>> -1044.014 Submodel fit* >>>>>>>>>>> 3176.157 Submodel df* >>>>>>>>>>> 1194 Difference Chi-squared >>>>>>>> 2.858 Difference d.f. >>>>>>>>>>>>>>> 1 Probability >>>>>>>>>>>>>>>>>>>> 0.091 Akaike's Information Criterion > 0.858 * Submodel statistic supplied by user This problem used 0.0% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.16 Execution 0: 0: 0: 9.78 TOTAL 0: 0: 0: 9.94 Total number of warnings issued: 2 ______________________________________________________________________________ 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 L FULL NVAR 1 FREE ! QTL GENETIC PATH, Q H FULL 1 1 ! SCALAR, 0.5 I FULL 1 1 ! SCALAR, 0.25 M FULL 1 1 FREE END MATRICES; LABEL ROW X ADD_GEN LABEL ROW Y COMM_ENV LABEL ROW Z SPEC_ENV LABEL ROW W DOM_GEN LABEL ROW L QTL MATRIX H .5 MATRIX I .25 START .6 ALL ! START ALL FREE PARAMETERS AT .6 START 0 M 1 1 1 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 Q= L*L'; ! QTL VARIANCE, Q^2 END ALGEBRA; OPTION NO_OUTPUT END The following MX script lines were read for group 2 TITLE G2: MZ DATA DATA NINPUT=2 RECTANGULAR FILE=MZDAT.REC Rectangular continuous data read initiated Note: Maximum ordinal/rectangular record length is: 1000 NOTE: Rectangular file contained 125 records with data that contained a total of 250 observations LABELS ARC1 ARC2 SELECT ARC1 ARC2 ; BEGIN MATRICES = GROUP 1; ! COPY ALL MATRICES FROM GROUP 1, BOTH DECLARED AND COMPUTED NOTE: Selection yields 125 data vectors for analysis NOTE: Vectors contain a total of 250 observations MEANS M|M ; ! 1X2 VECTOR OF EXPECTED MEANS COVARIANCE ! 2X2 COVARIANCE MATRIX A+C+E+D+Q | A+C+D+Q _ A+C+D+Q | A+C+E+D+Q ; OPTION RSIDUALS; OPTION NO_OUTPUT END The following MX script lines were read for group 3 TITLE G3: DZ DATA DATA NINPUT=15 RECTANGULAR FILE=GRID55.REC Rectangular continuous data read initiated Note: Maximum ordinal/rectangular record length is: 1000 NOTE: Rectangular file contained 392 records with data that contained a total of 4868 observations LABELS PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ARC1 ARC2 ARC3 ARC4 ARC5 SELECT PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ARC1 ARC2 ARC3 ARC4 ARC5 ; DEFINITON PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ; NOTE: Selection yields 392 data vectors for analysis NOTE: Vectors contain a total of 4868 observations Note: Global variable previously defined. Updating PAB Note: Global variable previously defined. Updating PAC Note: Global variable previously defined. Updating PAD Note: Global variable previously defined. Updating PAE Note: Global variable previously defined. Updating PBC Note: Global variable previously defined. Updating PBD Note: Global variable previously defined. Updating PBE Note: Global variable previously defined. Updating PCD Note: Global variable previously defined. Updating PCE Note: Global variable previously defined. Updating PDE NOTE: Definition yields 392 data vectors for analysis NOTE: Vectors contain a total of 948 observations BEGIN MATRICES ; A COMPUTED NVAR NVAR =A1 C COMPUTED NVAR NVAR =C1 D COMPUTED NVAR NVAR =D1 E COMPUTED NVAR NVAR =E1 Q COMPUTED NVAR NVAR =Q1 M FULL 1 1 =M1 H FULL 1 1 =H1 I FULL 1 1 =I1 P FULL 1 1 R FULL 1 1 S FULL 1 1 T FULL 1 1 U FULL 1 1 V FULL 1 1 W FULL 1 1 X FULL 1 1 Y FULL 1 1 Z FULL 1 1 END MATRICES; SPECIFY P PAB SPECIFY R PAC SPECIFY S PAD SPECIFY T PAE SPECIFY U PBC SPECIFY V PBD SPECIFY W PBE SPECIFY X PCD SPECIFY Y PCE SPECIFY Z PDE MEANS M|M|M|M|M ; COVARIANCE A+C+E+D+Q | H@A+C+I@D+P@Q | H@A+C+I@D+R@Q | H@A+C+I@D+S@Q | H@A+C+I@D+T@Q _ H@A+C+I@D+P@Q | A+C+E+D+Q | H@A+C+I@D+U@Q | H@A+C+I@D+V@Q | H@A+C+I@D+W@Q _ H@A+C+I@D+R@Q | H@A+C+I@D+U@Q | A+C+E+D+Q | H@A+C+I@D+X@Q | H@A+C+I@D+Y@Q _ H@A+C+I@D+S@Q | H@A+C+I@D+V@Q | H@A+C+I@D+X@Q | A+C+E+D+Q | H@A+C+I@D+Z@Q _ H@A+C+I@D+T@Q | H@A+C+I@D+W@Q | H@A+C+I@D+Y@Q | H@A+C+I@D+Z@Q | A+C+E+D+Q ; OPTION NO_OUTPUT END The following MX script lines were read for group 4 TITLE G4: STANDARDIZATION CALCULATION BEGIN MATRICES = GROUP 1; END MATRICES; BEGIN ALGEBRA; V=A+C+E+D+Q; ! TOTAL VARIANCE P=A|C|E|D|Q; ! CONCATENATE PARAMETER ESTIMATES S=P@V~; ! STANDARDIZED PARAMETER ESTIMATES END ALGEBRA; OPTION NDECIMALS=4 OPTION SUB=3176.157,1194 ! LIKELIHOOD AND DF FROM SATURATED MODEL: OZBMIYFSAT.MXO !DROP L 1 1 1 END Summary of VL file data for group 2 ARC1 ARC2 Code 1.0000 2.0000 Number 125.0000 125.0000 Mean -0.0090 -0.0190 Variance 0.8861 0.8757 Minimum -2.3713 -2.9867 Maximum 2.6268 2.2563 Summary of VL file data for group 3 PDE PCE PCD PBE PBD PBC PAE Code -10.0000 -9.0000 -8.0000 -7.0000 -6.0000 -5.0000 -4.0000 Number 392.0000 392.0000 392.0000 392.0000 392.0000 392.0000 392.0000 Mean 0.5025 0.5000 0.4898 0.5013 0.4931 0.5060 0.4987 Variance 0.0013 0.0013 0.0113 0.0006 0.0090 0.0340 0.0006 Minimum 0.5000 0.0000 0.0000 0.4950 0.0000 0.0000 0.0030 Maximum 0.9990 0.9980 0.9990 0.9990 0.9990 0.9990 0.5010 PAD PAC PAB ARC1 ARC2 ARC3 ARC4 Code -3.0000 -2.0000 -1.0000 1.0000 2.0000 3.0000 4.0000 Number 392.0000 392.0000 392.0000 392.0000 392.0000 132.0000 29.0000 Mean 0.4984 0.4940 0.5111 0.0844 0.0470 0.2260 0.2899 Variance 0.0092 0.0335 0.1047 0.9724 1.0117 1.1402 1.0366 Minimum 0.0010 0.0000 0.0000 -2.0600 -2.2000 -1.9900 -1.8400 Maximum 0.9990 0.9990 0.9990 3.1560 3.4310 3.5340 2.7000 ARC5 Code 5.0000 Number 3.0000 Mean -0.0407 Variance 0.5965 Minimum -0.9500 Maximum 0.9380 PARAMETER SPECIFICATIONS 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 I This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a FULL matrix of order 1 by 1 1 QTL 4 MATRIX M This is a FULL matrix of order 1 by 1 1 1 5 MATRIX P This is a computed FULL matrix of order 1 by 5 It has no free parameters specified MATRIX Q This is a computed 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 5 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 = 5 MX PARAMETER ESTIMATES GROUP NUMBER: 4 TITLE G4: STANDARDIZATION MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.6638 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.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.2389 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX I This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX L This is a FULL matrix of order 1 by 1 1 QTL 0.2893 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0639 MATRIX P This is a computed FULL matrix of order 1 by 5 [=A|C|E|D|Q] 1 2 3 4 5 1 0.6638 0.0000 0.2389 0.0000 0.0837 MATRIX Q This is a computed FULL matrix of order 1 by 1 [=L*L'] 1 1 0.0837 MATRIX S This is a computed FULL matrix of order 1 by 5 [=P@V~] 1 2 3 4 5 1 0.6729 0.0000 0.2422 0.0000 0.0849 MATRIX V This is a computed FULL matrix of order 1 by 1 [=A+C+E+D+Q] 1 1 0.9865 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 DOM_GEN 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 ADD_GEN 0.8147 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.4888 Your model has 5 estimated parameters and 1198 Observed statistics -2 times log-likelihood of data >>> 3175.174 Degrees of freedom >>>>>>>>>>>>>>>> 1193 Akaike's Information Criterion >>>> 789.174 Bayesian Information Criterion >>>> -2139.370 Sample size Adjusted BIC >>>> -245.964 Deviance Information Criterion >>>> -1043.077 Submodel fit* >>>>>>>>>>> 3176.157 Submodel df* >>>>>>>>>>> 1194 Difference Chi-squared >>>>>>>> 0.983 Difference d.f. >>>>>>>>>>>>>>> 1 Probability >>>>>>>>>>>>>>>>>>>> 0.322 Akaike's Information Criterion > -1.017 * Submodel statistic supplied by user This problem used 0.0% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.14 Execution 0: 0: 0:10.31 TOTAL 0: 0: 0:10.45 Total number of warnings issued: 2 ______________________________________________________________________________ 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 L FULL NVAR 1 FREE ! QTL GENETIC PATH, Q H FULL 1 1 ! SCALAR, 0.5 I FULL 1 1 ! SCALAR, 0.25 M FULL 1 1 FREE END MATRICES; LABEL ROW X ADD_GEN LABEL ROW Y COMM_ENV LABEL ROW Z SPEC_ENV LABEL ROW W DOM_GEN LABEL ROW L QTL MATRIX H .5 MATRIX I .25 START .6 ALL ! START ALL FREE PARAMETERS AT .6 START 0 M 1 1 1 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 Q= L*L'; ! QTL VARIANCE, Q^2 END ALGEBRA; OPTION NO_OUTPUT END The following MX script lines were read for group 2 TITLE G2: MZ DATA DATA NINPUT=2 RECTANGULAR FILE=MZDAT.REC Rectangular continuous data read initiated Note: Maximum ordinal/rectangular record length is: 1000 NOTE: Rectangular file contained 125 records with data that contained a total of 250 observations LABELS ARC1 ARC2 SELECT ARC1 ARC2 ; BEGIN MATRICES = GROUP 1; ! COPY ALL MATRICES FROM GROUP 1, BOTH DECLARED AND COMPUTED NOTE: Selection yields 125 data vectors for analysis NOTE: Vectors contain a total of 250 observations MEANS M|M ; ! 1X2 VECTOR OF EXPECTED MEANS COVARIANCE ! 2X2 COVARIANCE MATRIX A+C+E+D+Q | A+C+D+Q _ A+C+D+Q | A+C+E+D+Q ; OPTION RSIDUALS; OPTION NO_OUTPUT END The following MX script lines were read for group 3 TITLE G3: DZ DATA DATA NINPUT=15 RECTANGULAR FILE=GRID56.REC Rectangular continuous data read initiated Note: Maximum ordinal/rectangular record length is: 1000 NOTE: Rectangular file contained 392 records with data that contained a total of 4868 observations LABELS PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ARC1 ARC2 ARC3 ARC4 ARC5 SELECT PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ARC1 ARC2 ARC3 ARC4 ARC5 ; DEFINITON PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ; NOTE: Selection yields 392 data vectors for analysis NOTE: Vectors contain a total of 4868 observations Note: Global variable previously defined. Updating PAB Note: Global variable previously defined. Updating PAC Note: Global variable previously defined. Updating PAD Note: Global variable previously defined. Updating PAE Note: Global variable previously defined. Updating PBC Note: Global variable previously defined. Updating PBD Note: Global variable previously defined. Updating PBE Note: Global variable previously defined. Updating PCD Note: Global variable previously defined. Updating PCE Note: Global variable previously defined. Updating PDE NOTE: Definition yields 392 data vectors for analysis NOTE: Vectors contain a total of 948 observations BEGIN MATRICES ; A COMPUTED NVAR NVAR =A1 C COMPUTED NVAR NVAR =C1 D COMPUTED NVAR NVAR =D1 E COMPUTED NVAR NVAR =E1 Q COMPUTED NVAR NVAR =Q1 M FULL 1 1 =M1 H FULL 1 1 =H1 I FULL 1 1 =I1 P FULL 1 1 R FULL 1 1 S FULL 1 1 T FULL 1 1 U FULL 1 1 V FULL 1 1 W FULL 1 1 X FULL 1 1 Y FULL 1 1 Z FULL 1 1 END MATRICES; SPECIFY P PAB SPECIFY R PAC SPECIFY S PAD SPECIFY T PAE SPECIFY U PBC SPECIFY V PBD SPECIFY W PBE SPECIFY X PCD SPECIFY Y PCE SPECIFY Z PDE MEANS M|M|M|M|M ; COVARIANCE A+C+E+D+Q | H@A+C+I@D+P@Q | H@A+C+I@D+R@Q | H@A+C+I@D+S@Q | H@A+C+I@D+T@Q _ H@A+C+I@D+P@Q | A+C+E+D+Q | H@A+C+I@D+U@Q | H@A+C+I@D+V@Q | H@A+C+I@D+W@Q _ H@A+C+I@D+R@Q | H@A+C+I@D+U@Q | A+C+E+D+Q | H@A+C+I@D+X@Q | H@A+C+I@D+Y@Q _ H@A+C+I@D+S@Q | H@A+C+I@D+V@Q | H@A+C+I@D+X@Q | A+C+E+D+Q | H@A+C+I@D+Z@Q _ H@A+C+I@D+T@Q | H@A+C+I@D+W@Q | H@A+C+I@D+Y@Q | H@A+C+I@D+Z@Q | A+C+E+D+Q ; OPTION NO_OUTPUT END The following MX script lines were read for group 4 TITLE G4: STANDARDIZATION CALCULATION BEGIN MATRICES = GROUP 1; END MATRICES; BEGIN ALGEBRA; V=A+C+E+D+Q; ! TOTAL VARIANCE P=A|C|E|D|Q; ! CONCATENATE PARAMETER ESTIMATES S=P@V~; ! STANDARDIZED PARAMETER ESTIMATES END ALGEBRA; OPTION NDECIMALS=4 OPTION SUB=3176.157,1194 ! LIKELIHOOD AND DF FROM SATURATED MODEL: OZBMIYFSAT.MXO !DROP L 1 1 1 END Summary of VL file data for group 2 ARC1 ARC2 Code 1.0000 2.0000 Number 125.0000 125.0000 Mean -0.0090 -0.0190 Variance 0.8861 0.8757 Minimum -2.3713 -2.9867 Maximum 2.6268 2.2563 Summary of VL file data for group 3 PDE PCE PCD PBE PBD PBC PAE Code -10.0000 -9.0000 -8.0000 -7.0000 -6.0000 -5.0000 -4.0000 Number 392.0000 392.0000 392.0000 392.0000 392.0000 392.0000 392.0000 Mean 0.5025 0.5000 0.4895 0.5012 0.4928 0.5082 0.4988 Variance 0.0012 0.0013 0.0103 0.0006 0.0088 0.0339 0.0006 Minimum 0.5000 0.0030 0.0000 0.4960 0.0000 0.0000 0.0080 Maximum 0.9960 0.9930 0.9960 0.9950 0.9990 0.9990 0.5190 PAD PAC PAB ARC1 ARC2 ARC3 ARC4 Code -3.0000 -2.0000 -1.0000 1.0000 2.0000 3.0000 4.0000 Number 392.0000 392.0000 392.0000 392.0000 392.0000 132.0000 29.0000 Mean 0.4993 0.4985 0.5084 0.0844 0.0470 0.2260 0.2899 Variance 0.0092 0.0333 0.1005 0.9724 1.0117 1.1402 1.0366 Minimum 0.0040 0.0000 0.0000 -2.0600 -2.2000 -1.9900 -1.8400 Maximum 0.9990 0.9980 0.9990 3.1560 3.4310 3.5340 2.7000 ARC5 Code 5.0000 Number 3.0000 Mean -0.0407 Variance 0.5965 Minimum -0.9500 Maximum 0.9380 PARAMETER SPECIFICATIONS 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 I This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a FULL matrix of order 1 by 1 1 QTL 4 MATRIX M This is a FULL matrix of order 1 by 1 1 1 5 MATRIX P This is a computed FULL matrix of order 1 by 5 It has no free parameters specified MATRIX Q This is a computed 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 5 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 = 5 MX PARAMETER ESTIMATES GROUP NUMBER: 4 TITLE G4: STANDARDIZATION MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.6591 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.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.2392 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX I This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX L This is a FULL matrix of order 1 by 1 1 QTL 0.2968 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0638 MATRIX P This is a computed FULL matrix of order 1 by 5 [=A|C|E|D|Q] 1 2 3 4 5 1 0.6591 0.0000 0.2392 0.0000 0.0881 MATRIX Q This is a computed FULL matrix of order 1 by 1 [=L*L'] 1 1 0.0881 MATRIX S This is a computed FULL matrix of order 1 by 5 [=P@V~] 1 2 3 4 5 1 0.6682 0.0000 0.2425 0.0000 0.0893 MATRIX V This is a computed FULL matrix of order 1 by 1 [=A+C+E+D+Q] 1 1 0.9863 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 DOM_GEN 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 ADD_GEN 0.8118 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.4890 Your model has 5 estimated parameters and 1198 Observed statistics -2 times log-likelihood of data >>> 3175.079 Degrees of freedom >>>>>>>>>>>>>>>> 1193 Akaike's Information Criterion >>>> 789.079 Bayesian Information Criterion >>>> -2139.418 Sample size Adjusted BIC >>>> -246.012 Deviance Information Criterion >>>> -1043.124 Submodel fit* >>>>>>>>>>> 3176.157 Submodel df* >>>>>>>>>>> 1194 Difference Chi-squared >>>>>>>> 1.078 Difference d.f. >>>>>>>>>>>>>>> 1 Probability >>>>>>>>>>>>>>>>>>>> 0.299 Akaike's Information Criterion > -0.922 * Submodel statistic supplied by user This problem used 0.0% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.14 Execution 0: 0: 0:10.42 TOTAL 0: 0: 0:10.56 Total number of warnings issued: 2 ______________________________________________________________________________ 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 L FULL NVAR 1 FREE ! QTL GENETIC PATH, Q H FULL 1 1 ! SCALAR, 0.5 I FULL 1 1 ! SCALAR, 0.25 M FULL 1 1 FREE END MATRICES; LABEL ROW X ADD_GEN LABEL ROW Y COMM_ENV LABEL ROW Z SPEC_ENV LABEL ROW W DOM_GEN LABEL ROW L QTL MATRIX H .5 MATRIX I .25 START .6 ALL ! START ALL FREE PARAMETERS AT .6 START 0 M 1 1 1 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 Q= L*L'; ! QTL VARIANCE, Q^2 END ALGEBRA; OPTION NO_OUTPUT END The following MX script lines were read for group 2 TITLE G2: MZ DATA DATA NINPUT=2 RECTANGULAR FILE=MZDAT.REC Rectangular continuous data read initiated Note: Maximum ordinal/rectangular record length is: 1000 NOTE: Rectangular file contained 125 records with data that contained a total of 250 observations LABELS ARC1 ARC2 SELECT ARC1 ARC2 ; BEGIN MATRICES = GROUP 1; ! COPY ALL MATRICES FROM GROUP 1, BOTH DECLARED AND COMPUTED NOTE: Selection yields 125 data vectors for analysis NOTE: Vectors contain a total of 250 observations MEANS M|M ; ! 1X2 VECTOR OF EXPECTED MEANS COVARIANCE ! 2X2 COVARIANCE MATRIX A+C+E+D+Q | A+C+D+Q _ A+C+D+Q | A+C+E+D+Q ; OPTION RSIDUALS; OPTION NO_OUTPUT END The following MX script lines were read for group 3 TITLE G3: DZ DATA DATA NINPUT=15 RECTANGULAR FILE=GRID57.REC Rectangular continuous data read initiated Note: Maximum ordinal/rectangular record length is: 1000 NOTE: Rectangular file contained 392 records with data that contained a total of 4868 observations LABELS PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ARC1 ARC2 ARC3 ARC4 ARC5 SELECT PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ARC1 ARC2 ARC3 ARC4 ARC5 ; DEFINITON PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ; NOTE: Selection yields 392 data vectors for analysis NOTE: Vectors contain a total of 4868 observations Note: Global variable previously defined. Updating PAB Note: Global variable previously defined. Updating PAC Note: Global variable previously defined. Updating PAD Note: Global variable previously defined. Updating PAE Note: Global variable previously defined. Updating PBC Note: Global variable previously defined. Updating PBD Note: Global variable previously defined. Updating PBE Note: Global variable previously defined. Updating PCD Note: Global variable previously defined. Updating PCE Note: Global variable previously defined. Updating PDE NOTE: Definition yields 392 data vectors for analysis NOTE: Vectors contain a total of 948 observations BEGIN MATRICES ; A COMPUTED NVAR NVAR =A1 C COMPUTED NVAR NVAR =C1 D COMPUTED NVAR NVAR =D1 E COMPUTED NVAR NVAR =E1 Q COMPUTED NVAR NVAR =Q1 M FULL 1 1 =M1 H FULL 1 1 =H1 I FULL 1 1 =I1 P FULL 1 1 R FULL 1 1 S FULL 1 1 T FULL 1 1 U FULL 1 1 V FULL 1 1 W FULL 1 1 X FULL 1 1 Y FULL 1 1 Z FULL 1 1 END MATRICES; SPECIFY P PAB SPECIFY R PAC SPECIFY S PAD SPECIFY T PAE SPECIFY U PBC SPECIFY V PBD SPECIFY W PBE SPECIFY X PCD SPECIFY Y PCE SPECIFY Z PDE MEANS M|M|M|M|M ; COVARIANCE A+C+E+D+Q | H@A+C+I@D+P@Q | H@A+C+I@D+R@Q | H@A+C+I@D+S@Q | H@A+C+I@D+T@Q _ H@A+C+I@D+P@Q | A+C+E+D+Q | H@A+C+I@D+U@Q | H@A+C+I@D+V@Q | H@A+C+I@D+W@Q _ H@A+C+I@D+R@Q | H@A+C+I@D+U@Q | A+C+E+D+Q | H@A+C+I@D+X@Q | H@A+C+I@D+Y@Q _ H@A+C+I@D+S@Q | H@A+C+I@D+V@Q | H@A+C+I@D+X@Q | A+C+E+D+Q | H@A+C+I@D+Z@Q _ H@A+C+I@D+T@Q | H@A+C+I@D+W@Q | H@A+C+I@D+Y@Q | H@A+C+I@D+Z@Q | A+C+E+D+Q ; OPTION NO_OUTPUT END The following MX script lines were read for group 4 TITLE G4: STANDARDIZATION CALCULATION BEGIN MATRICES = GROUP 1; END MATRICES; BEGIN ALGEBRA; V=A+C+E+D+Q; ! TOTAL VARIANCE P=A|C|E|D|Q; ! CONCATENATE PARAMETER ESTIMATES S=P@V~; ! STANDARDIZED PARAMETER ESTIMATES END ALGEBRA; OPTION NDECIMALS=4 OPTION SUB=3176.157,1194 ! LIKELIHOOD AND DF FROM SATURATED MODEL: OZBMIYFSAT.MXO !DROP L 1 1 1 END Summary of VL file data for group 2 ARC1 ARC2 Code 1.0000 2.0000 Number 125.0000 125.0000 Mean -0.0090 -0.0190 Variance 0.8861 0.8757 Minimum -2.3713 -2.9867 Maximum 2.6268 2.2563 Summary of VL file data for group 3 PDE PCE PCD PBE PBD PBC PAE Code -10.0000 -9.0000 -8.0000 -7.0000 -6.0000 -5.0000 -4.0000 Number 392.0000 392.0000 392.0000 392.0000 392.0000 392.0000 392.0000 Mean 0.5022 0.4999 0.4905 0.5008 0.4931 0.5074 0.4993 Variance 0.0009 0.0009 0.0076 0.0005 0.0082 0.0253 0.0004 Minimum 0.5000 0.0540 0.0030 0.4420 0.0030 0.0030 0.1250 Maximum 0.9500 0.9040 0.9530 0.9340 0.9960 0.9960 0.5650 PAD PAC PAB ARC1 ARC2 ARC3 ARC4 Code -3.0000 -2.0000 -1.0000 1.0000 2.0000 3.0000 4.0000 Number 392.0000 392.0000 392.0000 392.0000 392.0000 132.0000 29.0000 Mean 0.5003 0.5007 0.5008 0.0844 0.0470 0.2260 0.2899 Variance 0.0077 0.0256 0.0763 0.9724 1.0117 1.1402 1.0366 Minimum 0.0580 0.0030 0.0030 -2.0600 -2.2000 -1.9900 -1.8400 Maximum 0.9960 0.9900 0.9960 3.1560 3.4310 3.5340 2.7000 ARC5 Code 5.0000 Number 3.0000 Mean -0.0407 Variance 0.5965 Minimum -0.9500 Maximum 0.9380 PARAMETER SPECIFICATIONS 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 I This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a FULL matrix of order 1 by 1 1 QTL 4 MATRIX M This is a FULL matrix of order 1 by 1 1 1 5 MATRIX P This is a computed FULL matrix of order 1 by 5 It has no free parameters specified MATRIX Q This is a computed 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 5 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 = 5 MX PARAMETER ESTIMATES GROUP NUMBER: 4 TITLE G4: STANDARDIZATION MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.6158 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.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.2398 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX I This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX L This is a FULL matrix of order 1 by 1 1 QTL 0.3609 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0627 MATRIX P This is a computed FULL matrix of order 1 by 5 [=A|C|E|D|Q] 1 2 3 4 5 1 0.6158 0.0000 0.2398 0.0000 0.1302 MATRIX Q This is a computed FULL matrix of order 1 by 1 [=L*L'] 1 1 0.1302 MATRIX S This is a computed FULL matrix of order 1 by 5 [=P@V~] 1 2 3 4 5 1 0.6246 0.0000 0.2433 0.0000 0.1321 MATRIX V This is a computed FULL matrix of order 1 by 1 [=A+C+E+D+Q] 1 1 0.9858 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 DOM_GEN 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 ADD_GEN 0.7847 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.4897 Your model has 5 estimated parameters and 1198 Observed statistics -2 times log-likelihood of data >>> 3174.419 Degrees of freedom >>>>>>>>>>>>>>>> 1193 Akaike's Information Criterion >>>> 788.419 Bayesian Information Criterion >>>> -2139.748 Sample size Adjusted BIC >>>> -246.342 Deviance Information Criterion >>>> -1043.454 Submodel fit* >>>>>>>>>>> 3176.157 Submodel df* >>>>>>>>>>> 1194 Difference Chi-squared >>>>>>>> 1.738 Difference d.f. >>>>>>>>>>>>>>> 1 Probability >>>>>>>>>>>>>>>>>>>> 0.187 Akaike's Information Criterion > -0.262 * Submodel statistic supplied by user This problem used 0.0% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.15 Execution 0: 0: 0: 9.42 TOTAL 0: 0: 0: 9.57 Total number of warnings issued: 2 ______________________________________________________________________________ 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 L FULL NVAR 1 FREE ! QTL GENETIC PATH, Q H FULL 1 1 ! SCALAR, 0.5 I FULL 1 1 ! SCALAR, 0.25 M FULL 1 1 FREE END MATRICES; LABEL ROW X ADD_GEN LABEL ROW Y COMM_ENV LABEL ROW Z SPEC_ENV LABEL ROW W DOM_GEN LABEL ROW L QTL MATRIX H .5 MATRIX I .25 START .6 ALL ! START ALL FREE PARAMETERS AT .6 START 0 M 1 1 1 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 Q= L*L'; ! QTL VARIANCE, Q^2 END ALGEBRA; OPTION NO_OUTPUT END The following MX script lines were read for group 2 TITLE G2: MZ DATA DATA NINPUT=2 RECTANGULAR FILE=MZDAT.REC Rectangular continuous data read initiated Note: Maximum ordinal/rectangular record length is: 1000 NOTE: Rectangular file contained 125 records with data that contained a total of 250 observations LABELS ARC1 ARC2 SELECT ARC1 ARC2 ; BEGIN MATRICES = GROUP 1; ! COPY ALL MATRICES FROM GROUP 1, BOTH DECLARED AND COMPUTED NOTE: Selection yields 125 data vectors for analysis NOTE: Vectors contain a total of 250 observations MEANS M|M ; ! 1X2 VECTOR OF EXPECTED MEANS COVARIANCE ! 2X2 COVARIANCE MATRIX A+C+E+D+Q | A+C+D+Q _ A+C+D+Q | A+C+E+D+Q ; OPTION RSIDUALS; OPTION NO_OUTPUT END The following MX script lines were read for group 3 TITLE G3: DZ DATA DATA NINPUT=15 RECTANGULAR FILE=GRID58.REC Rectangular continuous data read initiated Note: Maximum ordinal/rectangular record length is: 1000 NOTE: Rectangular file contained 392 records with data that contained a total of 4868 observations LABELS PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ARC1 ARC2 ARC3 ARC4 ARC5 SELECT PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ARC1 ARC2 ARC3 ARC4 ARC5 ; DEFINITON PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ; NOTE: Selection yields 392 data vectors for analysis NOTE: Vectors contain a total of 4868 observations Note: Global variable previously defined. Updating PAB Note: Global variable previously defined. Updating PAC Note: Global variable previously defined. Updating PAD Note: Global variable previously defined. Updating PAE Note: Global variable previously defined. Updating PBC Note: Global variable previously defined. Updating PBD Note: Global variable previously defined. Updating PBE Note: Global variable previously defined. Updating PCD Note: Global variable previously defined. Updating PCE Note: Global variable previously defined. Updating PDE NOTE: Definition yields 392 data vectors for analysis NOTE: Vectors contain a total of 948 observations BEGIN MATRICES ; A COMPUTED NVAR NVAR =A1 C COMPUTED NVAR NVAR =C1 D COMPUTED NVAR NVAR =D1 E COMPUTED NVAR NVAR =E1 Q COMPUTED NVAR NVAR =Q1 M FULL 1 1 =M1 H FULL 1 1 =H1 I FULL 1 1 =I1 P FULL 1 1 R FULL 1 1 S FULL 1 1 T FULL 1 1 U FULL 1 1 V FULL 1 1 W FULL 1 1 X FULL 1 1 Y FULL 1 1 Z FULL 1 1 END MATRICES; SPECIFY P PAB SPECIFY R PAC SPECIFY S PAD SPECIFY T PAE SPECIFY U PBC SPECIFY V PBD SPECIFY W PBE SPECIFY X PCD SPECIFY Y PCE SPECIFY Z PDE MEANS M|M|M|M|M ; COVARIANCE A+C+E+D+Q | H@A+C+I@D+P@Q | H@A+C+I@D+R@Q | H@A+C+I@D+S@Q | H@A+C+I@D+T@Q _ H@A+C+I@D+P@Q | A+C+E+D+Q | H@A+C+I@D+U@Q | H@A+C+I@D+V@Q | H@A+C+I@D+W@Q _ H@A+C+I@D+R@Q | H@A+C+I@D+U@Q | A+C+E+D+Q | H@A+C+I@D+X@Q | H@A+C+I@D+Y@Q _ H@A+C+I@D+S@Q | H@A+C+I@D+V@Q | H@A+C+I@D+X@Q | A+C+E+D+Q | H@A+C+I@D+Z@Q _ H@A+C+I@D+T@Q | H@A+C+I@D+W@Q | H@A+C+I@D+Y@Q | H@A+C+I@D+Z@Q | A+C+E+D+Q ; OPTION NO_OUTPUT END The following MX script lines were read for group 4 TITLE G4: STANDARDIZATION CALCULATION BEGIN MATRICES = GROUP 1; END MATRICES; BEGIN ALGEBRA; V=A+C+E+D+Q; ! TOTAL VARIANCE P=A|C|E|D|Q; ! CONCATENATE PARAMETER ESTIMATES S=P@V~; ! STANDARDIZED PARAMETER ESTIMATES END ALGEBRA; OPTION NDECIMALS=4 OPTION SUB=3176.157,1194 ! LIKELIHOOD AND DF FROM SATURATED MODEL: OZBMIYFSAT.MXO !DROP L 1 1 1 END Summary of VL file data for group 2 ARC1 ARC2 Code 1.0000 2.0000 Number 125.0000 125.0000 Mean -0.0090 -0.0190 Variance 0.8861 0.8757 Minimum -2.3713 -2.9867 Maximum 2.6268 2.2563 Summary of VL file data for group 3 PDE PCE PCD PBE PBD PBC PAE Code -10.0000 -9.0000 -8.0000 -7.0000 -6.0000 -5.0000 -4.0000 Number 392.0000 392.0000 392.0000 392.0000 392.0000 392.0000 392.0000 Mean 0.5019 0.4998 0.4917 0.5006 0.4936 0.5064 0.4997 Variance 0.0007 0.0007 0.0060 0.0004 0.0075 0.0204 0.0003 Minimum 0.5000 0.1090 0.0360 0.4190 0.0360 0.0360 0.2090 Maximum 0.8980 0.8300 0.9020 0.8780 0.9630 0.9630 0.5860 PAD PAC PAB ARC1 ARC2 ARC3 ARC4 Code -3.0000 -2.0000 -1.0000 1.0000 2.0000 3.0000 4.0000 Number 392.0000 392.0000 392.0000 392.0000 392.0000 132.0000 29.0000 Mean 0.5012 0.5019 0.4964 0.0844 0.0470 0.2260 0.2899 Variance 0.0067 0.0214 0.0639 0.9724 1.0117 1.1402 1.0366 Minimum 0.0360 0.0360 0.0360 -2.0600 -2.2000 -1.9900 -1.8400 Maximum 0.9630 0.9630 0.9630 3.1560 3.4310 3.5340 2.7000 ARC5 Code 5.0000 Number 3.0000 Mean -0.0407 Variance 0.5965 Minimum -0.9500 Maximum 0.9380 PARAMETER SPECIFICATIONS 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 I This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a FULL matrix of order 1 by 1 1 QTL 4 MATRIX M This is a FULL matrix of order 1 by 1 1 1 5 MATRIX P This is a computed FULL matrix of order 1 by 5 It has no free parameters specified MATRIX Q This is a computed 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 5 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 = 5 MX PARAMETER ESTIMATES GROUP NUMBER: 4 TITLE G4: STANDARDIZATION MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.5859 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.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.2402 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX I This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX L This is a FULL matrix of order 1 by 1 1 QTL 0.3992 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0620 MATRIX P This is a computed FULL matrix of order 1 by 5 [=A|C|E|D|Q] 1 2 3 4 5 1 0.5859 0.0000 0.2402 0.0000 0.1594 MATRIX Q This is a computed FULL matrix of order 1 by 1 [=L*L'] 1 1 0.1594 MATRIX S This is a computed FULL matrix of order 1 by 5 [=P@V~] 1 2 3 4 5 1 0.5945 0.0000 0.2438 0.0000 0.1617 MATRIX V This is a computed FULL matrix of order 1 by 1 [=A+C+E+D+Q] 1 1 0.9855 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 DOM_GEN 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 ADD_GEN 0.7654 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.4901 Your model has 5 estimated parameters and 1198 Observed statistics -2 times log-likelihood of data >>> 3174.064 Degrees of freedom >>>>>>>>>>>>>>>> 1193 Akaike's Information Criterion >>>> 788.064 Bayesian Information Criterion >>>> -2139.926 Sample size Adjusted BIC >>>> -246.520 Deviance Information Criterion >>>> -1043.632 Submodel fit* >>>>>>>>>>> 3176.157 Submodel df* >>>>>>>>>>> 1194 Difference Chi-squared >>>>>>>> 2.093 Difference d.f. >>>>>>>>>>>>>>> 1 Probability >>>>>>>>>>>>>>>>>>>> 0.148 Akaike's Information Criterion > 0.093 * Submodel statistic supplied by user This problem used 0.0% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.13 Execution 0: 0: 0: 9.57 TOTAL 0: 0: 0: 9.70 Total number of warnings issued: 2 ______________________________________________________________________________ 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 L FULL NVAR 1 FREE ! QTL GENETIC PATH, Q H FULL 1 1 ! SCALAR, 0.5 I FULL 1 1 ! SCALAR, 0.25 M FULL 1 1 FREE END MATRICES; LABEL ROW X ADD_GEN LABEL ROW Y COMM_ENV LABEL ROW Z SPEC_ENV LABEL ROW W DOM_GEN LABEL ROW L QTL MATRIX H .5 MATRIX I .25 START .6 ALL ! START ALL FREE PARAMETERS AT .6 START 0 M 1 1 1 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 Q= L*L'; ! QTL VARIANCE, Q^2 END ALGEBRA; OPTION NO_OUTPUT END The following MX script lines were read for group 2 TITLE G2: MZ DATA DATA NINPUT=2 RECTANGULAR FILE=MZDAT.REC Rectangular continuous data read initiated Note: Maximum ordinal/rectangular record length is: 1000 NOTE: Rectangular file contained 125 records with data that contained a total of 250 observations LABELS ARC1 ARC2 SELECT ARC1 ARC2 ; BEGIN MATRICES = GROUP 1; ! COPY ALL MATRICES FROM GROUP 1, BOTH DECLARED AND COMPUTED NOTE: Selection yields 125 data vectors for analysis NOTE: Vectors contain a total of 250 observations MEANS M|M ; ! 1X2 VECTOR OF EXPECTED MEANS COVARIANCE ! 2X2 COVARIANCE MATRIX A+C+E+D+Q | A+C+D+Q _ A+C+D+Q | A+C+E+D+Q ; OPTION RSIDUALS; OPTION NO_OUTPUT END The following MX script lines were read for group 3 TITLE G3: DZ DATA DATA NINPUT=15 RECTANGULAR FILE=GRID59.REC Rectangular continuous data read initiated Note: Maximum ordinal/rectangular record length is: 1000 NOTE: Rectangular file contained 392 records with data that contained a total of 4868 observations LABELS PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ARC1 ARC2 ARC3 ARC4 ARC5 SELECT PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ARC1 ARC2 ARC3 ARC4 ARC5 ; DEFINITON PAB PAC PAD PAE PBC PBD PBE PCD PCE PDE ; NOTE: Selection yields 392 data vectors for analysis NOTE: Vectors contain a total of 4868 observations Note: Global variable previously defined. Updating PAB Note: Global variable previously defined. Updating PAC Note: Global variable previously defined. Updating PAD Note: Global variable previously defined. Updating PAE Note: Global variable previously defined. Updating PBC Note: Global variable previously defined. Updating PBD Note: Global variable previously defined. Updating PBE Note: Global variable previously defined. Updating PCD Note: Global variable previously defined. Updating PCE Note: Global variable previously defined. Updating PDE NOTE: Definition yields 392 data vectors for analysis NOTE: Vectors contain a total of 948 observations BEGIN MATRICES ; A COMPUTED NVAR NVAR =A1 C COMPUTED NVAR NVAR =C1 D COMPUTED NVAR NVAR =D1 E COMPUTED NVAR NVAR =E1 Q COMPUTED NVAR NVAR =Q1 M FULL 1 1 =M1 H FULL 1 1 =H1 I FULL 1 1 =I1 P FULL 1 1 R FULL 1 1 S FULL 1 1 T FULL 1 1 U FULL 1 1 V FULL 1 1 W FULL 1 1 X FULL 1 1 Y FULL 1 1 Z FULL 1 1 END MATRICES; SPECIFY P PAB SPECIFY R PAC SPECIFY S PAD SPECIFY T PAE SPECIFY U PBC SPECIFY V PBD SPECIFY W PBE SPECIFY X PCD SPECIFY Y PCE SPECIFY Z PDE MEANS M|M|M|M|M ; COVARIANCE A+C+E+D+Q | H@A+C+I@D+P@Q | H@A+C+I@D+R@Q | H@A+C+I@D+S@Q | H@A+C+I@D+T@Q _ H@A+C+I@D+P@Q | A+C+E+D+Q | H@A+C+I@D+U@Q | H@A+C+I@D+V@Q | H@A+C+I@D+W@Q _ H@A+C+I@D+R@Q | H@A+C+I@D+U@Q | A+C+E+D+Q | H@A+C+I@D+X@Q | H@A+C+I@D+Y@Q _ H@A+C+I@D+S@Q | H@A+C+I@D+V@Q | H@A+C+I@D+X@Q | A+C+E+D+Q | H@A+C+I@D+Z@Q _ H@A+C+I@D+T@Q | H@A+C+I@D+W@Q | H@A+C+I@D+Y@Q | H@A+C+I@D+Z@Q | A+C+E+D+Q ; OPTION NO_OUTPUT END The following MX script lines were read for group 4 TITLE G4: STANDARDIZATION CALCULATION BEGIN MATRICES = GROUP 1; END MATRICES; BEGIN ALGEBRA; V=A+C+E+D+Q; ! TOTAL VARIANCE P=A|C|E|D|Q; ! CONCATENATE PARAMETER ESTIMATES S=P@V~; ! STANDARDIZED PARAMETER ESTIMATES END ALGEBRA; OPTION NDECIMALS=4 OPTION SUB=3176.157,1194 ! LIKELIHOOD AND DF FROM SATURATED MODEL: OZBMIYFSAT.MXO !DROP L 1 1 1 END Summary of VL file data for group 2 ARC1 ARC2 Code 1.0000 2.0000 Number 125.0000 125.0000 Mean -0.0090 -0.0190 Variance 0.8861 0.8757 Minimum -2.3713 -2.9867 Maximum 2.6268 2.2563 Summary of VL file data for group 3 PDE PCE PCD PBE PBD PBC PAE Code -10.0000 -9.0000 -8.0000 -7.0000 -6.0000 -5.0000 -4.0000 Number 392.0000 392.0000 392.0000 392.0000 392.0000 392.0000 392.0000 Mean 0.5012 0.5000 0.5007 0.4988 0.5030 0.4910 0.5000 Variance 0.0005 0.0011 0.0047 0.0005 0.0066 0.0327 0.0011 Minimum 0.5000 0.0400 0.0400 0.0400 0.0400 0.0400 0.0400 Maximum 0.9590 0.9590 0.9590 0.5000 0.9590 0.9590 0.9590 PAD PAC PAB ARC1 ARC2 ARC3 ARC4 Code -3.0000 -2.0000 -1.0000 1.0000 2.0000 3.0000 4.0000 Number 392.0000 392.0000 392.0000 392.0000 392.0000 132.0000 29.0000 Mean 0.4980 0.5162 0.4812 0.0844 0.0470 0.2260 0.2899 Variance 0.0072 0.0292 0.0813 0.9724 1.0117 1.1402 1.0366 Minimum 0.0400 0.0400 0.0400 -2.0600 -2.2000 -1.9900 -1.8400 Maximum 0.9590 0.9590 0.9590 3.1560 3.4310 3.5340 2.7000 ARC5 Code 5.0000 Number 3.0000 Mean -0.0407 Variance 0.5965 Minimum -0.9500 Maximum 0.9380 PARAMETER SPECIFICATIONS 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 I This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a FULL matrix of order 1 by 1 1 QTL 4 MATRIX M This is a FULL matrix of order 1 by 1 1 1 5 MATRIX P This is a computed FULL matrix of order 1 by 5 It has no free parameters specified MATRIX Q This is a computed 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 5 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 = 5 MX PARAMETER ESTIMATES GROUP NUMBER: 4 TITLE G4: STANDARDIZATION MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.4889 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.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.2405 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX I This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX L This is a FULL matrix of order 1 by 1 1 QTL 0.5043 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0625 MATRIX P This is a computed FULL matrix of order 1 by 5 [=A|C|E|D|Q] 1 2 3 4 5 1 0.4889 0.0000 0.2405 0.0000 0.2543 MATRIX Q This is a computed FULL matrix of order 1 by 1 [=L*L'] 1 1 0.2543 MATRIX S This is a computed FULL matrix of order 1 by 5 [=P@V~] 1 2 3 4 5 1 0.4970 0.0000 0.2445 0.0000 0.2585 MATRIX V This is a computed FULL matrix of order 1 by 1 [=A+C+E+D+Q] 1 1 0.9837 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 DOM_GEN 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 ADD_GEN 0.6992 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.4904 Your model has 5 estimated parameters and 1198 Observed statistics -2 times log-likelihood of data >>> 3168.059 Degrees of freedom >>>>>>>>>>>>>>>> 1193 Akaike's Information Criterion >>>> 782.059 Bayesian Information Criterion >>>> -2142.928 Sample size Adjusted BIC >>>> -249.522 Deviance Information Criterion >>>> -1046.634 Submodel fit* >>>>>>>>>>> 3176.157 Submodel df* >>>>>>>>>>> 1194 Difference Chi-squared >>>>>>>> 8.098 Difference d.f. >>>>>>>>>>>>>>> 1 Probability >>>>>>>>>>>>>>>>>>>> 0.004 Akaike's Information Criterion > 6.098 * Submodel statistic supplied by user This problem used 0.0% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.14 Execution 0: 0: 0: 9.28 TOTAL 0: 0: 0: 9.42 Total number of warnings issued: 2 ______________________________________________________________________________ Note: #Loop section ends here ______________________________________________________________________________