! SCRIPT NAME : Cholesky_M+F.mx (mb) ! GOAL : ACE model with sex limitation for a trait rated by mother and father of the subject ! DATA : continuous ! INPUT : raw data ! UNI/BI/MULTI : univariate ! DATA-GROUPS : MZM, DZM, MZF, DZF, DOSm-f ! MEANS MODEL : means ! VARIANCE COVARIANCE MODEL(S) : 1. ACE ! Downloading Mx software: http://www.vcu.edu/mx ! Mx script's library: http://www.psy.vu.nl/mxbib ! ! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! The following MX script lines were read for group 1 #NGROUPS 9 Note: #NGroup set number of groups to 9 #DEFINE NVAR 2 ! NUMBER OF VARIABLES (PHENOTYPES) UNDER STUDY #DEFINE NRATER 2 ! NUMBER OF RATERS (E.G. MOTHER AND FATHER) #DEFINE NRATNSIB 4 ! NUMBER OF RATER * NUMBER OF SIBS (IN THIS CASE TWINPAIRS) * NVAR G1: BOYS DATA CALC BEGIN MATRICES; X LOWER NVAR NVAR FREE ! GENETIC STRUCTURE FOR BEHAVIOR SIMILAR ASSESSED BY BOTH PARENTS Y LOWER NVAR NVAR FREE ! SHARED ENVIRONMENTAL STRUCTURE FOR BEHAVIOR SIMILAR ASSESSED BY BOTH PARENTS Z LOWER NVAR NVAR FREE ! NONSHARED ENVIRONMENTAL STRUCTURE FOR BEHAVIOR SIMILAR ASSESSED BY BOTH PARENTS R FULL NRATNSIB 1 FREE ! MEANS H FULL 1 1 FIXED END MATRICES; MATRIX H .5 ! MATRIX H IS FIXED AND CONTAINS .5 ST .5 ALL START 1 R 1 1 - R 4 1 BEGIN ALGEBRA; A= X*X' ; ! ADDITIVE GENETIC VARIANCE OF BEHAVIOR SIMILAR ASSESSED BY BOTH PARENTS C= Y*Y' ; ! SHARED ENVIRONMENTAL VARIANCE OF BEHAVIOR SIMILAR ASSESED BY BOTH PARENTS E= Z*Z' ; ! NONSHARED ENVIRONMENTAL VARIANCE OF BEHAVIOR SIMILAR ASSESSED BY BOTH PARENTS END ALGEBRA; END The following MX script lines were read for group 2 G2: GIRLS DATA CALC BEGIN MATRICES; X LOWER NVAR NVAR FREE ! GENETIC STRUCTURE FOR BEHAVIOR SIMILAR ASSESSED BY BOTH PARENTS Y LOWER NVAR NVAR FREE ! SHARED ENVIRONMENTAL STRUCTURE FOR BEHAVIOR SIMILAR ASSESSED BY BOTH PARENTS Z LOWER NVAR NVAR FREE ! NONSHARED ENVIRONMENTAL STRUCTURE FOR BEHAVIOR SIMILAR ASSESSED BY BOTH PARENTS R FULL NRATNSIB 1 FREE ! MEANS H FULL 1 1 FIXED END MATRICES; MATRIX H .5 ! MATRIX H IS FIXED AND CONTAINS .5 ST .5 ALL START 1 R 1 1 - R 4 1 BEGIN ALGEBRA; A= X*X' ; ! ADDITIVE GENETIC VARIANCE OF BEHAVIOR SIMILAR ASSESSED BY BOTH PARENTS C= Y*Y' ; ! SHARED ENVIRONMENTAL VARIANCE OF BEHAVIOR SIMILAR ASSESED BY BOTH PARENTS E= Z*Z' ; ! NONSHARED ENVIRONMENTAL VARIANCE OF BEHAVIOR SIMILAR ASSESSED BY BOTH PARENTS END ALGEBRA; END The following MX script lines were read for group 3 G3: MALE MZ TWIN PAIRS DATA NINPUT_VARS=29 ! NUMBER OF VARIABLES IN DATA FILE MISSING=-1.00 ! MISSING VALUE = -1.00 NOTE: Missing value *string* set to '-1.00' RECTANGULAR FILE=TAD.DAT ! READ IN RAW DATA FILE Rectangular continuous data read initiated NOTE: Rectangular file contained 6438 records with data that contained a total of 99821 observations LABELS NTRID SEXA SEXB ZYG COHORT TDELM7A TAGGM7A TDELV7A TAGGV7A TDELM7B TAGGM7B TDELV7B TAGGV7B TDELM10A TAGGM10A TDELV10A TAGGV10A TDELM10B TAGGM10B TDELV10B TAGGV10B TDELM12A TAGGM12A TDELV12A TAGGV12A TDELM12B TAGGM12B TDELV12B TAGGV12B ! TELL MX WHAT IS FOUND IN THE DATAFILE SELECT IF ZYG = 1; ! SELECT MZM TWINS NOTE: Select if yields 1066 data vectors for analysis NOTE: Vectors contain a total of 17248 observations SELECT TAGGM10A TAGGV10A TAGGM10B TAGGV10B ; ! SELECT A PHENOTYPE (RATED BY MOTHER AND RATED BY FATHER) MATRICES= GROUP 1 NOTE: Selection yields 525 data vectors for analysis NOTE: Vectors contain a total of 1841 observations MEANS R ; ! MODEL FOR THE MEANS COVARIANCES A+C+E | A+C_ A+C | A+C+E ; ! MODEL FOR MZ VARIANCE/COVARIANCE OPTION RSIDUAL END The following MX script lines were read for group 4 G4: MALE DZ TWIN PAIRS DATA NINPUT_VARS=29 ! NUMBER OF VARIABLES IN DATA FILE MISSING=-1.00 ! MISSING VALUE = -1.00 NOTE: Missing value *string* set to '-1.00' RECTANGULAR FILE=TAD.DAT ! READ IN RAW DATA FILE Rectangular continuous data read initiated NOTE: Rectangular file contained 6438 records with data that contained a total of 99821 observations LABELS NTRID SEXA SEXB ZYG COHORT TDELM7A TAGGM7A TDELV7A TAGGV7A TDELM7B TAGGM7B TDELV7B TAGGV7B TDELM10A TAGGM10A TDELV10A TAGGV10A TDELM10B TAGGM10B TDELV10B TAGGV10B TDELM12A TAGGM12A TDELV12A TAGGV12A TDELM12B TAGGM12B TDELV12B TAGGV12B ! TELL MX WHAT IS FOUND IN THE DATAFILE SELECT IF ZYG = 2; ! SELECT DZM TWINS NOTE: Select if yields 1040 data vectors for analysis NOTE: Vectors contain a total of 16025 observations SELECT TAGGM10A TAGGV10A TAGGM10B TAGGV10B ; ! SELECT A PHENOTYPE (RATED BY MOTHER AND RATED BY FATHER) MATRICES= GROUP 1 NOTE: Selection yields 472 data vectors for analysis NOTE: Vectors contain a total of 1631 observations MEANS R ; ! MODEL FOR THE MEANS COVARIANCES A+C+E | H@A+C_ H@A+C | A+C+E ; ! MODEL FOR DZ VARIANCE/COVARIANCE OPTION RSIDUAL END The following MX script lines were read for group 5 G5: FEMALE MZ TWIN PAIRS DATA NINPUT_VARS=29 ! NUMBER OF VARIABLES IN DATA FILE MISSING=-1.00 ! MISSING VALUE = -1.00 NOTE: Missing value *string* set to '-1.00' RECTANGULAR FILE=TAD.DAT ! READ IN RAW DATA FILE Rectangular continuous data read initiated NOTE: Rectangular file contained 6438 records with data that contained a total of 99821 observations LABELS NTRID SEXA SEXB ZYG COHORT TDELM7A TAGGM7A TDELV7A TAGGV7A TDELM7B TAGGM7B TDELV7B TAGGV7B TDELM10A TAGGM10A TDELV10A TAGGV10A TDELM10B TAGGM10B TDELV10B TAGGV10B TDELM12A TAGGM12A TDELV12A TAGGV12A TDELM12B TAGGM12B TDELV12B TAGGV12B ! TELL MX WHAT IS FOUND IN THE DATAFILE SELECT IF ZYG = 3; ! SELECT MZF TWINS NOTE: Select if yields 1222 data vectors for analysis NOTE: Vectors contain a total of 19750 observations SELECT TAGGM10A TAGGV10A TAGGM10B TAGGV10B ; ! SELECT A PHENOTYPE (RATED BY MOTHER AND RATED BY FATHER) MATRICES= GROUP 2 NOTE: Selection yields 619 data vectors for analysis NOTE: Vectors contain a total of 2170 observations MEANS R ; ! MODEL FOR THE MEANS COVARIANCES A+C+E | A+C_ A+C | A+C+E ; ! MODEL FOR MZ VARIANCE/COVARIANCE OPTION RSIDUAL END The following MX script lines were read for group 6 G6: FEMALE DZ TWIN PAIRS DATA NINPUT_VARS=29 ! NUMBER OF VARIABLES IN DATA FILE MISSING=-1.00 ! MISSING VALUE = -1.00 NOTE: Missing value *string* set to '-1.00' RECTANGULAR FILE=TAD.DAT ! READ IN RAW DATA FILE Rectangular continuous data read initiated NOTE: Rectangular file contained 6438 records with data that contained a total of 99821 observations LABELS NTRID SEXA SEXB ZYG COHORT TDELM7A TAGGM7A TDELV7A TAGGV7A TDELM7B TAGGM7B TDELV7B TAGGV7B TDELM10A TAGGM10A TDELV10A TAGGV10A TDELM10B TAGGM10B TDELV10B TAGGV10B TDELM12A TAGGM12A TDELV12A TAGGV12A TDELM12B TAGGM12B TDELV12B TAGGV12B ! TELL MX WHAT IS FOUND IN THE DATAFILE SELECT IF ZYG = 4; ! SELECT DZF TWINS NOTE: Select if yields 997 data vectors for analysis NOTE: Vectors contain a total of 15592 observations SELECT TAGGM10A TAGGV10A TAGGM10B TAGGV10B ; ! SELECT A PHENOTYPE (RATED BY MOTHER AND RATED BY FATHER) MATRICES= GROUP 2 NOTE: Selection yields 464 data vectors for analysis NOTE: Vectors contain a total of 1617 observations MEANS R ; ! MODEL FOR THE MEANS COVARIANCES A+C+E | H@A+C_ H@A+C | A+C+E ; ! MODEL FOR DZ VARIANCE/COVARIANCE OPTION RSIDUAL END The following MX script lines were read for group 7 G7: MALE-FEMALE DZ TWIN PAIRS DATA NINPUT_VARS=29 ! NUMBER OF VARIABLES IN DATA FILE MISSING=-1.00 ! MISSING VALUE = -1.00 NOTE: Missing value *string* set to '-1.00' RECTANGULAR FILE=TAD.DAT ! READ IN RAW DATA FILE Rectangular continuous data read initiated NOTE: Rectangular file contained 6438 records with data that contained a total of 99821 observations LABELS NTRID SEXA SEXB ZYG COHORT TDELM7A TAGGM7A TDELV7A TAGGV7A TDELM7B TAGGM7B TDELV7B TAGGV7B TDELM10A TAGGM10A TDELV10A TAGGV10A TDELM10B TAGGM10B TDELV10B TAGGV10B TDELM12A TAGGM12A TDELV12A TAGGV12A TDELM12B TAGGM12B TDELV12B TAGGV12B ! TELL MX WHAT IS FOUND IN THE DATAFILE SELECT IF ZYG = 5; ! SELECT DOSMF TWINS NOTE: Select if yields 2113 data vectors for analysis NOTE: Vectors contain a total of 31206 observations SELECT TAGGM10A TAGGV10A TAGGM10B TAGGV10B ; ! SELECT A PHENOTYPE (RATED BY MOTHER AND RATED BY FATHER) MATRICES NOTE: Selection yields 899 data vectors for analysis NOTE: Vectors contain a total of 3127 observations X LOWER NVAR NVAR =X1 ! GENETIC STRUCTURE BOYS Y LOWER NVAR NVAR =Y1 ! COMMON ENVIRONMENTAL STRUCTURE BOYS Z LOWER NVAR NVAR =Z1 ! SPECIFIC ENVIRONMENTAL STRUCTURE BOYS A LOWER NVAR NVAR =X2 ! GENETIC STRUCTURE GIRLS C LOWER NVAR NVAR =Y2 ! COMMON ENVIRONMENTAL STRUCTURE GIRLS E LOWER NVAR NVAR =Z2 ! SPECIFIC ENVIRONMENTAL STRUCTURE GIRLS H FULL 1 1 = H5 R FULL NRATNSIB 1 FREE ! MEANS END MATRICES MEANS R ; COVARIANCES (X*X')+(Y*Y')+(Z*Z') |H@(X*A')+(Y*C')_ H@(A*X')+(C*Y') | (A*A') + (C*C') + (E*E') ; ! MODEL FOR DOSM-F VARIANCE/COVARIANCE OPTION RSIDUAL END The following MX script lines were read for group 8 GR 8: CALCULATE STANDARDIZED ESTIMATES FOR BOYS DATA CALCULATION MATRICES A COMP NVAR NVAR =A1 C COMP NVAR NVAR =C1 E COMP NVAR NVAR =E1 END MATRICES BEGIN ALGEBRA; X= A%(A+C+E) ; Y= C%(A+C+E) ; Z= E%(A+C+E) ; END ALGEBRA; INTERVAL @95 X 1 1 - X 2 2 INTERVAL @95 Y 1 1 - Y 2 2 INTERVAL @95 Z 1 1 - Z 2 2 OPTION ND=4 END The following MX script lines were read for group 9 GR 9: CALCULATE STANDARDIZED ESTIMATES FOR GIRLS DATA CALCULATION MATRICES A COMP NVAR NVAR =A2 C COMP NVAR NVAR =C2 E COMP NVAR NVAR =E2 END MATRICES BEGIN ALGEBRA; X= A%(A+C+E) ; Y= C%(A+C+E) ; Z= E%(A+C+E) ; END ALGEBRA; INTERVAL @95 X 1 1 - X 2 2 INTERVAL @95 Y 1 1 - Y 2 2 INTERVAL @95 Z 1 1 - Z 2 2 OPTION ND=4 END Summary of VL file data for group 3 TAGGM10A TAGGV10A TAGGM10B TAGGV10B Code 1.0000 2.0000 3.0000 4.0000 Number 515.0000 406.0000 514.0000 406.0000 Mean 2.4110 2.1841 2.3842 2.1265 Variance 1.6183 1.5048 1.5618 1.5028 Minimum 0.0000 0.0000 0.0000 0.0000 Maximum 5.2900 5.3900 5.6600 5.2000 Summary of VL file data for group 4 TAGGM10A TAGGV10A TAGGM10B TAGGV10B Code 1.0000 2.0000 3.0000 4.0000 Number 468.0000 347.0000 469.0000 347.0000 Mean 2.2404 2.0840 2.3056 2.1064 Variance 1.6156 1.3550 1.4906 1.3829 Minimum 0.0000 0.0000 0.0000 0.0000 Maximum 6.0800 5.2900 5.7400 5.4800 Summary of VL file data for group 5 TAGGM10A TAGGV10A TAGGM10B TAGGV10B Code 1.0000 2.0000 3.0000 4.0000 Number 613.0000 474.0000 610.0000 473.0000 Mean 1.9859 1.7216 1.8960 1.6630 Variance 1.3027 1.4009 1.2773 1.2452 Minimum 0.0000 0.0000 0.0000 0.0000 Maximum 4.6900 4.6900 5.5700 5.2900 Summary of VL file data for group 6 TAGGM10A TAGGV10A TAGGM10B TAGGV10B Code 1.0000 2.0000 3.0000 4.0000 Number 455.0000 355.0000 452.0000 355.0000 Mean 1.9569 1.7613 1.9102 1.7538 Variance 1.3630 1.3092 1.2583 1.1860 Minimum 0.0000 0.0000 0.0000 0.0000 Maximum 5.2000 5.2000 5.0000 4.6900 Summary of VL file data for group 7 TAGGM10A TAGGV10A TAGGM10B TAGGV10B Code 1.0000 2.0000 3.0000 4.0000 Number 887.0000 676.0000 887.0000 677.0000 Mean 2.2269 2.0108 1.8164 1.6126 Variance 1.5640 1.5070 1.2925 1.2329 Minimum 0.0000 0.0000 0.0000 0.0000 Maximum 5.4800 5.5700 5.1000 5.2000 PARAMETER SPECIFICATIONS GROUP NUMBER: 1 G1: boys MATRIX A This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 2 by 2 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 R This is a FULL matrix of order 4 by 1 1 1 10 2 11 3 12 4 13 MATRIX X This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 1 2 2 3 MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 4 2 5 6 MATRIX Z This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 7 2 8 9 GROUP NUMBER: 2 G2: Girls MATRIX A This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 2 by 2 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 R This is a FULL matrix of order 4 by 1 1 1 23 2 24 3 25 4 26 MATRIX X This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 14 2 15 16 MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 17 2 18 19 MATRIX Z This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 20 2 21 22 GROUP NUMBER: 3 G3: male MZ twin pairs MATRIX A This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 2 by 2 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 R This is a FULL matrix of order 4 by 1 1 1 10 2 11 3 12 4 13 MATRIX X This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 1 2 2 3 MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 4 2 5 6 MATRIX Z This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 7 2 8 9 GROUP NUMBER: 4 G4: male DZ twin pairs MATRIX A This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 2 by 2 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 R This is a FULL matrix of order 4 by 1 1 1 10 2 11 3 12 4 13 MATRIX X This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 1 2 2 3 MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 4 2 5 6 MATRIX Z This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 7 2 8 9 GROUP NUMBER: 5 G5: female MZ twin pairs MATRIX A This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 2 by 2 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 R This is a FULL matrix of order 4 by 1 1 1 23 2 24 3 25 4 26 MATRIX X This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 14 2 15 16 MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 17 2 18 19 MATRIX Z This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 20 2 21 22 GROUP NUMBER: 6 G6: female DZ twin pairs MATRIX A This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 2 by 2 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 R This is a FULL matrix of order 4 by 1 1 1 23 2 24 3 25 4 26 MATRIX X This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 14 2 15 16 MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 17 2 18 19 MATRIX Z This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 20 2 21 22 GROUP NUMBER: 7 G7: Male-Female DZ twin pairs MATRIX A This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 14 2 15 16 MATRIX C This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 17 2 18 19 MATRIX E This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 20 2 21 22 MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX R This is a FULL matrix of order 4 by 1 1 1 27 2 28 3 29 4 30 MATRIX X This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 1 2 2 3 MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 4 2 5 6 MATRIX Z This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 7 2 8 9 GROUP NUMBER: 8 Gr 8: calculate standardized estimates for BOYS MATRIX A This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX X This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX Y This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX Z This is a computed FULL matrix of order 2 by 2 It has no free parameters specified GROUP NUMBER: 9 Gr 9: calculate standardized estimates for GIRLS MATRIX A This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX X This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX Y This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX Z This is a computed FULL matrix of order 2 by 2 It has no free parameters specified Mx starting optimization; number of parameters = 30 MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1: boys MATRIX A This is a computed FULL matrix of order 2 by 2 [=X*X'] 1 2 1 1.1395 0.8811 2 0.8811 0.8849 MATRIX C This is a computed FULL matrix of order 2 by 2 [=Y*Y'] 1 2 1 0.1754 0.1241 2 0.1241 0.3734 MATRIX E This is a computed FULL matrix of order 2 by 2 [=Z*Z'] 1 2 1 0.2561 0.1495 2 0.1495 0.2389 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX R This is a FULL matrix of order 4 by 1 1 1 2.3398 2 2.1694 3 2.3335 4 2.1489 MATRIX X This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 1.0675 2 0.8254 0.4513 MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.4188 2 0.2962 0.5345 MATRIX Z This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.5060 2 0.2955 0.3894 GROUP NUMBER: 2 G2: Girls MATRIX A This is a computed FULL matrix of order 2 by 2 [=X*X'] 1 2 1 0.6725 0.6322 2 0.6322 0.7327 MATRIX C This is a computed FULL matrix of order 2 by 2 [=Y*Y'] 1 2 1 0.3630 0.1387 2 0.1387 0.3024 MATRIX E This is a computed FULL matrix of order 2 by 2 [=Z*Z'] 1 2 1 0.2686 0.1470 2 0.1470 0.2427 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX R This is a FULL matrix of order 4 by 1 1 1 1.9797 2 1.7689 3 1.9037 4 1.7239 MATRIX X This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.8200 2 0.7709 0.3721 MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.6025 2 0.2302 0.4994 MATRIX Z This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.5183 2 0.2837 0.4028 GROUP NUMBER: 3 G3: male MZ twin pairs MATRIX A This is a computed FULL matrix of order 2 by 2 [=X*X'] 1 2 1 1.1395 0.8811 2 0.8811 0.8849 MATRIX C This is a computed FULL matrix of order 2 by 2 [=Y*Y'] 1 2 1 0.1754 0.1241 2 0.1241 0.3734 MATRIX E This is a computed FULL matrix of order 2 by 2 [=Z*Z'] 1 2 1 0.2561 0.1495 2 0.1495 0.2389 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX R This is a FULL matrix of order 4 by 1 1 1 2.3398 2 2.1694 3 2.3335 4 2.1489 MATRIX X This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 1.0675 2 0.8254 0.4513 MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.4188 2 0.2962 0.5345 MATRIX Z This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.5060 2 0.2955 0.3894 Vector of OBSERVED means TAGGM10A TAGGV10A TAGGM10B TAGGV10B Mean 2.4110 2.1841 2.3842 2.1265 Vector of EXPECTED means TAGGM10A TAGGV10A TAGGM10B TAGGV10B Mean 2.3398 2.1694 2.3335 2.1489 EXPECTED COVARIANCE MATRIX TAGGM10A TAGGV10A TAGGM10B TAGGV10B TAGGM10A 1.5710 TAGGV10A 1.1547 1.4973 TAGGM10B 1.3149 1.0052 1.5710 TAGGV10B 1.0052 1.2584 1.1547 1.4973 Function value of this group: 4406.7647 Where the fit function is -2 * Log-likelihood of raw data GROUP NUMBER: 4 G4: male DZ twin pairs MATRIX A This is a computed FULL matrix of order 2 by 2 [=X*X'] 1 2 1 1.1395 0.8811 2 0.8811 0.8849 MATRIX C This is a computed FULL matrix of order 2 by 2 [=Y*Y'] 1 2 1 0.1754 0.1241 2 0.1241 0.3734 MATRIX E This is a computed FULL matrix of order 2 by 2 [=Z*Z'] 1 2 1 0.2561 0.1495 2 0.1495 0.2389 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX R This is a FULL matrix of order 4 by 1 1 1 2.3398 2 2.1694 3 2.3335 4 2.1489 MATRIX X This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 1.0675 2 0.8254 0.4513 MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.4188 2 0.2962 0.5345 MATRIX Z This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.5060 2 0.2955 0.3894 Vector of OBSERVED means TAGGM10A TAGGV10A TAGGM10B TAGGV10B Mean 2.2404 2.0840 2.3056 2.1064 Vector of EXPECTED means TAGGM10A TAGGV10A TAGGM10B TAGGV10B Mean 2.3398 2.1694 2.3335 2.1489 EXPECTED COVARIANCE MATRIX TAGGM10A TAGGV10A TAGGM10B TAGGV10B TAGGM10A 1.5710 TAGGV10A 1.1547 1.4973 TAGGM10B 0.7452 0.5646 1.5710 TAGGV10B 0.5646 0.8159 1.1547 1.4973 Function value of this group: 4489.1227 Where the fit function is -2 * Log-likelihood of raw data GROUP NUMBER: 5 G5: female MZ twin pairs MATRIX A This is a computed FULL matrix of order 2 by 2 [=X*X'] 1 2 1 0.6725 0.6322 2 0.6322 0.7327 MATRIX C This is a computed FULL matrix of order 2 by 2 [=Y*Y'] 1 2 1 0.3630 0.1387 2 0.1387 0.3024 MATRIX E This is a computed FULL matrix of order 2 by 2 [=Z*Z'] 1 2 1 0.2686 0.1470 2 0.1470 0.2427 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX R This is a FULL matrix of order 4 by 1 1 1 1.9797 2 1.7689 3 1.9037 4 1.7239 MATRIX X This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.8200 2 0.7709 0.3721 MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.6025 2 0.2302 0.4994 MATRIX Z This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.5183 2 0.2837 0.4028 Vector of OBSERVED means TAGGM10A TAGGV10A TAGGM10B TAGGV10B Mean 1.9859 1.7216 1.8960 1.6630 Vector of EXPECTED means TAGGM10A TAGGV10A TAGGM10B TAGGV10B Mean 1.9797 1.7689 1.9037 1.7239 EXPECTED COVARIANCE MATRIX TAGGM10A TAGGV10A TAGGM10B TAGGV10B TAGGM10A 1.3041 TAGGV10A 0.9179 1.2779 TAGGM10B 1.0355 0.7709 1.3041 TAGGV10B 0.7709 1.0352 0.9179 1.2779 Function value of this group: 5131.1353 Where the fit function is -2 * Log-likelihood of raw data GROUP NUMBER: 6 G6: female DZ twin pairs MATRIX A This is a computed FULL matrix of order 2 by 2 [=X*X'] 1 2 1 0.6725 0.6322 2 0.6322 0.7327 MATRIX C This is a computed FULL matrix of order 2 by 2 [=Y*Y'] 1 2 1 0.3630 0.1387 2 0.1387 0.3024 MATRIX E This is a computed FULL matrix of order 2 by 2 [=Z*Z'] 1 2 1 0.2686 0.1470 2 0.1470 0.2427 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX R This is a FULL matrix of order 4 by 1 1 1 1.9797 2 1.7689 3 1.9037 4 1.7239 MATRIX X This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.8200 2 0.7709 0.3721 MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.6025 2 0.2302 0.4994 MATRIX Z This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.5183 2 0.2837 0.4028 Vector of OBSERVED means TAGGM10A TAGGV10A TAGGM10B TAGGV10B Mean 1.9569 1.7613 1.9102 1.7538 Vector of EXPECTED means TAGGM10A TAGGV10A TAGGM10B TAGGV10B Mean 1.9797 1.7689 1.9037 1.7239 EXPECTED COVARIANCE MATRIX TAGGM10A TAGGV10A TAGGM10B TAGGV10B TAGGM10A 1.3041 TAGGV10A 0.9179 1.2779 TAGGM10B 0.6993 0.4548 1.3041 TAGGV10B 0.4548 0.6688 0.9179 1.2779 Function value of this group: 4193.4375 Where the fit function is -2 * Log-likelihood of raw data GROUP NUMBER: 7 G7: Male-Female DZ twin pairs MATRIX A This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.8200 2 0.7709 0.3721 MATRIX C This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.6025 2 0.2302 0.4994 MATRIX E This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.5183 2 0.2837 0.4028 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX R This is a FULL matrix of order 4 by 1 1 1 2.2289 2 2.0446 3 1.8140 4 1.6314 MATRIX X This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 1.0675 2 0.8254 0.4513 MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.4188 2 0.2962 0.5345 MATRIX Z This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.5060 2 0.2955 0.3894 Vector of OBSERVED means TAGGM10A TAGGV10A TAGGM10B TAGGV10B Mean 2.2269 2.0108 1.8164 1.6126 Vector of EXPECTED means TAGGM10A TAGGV10A TAGGM10B TAGGV10B Mean 2.2289 2.0446 1.8140 1.6314 EXPECTED COVARIANCE MATRIX TAGGM10A TAGGV10A TAGGM10B TAGGV10B TAGGM10A 1.5710 TAGGV10A 1.1547 1.4973 TAGGM10B 0.6900 0.5169 1.3041 TAGGV10B 0.5079 0.7373 0.9179 1.2779 Function value of this group: 8344.0194 Where the fit function is -2 * Log-likelihood of raw data GROUP NUMBER: 8 Gr 8: calculate standardized estimates for BOYS MATRIX A This is a computed FULL matrix of order 2 by 2 [=X*X'] 1 2 1 1.1395 0.8811 2 0.8811 0.8849 MATRIX C This is a computed FULL matrix of order 2 by 2 [=Y*Y'] 1 2 1 0.1754 0.1241 2 0.1241 0.3734 MATRIX E This is a computed FULL matrix of order 2 by 2 1 2 1 0.2561 0.1495 2 0.1495 0.2389 MATRIX X This is a computed FULL matrix of order 2 by 2 [=A%(A+C+E)] 1 2 1 0.7254 0.7631 2 0.7631 0.5910 MATRIX Y This is a computed FULL matrix of order 2 by 2 [=C%(A+C+E)] 1 2 1 0.1116 0.1074 2 0.1074 0.2494 MATRIX Z This is a computed FULL matrix of order 2 by 2 [=E%(A+C+E)] 1 2 1 0.1630 0.1295 2 0.1295 0.1596 GROUP NUMBER: 9 Gr 9: calculate standardized estimates for GIRLS MATRIX A This is a computed FULL matrix of order 2 by 2 [=X*X'] 1 2 1 0.6725 0.6322 2 0.6322 0.7327 MATRIX C This is a computed FULL matrix of order 2 by 2 [=Y*Y'] 1 2 1 0.3630 0.1387 2 0.1387 0.3024 MATRIX E This is a computed FULL matrix of order 2 by 2 1 2 1 0.2686 0.1470 2 0.1470 0.2427 MATRIX X This is a computed FULL matrix of order 2 by 2 [=A%(A+C+E)] 1 2 1 0.5157 0.6887 2 0.6887 0.5734 MATRIX Y This is a computed FULL matrix of order 2 by 2 [=C%(A+C+E)] 1 2 1 0.2784 0.1511 2 0.1511 0.2367 MATRIX Z This is a computed FULL matrix of order 2 by 2 [=E%(A+C+E)] 1 2 1 0.2060 0.1602 2 0.1602 0.1900 Your model has 30 estimated parameters and 10386 Observed statistics -2 times log-likelihood of data >>> 26564.480 Degrees of freedom >>>>>>>>>>>>>>>> 10356 12 Confidence intervals requested in group 8 Matrix Element Int. Estimate Lower Upper Lfail Ufail X 8 1 1 95.0 0.7254 0.6468 0.7904 0 1 0 1 X 8 1 2 95.0 0.7631 0.6722 0.8352 0 1 0 1 X 8 2 1 95.0 0.7631 0.6722 0.8352 0 1 0 1 X 8 2 2 95.0 0.5910 0.5049 0.6716 0 1 0 1 Y 8 1 1 95.0 0.1116 0.0511 0.1863 0 1 0 1 Y 8 1 2 95.0 0.1074 0.0412 0.1931 0 1 0 1 Y 8 2 1 95.0 0.1074 0.0412 0.1931 0 1 0 1 Y 8 2 2 95.0 0.2494 0.1745 0.3288 0 1 0 1 Z 8 1 1 95.0 0.1630 0.1432 0.1861 0 1 6 1 Z 8 1 2 95.0 0.1295 0.1068 0.1559 0 1 0 1 Z 8 2 1 95.0 0.1295 0.1068 0.1559 0 1 0 1 Z 8 2 2 95.0 0.1596 0.1379 0.1853 0 1 0 1 12 Confidence intervals requested in group 9 Matrix Element Int. Estimate Lower Upper Lfail Ufail X 9 1 1 95.0 0.5157 0.4048 0.6313 0 1 0 1 X 9 1 2 95.0 0.6887 0.5527 0.8134 0 1 0 1 X 9 2 1 95.0 0.6887 0.5527 0.8134 0 1 0 1 X 9 2 2 95.0 0.5734 0.4659 0.6671 0 1 0 1 Y 9 1 1 95.0 0.2784 0.1671 0.3837 0 1 0 1 Y 9 1 2 95.0 0.1511 0.0314 0.2811 0 1 0 1 Y 9 2 1 95.0 0.1511 0.0314 0.2811 0 1 0 1 Y 9 2 2 95.0 0.2367 0.1478 0.3382 0 1 0 1 Z 9 1 1 95.0 0.2060 0.1825 0.2329 6 1 6 1 Z 9 1 2 95.0 0.1602 0.1336 0.1908 0 1 0 1 Z 9 2 1 95.0 0.1602 0.1336 0.1908 0 1 0 1 Z 9 2 2 95.0 0.1900 0.1664 0.2173 0 1 0 1 This problem used 11.0% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 2:25.16 Execution 0: 0:35:43.35 TOTAL 0: 0:38: 8.51 Total number of warnings issued: 0 ______________________________________________________________________________