! Fitting univariate models to raw data - rectangular data files ! ! univariateraw_withcovariates.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 FR ! ADDITIVE GENETIC PATH X LO 1 1 ! NON-ADDITIVE GENETIC PATH (FIXED IN THIS APPLICATION) Y LO 1 1 FI ! 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.6 W(1,1) ST 0.0 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 1 1 1 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 It has no free parameters specified 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 1 1 1 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 It has no free parameters specified 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 1 1 1 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 It has no free parameters specified MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 2 Mx starting optimization; number of parameters = 6 *** 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: 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.0196 MATRIX B This is a FULL matrix of order 1 by 6 1 2 3 4 5 6 1 -0.1011 -1.4016 1.2690 -0.1011 -1.4016 1.2690 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 [=X*X'] 1 1 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.0191 MATRIX I This is a computed FULL matrix of order 1 by 1 [=N~*A] 1 1 0.5070 MATRIX J This is a computed FULL matrix of order 1 by 1 [=N~*C] 1 1 0.0000 MATRIX L This is a computed FULL matrix of order 1 by 1 [=N~*E] 1 1 0.4930 MATRIX M This is a FULL matrix of order 1 by 2 1 2 1 0.6701 0.6701 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.1400 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.0000 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.1381 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.0196 MATRIX B This is a FULL matrix of order 1 by 6 1 2 3 4 5 6 1 -0.1011 -1.4016 1.2690 -0.1011 -1.4016 1.2690 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 [=X*X'] 1 1 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.0191 MATRIX I This is a computed FULL matrix of order 1 by 1 [=N~*A] 1 1 0.5070 MATRIX J This is a computed FULL matrix of order 1 by 1 [=N~*C] 1 1 0.0000 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.4930 MATRIX M This is a FULL matrix of order 1 by 2 1 2 1 0.6701 0.6701 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.1400 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.0000 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.1381 Vector of OBSERVED means EPQLT1 EPQLT2 Mean 0.5119 0.5117 Vector of EXPECTED means EPQLT1 EPQLT2 Mean 0.6040 0.6040 EXPECTED COVARIANCE MATRIX EPQLT1 EPQLT2 EPQLT1 0.0387 EPQLT2 0.0196 0.0387 Function value of this group: -1397.7346 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.0196 MATRIX B This is a FULL matrix of order 1 by 6 1 2 3 4 5 6 1 -0.1011 -1.4016 1.2690 -0.1011 -1.4016 1.2690 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 [=X*X'] 1 1 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.0191 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.5070 MATRIX J This is a computed FULL matrix of order 1 by 1 [=N~*C] 1 1 0.0000 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.4930 MATRIX M This is a FULL matrix of order 1 by 2 1 2 1 0.6701 0.6701 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.1400 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.0000 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.1381 Vector of OBSERVED means EPQLT1 EPQLT2 Mean 0.5301 0.5147 Vector of EXPECTED means EPQLT1 EPQLT2 Mean 0.5935 0.5935 EXPECTED COVARIANCE MATRIX EPQLT1 EPQLT2 EPQLT1 0.0387 EPQLT2 0.0098 0.0387 Function value of this group: -716.6877 Where the fit function is -2 * Log-likelihood of raw data *** WARNING! *** Minimization may not be successful. See above CODE GREEN - it probably was OK Your model has 6 estimated parameters and 4093 Observed statistics -2 times log-likelihood of data >>> -2114.422 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.00 Execution 0: 0: 0:28.96 TOTAL 0: 0: 0:44.96 Total number of warnings issued: 2 ______________________________________________________________________________