! My first Mx Model ! Univariate example testing ADE/ACE model ! Australian Females - Neuroticism (arcsin sqrt transformation) ! ! A - additive genetic variance (VA or a^2) ! C - shared environmental variance (VC or c^2) ! E - unique environmental variance (VE or e^2) ! D - dominance genetic variance (VD or d^2) ! ! Phenotypic variance = A + C + E + D ! MZ covariance = A + C + D ! DZ covariance = 0.5 A + C + 0.25 D ! ! Path Coefficients Model ! ! Group 1 calculates variance components for females ! Group 2 fits model to MZ female data ! Group 3 fits model to DZ female data The following MX script lines were read for group 1 #DEFINE NVAR 1 ! Mx replaces every occurrence of nvar with 1 #NGROUP 3 Note: #NGroup set number of groups to 3 TITLE GROUP 1: PARAMETERS OF MODEL CALCULATION BEGIN MATRICES; X FULL NVAR NVAR FREE ! A: ADDITIVE GENETIC PARAMETER Y FULL NVAR NVAR FREE ! C: SHARED ENVIRONMENTAL PARAMETER Z FULL NVAR NVAR FREE ! E: UNIQUE ENVIRONMENTAL PARAMETER W FULL NVAR NVAR ! D: DOMINANCE PARAMETER H FULL 1 1 ! SCALAR, .5 Q FULL 1 1 ! SCALAR, .25 END MATRICES; MATRIX H .5 MATRIX Q .25 START .1 ALL BEGIN ALGEBRA; A=X*X'; ! A^2: ADDITIVE GENETIC VARIANCE C=Y*Y'; ! C^2: SHARED ENVIRONMENTAL VARIANCE E=Z*Z'; ! E^2: UNIQUE ENVIRONMENTAL VARIANCE D=W*W'; ! D^2: DOMINANCE VARIANCE V=A+C+E+D; ! TOTAL VARIANCE P=A|C|E|D; ! PUT PARAMETER ESTIMATES IN ONE MATRIX ! the | operator concatenates matrices with same number of rows S=P@V~; ! STANDARDIZED PARAMETER ESTIMATES ! the @ kronecker product multiplies every element of the first matrix ! with the second matrix, ~ takes the inverse a matrix END ALGEBRA; LABELS ROW X PAREST_A LABELS ROW Y PAREST_C LABELS ROW Z PAREST_E LABELS ROW W PAREST_D LABELS ROW A A^2 LABELS ROW C C^2 LABELS ROW E E^2 LABELS ROW D D^2 LABELS ROW V VARIANCE LABELS ROW P ESTIMATE LABELS COL P A C E D LABELS ROW S STANDEST LABELS COL S A^2 C^2 E^2 D^2 END The following MX script lines were read for group 2 TITLE G2: FEMALE MZ TWIN PAIRS Note: Opening #include file 1 neurmzf.dat ! Neuroticism (arcsin sqrt) covariances - MZ FEMALES DATA NINPUT_VARS=2 NOBSERVATIONS=1233 CMATRIX .691992D-01 .343410D-01 .637396D-01 LABELS TW1-N TW2-N Note: Closing #include file 1 MATRICES = GROUP 1 COVARIANCES A+C+E+D | A+C+D _ A+C+D | A+C+E+D ; ! the _ operator concatenates matrices with same number of columns OPTION RSIDUALS ! PRINTS OBSERVED, EXPECTED & RESIDUAL MATRIX END The following MX script lines were read for group 3 TITLE G3: FEMALE DZ TWIN PAIRS Note: Opening #include file 1 neurdzf.dat ! Neuroticism (arcsin sqrt) covariances - DZ FEMALES DATA NINPUT_VARS=2 NOBSERVATIONS=750 CMATRIX .651865D-01 .183778D-01 .695700D-01 LABELS TW1-N TW2-N Note: Closing #include file 1 MATRICES= GROUP 1 COVARIANCES A+C+E+D | H@A+C+Q@D _ H@A+C+Q@D | A+C+E+D ; OPTION RSIDUALS OPTION NDECIMALS=4 OPTION MULTIPLE END PARAMETER SPECIFICATIONS GROUP NUMBER: 1 Title Group 1: parameters of model MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified 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 H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX P This is a computed FULL matrix of order 1 by 4 It has no free parameters specified MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX S This is a computed FULL matrix of order 1 by 4 It has no free parameters specified MATRIX V This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX W This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a FULL matrix of order 1 by 1 1 PAREST_A 1 MATRIX Y This is a FULL matrix of order 1 by 1 1 PAREST_C 2 MATRIX Z This is a FULL matrix of order 1 by 1 1 PAREST_E 3 GROUP NUMBER: 2 Title G2: female MZ twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified 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 H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX P This is a computed FULL matrix of order 1 by 4 It has no free parameters specified MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX S This is a computed FULL matrix of order 1 by 4 It has no free parameters specified MATRIX V This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX W This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a FULL matrix of order 1 by 1 1 PAREST_A 1 MATRIX Y This is a FULL matrix of order 1 by 1 1 PAREST_C 2 MATRIX Z This is a FULL matrix of order 1 by 1 1 PAREST_E 3 GROUP NUMBER: 3 Title G3: female DZ twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified 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 H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX P This is a computed FULL matrix of order 1 by 4 It has no free parameters specified MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX S This is a computed FULL matrix of order 1 by 4 It has no free parameters specified MATRIX V This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX W This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a FULL matrix of order 1 by 1 1 PAREST_A 1 MATRIX Y This is a FULL matrix of order 1 by 1 1 PAREST_C 2 MATRIX Z This is a FULL matrix of order 1 by 1 1 PAREST_E 3 MX PARAMETER ESTIMATES GROUP NUMBER: 1 Title Group 1: parameters of model MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 A^2 0.0331 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 C^2 0.0016 MATRIX D This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 D^2 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 E^2 0.0322 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX P This is a computed FULL matrix of order 1 by 4 [=A|C|E|D] A C E D ESTIMATE 0.0331 0.0016 0.0322 0.0000 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX S This is a computed FULL matrix of order 1 by 4 [=P@V~] A^2 C^2 E^2 D^2 STANDEST 0.4949 0.0235 0.4816 0.0000 MATRIX V This is a computed FULL matrix of order 1 by 1 [=A+C+E+D] 1 VARIANCE 0.0668 MATRIX W This is a FULL matrix of order 1 by 1 1 PAREST_D 0.0000 MATRIX X This is a FULL matrix of order 1 by 1 1 PAREST_A 0.1819 MATRIX Y This is a FULL matrix of order 1 by 1 1 PAREST_C 0.0396 MATRIX Z This is a FULL matrix of order 1 by 1 1 PAREST_E 0.1794 GROUP NUMBER: 2 Title G2: female MZ twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 A^2 0.0331 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 C^2 0.0016 MATRIX D This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 D^2 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 E^2 0.0322 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX P This is a computed FULL matrix of order 1 by 4 [=A|C|E|D] A C E D ESTIMATE 0.0331 0.0016 0.0322 0.0000 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX S This is a computed FULL matrix of order 1 by 4 [=P@V~] A^2 C^2 E^2 D^2 STANDEST 0.4949 0.0235 0.4816 0.0000 MATRIX V This is a computed FULL matrix of order 1 by 1 [=A+C+E+D] 1 VARIANCE 0.0668 MATRIX W This is a FULL matrix of order 1 by 1 1 PAREST_D 0.0000 MATRIX X This is a FULL matrix of order 1 by 1 1 PAREST_A 0.1819 MATRIX Y This is a FULL matrix of order 1 by 1 1 PAREST_C 0.0396 MATRIX Z This is a FULL matrix of order 1 by 1 1 PAREST_E 0.1794 OBSERVED COVARIANCE MATRIX TW1-N TW2-N TW1-N 0.0692 TW2-N 0.0343 0.0637 EXPECTED COVARIANCE MATRIX TW1-N TW2-N TW1-N 0.0668 TW2-N 0.0347 0.0668 RESIDUAL MATRIX TW1-N TW2-N TW1-N 0.0023 TW2-N -0.0003 -0.0031 Function value of this group: 2.8693 Where the fit function is Maximum Likelihood GROUP NUMBER: 3 Title G3: female DZ twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 A^2 0.0331 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 C^2 0.0016 MATRIX D This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 D^2 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 E^2 0.0322 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX P This is a computed FULL matrix of order 1 by 4 [=A|C|E|D] A C E D ESTIMATE 0.0331 0.0016 0.0322 0.0000 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX S This is a computed FULL matrix of order 1 by 4 [=P@V~] A^2 C^2 E^2 D^2 STANDEST 0.4949 0.0235 0.4816 0.0000 MATRIX V This is a computed FULL matrix of order 1 by 1 [=A+C+E+D] 1 VARIANCE 0.0668 MATRIX W This is a FULL matrix of order 1 by 1 1 PAREST_D 0.0000 MATRIX X This is a FULL matrix of order 1 by 1 1 PAREST_A 0.1819 MATRIX Y This is a FULL matrix of order 1 by 1 1 PAREST_C 0.0396 MATRIX Z This is a FULL matrix of order 1 by 1 1 PAREST_E 0.1794 OBSERVED COVARIANCE MATRIX TW1-N TW2-N TW1-N 0.0652 TW2-N 0.0184 0.0696 EXPECTED COVARIANCE MATRIX TW1-N TW2-N TW1-N 0.0668 TW2-N 0.0181 0.0668 RESIDUAL MATRIX TW1-N TW2-N TW1-N -0.0017 TW2-N 0.0003 0.0027 Function value of this group: 0.9002 Where the fit function is Maximum Likelihood Your model has 3 estimated parameters and 6 Observed statistics Chi-squared fit of model >>>>>>> 3.769 Degrees of freedom >>>>>>>>>>>>> 3 Probability >>>>>>>>>>>>>>>>>>>> 0.287 Akaike's Information Criterion > -2.231 RMSEA >>>>>>>>>>>>>>>>>>>>>>>>>> 0.014 This problem used 0.5% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.50 Execution 0: 0: 0: 0.11 TOTAL 0: 0: 0: 0.61 Total number of warnings issued: 0 ______________________________________________________________________________ Multiple fit option in effect. The following MX script lines have been read: SAVE ACE.MXS DROP Y 1 1 1 ! DROP C END PARAMETER SPECIFICATIONS GROUP NUMBER: 1 Title Group 1: parameters of model MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified 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 H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX P This is a computed FULL matrix of order 1 by 4 It has no free parameters specified MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX S This is a computed FULL matrix of order 1 by 4 It has no free parameters specified MATRIX V This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX W This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a FULL matrix of order 1 by 1 1 PAREST_A 1 MATRIX Y This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX Z This is a FULL matrix of order 1 by 1 1 PAREST_E 3 GROUP NUMBER: 2 Title G2: female MZ twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified 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 H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX P This is a computed FULL matrix of order 1 by 4 It has no free parameters specified MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX S This is a computed FULL matrix of order 1 by 4 It has no free parameters specified MATRIX V This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX W This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a FULL matrix of order 1 by 1 1 PAREST_A 1 MATRIX Y This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX Z This is a FULL matrix of order 1 by 1 1 PAREST_E 3 GROUP NUMBER: 3 Title G3: female DZ twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified 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 H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX P This is a computed FULL matrix of order 1 by 4 It has no free parameters specified MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX S This is a computed FULL matrix of order 1 by 4 It has no free parameters specified MATRIX V This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX W This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a FULL matrix of order 1 by 1 1 PAREST_A 1 MATRIX Y This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX Z This is a FULL matrix of order 1 by 1 1 PAREST_E 3 *** WARNING! *** I am not sure I have found a solution that satisfies Kuhn-Tucker conditions for a minimum. NAG's IFAIL parameter is 6 Looks like I got stuck here. Check the following: 1. The model is correctly specified 2. Starting values are good 3. You are not already at the solution The error can arise if the Hessian is ill-conditioned You can try resetting it to an identity matrix and fit from the solution by putting TH=-n on the OU line where n is the number of refits that you want to do If all else fails try putting NAG=30 on the OU line and examine the file NAGDUMP.OUT and the NAG manual MX PARAMETER ESTIMATES GROUP NUMBER: 1 Title Group 1: parameters of model MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 A^2 0.0347 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 C^2 0.0000 MATRIX D This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 D^2 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 E^2 0.0321 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX P This is a computed FULL matrix of order 1 by 4 [=A|C|E|D] A C E D ESTIMATE 0.0347 0.0000 0.0321 0.0000 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX S This is a computed FULL matrix of order 1 by 4 [=P@V~] A^2 C^2 E^2 D^2 STANDEST 0.5199 0.0000 0.4801 0.0000 MATRIX V This is a computed FULL matrix of order 1 by 1 [=A+C+E+D] 1 VARIANCE 0.0668 MATRIX W This is a FULL matrix of order 1 by 1 1 PAREST_D 0.0000 MATRIX X This is a FULL matrix of order 1 by 1 1 PAREST_A 0.1864 MATRIX Y This is a FULL matrix of order 1 by 1 1 PAREST_C 0.0000 MATRIX Z This is a FULL matrix of order 1 by 1 1 PAREST_E 0.1791 GROUP NUMBER: 2 Title G2: female MZ twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 A^2 0.0347 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 C^2 0.0000 MATRIX D This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 D^2 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 E^2 0.0321 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX P This is a computed FULL matrix of order 1 by 4 [=A|C|E|D] A C E D ESTIMATE 0.0347 0.0000 0.0321 0.0000 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX S This is a computed FULL matrix of order 1 by 4 [=P@V~] A^2 C^2 E^2 D^2 STANDEST 0.5199 0.0000 0.4801 0.0000 MATRIX V This is a computed FULL matrix of order 1 by 1 [=A+C+E+D] 1 VARIANCE 0.0668 MATRIX W This is a FULL matrix of order 1 by 1 1 PAREST_D 0.0000 MATRIX X This is a FULL matrix of order 1 by 1 1 PAREST_A 0.1864 MATRIX Y This is a FULL matrix of order 1 by 1 1 PAREST_C 0.0000 MATRIX Z This is a FULL matrix of order 1 by 1 1 PAREST_E 0.1791 OBSERVED COVARIANCE MATRIX TW1-N TW2-N TW1-N 0.0692 TW2-N 0.0343 0.0637 EXPECTED COVARIANCE MATRIX TW1-N TW2-N TW1-N 0.0668 TW2-N 0.0347 0.0668 RESIDUAL MATRIX TW1-N TW2-N TW1-N 0.0024 TW2-N -0.0004 -0.0031 Function value of this group: 2.8733 Where the fit function is Maximum Likelihood GROUP NUMBER: 3 Title G3: female DZ twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 A^2 0.0347 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 C^2 0.0000 MATRIX D This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 D^2 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 E^2 0.0321 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX P This is a computed FULL matrix of order 1 by 4 [=A|C|E|D] A C E D ESTIMATE 0.0347 0.0000 0.0321 0.0000 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX S This is a computed FULL matrix of order 1 by 4 [=P@V~] A^2 C^2 E^2 D^2 STANDEST 0.5199 0.0000 0.4801 0.0000 MATRIX V This is a computed FULL matrix of order 1 by 1 [=A+C+E+D] 1 VARIANCE 0.0668 MATRIX W This is a FULL matrix of order 1 by 1 1 PAREST_D 0.0000 MATRIX X This is a FULL matrix of order 1 by 1 1 PAREST_A 0.1864 MATRIX Y This is a FULL matrix of order 1 by 1 1 PAREST_C 0.0000 MATRIX Z This is a FULL matrix of order 1 by 1 1 PAREST_E 0.1791 OBSERVED COVARIANCE MATRIX TW1-N TW2-N TW1-N 0.0652 TW2-N 0.0184 0.0696 EXPECTED COVARIANCE MATRIX TW1-N TW2-N TW1-N 0.0668 TW2-N 0.0174 0.0668 RESIDUAL MATRIX TW1-N TW2-N TW1-N -0.0016 TW2-N 0.0010 0.0028 Function value of this group: 1.0155 Where the fit function is Maximum Likelihood *** WARNING! *** Minimization may not be successful. See above CODE RED - Hessian/precision problem Your model has 2 estimated parameters and 6 Observed statistics Chi-squared fit of model >>>>>>> 3.889 Degrees of freedom >>>>>>>>>>>>> 4 Probability >>>>>>>>>>>>>>>>>>>> 0.421 Akaike's Information Criterion > -4.111 RMSEA >>>>>>>>>>>>>>>>>>>>>>>>>> 0.009 This problem used 0.5% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.05 Execution 0: 0: 0: 0.28 TOTAL 0: 0: 0: 0.33 Total number of warnings issued: 2 ______________________________________________________________________________ Multiple fit option in effect. The following MX script lines have been read: GET ACE.MXS Note: 1351 binary records read DROP X 1 1 1 ! DROP A END PARAMETER SPECIFICATIONS GROUP NUMBER: 1 Title Group 1: parameters of model MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified 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 H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX P This is a computed FULL matrix of order 1 by 4 It has no free parameters specified MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX S This is a computed FULL matrix of order 1 by 4 It has no free parameters specified MATRIX V This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX W This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX Y This is a FULL matrix of order 1 by 1 1 PAREST_C 2 MATRIX Z This is a FULL matrix of order 1 by 1 1 PAREST_E 3 GROUP NUMBER: 2 Title G2: female MZ twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified 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 H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX P This is a computed FULL matrix of order 1 by 4 It has no free parameters specified MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX S This is a computed FULL matrix of order 1 by 4 It has no free parameters specified MATRIX V This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX W This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX Y This is a FULL matrix of order 1 by 1 1 PAREST_C 2 MATRIX Z This is a FULL matrix of order 1 by 1 1 PAREST_E 3 GROUP NUMBER: 3 Title G3: female DZ twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified 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 H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX P This is a computed FULL matrix of order 1 by 4 It has no free parameters specified MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX S This is a computed FULL matrix of order 1 by 4 It has no free parameters specified MATRIX V This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX W This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX Y This is a FULL matrix of order 1 by 1 1 PAREST_C 2 MATRIX Z This is a FULL matrix of order 1 by 1 1 PAREST_E 3 MX PARAMETER ESTIMATES GROUP NUMBER: 1 Title Group 1: parameters of model MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 A^2 0.0000 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 C^2 0.0283 MATRIX D This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 D^2 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 E^2 0.0385 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX P This is a computed FULL matrix of order 1 by 4 [=A|C|E|D] A C E D ESTIMATE 0.0000 0.0283 0.0385 0.0000 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX S This is a computed FULL matrix of order 1 by 4 [=P@V~] A^2 C^2 E^2 D^2 STANDEST 0.0000 0.4237 0.5763 0.0000 MATRIX V This is a computed FULL matrix of order 1 by 1 [=A+C+E+D] 1 VARIANCE 0.0668 MATRIX W This is a FULL matrix of order 1 by 1 1 PAREST_D 0.0000 MATRIX X This is a FULL matrix of order 1 by 1 1 PAREST_A 0.0000 MATRIX Y This is a FULL matrix of order 1 by 1 1 PAREST_C 0.1682 MATRIX Z This is a FULL matrix of order 1 by 1 1 PAREST_E 0.1962 GROUP NUMBER: 2 Title G2: female MZ twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 A^2 0.0000 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 C^2 0.0283 MATRIX D This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 D^2 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 E^2 0.0385 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX P This is a computed FULL matrix of order 1 by 4 [=A|C|E|D] A C E D ESTIMATE 0.0000 0.0283 0.0385 0.0000 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX S This is a computed FULL matrix of order 1 by 4 [=P@V~] A^2 C^2 E^2 D^2 STANDEST 0.0000 0.4237 0.5763 0.0000 MATRIX V This is a computed FULL matrix of order 1 by 1 [=A+C+E+D] 1 VARIANCE 0.0668 MATRIX W This is a FULL matrix of order 1 by 1 1 PAREST_D 0.0000 MATRIX X This is a FULL matrix of order 1 by 1 1 PAREST_A 0.0000 MATRIX Y This is a FULL matrix of order 1 by 1 1 PAREST_C 0.1682 MATRIX Z This is a FULL matrix of order 1 by 1 1 PAREST_E 0.1962 OBSERVED COVARIANCE MATRIX TW1-N TW2-N TW1-N 0.0692 TW2-N 0.0343 0.0637 EXPECTED COVARIANCE MATRIX TW1-N TW2-N TW1-N 0.0668 TW2-N 0.0283 0.0668 RESIDUAL MATRIX TW1-N TW2-N TW1-N 0.0024 TW2-N 0.0060 -0.0031 Function value of this group: 23.9991 Where the fit function is Maximum Likelihood GROUP NUMBER: 3 Title G3: female DZ twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 A^2 0.0000 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 C^2 0.0283 MATRIX D This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 D^2 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 E^2 0.0385 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX P This is a computed FULL matrix of order 1 by 4 [=A|C|E|D] A C E D ESTIMATE 0.0000 0.0283 0.0385 0.0000 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX S This is a computed FULL matrix of order 1 by 4 [=P@V~] A^2 C^2 E^2 D^2 STANDEST 0.0000 0.4237 0.5763 0.0000 MATRIX V This is a computed FULL matrix of order 1 by 1 [=A+C+E+D] 1 VARIANCE 0.0668 MATRIX W This is a FULL matrix of order 1 by 1 1 PAREST_D 0.0000 MATRIX X This is a FULL matrix of order 1 by 1 1 PAREST_A 0.0000 MATRIX Y This is a FULL matrix of order 1 by 1 1 PAREST_C 0.1682 MATRIX Z This is a FULL matrix of order 1 by 1 1 PAREST_E 0.1962 OBSERVED COVARIANCE MATRIX TW1-N TW2-N TW1-N 0.0652 TW2-N 0.0184 0.0696 EXPECTED COVARIANCE MATRIX TW1-N TW2-N TW1-N 0.0668 TW2-N 0.0283 0.0668 RESIDUAL MATRIX TW1-N TW2-N TW1-N -0.0016 TW2-N -0.0099 0.0028 Function value of this group: 28.3462 Where the fit function is Maximum Likelihood Your model has 2 estimated parameters and 6 Observed statistics Chi-squared fit of model >>>>>>> 52.345 Degrees of freedom >>>>>>>>>>>>> 4 Probability >>>>>>>>>>>>>>>>>>>> 0.000 Akaike's Information Criterion > 44.345 RMSEA >>>>>>>>>>>>>>>>>>>>>>>>>> 0.113 This problem used 0.5% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.00 Execution 0: 0: 0: 0.11 TOTAL 0: 0: 0: 0.11 Total number of warnings issued: 2 ______________________________________________________________________________