! Fitting univariate models to raw bivariate data - rectangular data files ! ! ordinalraw_withcovariates.mx ! ! The following MX script lines were read for group 1 ANALYSIS OF SOCIAL NON-CONFORMITY DATA - UNIVARIATE GENETIC MODEL DA CALC NG=4 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 FR ! SHARED ENVIRONMENTAL PATH Z LO 1 1 FR ! WITHIN-FAMILY ENVIRONMENTAL PATH M FU 1 2 FR ! VECTOR OF THRESHOLDS (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 0 0 102 0 0 ST 1.75 M(1,1) ST 0.6 Z(1,1) ST 0.8 W(1,1) ST 0.1 Y(1,1) ST 0.05 B(1,1) ST -0.00 B(1,2) ST -0.00 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 ORDINAL FI=ALCMZF.REC Ordinal data read initiated NOTE: Rectangular file contained 1011 records with data that contained a total of 7001 observations LABELS OBS ALCT1 ALCT2 AGE AGESQ AGECB ZYG DEFINITION_VARIABLES OBS AGE AGESQ AGECB ZYG / NOTE: Definition yields 1011 data vectors for analysis NOTE: Vectors contain a total of 1946 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) / TH M-B*K/ OPTIONS RS END The following MX script lines were read for group 3 FEMALE DZ PAIRS DA NI=7 ORDINAL FI=ALCDZF.REC Ordinal data read initiated NOTE: Rectangular file contained 592 records with data that contained a total of 4086 observations LABELS OBS ALCT1 ALCT2 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 592 data vectors for analysis NOTE: Vectors contain a total of 1126 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) / TH M-B*K/ MAT G 0.5 MAT H 0.25 OPTIONS RS END The following MX script lines were read for group 4 CONSTRAINT FUNCTION - CONSTRAIN PHENOTYPE VARIANCES TO UNITY CO NI=1 BEGIN MATRICES=GROUP 1; U UNIT 1 1 END MATRICES; CO \D2V(N) = U; 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.0110E+03 1.0110E+03 1.0110E+03 1.0110E+03 1.0110E+03 Mean 1.0000E+00 5.9503E+01 1.3445E+01 3.4314E+00 5.0600E+02 Variance 1.1245E-16 5.1042E+03 1.1048E+02 1.6706E+00 8.5177E+04 Minimum 1.0000E+00 4.9130E+00 2.8900E+00 1.7000E+00 1.0000E+00 Maximum 1.0000E+00 4.9304E+02 6.2410E+01 7.9000E+00 1.0110E+03 ALCT1 ALCT2 Code 1.0000E+00 2.0000E+00 Number 1.0110E+03 9.3500E+02 Mean 5.8358E-02 4.7059E-02 Variance 5.4952E-02 4.4844E-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 5.9200E+02 5.9200E+02 5.9200E+02 5.9200E+02 5.9200E+02 Mean 3.0000E+00 5.6704E+01 1.3072E+01 3.3905E+00 2.9650E+02 Variance 0.0000E+00 4.4793E+03 1.0050E+02 1.5759E+00 2.9205E+04 Minimum 3.0000E+00 4.9130E+00 2.8900E+00 1.7000E+00 1.0000E+00 Maximum 3.0000E+00 4.7455E+02 6.0840E+01 7.8000E+00 5.9200E+02 ALCT1 ALCT2 Code 1.0000E+00 2.0000E+00 Number 5.9200E+02 5.3400E+02 Mean 5.0676E-02 8.8015E-02 Variance 4.8108E-02 8.0268E-02 Minimum 0.0000E+00 0.0000E+00 Maximum 1.0000E+00 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 0 0 102 0 0 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 1 1 2 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 3 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 0 0 102 0 0 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 1 1 2 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 3 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 0 0 102 0 0 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 1 1 2 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 3 GROUP NUMBER: 4 CONSTRAINT FUNCTION - constrain phenotype variances to unity 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 0 0 102 0 0 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 U This is a UNIT matrix of order 1 by 1 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 1 1 2 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 3 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: 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.5382 MATRIX B This is a FULL matrix of order 1 by 6 1 2 3 4 5 6 1 -0.2407 0.0000 0.0000 -0.2407 0.0000 0.0000 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 2.0878E-13 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.4618 MATRIX I This is a computed FULL matrix of order 1 by 1 [=N~*A] 1 1 0.5382 MATRIX J This is a computed FULL matrix of order 1 by 1 [=N~*C] 1 1 2.0878E-13 MATRIX L This is a computed FULL matrix of order 1 by 1 [=N~*E] 1 1 0.4618 MATRIX M This is a FULL matrix of order 1 by 2 1 2 1 0.8092 0.8092 MATRIX N This is a computed FULL matrix of order 1 by 1 [=A+D+C+E] 1 1 1.0000 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.7336 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 -4.5693E-07 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.6796 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.5382 MATRIX B This is a FULL matrix of order 1 by 6 1 2 3 4 5 6 1 -0.2407 0.0000 0.0000 -0.2407 0.0000 0.0000 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 2.0878E-13 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.4618 MATRIX I This is a computed FULL matrix of order 1 by 1 [=N~*A] 1 1 0.5382 MATRIX J This is a computed FULL matrix of order 1 by 1 [=N~*C] 1 1 2.0878E-13 MATRIX K This is a FULL matrix of order 6 by 2 1 2 1 2.0000 0.0000 2 4.0000 0.0000 3 8.0000 0.0000 4 0.0000 2.0000 5 0.0000 4.0000 6 0.0000 8.0000 MATRIX L This is a computed FULL matrix of order 1 by 1 [=N~*E] 1 1 0.4618 MATRIX M This is a FULL matrix of order 1 by 2 1 2 1 0.8092 0.8092 MATRIX N This is a computed FULL matrix of order 1 by 1 [=A+D+C+E] 1 1 1.0000 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.7336 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 -4.5693E-07 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.6796 Matrix of EXPECTED thresholds ALCT1 ALCT2 Threshold 1 1.2906 1.2906 EXPECTED COVARIANCE MATRIX ALCT1 ALCT2 ALCT1 1.0000 ALCT2 0.5382 1.0000 Function value of this group: 749.6041 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 3 FEMALE DZ PAIRS MATRIX A This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 1 0.5382 MATRIX B This is a FULL matrix of order 1 by 6 1 2 3 4 5 6 1 -0.2407 0.0000 0.0000 -0.2407 0.0000 0.0000 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 2.0878E-13 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.4618 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.5382 MATRIX J This is a computed FULL matrix of order 1 by 1 [=N~*C] 1 1 2.0878E-13 MATRIX K This is a FULL matrix of order 6 by 2 1 2 1 2.8000 0.0000 2 7.8400 0.0000 3 21.9520 0.0000 4 0.0000 2.8000 5 0.0000 7.8400 6 0.0000 21.9520 MATRIX L This is a computed FULL matrix of order 1 by 1 [=N~*E] 1 1 0.4618 MATRIX M This is a FULL matrix of order 1 by 2 1 2 1 0.8092 0.8092 MATRIX N This is a computed FULL matrix of order 1 by 1 [=A+D+C+E] 1 1 1.0000 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.7336 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 -4.5693E-07 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.6796 Matrix of EXPECTED thresholds ALCT1 ALCT2 Threshold 1 1.4831 1.4831 EXPECTED COVARIANCE MATRIX ALCT1 ALCT2 ALCT1 1.0000 ALCT2 0.2691 1.0000 Function value of this group: 542.0682 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 4 CONSTRAINT FUNCTION - constrain phenotype variances to unity MATRIX A This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 1 0.5382 MATRIX B This is a FULL matrix of order 1 by 6 1 2 3 4 5 6 1 -0.2407 0.0000 0.0000 -0.2407 0.0000 0.0000 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 2.0878E-13 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.4618 MATRIX I This is a computed FULL matrix of order 1 by 1 [=N~*A] 1 1 0.5382 MATRIX J This is a computed FULL matrix of order 1 by 1 [=N~*C] 1 1 2.0878E-13 MATRIX L This is a computed FULL matrix of order 1 by 1 [=N~*E] 1 1 0.4618 MATRIX M This is a FULL matrix of order 1 by 2 1 2 1 0.8092 0.8092 MATRIX N This is a computed FULL matrix of order 1 by 1 [=A+D+C+E] 1 1 1.0000 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.7336 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 -4.5693E-07 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.6796 *** WARNING! *** Minimization may not be successful. See above CODE GREEN - it probably was OK Your model has 5 estimated parameters and 3073 Observed statistics Observed statistics include 1 constraints. -2 times log-likelihood of data >>> 1291.672 Degrees of freedom >>>>>>>>>>>>>>>> 3068 This problem used 14.4% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 9.32 Execution 0: 0: 0: 4.12 TOTAL 0: 0: 0:13.44 Total number of warnings issued: 2 ______________________________________________________________________________