!job for Falconer GxE model ! Analyse the simulated data The following MX script lines were read for group 1 G1: MODEL PARAMETERS DATA CALC NGROUPS=3 BEGIN MATRICES; X LOWER 2 2 FR ! GENETIC STRUCTURE Y LOWER 2 2 FR ! COMMON ENVIRONMENTAL STRUCTURE Z LOWER 2 2 FR ! SPECIFIC ENVIRONMENTAL STRUCTURE END MATRICES; FIX X 1 2 1 START .3 X 1 1 1 X 1 2 2 START .1 Y 1 1 1 Y 1 2 2 START .5 Z 1 1 1 Z 1 2 2 BEGIN ALGEBRA; A= X*X' ; C= Y*Y' ; E= Z*Z' ; K=\STND(X*X'); ! CALCULATE THE GENETIC CORRELATION L=\STND(C); ! COMMON E CORRELATION M=\STND(E); ! UNIQUE E CORRELATION END ALGEBRA; END The following MX script lines were read for group 2 G2: MZ TWIN PAIRS DATA NINPUT_VARS=4 NOBSERVATIONS=100 CMATRIX FULL FILE=MZSIM.COV MATRICES= GROUP 1 COVARIANCES A+C+E | A+C _ A+C | A+C+E / OPTION RSIDUALS END The following MX script lines were read for group 3 G3: DZ TWIN PAIRS DATA NINPUT_VARS=4 NOBSERVATIONS=100 CMATRIX FULL FILE=DZSIM.COV MATRICES= GROUP 1 H FULL 1 1 COVARIANCES A+C+E | H@A+C _ H@A+C | A+C+E / MATRIX H .5 END PARAMETER SPECIFICATIONS GROUP NUMBER: 1 G1: model parameters 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 K This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX L This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 1 2 0 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: 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 K This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX L This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 1 2 0 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: 3 G3: 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 K This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX L This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 1 2 0 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 Mx starting optimization; number of parameters = 8 MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1: model parameters MATRIX A This is a computed FULL matrix of order 2 by 2 [=X*X'] 1 2 1 0.1662 0.0000 2 0.0000 0.1594 MATRIX C This is a computed FULL matrix of order 2 by 2 [=Y*Y'] 1 2 1 0.0739 0.0735 2 0.0735 0.0731 MATRIX E This is a computed FULL matrix of order 2 by 2 [=Z*Z'] 1 2 1 0.2151 0.0200 2 0.0200 0.2151 MATRIX K This is a computed FULL matrix of order 2 by 2 [=\STND(X*X')] 1 2 1 1.0000 0.0000 2 0.0000 1.0000 MATRIX L This is a computed FULL matrix of order 2 by 2 [=\STND(C)] 1 2 1 1.0000 1.0000 2 1.0000 1.0000 MATRIX M This is a computed FULL matrix of order 2 by 2 [=\STND(E)] 1 2 1 1.0000 0.0931 2 0.0931 1.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.4077 2 0.0000 0.3992 MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 2.7179E-01 2 2.7031E-01 1.0688E-08 MATRIX Z This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.4637 2 0.0432 0.4618 GROUP NUMBER: 2 G2: MZ twin pairs MATRIX A This is a computed FULL matrix of order 2 by 2 [=X*X'] 1 2 1 0.1662 0.0000 2 0.0000 0.1594 MATRIX C This is a computed FULL matrix of order 2 by 2 [=Y*Y'] 1 2 1 0.0739 0.0735 2 0.0735 0.0731 MATRIX E This is a computed FULL matrix of order 2 by 2 [=Z*Z'] 1 2 1 0.2151 0.0200 2 0.0200 0.2151 MATRIX K This is a computed FULL matrix of order 2 by 2 [=\STND(X*X')] 1 2 1 1.0000 0.0000 2 0.0000 1.0000 MATRIX L This is a computed FULL matrix of order 2 by 2 [=\STND(C)] 1 2 1 1.0000 1.0000 2 1.0000 1.0000 MATRIX M This is a computed FULL matrix of order 2 by 2 [=\STND(E)] 1 2 1 1.0000 0.0931 2 0.0931 1.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.4077 2 0.0000 0.3992 MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 2.7179E-01 2 2.7031E-01 1.0688E-08 MATRIX Z This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.4637 2 0.0432 0.4618 OBSERVED COVARIANCE MATRIX 1 2 3 4 1 0.4550 2 0.1025 0.4475 3 0.2525 0.1025 0.4550 4 0.1025 0.2450 0.1025 0.4475 EXPECTED COVARIANCE MATRIX 1 2 3 4 1 0.4551 2 0.0935 0.4476 3 0.2401 0.0735 0.4551 4 0.0735 0.2324 0.0935 0.4476 RESIDUAL MATRIX 1 2 3 4 1 -1.1860E-04 2 9.0034E-03 -6.8958E-05 3 1.2440E-02 2.9032E-02 -1.1860E-04 4 2.9032E-02 1.2557E-02 9.0034E-03 -6.8958E-05 Function value of this group: 1.4115 Where the fit function is Maximum Likelihood GROUP NUMBER: 3 G3: DZ twin pairs MATRIX A This is a computed FULL matrix of order 2 by 2 [=X*X'] 1 2 1 0.1662 0.0000 2 0.0000 0.1594 MATRIX C This is a computed FULL matrix of order 2 by 2 [=Y*Y'] 1 2 1 0.0739 0.0735 2 0.0735 0.0731 MATRIX E This is a computed FULL matrix of order 2 by 2 [=Z*Z'] 1 2 1 0.2151 0.0200 2 0.0200 0.2151 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 2 by 2 [=\STND(X*X')] 1 2 1 1.0000 0.0000 2 0.0000 1.0000 MATRIX L This is a computed FULL matrix of order 2 by 2 [=\STND(C)] 1 2 1 1.0000 1.0000 2 1.0000 1.0000 MATRIX M This is a computed FULL matrix of order 2 by 2 [=\STND(E)] 1 2 1 1.0000 0.0931 2 0.0931 1.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.4077 2 0.0000 0.3992 MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 2.7179E-01 2 2.7031E-01 1.0688E-08 MATRIX Z This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.4637 2 0.0432 0.4618 Your model has 8 estimated parameters and 20 Observed statistics Chi-squared fit of model >>>>>>> 3.176 Degrees of freedom >>>>>>>>>>>>> 12 Probability >>>>>>>>>>>>>>>>>>>> 0.994 Akaike's Information Criterion > -20.824 RMSEA >>>>>>>>>>>>>>>>>>>>>>>>>> 0.000 This problem used 0.7% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.00 Execution 0: 0: 0: 0.05 TOTAL 0: 0: 0: 0.05 Total number of warnings issued: 0 ______________________________________________________________________________