! Fitting univariate models to raw data - rectangular data files ! ! univariateraw_withcovariates.mx ! ! Whitfield alcohol challenge male data, ADH1B genotype ! and alcohol dependence symptom count ! The following MX script lines were read for group 1 ANALYSIS OF ALCOHOL DEPENDENCE SYMPTOM COUNT - UNIVARIATE GENETIC MODEL ! Incorporates ADH1B genotype information, but NO control for stratification 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 MEANS (INTERCEPTS) B FU 1 2 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 102 ST 0.88 M(1,1) ST 0.25 Z(1,1) ST 0.25 W(1,1) ST 0.25 Y(1,1) ST -0.4 B(1,1) BO 0.01 10 Z(1,1) INTERVAL B(1,1) M(1,1) I(1,1) J(1,1) L(1,1) END; The following MX script lines were read for group 2 MZ MALE TWIN PAIRS DA NI=6 RECTANGULAR FI=WHITMZ.REC Rectangular continuous data read initiated NOTE: Rectangular file contained 39 records with data that contained a total of 228 observations LABELS OBS ALCSXT1 ALCSXT2 GTYPT1 GTYPT2 ZYG DEFINITION_VARIABLES OBS GTYPT1 GTYPT2 ZYG / NOTE: Definition yields 39 data vectors for analysis NOTE: Vectors contain a total of 72 observations BEGIN MATRICES = GROUP 1 ; K FU 2 2 ! MATRIX FOR DEFINITION VARIABLES END MATRICES; SP K -2 0 0 -3 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 MALE DZ PAIRS (BOTH TWINS GENOTYPED) DA NI=6 RECTANGULAR FI=WHITDZ2.REC Rectangular continuous data read initiated NOTE: Rectangular file contained 28 records with data that contained a total of 163 observations LABELS OBS ALCSXT1 ALCSXT2 GTYPT1 GTYPT2 ZYG DEFINITION_VARIABLES OBS GTYPT1 GTYPT2 ZYG / Note: Global variable previously defined. Updating OBS Note: Global variable previously defined. Updating GTYPT1 Note: Global variable previously defined. Updating GTYPT2 Note: Global variable previously defined. Updating ZYG NOTE: Definition yields 28 data vectors for analysis NOTE: Vectors contain a total of 51 observations BEGIN MATRICES = GROUP 1 ; K FU 2 2 G DI 1 1 H DI 1 1 END MATRICES; SP K -2 0 0 -3 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 The following MX script lines were read for group 4 DZ MALES WITH COTWIN NOT GENOTYPED, PLUS MALES FROM UNLIKE-SEX TWIN PAIRS DA NI=4 RECTANGULAR FI=WHITDZ1.REC Rectangular continuous data read initiated NOTE: Rectangular file contained 26 records with data that contained a total of 101 observations LABELS OBS ALCSX GTYP ZYG DEFINITION_VARIABLES OBS GTYP ZYG / Note: Global variable previously defined. Updating OBS Note: Global variable previously defined. Updating ZYG NOTE: Definition yields 23 data vectors for analysis NOTE: Vectors contain a total of 23 observations BEGIN MATRICES; V CO 1 1 = N1 ! VARIANCE (FROM GROUP 1) M FU 1 1 FR ! MEAN (INTERCEPTS) B FU 1 1 FR ! VECTOR OF REGRESSION COEFFICIENTS K FU 1 1 ! MATRIX FOR DEFINITION VARIABLES END MATRICES; SP M 101 SP B 102 SP K -2 ! st 0 m(1,1) ! st 0.20 b(1,1) CO V / ME M+B*K/ OPTIONS RS END Summary of VL file data for group 2 ZYG GTYPT2 GTYPT1 OBS ALCSXT1 ALCSXT2 Code -4.0000 -3.0000 -2.0000 -1.0000 1.0000 2.0000 Number 39.0000 39.0000 39.0000 39.0000 36.0000 36.0000 Mean 2.0000 0.1282 0.1282 20.0000 0.8744 0.7520 Variance 0.0000 0.1118 0.1118 126.6667 0.3426 0.3253 Minimum 2.0000 0.0000 0.0000 1.0000 0.0000 0.0000 Maximum 2.0000 1.0000 1.0000 39.0000 1.7917 1.6094 Summary of VL file data for group 3 ZYG GTYPT2 GTYPT1 OBS ALCSXT1 Code -4.0000E+00 -3.0000E+00 -2.0000E+00 -1.0000E+00 1.0000E+00 Number 2.8000E+01 2.8000E+01 2.8000E+01 2.8000E+01 2.7000E+01 Mean 4.0000E+00 1.0714E-01 1.0714E-01 1.4500E+01 9.5810E-01 Variance 4.0602E-15 9.5663E-02 9.5663E-02 6.5250E+01 2.7801E-01 Minimum 4.0000E+00 0.0000E+00 0.0000E+00 1.0000E+00 0.0000E+00 Maximum 4.0000E+00 1.0000E+00 1.0000E+00 2.8000E+01 1.7917E+00 ALCSXT2 Code 2.0000E+00 Number 2.4000E+01 Mean 8.1523E-01 Variance 3.3178E-01 Minimum 0.0000E+00 Maximum 1.6094E+00 Summary of VL file data for group 4 ZYG GTYP OBS ALCSX Code -3.0000 -2.0000 -1.0000 1.0000 Number 23.0000 23.0000 23.0000 23.0000 Mean 4.7826 0.0435 28.5652 0.7273 Variance 0.1701 0.0416 1184.7675 0.4812 Minimum 4.0000 0.0000 1.0000 0.0000 Maximum 5.0000 1.0000 95.0000 2.0791 PARAMETER SPECIFICATIONS GROUP NUMBER: 1 Analysis of alcohol dependence symptom count - 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 2 1 2 1 102 102 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 male 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 2 1 2 1 102 102 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 2 by 2 1 2 1 -2 0 2 0 -3 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 MALE DZ PAIRS (both twins genotyped) 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 2 1 2 1 102 102 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 2 by 2 1 2 1 -2 0 2 0 -3 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 DZ males with cotwin not genotyped, plus males from unlike-sex twin pairs MATRIX B This is a FULL matrix of order 1 by 1 1 1 102 MATRIX K This is a FULL matrix of order 1 by 1 1 1 -2 MATRIX M This is a FULL matrix of order 1 by 1 1 1 101 MATRIX V This is a computed FULL matrix of order 1 by 1 It has no free parameters specified Mx starting optimization; number of parameters = 5 MX PARAMETER ESTIMATES GROUP NUMBER: 1 Analysis of alcohol dependence symptom count - univariate genetic model MATRIX A This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 1 0.1761 MATRIX B This is a FULL matrix of order 1 by 2 1 2 1 -0.4812 -0.4812 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 1.1183E-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.1559 MATRIX I This is a computed FULL matrix of order 1 by 1 [=N~*A] 1 1 0.5304 MATRIX J This is a computed FULL matrix of order 1 by 1 [=N~*C] 1 1 3.3677E-13 MATRIX L This is a computed FULL matrix of order 1 by 1 [=N~*E] 1 1 0.4696 MATRIX M This is a FULL matrix of order 1 by 2 1 2 1 0.8873 0.8873 MATRIX N This is a computed FULL matrix of order 1 by 1 [=A+D+C+E] 1 1 0.3321 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.4197 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 3.3441E-07 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.3949 GROUP NUMBER: 2 MZ male twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 1 0.1761 MATRIX B This is a FULL matrix of order 1 by 2 1 2 1 -0.4812 -0.4812 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 1.1183E-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.1559 MATRIX I This is a computed FULL matrix of order 1 by 1 [=N~*A] 1 1 0.5304 MATRIX J This is a computed FULL matrix of order 1 by 1 [=N~*C] 1 1 3.3677E-13 MATRIX K This is a FULL matrix of order 2 by 2 1 2 1 0.0000 0.0000 2 0.0000 0.0000 MATRIX L This is a computed FULL matrix of order 1 by 1 [=N~*E] 1 1 0.4696 MATRIX M This is a FULL matrix of order 1 by 2 1 2 1 0.8873 0.8873 MATRIX N This is a computed FULL matrix of order 1 by 1 [=A+D+C+E] 1 1 0.3321 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.4197 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 3.3441E-07 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.3949 Vector of OBSERVED means ALCSXT1 ALCSXT2 Mean 0.8744 0.7520 Vector of EXPECTED means ALCSXT1 ALCSXT2 Mean 0.8873 0.8873 EXPECTED COVARIANCE MATRIX ALCSXT1 ALCSXT2 ALCSXT1 0.3321 ALCSXT2 0.1761 0.3321 Function value of this group: 107.3958 Where the fit function is -2 * Log-likelihood of raw data GROUP NUMBER: 3 MALE DZ PAIRS (both twins genotyped) MATRIX A This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 1 0.1761 MATRIX B This is a FULL matrix of order 1 by 2 1 2 1 -0.4812 -0.4812 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 1.1183E-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.1559 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.5304 MATRIX J This is a computed FULL matrix of order 1 by 1 [=N~*C] 1 1 3.3677E-13 MATRIX K This is a FULL matrix of order 2 by 2 1 2 1 0.0000 0.0000 2 0.0000 0.0000 MATRIX L This is a computed FULL matrix of order 1 by 1 [=N~*E] 1 1 0.4696 MATRIX M This is a FULL matrix of order 1 by 2 1 2 1 0.8873 0.8873 MATRIX N This is a computed FULL matrix of order 1 by 1 [=A+D+C+E] 1 1 0.3321 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.4197 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 3.3441E-07 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.3949 Vector of OBSERVED means ALCSXT1 ALCSXT2 Mean 0.9581 0.8152 Vector of EXPECTED means ALCSXT1 ALCSXT2 Mean 0.8873 0.8873 EXPECTED COVARIANCE MATRIX ALCSXT1 ALCSXT2 ALCSXT1 0.3321 ALCSXT2 0.0881 0.3321 Function value of this group: 83.2562 Where the fit function is -2 * Log-likelihood of raw data GROUP NUMBER: 4 DZ males with cotwin not genotyped, plus males from unlike-sex twin pairs MATRIX B This is a FULL matrix of order 1 by 1 1 1 -0.4812 MATRIX K This is a FULL matrix of order 1 by 1 1 1 0.0000 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.8873 MATRIX V This is a computed FULL matrix of order 1 by 1 1 1 0.3321 Vector of OBSERVED means ALCSX Mean 0.7273 Vector of EXPECTED means ALCSX Mean 0.8873 EXPECTED COVARIANCE MATRIX ALCSX ALCSX 0.3321 Function value of this group: 50.1449 Where the fit function is -2 * Log-likelihood of raw data Your model has 5 estimated parameters and 146 Observed statistics -2 times log-likelihood of data >>> 240.797 Degrees of freedom >>>>>>>>>>>>>>>> 141 5 Confidence intervals requested in group 1 Matrix Element Int. Estimate Lower Upper Lfail Ufail B 1 1 1 95.0 -0.4812 -0.8182 -0.1438 0 0 0 0 M 1 1 1 95.0 0.8873 0.7741 1.0006 0 0 0 0 I 1 1 1 95.0 0.5304 0.0000 0.7153 0 0 0 0 J 1 1 1 95.0 0.0000 0.0000 0.5498 0 0 0 0 L 1 1 1 95.0 0.4696 0.2847 0.7603 0 0 0 0 This problem used 2.0% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.59 Execution 0: 0: 0:10.71 TOTAL 0: 0: 0:11.30 Total number of warnings issued: 0 ______________________________________________________________________________