! Genetic Simplex Neale & Cardon (1992) Example #define NVAR 6 G1: Genetic structure Data Calc NGroups=3 Begin Matrices; X Diag NVAR NVAR Free ! genetic innovations Y Diag NVAR NVAR Free ! common environmental innovations Z Diag NVAR NVAR Free ! unique environmental innovations U Diag NVAR NVAR Free ! measurement error variances N Lower NVAR NVAR Free ! genetic (beta) transmissions O Lower NVAR NVAR Free ! common environmental (beta) transmissions P Lower NVAR NVAR Free ! unique environmental (beta) transmissions I Iden NVAR NVAR ! identity matrix H Full 1 1 End Matrices; Matrix H .5 Specify N 0 101 0 0 102 0 0 0 103 0 0 0 0 104 0 0 0 0 0 105 0 Specify O 0 106 0 0 107 0 0 0 108 0 0 0 0 109 0 0 0 0 0 110 0 Specify P 0 111 0 0 112 0 0 0 113 0 0 0 0 114 0 0 0 0 0 115 0 EQ u 1 1 u 6 6 !Measurement error variance the same at first and last occasion Begin Algebra; A = (I-N)~ & (X*X'); C = (I-O)~ & (Y*Y'); E = (I-P)~ & (Z*Z') + U*U'; End Algebra; End G2: Female MZ twin pairs Data NInput_vars=12 NObservations=32 Labels t1w1 t1w2 t1w3 t1w4 t1w5 t1w6 t2w1 t2w2 t2w3 t2w4 t2w5 t2w6 CMatrix 30.347 33.096 37.443 33.993 38.210 42.319 34.778 39.356 44.020 47.606 34.703 39.403 44.038 48.239 50.931 33.202 37.419 41.574 45.373 47.971 48.154 33.469 36.124 36.005 36.979 38.484 37.885 44.188 35.551 39.813 39.439 41.009 42.988 41.554 45.593 49.280 37.010 40.972 42.773 44.962 46.657 45.031 46.284 49.359 52.083 38.504 42.658 45.377 48.841 50.994 49.015 47.113 50.631 54.266 58.924 38.829 42.866 46.065 49.803 53.449 51.444 48.183 51.669 55.154 59.748 63.154 38.266 42.058 44.320 47.550 50.958 51.082 48.319 51.284 53.965 57.913 61.290 62.310 Matrices = Group 1; Covariances A+E+C | A+C _ A+C | A+E+C / Options rs End G3: Female DZ twin pairs Data NInput_vars=12 NObservations=51 Labels t1w1 t1w2 t1w3 t1w4 t1w5 t1w6 t2w1 t2w2 t2w3 t2w4 t2w5 t2w6 CMatrix 28.873 29.423 31.546 31.087 33.341 38.996 30.395 32.557 37.267 38.366 30.987 32.649 37.896 38.620 41.418 29.880 31.727 36.752 37.826 40.317 42.379 16.157 15.524 16.402 16.393 16.040 16.231 20.032 16.319 15.992 16.393 16.431 15.513 15.809 20.140 21.577 15.946 15.356 16.315 15.873 15.253 15.175 21.600 22.822 26.738 16.367 15.565 16.465 17.038 15.913 14.977 21.730 23.331 26.453 28.948 15.763 14.633 16.289 15.970 16.144 14.905 21.520 22.878 26.361 27.891 30.103 13.692 12.443 13.782 13.592 13.583 12.886 20.940 22.117 25.784 27.361 29.373 31.204 Matrices = Group 1; Covariances A+E+C | H@A+C _ H@A+C | A+E+C / Options rs Options Multiple Option issat Start 1 All Start 4 X 1 1 1 End !Equate measurement errors EQ U 1 1 1 U 1 2 2 U 1 3 3 U 1 4 4 U 1 5 5 U 1 6 6 End !Drop measurement error term Drop U 1 1 1 End !Drop Common environmental structure Drop Y 1 1 1 Y 1 2 2 Y 1 3 3 Y 1 4 4 Y 1 5 5 Y 1 6 6 Drop 106 to 110 End