** Mx startup successful ** **MX-PC 1.52** Job started on 03/06/02 at 12:15:00 !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 !Simulate the data The following MX script lines were read for group 1 TITLE CARD: GROUP 1: MODEL PARAMETERS DATA CALCULATION NGROUPS=3 BEGIN MATRICES; X LOWER 2 2 FIXED ! GENETIC STRUCTURE Y LOWER 2 2 FIXED ! COMMON ENVIRONMENTAL STRUCTURE Z LOWER 2 2 FIXED ! SPECIFIC ENVIRONMENTAL STRUCTURE END MATRICES; MATRIX X .50 !GENETIC CORRELATION IS DETERMINED BY OFF-DIAGONAL ELEMENT .20 .45 !AND THE GENETIC VARIANCES OF THE TWO TRAITS MATRIX Y .05 .05 0 ! RC IS 1 IN THIS MODEL MATRIX Z .45 0 .45 ! RE IS 0 IN THIS MODEL BEGIN ALGEBRA; A= X*X' ; !GENETIC VARIANCE AND COVARIANCE C= Y*Y' ; E= Z*Z' ; K=\STND(A); ! COMPUTE 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: SIMUTION FOR MZ TWIN PAIRS DATA CALC NINPUT_VARS=4 MATRICES= GROUP 1 COVARIANCES A+C+E | A+C _ A+C | A+C+E / OPTIONS MX%E=MZSIM.COV ! WRITE THE MZ MATRIX TO A FILE OPTIONS RS ! SHOW THE MZ MATRIX END The following MX script lines were read for group 3 G3: SIMULATION FOR DZ TWIN PAIRS DATA CALC NINPUT_VARS=4 MATRICES= GROUP 1 H FULL 1 1 COVARIANCES A+C+E | H@A+C _ H@A+C | A+C+E / MATRIX H .5 OPTIONS MX%E=DZSIM.COV OPTIONS RS END PARAMETER SPECIFICATIONS No free parameters specified MX PARAMETER ESTIMATES GROUP NUMBER: 1 title card: Group 1: model parameters MATRIX A This is a computed FULL matrix of order 2 by 2 [=X*X'] 1 2 1 0.2500 0.1000 2 0.1000 0.2425 MATRIX C This is a computed FULL matrix of order 2 by 2 [=Y*Y'] 1 2 1 0.0025 0.0025 2 0.0025 0.0025 MATRIX E This is a computed FULL matrix of order 2 by 2 [=Z*Z'] 1 2 1 0.2025 0.0000 2 0.0000 0.2025 MATRIX K This is a computed FULL matrix of order 2 by 2 [=\STND(A)] 1 2 1 1.0000 0.4061 2 0.4061 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.0000 2 0.0000 1.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.5000 2 0.2000 0.4500 MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.0500 2 0.0500 0.0000 MATRIX Z This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.4500 2 0.0000 0.4500 GROUP NUMBER: 2 G2: simution for MZ twin pairs Estimates of matrix % will be written to file: mzsim.cov MATRIX A This is a computed FULL matrix of order 2 by 2 [=X*X'] 1 2 1 0.2500 0.1000 2 0.1000 0.2425 MATRIX C This is a computed FULL matrix of order 2 by 2 [=Y*Y'] 1 2 1 0.0025 0.0025 2 0.0025 0.0025 MATRIX E This is a computed FULL matrix of order 2 by 2 [=Z*Z'] 1 2 1 0.2025 0.0000 2 0.0000 0.2025 MATRIX K This is a computed FULL matrix of order 2 by 2 [=\STND(A)] 1 2 1 1.0000 0.4061 2 0.4061 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.0000 2 0.0000 1.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.5000 2 0.2000 0.4500 MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.0500 2 0.0500 0.0000 MATRIX Z This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.4500 2 0.0000 0.4500 EXPECTED MATRIX of this CALCULATION group 1 2 3 4 1 0.4550 0.1025 0.2525 0.1025 2 0.1025 0.4475 0.1025 0.2450 3 0.2525 0.1025 0.4550 0.1025 4 0.1025 0.2450 0.1025 0.4475 GROUP NUMBER: 3 G3: simulation for DZ twin pairs Estimates of matrix % will be written to file: dzsim.cov MATRIX A This is a computed FULL matrix of order 2 by 2 [=X*X'] 1 2 1 0.2500 0.1000 2 0.1000 0.2425 MATRIX C This is a computed FULL matrix of order 2 by 2 [=Y*Y'] 1 2 1 0.0025 0.0025 2 0.0025 0.0025 MATRIX E This is a computed FULL matrix of order 2 by 2 [=Z*Z'] 1 2 1 0.2025 0.0000 2 0.0000 0.2025 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(A)] 1 2 1 1.0000 0.4061 2 0.4061 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.0000 2 0.0000 1.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.5000 2 0.2000 0.4500 MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.0500 2 0.0500 0.0000 MATRIX Z This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.4500 2 0.0000 0.4500 EXPECTED MATRIX of this CALCULATION group 1 2 3 4 1 0.4550 0.1025 0.1275 0.0525 2 0.1025 0.4475 0.0525 0.1238 3 0.1275 0.0525 0.4550 0.1025 4 0.0525 0.1238 0.1025 0.4475 Your model has 0 estimated parameters and 0 Observed statistics Chi-squared fit of model >>>>>>> 0.000 Degrees of freedom >>>>>>>>>>>>> 0 Probability incalculable Akaike's Information Criterion > 0.000 RMSEA >>>>>>>>>>>>>>>>>>>>>>>>>> 0.000 This problem used 0.6% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.55 Execution 0: 0: 0:-0.27 TOTAL 0: 0: 0: 0.28 Total number of warnings issued: 0 ______________________________________________________________________________