** Mx startup successful ** **MX-PC 1.54** Job started on 03/03/03 at 18:47:43 ! Boulder 2003 ! multivariate Mx script (cholesky decomposition) for the analysis of three continuous phenotypes ! using raw data input file ! MZM, DZM, MZF, DZF, DOSMF, DOSFM twins ! Evaluated models: ! ACE males, ACE females, DOS A-correlation = DZ A-correlation ! Mx homepage: http://views.vcu.edu/mx The following MX script lines were read for group 1 #DEFINE NVAR 3 ! MX READS 1. THIS CORRESPONDS TO ONE PHENOTYPE TO BE ANALYSED G1: CALCULATION GROUP DATA CALC NGROUPS=8 BEGIN MATRICES; X LOWER NVAR NVAR FREE ! MALES ADDITIVE GENETIC PATH COEFFICIENT (STANDARD DEVIATION) Y LOWER NVAR NVAR FIX ! FREE ! ASSUME NO COMMON ENVIRONMENT ! MALES COMMON ENVIRONMENTAL PATH COEFFICIENT (STANDARD DEVIATION) Z LOWER NVAR NVAR FREE ! MALES UNIQUE ENVIRONMENTAL PATH COEFFICIENT (STANDARD DEVIATION) S LOWER NVAR NVAR FREE ! FEMALES ADDITIVE GENETIC PATH COEFFICIENT (STANDARD DEVIATION) T LOWER NVAR NVAR FIX ! FREE ! ASSUME NO COMMON ENVIRONMENT ! FEMALES COMMON ENVIRONMENTAL PATH COEFFICIENT (STANDARD DEVIATION) U LOWER NVAR NVAR FREE ! FEMALES UNIQUE ENVIRONMENTAL PATH COEFFICIENT (STANDARD DEVIATION) H FULL 1 1 FIX ! .5 FOR DZ COV A F FULL 1 1 FIX !FREE ! .5 DOS COV A END MATRICES; MATRIX H .5 ! MATRIX H IS FIXED AND CONTAINS .5 MATRIX F .5 ! MATRIX H IS FIXED AND CONTAINS .5 BEGIN ALGEBRA; A= X*X' ; ! MALES ADDITIVE GENETIC VARIANCE C= Y*Y' ; ! MALES COMMON ENVIRONMENTAL VARIANCE E= Z*Z' ; ! MALES UNIQUE ENVIRONMENTAL VARIANCE K= S*S' ; ! FEMALES ADDITIVE GENETIC VARIANCE L= T*T' ; ! FEMALES COMMON ENVIRONMENTAL VARIANCE M= U*U' ; ! FEMALES UNIQUE ENVIRONMENTAL VARIANCE END ALGEBRA; ST 0.1 ALL ST 1.5 X 1 1 1 X 1 2 1 X 1 3 1 ! PROVIDE STARTING VALUES ! st 1.5 Y 1 1 1 Y 1 2 1 Y 1 3 1 ST 1.5 Z 1 1 1 Z 1 2 2 Z 1 3 3 ST 1.5 S 1 1 1 S 1 2 1 S 1 3 1 ! st 1.5 T 1 1 1 T 1 2 1 T 1 3 1 ST 1.5 U 1 1 1 U 1 2 2 U 1 3 3 ! st 0.5 F 1 1 1 ! Bound -0.5 0.5 F 1 1 1 END The following MX script lines were read for group 2 G2: MZM TWINS, DATAGROUP DATA NINPUT_VARS= 14 ! NR OF INPUTVARS PER FAMILY MISSING=-1.00 ! MISSING VALUES = -1.00 NOTE: Missing value *string* set to '-1.00' RECTANGULAR FILE=ALL.DAT ! READ IN RAW DATAFILE, MUST BE IN SAME DIRECTORY Rectangular continuous data read initiated NOTE: Rectangular file contained 2028 records with data that contained a total of 27296 observations LABELS ! TELL MX WHAT IS FOUND IN THE DATAFILE NIH_ID COUNTRY ZYGOT6 REAR AGEA SEXA LDLLNA APOBLNA APOELNA AGEB SEXB LDLLNB APOBLNB APOELNB SELECT IF COUNTRY = 1 ; ! SELECT DUTCH TWINS, YOUNGER SAMPLE (DORRET) NOTE: Select if yields 160 data vectors for analysis NOTE: Vectors contain a total of 2236 observations SELECT IF ZYGOT6 = 1 ; ! SELECT MZM'S NOTE: Select if yields 35 data vectors for analysis NOTE: Vectors contain a total of 486 observations SELECT LDLLNA APOBLNA APOELNA LDLLNB APOBLNB APOELNB ; ! SELECT NVAR PHENOTYPES PER TWIN MEMBER MATRICES = GROUP 1 NOTE: Selection yields 35 data vectors for analysis NOTE: Vectors contain a total of 206 observations G FULL 1 6 FREE ! MEAN PHENOTYPES END MATRICES MEANS G ; ! MODEL FOR MEANS COVARIANCES ! MODEL FOR MZM VARIANCE/COVARIANCES A+C+E | A+C _ A+C | A+C+E ; ST 15 G 1 1 - G 1 NVAR OPTIONS RSIDUALS END The following MX script lines were read for group 3 G3: DZM TWINS, DATAGROUP DATA NINPUT_VARS= 14 ! NR OF INPUTVARS PER FAMILY MISSING=-1.00 ! MISSING VALUES = -1.00 NOTE: Missing value *string* set to '-1.00' RECTANGULAR FILE=ALL.DAT ! READ IN RAW DATAFILE, MUST BE IN SAME DIRECTORY Rectangular continuous data read initiated NOTE: Rectangular file contained 2028 records with data that contained a total of 27296 observations LABELS ! TELL MX WHAT IS FOUND IN THE DATAFILE NIH_ID COUNTRY ZYGOT6 REAR AGEA SEXA LDLLNA APOBLNA APOELNA AGEB SEXB LDLLNB APOBLNB APOELNB SELECT IF COUNTRY = 1 ; ! SELECT DUTCH TWINS, YOUNGER SAMPLE (DORRET) NOTE: Select if yields 160 data vectors for analysis NOTE: Vectors contain a total of 2236 observations SELECT IF ZYGOT6 = 2 ; ! SELECT DZM'S NOTE: Select if yields 31 data vectors for analysis NOTE: Vectors contain a total of 434 observations SELECT LDLLNA APOBLNA APOELNA LDLLNB APOBLNB APOELNB ; ! SELECT NVAR PHENOTYPES PER TWIN MEMBER MATRICES = GROUP 1 NOTE: Selection yields 31 data vectors for analysis NOTE: Vectors contain a total of 186 observations G FULL 1 6 FREE ! MEAN PHENOTYPES END MATRICES; MEANS G ; ! MODEL FOR MEANS COVARIANCES ! MODEL FOR DZM VARIANCE/COVARIANCES A+C+E | H@A+C _ H@A+C | A+C+E ; ST 15 G 1 1 - G 1 NVAR OPTIONS RSIDUALS END The following MX script lines were read for group 4 G4: MZF TWINS, DATAGROUP DATA NINPUT_VARS= 14 ! NR OF INPUTVARS PER FAMILY MISSING=-1.00 ! MISSING VALUES = -1.00 NOTE: Missing value *string* set to '-1.00' RECTANGULAR FILE=ALL.DAT ! READ IN RAW DATAFILE, MUST BE IN SAME DIRECTORY Rectangular continuous data read initiated NOTE: Rectangular file contained 2028 records with data that contained a total of 27296 observations LABELS ! TELL MX WHAT IS FOUND IN THE DATAFILE NIH_ID COUNTRY ZYGOT6 REAR AGEA SEXA LDLLNA APOBLNA APOELNA AGEB SEXB LDLLNB APOBLNB APOELNB SELECT IF COUNTRY = 1 ; ! SELECT DUTCH TWINS, YOUNGER SAMPLE (DORRET) NOTE: Select if yields 160 data vectors for analysis NOTE: Vectors contain a total of 2236 observations SELECT IF ZYGOT6 = 3 ; ! SELECT MZF'S NOTE: Select if yields 35 data vectors for analysis NOTE: Vectors contain a total of 490 observations SELECT LDLLNA APOBLNA APOELNA LDLLNB APOBLNB APOELNB ; ! SELECT NVAR PHENOTYPES PER TWIN MEMBER MATRICES = GROUP 1 NOTE: Selection yields 35 data vectors for analysis NOTE: Vectors contain a total of 210 observations G FULL 1 6 FREE ! MEAN PHENOTYPES END MATRICES; MEANS G ; ! MODEL FOR MEANS, ASSUMING MEAN TWIN 1 = MEAN TWIN 2, AND DZ MEAN=MZ MEAN COVARIANCES ! MODEL FOR MZF VARIANCE/COVARIANCES K+L+M | K+L _ K+L | K+L+M ; ST 15 G 1 1 - G 1 NVAR OPTIONS RSIDUALS END The following MX script lines were read for group 5 G5: DZF TWINS, DATAGROUP DATA NINPUT_VARS= 14 ! NR OF INPUTVARS PER FAMILY MISSING=-1.00 ! MISSING VALUES = -1.00 NOTE: Missing value *string* set to '-1.00' RECTANGULAR FILE=ALL.DAT ! READ IN RAW DATAFILE, MUST BE IN SAME DIRECTORY Rectangular continuous data read initiated NOTE: Rectangular file contained 2028 records with data that contained a total of 27296 observations LABELS ! TELL MX WHAT IS FOUND IN THE DATAFILE NIH_ID COUNTRY ZYGOT6 REAR AGEA SEXA LDLLNA APOBLNA APOELNA AGEB SEXB LDLLNB APOBLNB APOELNB SELECT IF COUNTRY = 1 ; ! SELECT DUTCH TWINS, YOUNGER SAMPLE (DORRET) NOTE: Select if yields 160 data vectors for analysis NOTE: Vectors contain a total of 2236 observations SELECT IF ZYGOT6 = 4 ; ! SELECT DZF'S NOTE: Select if yields 30 data vectors for analysis NOTE: Vectors contain a total of 420 observations SELECT LDLLNA APOBLNA APOELNA LDLLNB APOBLNB APOELNB ; ! SELECT NVAR PHENOTYPES PER TWIN MEMBER, MATRICES = GROUP 1 NOTE: Selection yields 30 data vectors for analysis NOTE: Vectors contain a total of 180 observations G FULL 1 6 FREE ! MEAN PHENOTYPES END MATRICES; MEANS G ; ! MODEL FOR MEANS, COVARIANCES ! MODEL FOR DZF VARIANCE/COVARIANCES K+L+M | H@K+L _ H@K+L | K+L+M ; ST 15 G 1 1 - G 1 NVAR OPTIONS RSIDUALS END The following MX script lines were read for group 6 G6: DOS-MF TWINS, DATAGROUP DATA NINPUT_VARS= 14 ! NR OF INPUTVARS PER FAMILY MISSING=-1.00 ! MISSING VALUES = -1.00 NOTE: Missing value *string* set to '-1.00' RECTANGULAR FILE=ALL.DAT ! READ IN RAW DATAFILE, MUST BE IN SAME DIRECTORY Rectangular continuous data read initiated NOTE: Rectangular file contained 2028 records with data that contained a total of 27296 observations LABELS ! TELL MX WHAT IS FOUND IN THE DATAFILE NIH_ID COUNTRY ZYGOT6 REAR AGEA SEXA LDLLNA APOBLNA APOELNA AGEB SEXB LDLLNB APOBLNB APOELNB SELECT IF COUNTRY = 1 ; ! SELECT DUTCH TWINS, YOUNGER SAMPLE (DORRET) NOTE: Select if yields 160 data vectors for analysis NOTE: Vectors contain a total of 2236 observations SELECT IF ZYGOT6 = 5 ; ! SELECT DOSMF'S NOTE: Select if yields 16 data vectors for analysis NOTE: Vectors contain a total of 224 observations SELECT LDLLNA APOBLNA APOELNA LDLLNB APOBLNB APOELNB ; ! SELECT NVAR PHENOTYPES PER TWIN MEMBER MATRICES = GROUP 1 NOTE: Selection yields 16 data vectors for analysis NOTE: Vectors contain a total of 96 observations G FULL 1 6 FREE ! MEAN PHENOTYPES END MATRICES; MEANS G ; ! MODEL FOR MEANS, ASSUMING MEAN TWIN 1 = MEAN TWIN 2, AND DZ MEAN=MZ MEAN COVARIANCES ! MODEL FOR DOSMF VARIANCE/COVARIANCES A+C+E | F@(X*S')+(Y*T') _ F@(S*X')+(T*Y') | K+L+M ; ST 15 G 1 1 - G 1 NVAR OPTIONS RSIDUALS END The following MX script lines were read for group 7 G7: DOSFM TWINS, DATAGROUP DATA NINPUT_VARS= 14 ! NR OF INPUTVARS PER FAMILY MISSING=-1.00 ! MISSING VALUES = -1.00 NOTE: Missing value *string* set to '-1.00' RECTANGULAR FILE=ALL.DAT ! READ IN RAW DATAFILE, MUST BE IN SAME DIRECTORY Rectangular continuous data read initiated NOTE: Rectangular file contained 2028 records with data that contained a total of 27296 observations LABELS ! TELL MX WHAT IS FOUND IN THE DATAFILE NIH_ID COUNTRY ZYGOT6 REAR AGEA SEXA LDLLNA APOBLNA APOELNA AGEB SEXB LDLLNB APOBLNB APOELNB SELECT IF COUNTRY = 1 ; ! SELECT DUTCH TWINS, YOUNGER SAMPLE (DORRET) NOTE: Select if yields 160 data vectors for analysis NOTE: Vectors contain a total of 2236 observations SELECT IF ZYGOT6 = 6 ; ! SELECT DOSFM'S NOTE: Select if yields 13 data vectors for analysis NOTE: Vectors contain a total of 182 observations SELECT LDLLNA APOBLNA APOELNA LDLLNB APOBLNB APOELNB ; ! SELECT NVAR PHENOTYPES PER TWIN MEMBER MATRICES = GROUP 1 NOTE: Selection yields 13 data vectors for analysis NOTE: Vectors contain a total of 78 observations G FULL 1 6 FREE ! MEAN PHENOTYPES END MATRICES; MEANS G ; ! MODEL FOR MEANS, ASSUMING MEAN TWIN 1 = MEAN TWIN 2, AND DZ MEAN=MZ MEAN COVARIANCES ! MODEL FOR DZ VARIANCE/COVARIANCES K+L+M | F@(S*X')+(T*Y') _ F@(X*S')+(Y*T') | A+C+E ; ST 15 G 1 1 - G 1 NVAR OPTIONS RSIDUALS !Option Multiple issat ! allows comparison of submodels with present model (ACE KLM) !Interval @95 V 1 1 1 V 1 1 2 V 1 1 3 V 1 2 1 V 1 2 2 V 1 2 3 P 1 1 1 P 1 2 1 G 1 1 1 F 1 1 1 END The following MX script lines were read for group 8 G8: CALCULATE GENETIC AND ENVIRONMENTAL CORRELATIONS AND STANDARDIZED VARIANCE-COVARIANCE COMPONENTS CALCULATION BEGIN MATRICES A COMP NVAR NVAR = A1 C COMP NVAR NVAR = C1 E COMP NVAR NVAR = E1 K COMP NVAR NVAR = K1 L COMP NVAR NVAR = L1 M COMP NVAR NVAR = M1 I IDEN NVAR NVAR END MATRICES; BEGIN ALGEBRA; !men: calculation of correlations using more elaborated approach (no error messages in output) R=\SQRT(I.A)~*A*\SQRT(I.A)~ ; ! GENETIC CORRELATIONS S=\SQRT(I.C)~*C*\SQRT(I.C)~ ; ! SHARED ENVIRONMENTAL CORRELATIONS T=\SQRT(I.E)~*E*\SQRT(I.E)~ ; ! UNIQUE ENVIRONMENTAL CORRELATIONS ! women: calculation of correlations using \stnd function (mx will give error messages in output, but this is okay) J = \STND(K) ; ! STANDARDIZE ADDITIVE GENETIC VARIANCE/COVARIANCE MATRIX FEMALES TO OBTAIN GENETIC CORRELATIONS BETWEEN VARIABLES How am I supposed to take the square root of 0.000000000000000E+00 ? Diagonal elements are: 0.000000000000000E+00 0.000000000000000E+00 Matrix is 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 Matrix not fully standardized N = \STND(L) ; ! STANDARDIZE COMMON ENVIRONMENTAL VARIANCE/COVARIANCE MATRIX FEMALES TO OBTAIN COM. ENV. CORRELATIONS BETWEEN VARIABLES How am I supposed to take the square root of 0.000000000000000E+00 ? Diagonal elements are: 0.000000000000000E+00 0.000000000000000E+00 Matrix is 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 Matrix not fully standardized O = \STND(M) ; ! STANDARDIZE UNIQUE ENVIRONMENTAL VARIANCE/COVARIANCE MATRIX FEMALES TO OBTAIN UNIQUE ENV. CORRELATIONS BETWEEN VARIABLES How am I supposed to take the square root of 0.000000000000000E+00 ? Diagonal elements are: 0.000000000000000E+00 0.000000000000000E+00 Matrix is 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 Matrix not fully standardized !standardized variance and covariance components V = A%(A+C+E) | C%(A+C+E) | E%(A+C+E) _ !MALES FIRST ROWS STANDARDIZED VARIANCE/COVARIANCE COMPONENTS K%(K+L+M) | L%(K+L+M) | M%(K+L+M) ; !FEMALES SECOND ROWS STANDARDIZED VARIANCE/COVARIANCE COMPONENTS !all correlations in 1 matrix W = R | S | T _ J | N | O ; END ALGEBRA; LABEL COL W R_A R_A R_A R_C R_C R_C R_E R_E R_E LABEL ROW W LDLMEN APOBMEN APOEMEN LDLWOMEN APOBWOMEN APOEWOMEN LABEL COL V H2_LDL H2_APOB H2_APOE C2_LDL C2_APOB C2_APOE E2_LDL E2_APOB E2_APOE LABEL ROW V MEN MEN MEN WOMEN WOMEN WOMEN OPTION ITERATIONS=5000 OPTION NDECIMALS=4 ! THARD=-2 OPTION MULTIPLE OPTION ISSAT ! COMPARE FIT OF SUBMODELS TO FULL MODEL END How am I supposed to take the square root of 0.000000000000000E+00 ? Diagonal elements are: 0.000000000000000E+00 0.000000000000000E+00 Matrix is 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 Matrix not fully standardized Summary of VL file data for group 2 LDLLNA APOBLNA APOELNA LDLLNB APOBLNB APOELNB Code 1.0000 2.0000 3.0000 4.0000 5.0000 6.0000 Number 35.0000 35.0000 33.0000 35.0000 35.0000 33.0000 Mean 19.0071 16.6860 17.3400 18.8937 16.3049 17.3733 Variance 7.4983 4.5548 4.6655 6.0212 4.9993 4.1953 Minimum 13.5100 11.9200 11.2200 14.5700 11.4400 13.4000 Maximum 26.6400 22.9600 22.1400 25.9300 22.8700 22.0300 Summary of VL file data for group 3 LDLLNA APOBLNA APOELNA LDLLNB APOBLNB APOELNB Code 1.0000 2.0000 3.0000 4.0000 5.0000 6.0000 Number 31.0000 31.0000 31.0000 31.0000 31.0000 31.0000 Mean 18.8848 17.2861 16.7829 18.7997 16.9290 16.9206 Variance 2.7131 3.9938 6.7878 6.4321 4.6640 6.5336 Minimum 15.3100 13.6700 11.6900 12.0700 12.3800 9.1500 Maximum 21.6600 21.1400 22.5800 23.9900 20.8100 21.6600 Summary of VL file data for group 4 LDLLNA APOBLNA APOELNA LDLLNB APOBLNB APOELNB Code 1.0000 2.0000 3.0000 4.0000 5.0000 6.0000 Number 35.0000 35.0000 35.0000 35.0000 35.0000 35.0000 Mean 20.5186 17.4634 19.2240 20.2783 17.3709 19.2054 Variance 6.6600 7.7347 4.1145 7.4321 7.4772 4.1662 Minimum 16.1000 11.6800 14.6900 12.4700 8.7000 14.5300 Maximum 25.3000 23.3300 23.4500 24.9500 22.0900 24.8200 Summary of VL file data for group 5 LDLLNA APOBLNA APOELNA LDLLNB APOBLNB APOELNB Code 1.0000 2.0000 3.0000 4.0000 5.0000 6.0000 Number 30.0000 30.0000 30.0000 30.0000 30.0000 30.0000 Mean 19.8897 17.0177 18.8447 18.9287 16.0387 18.4547 Variance 8.5730 9.8391 3.8030 10.8826 15.1407 4.1787 Minimum 12.9300 9.5900 14.7200 9.9000 2.3600 13.9300 Maximum 24.9300 21.8800 23.2400 23.8900 22.9600 24.2300 Summary of VL file data for group 6 LDLLNA APOBLNA APOELNA LDLLNB APOBLNB APOELNB Code 1.0000 2.0000 3.0000 4.0000 5.0000 6.0000 Number 16.0000 16.0000 16.0000 16.0000 16.0000 16.0000 Mean 19.8394 17.6194 17.5769 19.3756 17.2925 17.5050 Variance 5.1732 6.2578 5.9011 4.6168 6.1375 5.7637 Minimum 16.2600 13.0400 13.2900 15.2500 12.8300 12.7100 Maximum 24.6800 21.3600 21.3000 23.1100 22.0900 20.9000 Summary of VL file data for group 7 LDLLNA APOBLNA APOELNA LDLLNB APOBLNB APOELNB Code 1.0000 2.0000 3.0000 4.0000 5.0000 6.0000 Number 13.0000 13.0000 13.0000 13.0000 13.0000 13.0000 Mean 18.3408 15.4908 19.3377 18.6462 16.2985 17.4277 Variance 9.2296 2.8979 5.7567 4.4996 3.9284 5.3906 Minimum 13.2200 12.8300 16.2100 15.2500 11.9200 13.8300 Maximum 22.9200 18.8700 23.7800 22.3000 20.1200 22.1600 PARAMETER SPECIFICATIONS GROUP NUMBER: 1 G1: calculation group MATRIX A This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX F This is a 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 K This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX L This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX S This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 13 2 14 15 3 16 17 18 MATRIX T This is a LOWER TRIANGULAR matrix of order 3 by 3 It has no free parameters specified MATRIX U This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 19 2 20 21 3 22 23 24 MATRIX X This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 1 2 2 3 3 4 5 6 MATRIX Y This is a LOWER TRIANGULAR matrix of order 3 by 3 It has no free parameters specified MATRIX Z This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 7 2 8 9 3 10 11 12 GROUP NUMBER: 2 G2: MZM twins, datagroup MATRIX A This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX F This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX G This is a FULL matrix of order 1 by 6 1 2 3 4 5 6 1 25 26 27 28 29 30 MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX K This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX L This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX S This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 13 2 14 15 3 16 17 18 MATRIX T This is a LOWER TRIANGULAR matrix of order 3 by 3 It has no free parameters specified MATRIX U This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 19 2 20 21 3 22 23 24 MATRIX X This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 1 2 2 3 3 4 5 6 MATRIX Y This is a LOWER TRIANGULAR matrix of order 3 by 3 It has no free parameters specified MATRIX Z This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 7 2 8 9 3 10 11 12 GROUP NUMBER: 3 G3: DZM twins, datagroup MATRIX A This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX F This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX G This is a FULL matrix of order 1 by 6 1 2 3 4 5 6 1 31 32 33 34 35 36 MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX K This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX L This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX S This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 13 2 14 15 3 16 17 18 MATRIX T This is a LOWER TRIANGULAR matrix of order 3 by 3 It has no free parameters specified MATRIX U This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 19 2 20 21 3 22 23 24 MATRIX X This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 1 2 2 3 3 4 5 6 MATRIX Y This is a LOWER TRIANGULAR matrix of order 3 by 3 It has no free parameters specified MATRIX Z This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 7 2 8 9 3 10 11 12 GROUP NUMBER: 4 G4: MZF twins, datagroup MATRIX A This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX F This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX G This is a FULL matrix of order 1 by 6 1 2 3 4 5 6 1 37 38 39 40 41 42 MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX K This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX L This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX S This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 13 2 14 15 3 16 17 18 MATRIX T This is a LOWER TRIANGULAR matrix of order 3 by 3 It has no free parameters specified MATRIX U This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 19 2 20 21 3 22 23 24 MATRIX X This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 1 2 2 3 3 4 5 6 MATRIX Y This is a LOWER TRIANGULAR matrix of order 3 by 3 It has no free parameters specified MATRIX Z This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 7 2 8 9 3 10 11 12 GROUP NUMBER: 5 G5: DZF twins, datagroup MATRIX A This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX F This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX G This is a FULL matrix of order 1 by 6 1 2 3 4 5 6 1 43 44 45 46 47 48 MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX K This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX L This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX S This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 13 2 14 15 3 16 17 18 MATRIX T This is a LOWER TRIANGULAR matrix of order 3 by 3 It has no free parameters specified MATRIX U This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 19 2 20 21 3 22 23 24 MATRIX X This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 1 2 2 3 3 4 5 6 MATRIX Y This is a LOWER TRIANGULAR matrix of order 3 by 3 It has no free parameters specified MATRIX Z This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 7 2 8 9 3 10 11 12 GROUP NUMBER: 6 G6: DOS-MF twins, datagroup MATRIX A This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX F This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX G This is a FULL matrix of order 1 by 6 1 2 3 4 5 6 1 49 50 51 52 53 54 MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX K This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX L This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX S This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 13 2 14 15 3 16 17 18 MATRIX T This is a LOWER TRIANGULAR matrix of order 3 by 3 It has no free parameters specified MATRIX U This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 19 2 20 21 3 22 23 24 MATRIX X This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 1 2 2 3 3 4 5 6 MATRIX Y This is a LOWER TRIANGULAR matrix of order 3 by 3 It has no free parameters specified MATRIX Z This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 7 2 8 9 3 10 11 12 GROUP NUMBER: 7 G7: DOSFM twins, datagroup MATRIX A This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX F This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX G This is a FULL matrix of order 1 by 6 1 2 3 4 5 6 1 55 56 57 58 59 60 MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX K This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX L This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX S This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 13 2 14 15 3 16 17 18 MATRIX T This is a LOWER TRIANGULAR matrix of order 3 by 3 It has no free parameters specified MATRIX U This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 19 2 20 21 3 22 23 24 MATRIX X This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 1 2 2 3 3 4 5 6 MATRIX Y This is a LOWER TRIANGULAR matrix of order 3 by 3 It has no free parameters specified MATRIX Z This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 7 2 8 9 3 10 11 12 GROUP NUMBER: 8 G8: Calculate genetic and environmental correlations and standardized variance-covariance components MATRIX A This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX I This is an IDENTITY matrix of order 3 by 3 MATRIX J This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX K This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX L This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX N This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX O This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX R This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX S This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX T This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX V This is a computed FULL matrix of order 6 by 9 It has no free parameters specified MATRIX W This is a computed FULL matrix of order 6 by 9 It has no free parameters specified Mx starting optimization; number of parameters = 60 How am I supposed to take the square root of 0.000000000000000E+00 ? Diagonal elements are: 0.000000000000000E+00 0.000000000000000E+00 Matrix is 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 Matrix not fully standardized How am I supposed to take the square root of 0.000000000000000E+00 ? Diagonal elements are: 0.000000000000000E+00 0.000000000000000E+00 Matrix is 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 Matrix not fully standardized MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1: calculation group MATRIX A This is a computed FULL matrix of order 3 by 3 [=X*X'] 1 2 3 1 4.5118 3.3173 0.5625 2 3.3173 3.7024 0.0101 3 0.5625 0.0101 5.6886 MATRIX C This is a computed FULL matrix of order 3 by 3 [=Y*Y'] 1 2 3 1 0.0000 0.0000 0.0000 2 0.0000 0.0000 0.0000 3 0.0000 0.0000 0.0000 MATRIX E This is a computed FULL matrix of order 3 by 3 [=Z*Z'] 1 2 3 1 1.0008 0.7430 0.2765 2 0.7430 0.9052 0.2429 3 0.2765 0.2429 0.8081 MATRIX F This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX K This is a computed FULL matrix of order 3 by 3 [=S*S'] 1 2 3 1 6.3672 5.6406 -0.6563 2 5.6406 7.0662 -0.3507 3 -0.6563 -0.3507 4.4656 MATRIX L This is a computed FULL matrix of order 3 by 3 [=T*T'] 1 2 3 1 0.0000 0.0000 0.0000 2 0.0000 0.0000 0.0000 3 0.0000 0.0000 0.0000 MATRIX M This is a computed FULL matrix of order 3 by 3 [=U*U'] 1 2 3 1 1.6032 1.0776 0.2917 2 1.0776 1.4860 0.0836 3 0.2917 0.0836 0.7028 MATRIX S This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 2.5233 2 2.2354 -1.4385 3 -0.2601 -0.1603 2.0910 MATRIX T This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.0000 2 0.0000 0.0000 3 0.0000 0.0000 0.0000 MATRIX U This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 1.2662 2 0.8510 0.8728 3 0.2304 -0.1289 0.7957 MATRIX X This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 2.1241 2 1.5617 -1.1240 3 0.2648 0.3590 -2.3430 MATRIX Y This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.0000 2 0.0000 0.0000 3 0.0000 0.0000 0.0000 MATRIX Z This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 1.0004 2 0.7427 0.5947 3 0.2764 0.0633 0.8531 GROUP NUMBER: 2 G2: MZM twins, datagroup MATRIX A This is a computed FULL matrix of order 3 by 3 [=X*X'] 1 2 3 1 4.5118 3.3173 0.5625 2 3.3173 3.7024 0.0101 3 0.5625 0.0101 5.6886 MATRIX C This is a computed FULL matrix of order 3 by 3 [=Y*Y'] 1 2 3 1 0.0000 0.0000 0.0000 2 0.0000 0.0000 0.0000 3 0.0000 0.0000 0.0000 MATRIX E This is a computed FULL matrix of order 3 by 3 [=Z*Z'] 1 2 3 1 1.0008 0.7430 0.2765 2 0.7430 0.9052 0.2429 3 0.2765 0.2429 0.8081 MATRIX F This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX G This is a FULL matrix of order 1 by 6 1 2 3 4 5 6 1 19.0082 16.6868 17.2824 18.8948 16.3057 17.3308 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX K This is a computed FULL matrix of order 3 by 3 [=S*S'] 1 2 3 1 6.3672 5.6406 -0.6563 2 5.6406 7.0662 -0.3507 3 -0.6563 -0.3507 4.4656 MATRIX L This is a computed FULL matrix of order 3 by 3 [=T*T'] 1 2 3 1 0.0000 0.0000 0.0000 2 0.0000 0.0000 0.0000 3 0.0000 0.0000 0.0000 MATRIX M This is a computed FULL matrix of order 3 by 3 [=U*U'] 1 2 3 1 1.6032 1.0776 0.2917 2 1.0776 1.4860 0.0836 3 0.2917 0.0836 0.7028 MATRIX S This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 2.5233 2 2.2354 -1.4385 3 -0.2601 -0.1603 2.0910 MATRIX T This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.0000 2 0.0000 0.0000 3 0.0000 0.0000 0.0000 MATRIX U This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 1.2662 2 0.8510 0.8728 3 0.2304 -0.1289 0.7957 MATRIX X This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 2.1241 2 1.5617 -1.1240 3 0.2648 0.3590 -2.3430 MATRIX Y This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.0000 2 0.0000 0.0000 3 0.0000 0.0000 0.0000 MATRIX Z This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 1.0004 2 0.7427 0.5947 3 0.2764 0.0633 0.8531 Vector of OBSERVED means LDLLNA APOBLNA APOELNA LDLLNB APOBLNB APOELNB Mean 19.0071 16.6860 17.3400 18.8937 16.3049 17.3733 Vector of EXPECTED means LDLLNA APOBLNA APOELNA LDLLNB APOBLNB APOELNB Mean 19.0082 16.6868 17.2824 18.8948 16.3057 17.3308 EXPECTED COVARIANCE MATRIX LDLLNA APOBLNA APOELNA LDLLNB APOBLNB APOELNB LDLLNA 5.5126 APOBLNA 4.0602 4.6076 APOELNA 0.8390 0.2530 6.4967 LDLLNB 4.5118 3.3173 0.5625 5.5126 APOBLNB 3.3173 3.7024 0.0101 4.0602 4.6076 APOELNB 0.5625 0.0101 5.6886 0.8390 0.2530 6.4967 Function value of this group: 733.5078 Where the fit function is -2 * Log-likelihood of raw data GROUP NUMBER: 3 G3: DZM twins, datagroup MATRIX A This is a computed FULL matrix of order 3 by 3 [=X*X'] 1 2 3 1 4.5118 3.3173 0.5625 2 3.3173 3.7024 0.0101 3 0.5625 0.0101 5.6886 MATRIX C This is a computed FULL matrix of order 3 by 3 [=Y*Y'] 1 2 3 1 0.0000 0.0000 0.0000 2 0.0000 0.0000 0.0000 3 0.0000 0.0000 0.0000 MATRIX E This is a computed FULL matrix of order 3 by 3 [=Z*Z'] 1 2 3 1 1.0008 0.7430 0.2765 2 0.7430 0.9052 0.2429 3 0.2765 0.2429 0.8081 MATRIX F This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX G This is a FULL matrix of order 1 by 6 1 2 3 4 5 6 1 18.8852 17.2864 16.7841 18.8000 16.9292 16.9220 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX K This is a computed FULL matrix of order 3 by 3 [=S*S'] 1 2 3 1 6.3672 5.6406 -0.6563 2 5.6406 7.0662 -0.3507 3 -0.6563 -0.3507 4.4656 MATRIX L This is a computed FULL matrix of order 3 by 3 [=T*T'] 1 2 3 1 0.0000 0.0000 0.0000 2 0.0000 0.0000 0.0000 3 0.0000 0.0000 0.0000 MATRIX M This is a computed FULL matrix of order 3 by 3 [=U*U'] 1 2 3 1 1.6032 1.0776 0.2917 2 1.0776 1.4860 0.0836 3 0.2917 0.0836 0.7028 MATRIX S This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 2.5233 2 2.2354 -1.4385 3 -0.2601 -0.1603 2.0910 MATRIX T This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.0000 2 0.0000 0.0000 3 0.0000 0.0000 0.0000 MATRIX U This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 1.2662 2 0.8510 0.8728 3 0.2304 -0.1289 0.7957 MATRIX X This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 2.1241 2 1.5617 -1.1240 3 0.2648 0.3590 -2.3430 MATRIX Y This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.0000 2 0.0000 0.0000 3 0.0000 0.0000 0.0000 MATRIX Z This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 1.0004 2 0.7427 0.5947 3 0.2764 0.0633 0.8531 Vector of OBSERVED means LDLLNA APOBLNA APOELNA LDLLNB APOBLNB APOELNB Mean 18.8848 17.2861 16.7829 18.7997 16.9290 16.9206 Vector of EXPECTED means LDLLNA APOBLNA APOELNA LDLLNB APOBLNB APOELNB Mean 18.8852 17.2864 16.7841 18.8000 16.9292 16.9220 EXPECTED COVARIANCE MATRIX LDLLNA APOBLNA APOELNA LDLLNB APOBLNB APOELNB LDLLNA 5.5126 APOBLNA 4.0602 4.6076 APOELNA 0.8390 0.2530 6.4967 LDLLNB 2.2559 1.6586 0.2812 5.5126 APOBLNB 1.6586 1.8512 0.0050 4.0602 4.6076 APOELNB 0.2812 0.0050 2.8443 0.8390 0.2530 6.4967 Function value of this group: 753.4575 Where the fit function is -2 * Log-likelihood of raw data GROUP NUMBER: 4 G4: MZF twins, datagroup MATRIX A This is a computed FULL matrix of order 3 by 3 [=X*X'] 1 2 3 1 4.5118 3.3173 0.5625 2 3.3173 3.7024 0.0101 3 0.5625 0.0101 5.6886 MATRIX C This is a computed FULL matrix of order 3 by 3 [=Y*Y'] 1 2 3 1 0.0000 0.0000 0.0000 2 0.0000 0.0000 0.0000 3 0.0000 0.0000 0.0000 MATRIX E This is a computed FULL matrix of order 3 by 3 [=Z*Z'] 1 2 3 1 1.0008 0.7430 0.2765 2 0.7430 0.9052 0.2429 3 0.2765 0.2429 0.8081 MATRIX F This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX G This is a FULL matrix of order 1 by 6 1 2 3 4 5 6 1 20.5175 17.4626 19.2246 20.2774 17.3702 19.2061 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX K This is a computed FULL matrix of order 3 by 3 [=S*S'] 1 2 3 1 6.3672 5.6406 -0.6563 2 5.6406 7.0662 -0.3507 3 -0.6563 -0.3507 4.4656 MATRIX L This is a computed FULL matrix of order 3 by 3 [=T*T'] 1 2 3 1 0.0000 0.0000 0.0000 2 0.0000 0.0000 0.0000 3 0.0000 0.0000 0.0000 MATRIX M This is a computed FULL matrix of order 3 by 3 [=U*U'] 1 2 3 1 1.6032 1.0776 0.2917 2 1.0776 1.4860 0.0836 3 0.2917 0.0836 0.7028 MATRIX S This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 2.5233 2 2.2354 -1.4385 3 -0.2601 -0.1603 2.0910 MATRIX T This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.0000 2 0.0000 0.0000 3 0.0000 0.0000 0.0000 MATRIX U This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 1.2662 2 0.8510 0.8728 3 0.2304 -0.1289 0.7957 MATRIX X This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 2.1241 2 1.5617 -1.1240 3 0.2648 0.3590 -2.3430 MATRIX Y This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.0000 2 0.0000 0.0000 3 0.0000 0.0000 0.0000 MATRIX Z This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 1.0004 2 0.7427 0.5947 3 0.2764 0.0633 0.8531 Vector of OBSERVED means LDLLNA APOBLNA APOELNA LDLLNB APOBLNB APOELNB Mean 20.5186 17.4634 19.2240 20.2783 17.3709 19.2054 Vector of EXPECTED means LDLLNA APOBLNA APOELNA LDLLNB APOBLNB APOELNB Mean 20.5175 17.4626 19.2246 20.2774 17.3702 19.2061 EXPECTED COVARIANCE MATRIX LDLLNA APOBLNA APOELNA LDLLNB APOBLNB APOELNB LDLLNA 7.9704 APOBLNA 6.7181 8.5521 APOELNA -0.3646 -0.2672 5.1684 LDLLNB 6.3672 5.6406 -0.6563 7.9704 APOBLNB 5.6406 7.0662 -0.3507 6.7181 8.5521 APOELNB -0.6563 -0.3507 4.4656 -0.3646 -0.2672 5.1684 Function value of this group: 808.0811 Where the fit function is -2 * Log-likelihood of raw data GROUP NUMBER: 5 G5: DZF twins, datagroup MATRIX A This is a computed FULL matrix of order 3 by 3 [=X*X'] 1 2 3 1 4.5118 3.3173 0.5625 2 3.3173 3.7024 0.0101 3 0.5625 0.0101 5.6886 MATRIX C This is a computed FULL matrix of order 3 by 3 [=Y*Y'] 1 2 3 1 0.0000 0.0000 0.0000 2 0.0000 0.0000 0.0000 3 0.0000 0.0000 0.0000 MATRIX E This is a computed FULL matrix of order 3 by 3 [=Z*Z'] 1 2 3 1 1.0008 0.7430 0.2765 2 0.7430 0.9052 0.2429 3 0.2765 0.2429 0.8081 MATRIX F This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX G This is a FULL matrix of order 1 by 6 1 2 3 4 5 6 1 19.8897 17.0179 18.8446 18.9295 16.0397 18.4547 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX K This is a computed FULL matrix of order 3 by 3 [=S*S'] 1 2 3 1 6.3672 5.6406 -0.6563 2 5.6406 7.0662 -0.3507 3 -0.6563 -0.3507 4.4656 MATRIX L This is a computed FULL matrix of order 3 by 3 [=T*T'] 1 2 3 1 0.0000 0.0000 0.0000 2 0.0000 0.0000 0.0000 3 0.0000 0.0000 0.0000 MATRIX M This is a computed FULL matrix of order 3 by 3 [=U*U'] 1 2 3 1 1.6032 1.0776 0.2917 2 1.0776 1.4860 0.0836 3 0.2917 0.0836 0.7028 MATRIX S This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 2.5233 2 2.2354 -1.4385 3 -0.2601 -0.1603 2.0910 MATRIX T This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.0000 2 0.0000 0.0000 3 0.0000 0.0000 0.0000 MATRIX U This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 1.2662 2 0.8510 0.8728 3 0.2304 -0.1289 0.7957 MATRIX X This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 2.1241 2 1.5617 -1.1240 3 0.2648 0.3590 -2.3430 MATRIX Y This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.0000 2 0.0000 0.0000 3 0.0000 0.0000 0.0000 MATRIX Z This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 1.0004 2 0.7427 0.5947 3 0.2764 0.0633 0.8531 Vector of OBSERVED means LDLLNA APOBLNA APOELNA LDLLNB APOBLNB APOELNB Mean 19.8897 17.0177 18.8447 18.9287 16.0387 18.4547 Vector of EXPECTED means LDLLNA APOBLNA APOELNA LDLLNB APOBLNB APOELNB Mean 19.8897 17.0179 18.8446 18.9295 16.0397 18.4547 EXPECTED COVARIANCE MATRIX LDLLNA APOBLNA APOELNA LDLLNB APOBLNB APOELNB LDLLNA 7.9704 APOBLNA 6.7181 8.5521 APOELNA -0.3646 -0.2672 5.1684 LDLLNB 3.1836 2.8203 -0.3281 7.9704 APOBLNB 2.8203 3.5331 -0.1754 6.7181 8.5521 APOELNB -0.3281 -0.1754 2.2328 -0.3646 -0.2672 5.1684 Function value of this group: 783.9809 Where the fit function is -2 * Log-likelihood of raw data GROUP NUMBER: 6 G6: DOS-MF twins, datagroup MATRIX A This is a computed FULL matrix of order 3 by 3 [=X*X'] 1 2 3 1 4.5118 3.3173 0.5625 2 3.3173 3.7024 0.0101 3 0.5625 0.0101 5.6886 MATRIX C This is a computed FULL matrix of order 3 by 3 [=Y*Y'] 1 2 3 1 0.0000 0.0000 0.0000 2 0.0000 0.0000 0.0000 3 0.0000 0.0000 0.0000 MATRIX E This is a computed FULL matrix of order 3 by 3 [=Z*Z'] 1 2 3 1 1.0008 0.7430 0.2765 2 0.7430 0.9052 0.2429 3 0.2765 0.2429 0.8081 MATRIX F This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX G This is a FULL matrix of order 1 by 6 1 2 3 4 5 6 1 19.8405 17.6204 17.5759 19.3768 17.2937 17.5047 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX K This is a computed FULL matrix of order 3 by 3 [=S*S'] 1 2 3 1 6.3672 5.6406 -0.6563 2 5.6406 7.0662 -0.3507 3 -0.6563 -0.3507 4.4656 MATRIX L This is a computed FULL matrix of order 3 by 3 [=T*T'] 1 2 3 1 0.0000 0.0000 0.0000 2 0.0000 0.0000 0.0000 3 0.0000 0.0000 0.0000 MATRIX M This is a computed FULL matrix of order 3 by 3 [=U*U'] 1 2 3 1 1.6032 1.0776 0.2917 2 1.0776 1.4860 0.0836 3 0.2917 0.0836 0.7028 MATRIX S This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 2.5233 2 2.2354 -1.4385 3 -0.2601 -0.1603 2.0910 MATRIX T This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.0000 2 0.0000 0.0000 3 0.0000 0.0000 0.0000 MATRIX U This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 1.2662 2 0.8510 0.8728 3 0.2304 -0.1289 0.7957 MATRIX X This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 2.1241 2 1.5617 -1.1240 3 0.2648 0.3590 -2.3430 MATRIX Y This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.0000 2 0.0000 0.0000 3 0.0000 0.0000 0.0000 MATRIX Z This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 1.0004 2 0.7427 0.5947 3 0.2764 0.0633 0.8531 Vector of OBSERVED means LDLLNA APOBLNA APOELNA LDLLNB APOBLNB APOELNB Mean 19.8394 17.6194 17.5769 19.3756 17.2925 17.5050 Vector of EXPECTED means LDLLNA APOBLNA APOELNA LDLLNB APOBLNB APOELNB Mean 19.8405 17.6204 17.5759 19.3768 17.2937 17.5047 EXPECTED COVARIANCE MATRIX LDLLNA APOBLNA APOELNA LDLLNB APOBLNB APOELNB LDLLNA 5.5126 APOBLNA 4.0602 4.6076 APOELNA 0.8390 0.2530 6.4967 LDLLNB 2.6799 1.9704 0.3341 7.9704 APOBLNB 2.3741 2.5540 0.0378 6.7181 8.5521 APOELNB -0.2762 -0.1130 -2.5128 -0.3646 -0.2672 5.1684 Function value of this group: 403.0762 Where the fit function is -2 * Log-likelihood of raw data GROUP NUMBER: 7 G7: DOSFM twins, datagroup MATRIX A This is a computed FULL matrix of order 3 by 3 [=X*X'] 1 2 3 1 4.5118 3.3173 0.5625 2 3.3173 3.7024 0.0101 3 0.5625 0.0101 5.6886 MATRIX C This is a computed FULL matrix of order 3 by 3 [=Y*Y'] 1 2 3 1 0.0000 0.0000 0.0000 2 0.0000 0.0000 0.0000 3 0.0000 0.0000 0.0000 MATRIX E This is a computed FULL matrix of order 3 by 3 [=Z*Z'] 1 2 3 1 1.0008 0.7430 0.2765 2 0.7430 0.9052 0.2429 3 0.2765 0.2429 0.8081 MATRIX F This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX G This is a FULL matrix of order 1 by 6 1 2 3 4 5 6 1 18.3409 15.4907 19.3360 18.6466 16.2983 17.4303 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX K This is a computed FULL matrix of order 3 by 3 [=S*S'] 1 2 3 1 6.3672 5.6406 -0.6563 2 5.6406 7.0662 -0.3507 3 -0.6563 -0.3507 4.4656 MATRIX L This is a computed FULL matrix of order 3 by 3 [=T*T'] 1 2 3 1 0.0000 0.0000 0.0000 2 0.0000 0.0000 0.0000 3 0.0000 0.0000 0.0000 MATRIX M This is a computed FULL matrix of order 3 by 3 [=U*U'] 1 2 3 1 1.6032 1.0776 0.2917 2 1.0776 1.4860 0.0836 3 0.2917 0.0836 0.7028 MATRIX S This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 2.5233 2 2.2354 -1.4385 3 -0.2601 -0.1603 2.0910 MATRIX T This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.0000 2 0.0000 0.0000 3 0.0000 0.0000 0.0000 MATRIX U This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 1.2662 2 0.8510 0.8728 3 0.2304 -0.1289 0.7957 MATRIX X This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 2.1241 2 1.5617 -1.1240 3 0.2648 0.3590 -2.3430 MATRIX Y This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.0000 2 0.0000 0.0000 3 0.0000 0.0000 0.0000 MATRIX Z This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 1.0004 2 0.7427 0.5947 3 0.2764 0.0633 0.8531 Vector of OBSERVED means LDLLNA APOBLNA APOELNA LDLLNB APOBLNB APOELNB Mean 18.3408 15.4908 19.3377 18.6462 16.2985 17.4277 Vector of EXPECTED means LDLLNA APOBLNA APOELNA LDLLNB APOBLNB APOELNB Mean 18.3409 15.4907 19.3360 18.6466 16.2983 17.4303 EXPECTED COVARIANCE MATRIX LDLLNA APOBLNA APOELNA LDLLNB APOBLNB APOELNB LDLLNA 7.9704 APOBLNA 6.7181 8.5521 APOELNA -0.3646 -0.2672 5.1684 LDLLNB 2.6799 2.3741 -0.2762 5.5126 APOBLNB 1.9704 2.5540 -0.1130 4.0602 4.6076 APOELNB 0.3341 0.0378 -2.5128 0.8390 0.2530 6.4967 Function value of this group: 345.1131 Where the fit function is -2 * Log-likelihood of raw data GROUP NUMBER: 8 G8: Calculate genetic and environmental correlations and standardized variance-covariance components MATRIX A This is a computed FULL matrix of order 3 by 3 [=X*X'] 1 2 3 1 4.5118 3.3173 0.5625 2 3.3173 3.7024 0.0101 3 0.5625 0.0101 5.6886 MATRIX C This is a computed FULL matrix of order 3 by 3 [=Y*Y'] 1 2 3 1 0.0000 0.0000 0.0000 2 0.0000 0.0000 0.0000 3 0.0000 0.0000 0.0000 MATRIX E This is a computed FULL matrix of order 3 by 3 1 2 3 1 1.0008 0.7430 0.2765 2 0.7430 0.9052 0.2429 3 0.2765 0.2429 0.8081 MATRIX I This is an IDENTITY matrix of order 3 by 3 MATRIX J This is a computed FULL matrix of order 3 by 3 [=\STND(K)] 1 2 3 1 1.0000 0.8409 -0.1231 2 0.8409 1.0000 -0.0624 3 -0.1231 -0.0624 1.0000 MATRIX K This is a computed FULL matrix of order 3 by 3 [=S*S'] 1 2 3 1 6.3672 5.6406 -0.6563 2 5.6406 7.0662 -0.3507 3 -0.6563 -0.3507 4.4656 MATRIX L This is a computed FULL matrix of order 3 by 3 1 2 3 1 0.0000 0.0000 0.0000 2 0.0000 0.0000 0.0000 3 0.0000 0.0000 0.0000 MATRIX M This is a computed FULL matrix of order 3 by 3 1 2 3 1 1.6032 1.0776 0.2917 2 1.0776 1.4860 0.0836 3 0.2917 0.0836 0.7028 MATRIX N This is a computed FULL matrix of order 3 by 3 [=\STND(L)] 1 2 3 1 0.0000 0.0000 0.0000 2 0.0000 0.0000 0.0000 3 0.0000 0.0000 0.0000 MATRIX O This is a computed FULL matrix of order 3 by 3 [=\STND(M)] 1 2 3 1 1.0000 0.6981 0.2748 2 0.6981 1.0000 0.0818 3 0.2748 0.0818 1.0000 MATRIX R This is a computed FULL matrix of order 3 by 3 [=\SQRT(I.A)~*A*\SQRT(I.A)~] 1 2 3 1 1.0000 0.8116 0.1110 2 0.8116 1.0000 0.0022 3 0.1110 0.0022 1.0000 MATRIX S This is a computed FULL matrix of order 3 by 3 [=\SQRT(I.C)~*C*\SQRT(I.C)~] 1 2 3 1 0.0000 0.0000 0.0000 2 0.0000 0.0000 0.0000 3 0.0000 0.0000 0.0000 MATRIX T This is a computed FULL matrix of order 3 by 3 [=\SQRT(I.E)~*E*\SQRT(I.E)~] 1 2 3 1 1.0000 0.7806 0.3075 2 0.7806 1.0000 0.2840 3 0.3075 0.2840 1.0000 MATRIX V This is a computed FULL matrix of order 6 by 9 [=A%(A+C+E)|C%(A+C+E)|E%(A+C+E)_K%(K+L+M)|L%(K+L+M)|M%(K+L+M)] H2_LDL H2_APOB H2_APOE C2_LDL C2_APOB C2_APOE E2_LDL MEN 0.8185 0.8170 0.6704 0.0000 0.0000 0.0000 0.1815 MEN 0.8170 0.8035 0.0398 0.0000 0.0000 0.0000 0.1830 MEN 0.6704 0.0398 0.8756 0.0000 0.0000 0.0000 0.3296 WOMEN 0.7989 0.8396 1.7999 0.0000 0.0000 0.0000 0.2011 WOMEN 0.8396 0.8262 1.3129 0.0000 0.0000 0.0000 0.1604 WOMEN 1.7999 1.3129 0.8640 0.0000 0.0000 0.0000 -0.7999 E2_APOB E2_APOE MEN 0.1830 0.3296 MEN 0.1965 0.9602 MEN 0.9602 0.1244 WOMEN 0.1604 -0.7999 WOMEN 0.1738 -0.3129 WOMEN -0.3129 0.1360 MATRIX W This is a computed FULL matrix of order 6 by 9 [=R|S|T_J|N|O] R_A R_A R_A R_C R_C R_C LDLMEN 1.0000 0.8116 0.1110 0.0000 0.0000 0.0000 APOBMEN 0.8116 1.0000 0.0022 0.0000 0.0000 0.0000 APOEMEN 0.1110 0.0022 1.0000 0.0000 0.0000 0.0000 LDLWOMEN 1.0000 0.8409 -0.1231 0.0000 0.0000 0.0000 APOBWOMEN 0.8409 1.0000 -0.0624 0.0000 0.0000 0.0000 APOEWOMEN -0.1231 -0.0624 1.0000 0.0000 0.0000 0.0000 R_E R_E R_E LDLMEN 1.0000 0.7806 0.3075 APOBMEN 0.7806 1.0000 0.2840 APOEMEN 0.3075 0.2840 1.0000 LDLWOMEN 1.0000 0.6981 0.2748 APOBWOMEN 0.6981 1.0000 0.0818 APOEWOMEN 0.2748 0.0818 1.0000 Your model has 60 estimated parameters and 956 Observed statistics -2 times log-likelihood of data >>> 3827.217 Degrees of freedom >>>>>>>>>>>>>>>> 896 This problem used 4.7% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 1: 4.87 Execution 0: 0: 0:***** TOTAL 0: 0: 0:14.98 Total number of warnings issued: 0 ______________________________________________________________________________ ! FIX DOS COVARIANCE FOR A AT 0.5, WHICH EQUALS THE DZ COVARIANCE ! THIS PROVIDES A TEST OF WHETHER THE SAME GENES ARE EXPRESSED IN MALES AND FEMALES !DROP @.5 F 1 1 1 !end