! Fitting univariate models to raw data - rectangular data files ! ! univariateraw_withcovariates_ce.mx ! ! The following MX script lines were read for group 1 ANALYSIS OF SOCIAL NON-CONFORMITY DATA - UNIVARIATE GENETIC MODEL DA CALC NG=3 BEGIN MATRICES; W LO 1 1 FI ! ADDITIVE GENETIC PATH X LO 1 1 ! NON-ADDITIVE GENETIC PATH (FIXED IN THIS APPLICATION) Y LO 1 1 FR ! SHARED ENVIRONMENTAL PATH Z LO 1 1 FR ! WITHIN-FAMILY ENVIRONMENTAL PATH M FU 1 2 FR ! VECTOR OF MEANS (INTERCEPTS) B FU 1 6 FR ! VECTOR OF REGRESSION COEFFICIENTS END MATRICES; BEGIN ALGEBRA; A=W*W'; D=X*X'; C=Y*Y'; E=Z*Z'; N=A+D+C+E; I=N~*A; J=N~*C; L=N~*E; END ALGEBRA; SP M 101 101 SP B 102 103 104 102 103 104 ST 0.75 M(1,1) ST 0.6 Z(1,1) ST 0.0 W(1,1) ST 0.6 Y(1,1) ST 0.5 B(1,1) ST -0.05 B(1,2) ST -0.05 B(1,3) BO 0.01 10 Z(1,1) END; The following MX script lines were read for group 2 MZ FEMALE TWIN PAIRS DA NI=7 RECTANGULAR FI=EPQLMZF.REC Rectangular continuous data read initiated NOTE: Rectangular file contained 1287 records with data that contained a total of 8954 observations LABELS OBS EPQLT1 EPQLT2 AGE AGESQ AGECB ZYG DEFINITION_VARIABLES OBS AGE AGESQ AGECB ZYG / NOTE: Definition yields 1287 data vectors for analysis NOTE: Vectors contain a total of 2519 observations BEGIN MATRICES = GROUP 1 ; K FU 6 2 ! MATRIX FOR DEFINITION VARIABLES END MATRICES; SP K -2 0 -3 0 -4 0 0 -2 0 -3 0 -4 CO (A+D+C+E | A+D+C_ A+D+C | A+D+C+E) / ME M+B*K/ OPTIONS RS END The following MX script lines were read for group 3 FEMALE DZ PAIRS DA NI=7 RECTANGULAR FI=EPQLDZF.REC Rectangular continuous data read initiated NOTE: Rectangular file contained 828 records with data that contained a total of 5714 observations LABELS OBS EPQLT1 EPQLT2 AGE AGESQ AGECB ZYG DEFINITION_VARIABLES OBS AGE AGESQ AGECB ZYG/ Note: Global variable previously defined. Updating OBS Note: Global variable previously defined. Updating AGE Note: Global variable previously defined. Updating AGESQ Note: Global variable previously defined. Updating AGECB Note: Global variable previously defined. Updating ZYG BEGIN MATRICES = GROUP 1 ; NOTE: Definition yields 828 data vectors for analysis NOTE: Vectors contain a total of 1574 observations K FU 6 2 G DI 1 1 H DI 1 1 END MATRICES; SP K -2 0 -3 0 -4 0 0 -2 0 -3 0 -4 CO (A+D+C+E | G*A+H*D+C_ G*A+H*D+C | A+D+C+E) / ME M+B*K/ MAT G 0.5 MAT H 0.25 OPTIONS RS END Summary of VL file data for group 2 ZYG AGECB AGESQ AGE OBS Code -5.0000E+00 -4.0000E+00 -3.0000E+00 -2.0000E+00 -1.0000E+00 Number 1.2870E+03 1.2870E+03 1.2870E+03 1.2870E+03 1.2870E+03 Mean 1.0000E+00 6.8901E-02 1.4668E-01 3.5595E-01 6.5600E+02 Variance 0.0000E+00 7.3719E-03 1.4418E-02 1.9980E-02 1.3803E+05 Minimum 1.0000E+00 4.9100E-03 2.8900E-02 1.7000E-01 1.3000E+01 Maximum 1.0000E+00 5.3144E-01 6.5610E-01 8.1000E-01 1.2990E+03 EPQLT1 EPQLT2 Code 1.0000E+00 2.0000E+00 Number 1.2870E+03 1.2320E+03 Mean 5.1192E-01 5.1168E-01 Variance 4.3913E-02 4.4467E-02 Minimum 0.0000E+00 0.0000E+00 Maximum 1.0000E+00 1.0000E+00 Summary of VL file data for group 3 ZYG AGECB AGESQ AGE OBS Code -5.0000E+00 -4.0000E+00 -3.0000E+00 -2.0000E+00 -1.0000E+00 Number 8.2800E+02 8.2800E+02 8.2800E+02 8.2800E+02 8.2800E+02 Mean 3.0000E+00 6.8228E-02 1.4479E-01 3.5295E-01 4.2850E+02 Variance 1.0984E-15 8.0761E-03 1.5057E-02 2.0215E-02 5.7132E+04 Minimum 3.0000E+00 4.9100E-03 2.8900E-02 1.7000E-01 1.5000E+01 Maximum 3.0000E+00 5.3144E-01 6.5610E-01 8.1000E-01 8.4200E+02 EPQLT1 EPQLT2 Code 1.0000E+00 2.0000E+00 Number 8.2700E+02 7.4700E+02 Mean 5.3008E-01 5.1466E-01 Variance 4.0900E-02 5.0174E-02 Minimum 4.7620E-02 0.0000E+00 Maximum 9.5227E-01 1.0000E+00 PARAMETER SPECIFICATIONS GROUP NUMBER: 1 Analysis of Social Non-conformity data - univariate genetic model MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX B This is a FULL matrix of order 1 by 6 1 2 3 4 5 6 1 102 103 104 102 103 104 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 I This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX J This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX M This is a FULL matrix of order 1 by 2 1 2 1 101 101 MATRIX N 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 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 1 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 2 GROUP NUMBER: 2 MZ female twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX B This is a FULL matrix of order 1 by 6 1 2 3 4 5 6 1 102 103 104 102 103 104 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 I This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX J This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX K This is a FULL matrix of order 6 by 2 1 2 1 -2 0 2 -3 0 3 -4 0 4 0 -2 5 0 -3 6 0 -4 MATRIX L This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX M This is a FULL matrix of order 1 by 2 1 2 1 101 101 MATRIX N 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 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 1 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 2 GROUP NUMBER: 3 FEMALE DZ PAIRS MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX B This is a FULL matrix of order 1 by 6 1 2 3 4 5 6 1 102 103 104 102 103 104 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 G This is a DIAGONAL matrix of order 1 by 1 It has no free parameters specified MATRIX H This is a DIAGONAL matrix of order 1 by 1 It has no free parameters specified MATRIX I This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX J This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX K This is a FULL matrix of order 6 by 2 1 2 1 -2 0 2 -3 0 3 -4 0 4 0 -2 5 0 -3 6 0 -4 MATRIX L This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX M This is a FULL matrix of order 1 by 2 1 2 1 101 101 MATRIX N 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 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 1 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 2 Mx starting optimization; number of parameters = 6 MX PARAMETER ESTIMATES GROUP NUMBER: 1 Analysis of Social Non-conformity data - univariate genetic model MATRIX A This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 1 0.0000 MATRIX B This is a FULL matrix of order 1 by 6 1 2 3 4 5 6 1 -0.0563 -1.5000 1.3401 -0.0563 -1.5000 1.3401 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.0162 MATRIX D This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.0225 MATRIX I This is a computed FULL matrix of order 1 by 1 [=N~*A] 1 1 0.0000 MATRIX J This is a computed FULL matrix of order 1 by 1 [=N~*C] 1 1 0.4188 MATRIX L This is a computed FULL matrix of order 1 by 1 [=N~*E] 1 1 0.5812 MATRIX M This is a FULL matrix of order 1 by 2 1 2 1 0.6634 0.6634 MATRIX N This is a computed FULL matrix of order 1 by 1 [=A+D+C+E] 1 1 0.0387 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 -0.1273 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.1500 GROUP NUMBER: 2 MZ female twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 1 0.0000 MATRIX B This is a FULL matrix of order 1 by 6 1 2 3 4 5 6 1 -0.0563 -1.5000 1.3401 -0.0563 -1.5000 1.3401 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.0162 MATRIX D This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.0225 MATRIX I This is a computed FULL matrix of order 1 by 1 [=N~*A] 1 1 0.0000 MATRIX J This is a computed FULL matrix of order 1 by 1 [=N~*C] 1 1 0.4188 MATRIX K This is a FULL matrix of order 6 by 2 1 2 1 0.2000 0.0000 2 0.0400 0.0000 3 0.0080 0.0000 4 0.0000 0.2000 5 0.0000 0.0400 6 0.0000 0.0080 MATRIX L This is a computed FULL matrix of order 1 by 1 [=N~*E] 1 1 0.5812 MATRIX M This is a FULL matrix of order 1 by 2 1 2 1 0.6634 0.6634 MATRIX N This is a computed FULL matrix of order 1 by 1 [=A+D+C+E] 1 1 0.0387 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 -0.1273 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.1500 Vector of OBSERVED means EPQLT1 EPQLT2 Mean 0.5119 0.5117 Vector of EXPECTED means EPQLT1 EPQLT2 Mean 0.6028 0.6028 EXPECTED COVARIANCE MATRIX EPQLT1 EPQLT2 EPQLT1 0.0387 EPQLT2 0.0162 0.0387 Function value of this group: -1381.5453 Where the fit function is -2 * Log-likelihood of raw data GROUP NUMBER: 3 FEMALE DZ PAIRS MATRIX A This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 1 0.0000 MATRIX B This is a FULL matrix of order 1 by 6 1 2 3 4 5 6 1 -0.0563 -1.5000 1.3401 -0.0563 -1.5000 1.3401 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.0162 MATRIX D This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.0225 MATRIX G This is a DIAGONAL matrix of order 1 by 1 1 1 0.5000 MATRIX H This is a DIAGONAL matrix of order 1 by 1 1 1 0.2500 MATRIX I This is a computed FULL matrix of order 1 by 1 [=N~*A] 1 1 0.0000 MATRIX J This is a computed FULL matrix of order 1 by 1 [=N~*C] 1 1 0.4188 MATRIX K This is a FULL matrix of order 6 by 2 1 2 1 0.2200 0.0000 2 0.0484 0.0000 3 0.0106 0.0000 4 0.0000 0.2200 5 0.0000 0.0484 6 0.0000 0.0106 MATRIX L This is a computed FULL matrix of order 1 by 1 [=N~*E] 1 1 0.5812 MATRIX M This is a FULL matrix of order 1 by 2 1 2 1 0.6634 0.6634 MATRIX N This is a computed FULL matrix of order 1 by 1 [=A+D+C+E] 1 1 0.0387 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 -0.1273 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.1500 Vector of OBSERVED means EPQLT1 EPQLT2 Mean 0.5301 0.5147 Vector of EXPECTED means EPQLT1 EPQLT2 Mean 0.5927 0.5927 EXPECTED COVARIANCE MATRIX EPQLT1 EPQLT2 EPQLT1 0.0387 EPQLT2 0.0162 0.0387 Function value of this group: -694.4331 Where the fit function is -2 * Log-likelihood of raw data Your model has 6 estimated parameters and 4093 Observed statistics -2 times log-likelihood of data >>> -2075.978 Degrees of freedom >>>>>>>>>>>>>>>> 4087 This problem used 18.6% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0:16.23 Execution 0: 0: 0:34.12 TOTAL 0: 0: 0:50.35 Total number of warnings issued: 0 ______________________________________________________________________________