!job for Falconer GxE model !first simulate the bivariate data and then analyze the data !next test if the genet corr is different from 1 or 0 ! Fit the wrong model to the simulated data: rg=1 ! Use the constraint option in Mx to constrain the genetic correlation at 1 The following MX script lines were read for group 1 G1: MODEL PARAMETERS DATA CALC NGROUPS=4 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 I FULL 2 2 FI ! TO BE USED IN CONSTRAINT GROUP END MATRICES; 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'); ! GENETIC CORRELATION (USED IN THE CONSTRAINT) L=\STND(C); ! COMMON E CORRELATION M=\STND(E); ! UNIQUE E CORRELATION END ALGEBRA; MATRIX I 1 0 0 1 ! I HAS '1' ON THE DIAGONAL AND OFF-DIAGONAL ELEMENTS 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 OPTIONS RSIDUALS POWER= .05,1 ! .05 SIG LEVEL & 1 DF END The following MX script lines were read for group 4 G4: CONSTRAINT GROUP TO ENSURE GENET CORR=1 CONSTRAINT_GROUP MATRICES = GROUP 1 CO I - K ; !CONSTRAINT VALUE SHOULD BE ZERO OPTION RS 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 I This is a 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 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: 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 I This is a 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 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: 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 I This is a 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 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: CONSTRAINT GROUP TO ENSURE genet corr=1 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 I This is a 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 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 Mx starting optimization; number of parameters = 9 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 1.6619E-01 5.9207E-27 2 5.9207E-27 1.5938E-01 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 I This is a FULL matrix of order 2 by 2 1 2 1 1.0000 0.0000 2 0.0000 1.0000 MATRIX K This is a computed FULL matrix of order 2 by 2 [=\STND(X*X')] 1 2 1 1.0000E+00 3.6380E-26 2 3.6380E-26 1.0000E+00 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 4.0766E-01 2 1.4523E-26 3.9922E-01 MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 2.7179E-01 2 2.7031E-01 -4.4687E-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 1.6619E-01 5.9207E-27 2 5.9207E-27 1.5938E-01 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 I This is a FULL matrix of order 2 by 2 1 2 1 1.0000 0.0000 2 0.0000 1.0000 MATRIX K This is a computed FULL matrix of order 2 by 2 [=\STND(X*X')] 1 2 1 1.0000E+00 3.6380E-26 2 3.6380E-26 1.0000E+00 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 4.0766E-01 2 1.4523E-26 3.9922E-01 MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 2.7179E-01 2 2.7031E-01 -4.4687E-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.1863E-04 2 9.0034E-03 -6.8959E-05 3 1.2440E-02 2.9032E-02 -1.1863E-04 4 2.9032E-02 1.2557E-02 9.0034E-03 -6.8959E-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 1.6619E-01 5.9207E-27 2 5.9207E-27 1.5938E-01 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 I This is a FULL matrix of order 2 by 2 1 2 1 1.0000 0.0000 2 0.0000 1.0000 MATRIX K This is a computed FULL matrix of order 2 by 2 [=\STND(X*X')] 1 2 1 1.0000E+00 3.6380E-26 2 3.6380E-26 1.0000E+00 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 4.0766E-01 2 1.4523E-26 3.9922E-01 MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 2.7179E-01 2 2.7031E-01 -4.4687E-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.1275 0.0525 0.4550 4 0.0525 0.1237 0.1025 0.4475 EXPECTED COVARIANCE MATRIX 1 2 3 4 1 0.4551 2 0.0935 0.4476 3 0.1570 0.0735 0.4551 4 0.0735 0.1528 0.0935 0.4476 RESIDUAL MATRIX 1 2 3 4 1 -1.1863E-04 2 9.0034E-03 -6.8959E-05 3 -2.9465E-02 -2.0968E-02 -1.1863E-04 4 -2.0968E-02 -2.9005E-02 9.0034E-03 -6.8959E-05 Function value of this group: 1.7642 Where the fit function is Maximum Likelihood GROUP NUMBER: 4 G4: CONSTRAINT GROUP TO ENSURE genet corr=1 MATRIX A This is a computed FULL matrix of order 2 by 2 [=X*X'] 1 2 1 1.6619E-01 5.9207E-27 2 5.9207E-27 1.5938E-01 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 I This is a FULL matrix of order 2 by 2 1 2 1 1.0000 0.0000 2 0.0000 1.0000 MATRIX K This is a computed FULL matrix of order 2 by 2 [=\STND(X*X')] 1 2 1 1.0000E+00 3.6380E-26 2 3.6380E-26 1.0000E+00 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 4.0766E-01 2 1.4523E-26 3.9922E-01 MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 2.7179E-01 2 2.7031E-01 -4.4687E-08 MATRIX Z This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.4637 2 0.0432 0.4618 CONSTRAINT VALUES (should be near zero) 1 2 1 0.0000E+00 -3.6380E-26 2 -3.6380E-26 0.0000E+00 Your model has 9 estimated parameters and 24 Observed statistics Observed statistics include 4 constraints. Chi-squared fit of model >>>>>>> 3.176 Degrees of freedom >>>>>>>>>>>>> 15 Probability >>>>>>>>>>>>>>>>>>>> 0.999 Akaike's Information Criterion > -26.824 RMSEA >>>>>>>>>>>>>>>>>>>>>>>>>> 0.000 Power of this test, at the 0.0500 significance level with 1. df is 0.429488 Based on your combined observed sample size of 200. The following sample sizes would be required to reject the hypothesis: Power Total N .25 104. (approximately) .50 242. (approximately) .75 437. (approximately) .80 494. (approximately) .90 662. (approximately) .95 818. (approximately) .99 1157. (approximately) This problem used 1.0% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.00 Execution 0: 0: 0: 0.50 TOTAL 0: 0: 0: 0.50 Total number of warnings issued: 0 ______________________________________________________________________________