! 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 | A+C _ A+C | P ; Q = \ALLINT (R_M_N_S); W = U@Q; END ALGEBRA; OPT FUNC=1.E-12 WRITE W OU2525MZ.REC ! opt mxu=ou2525n.rec ! opt mxw=ou2525.rec END PARAMETER SPECIFICATIONS No free parameters specified MX PARAMETER ESTIMATES GROUP NUMBER: 1 data simulation Estimates of matrix W will be written to file: u2525mz.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.5600 0.0336 0.0336 0.0604 0.0336 0.0058 0.0058 0.0172 0.0336 10 11 12 13 14 15 16 1 0.0058 0.0058 0.0172 0.0603 0.0172 0.0173 0.0927 MATRIX R This is a computed FULL matrix of order 4 by 4 [=P|A+C_A+C|P] 1 2 3 4 1 1.0000 0.8799 0.5900 0.5900 2 0.8799 1.0000 0.5900 0.5900 3 0.5900 0.5900 1.0000 0.8799 4 0.5900 0.5900 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 W This is a computed FULL matrix of order 1 by 16 [=U@Q] 1 2 3 4 5 6 7 1 559.9707 33.5915 33.5707 60.3566 33.5935 5.8234 5.8220 8 9 10 11 12 13 14 1 17.2429 33.5885 5.8245 5.8242 17.2356 60.3417 17.2447 15 16 1 17.2509 92.7041 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.06 Execution 0: 0: 0: 0.10 TOTAL 0: 0: 0: 0.16 Total number of warnings issued: 0 ______________________________________________________________________________