! simulate two-wave data on alcohol dependence for male MZ & DZ pairs The following MX script lines were read for group 1 #NG 1 Note: #NGroup set number of groups to 1 #DEFINE NVAR 2 #DEFINE NVAR2 4 #DEFINE MAXTHRES 1 ! MAX NUMBER OF THRESHOLDS FOR ANY DIMENSION DATA SIMULATION CALCULATION GROUP BEGIN MATRICES; X LO NVAR NVAR ! ADDITIVE GENETIC CHOLESKY Y LO NVAR NVAR ! SHARED ENVIRONMENTAL CHOLESKY Z LO NVAR NVAR ! NON-SHARED ENVIRONMENTAL CHOLESKY N FU 1 NVAR2 ! MATRIX TO STORE NUMBER OF THRESHOLDS FOR EACH DIMENSION M FU 1 NVAR2 ! MATRIX OF MEANS S FU MAXTHRES NVAR2 ! THRESHOLDS G DI 1 1 U DI 1 1 END MATRICES; MAT G 0.5 MAT U 1000 MAT X 0.7071 0.7071 0.0 MAT Y 0.3 0.3 0.0 MAT Z 0.6403 0.4528 0.4528 MAT N 1 1 1 1 MAT M 0 0 0 0 MAT S 0.67449 0.67449 0.67449 0.67449 BEGIN ALGEBRA; A=X*X'; C=Y*Y'; E=Z*Z'; P=A+C+E; R=P | G@A+C _ G@A'+C' | P ; Q = \ALLINT (R_M_N_S); V = U@Q; END ALGEBRA; OPT FUNC=1.E-12 WRITE Q XOUT2526DZ.REC ! opt mxq=out2526 END GROUP PARAMETER SPECIFICATIONS No free parameters specified MX PARAMETER ESTIMATES GROUP NUMBER: 1 data simulation Estimates of matrix Q will be written to file: out2526dz.rec MATRIX A This is a computed FULL matrix of order 2 by 2 [=X*X'] 1 2 1 0.5000 0.5000 2 0.5000 0.5000 MATRIX C This is a computed FULL matrix of order 2 by 2 [=Y*Y'] 1 2 1 0.0900 0.0900 2 0.0900 0.0900 MATRIX E This is a computed FULL matrix of order 2 by 2 [=Z*Z'] 1 2 1 0.4100 0.2899 2 0.2899 0.4101 MATRIX G This is a DIAGONAL matrix of order 1 by 1 1 1 0.5000 MATRIX M This is a FULL matrix of order 1 by 4 1 2 3 4 1 0.0000 0.0000 0.0000 0.0000 MATRIX N This is a FULL matrix of order 1 by 4 1 2 3 4 1 1.0000 1.0000 1.0000 1.0000 MATRIX P This is a computed FULL matrix of order 2 by 2 [=A+C+E] 1 2 1 1.0000 0.8799 2 0.8799 1.0000 MATRIX Q This is a computed FULL matrix of order 1 by 16 [=\ALLINT(R_M_N_S)] 1 2 3 4 5 6 7 8 9 1 0.5195 0.0378 0.0378 0.0924 0.0378 0.0047 0.0047 0.0153 0.0378 10 11 12 13 14 15 16 1 0.0047 0.0047 0.0153 0.0924 0.0153 0.0153 0.0646 MATRIX R This is a computed FULL matrix of order 4 by 4 [=P|G@A+C_G@A'+C'|P] 1 2 3 4 1 1.0000 0.8799 0.3400 0.3400 2 0.8799 1.0000 0.3400 0.3400 3 0.3400 0.3400 1.0000 0.8799 4 0.3400 0.3400 0.8799 1.0000 MATRIX S This is a FULL matrix of order 1 by 4 1 2 3 4 1 0.6745 0.6745 0.6745 0.6745 MATRIX U This is a DIAGONAL matrix of order 1 by 1 1 1 1000.0000 MATRIX V This is a computed FULL matrix of order 1 by 16 [=U@Q] 1 2 3 4 5 6 7 1 519.5187 37.8212 37.8055 92.3817 37.8179 4.6910 4.6907 8 9 10 11 12 13 14 1 15.2703 37.8153 4.6918 4.6914 15.2792 92.3601 15.2865 15 16 1 15.2787 64.6498 MATRIX X This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.7071 2 0.7071 0.0000 MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.3000 2 0.3000 0.0000 MATRIX Z This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.6403 2 0.4528 0.4528 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.3% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.05 Execution 0: 0: 0: 0.12 TOTAL 0: 0: 0: 0.17 Total number of warnings issued: 0 ______________________________________________________________________________