! 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 ! 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 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 2 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 2 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 2 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 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.11 TOTAL 0: 0: 0: 0.16 Total number of warnings issued: 0 ______________________________________________________________________________