LEM: log-linear and event history analysis with missing data. Developed by Jeroen K. Vermunt (c), Tilburg University, The Netherlands. Version 1.2 (July 10, 1998). *** INPUT *** * Wickens (1992) MLE of a multivariate guassian rating * model with excluded data. J Math Psych, 36, 213-234. * Subject A (Wickens & Olzak 1989) * * Log-linear models using 4 variables: (HLX,,HY,XY,) * man 4 dim 2 2 6 6 lab L H X Y mod {L H X Y HL HY LX XY HLX} nco dat[ 44 4 9 7 6 7 13 30 20 8 14 7 9 23 17 17 3 0 16 17 10 20 2 2 5 4 9 10 4 0 3 3 0 1 4 1 7 4 5 5 14 69 5 7 13 15 38 37 6 7 8 10 10 15 4 12 5 13 6 14 2 3 1 1 3 5 0 0 1 1 1 3 8 2 2 1 0 4 5 5 5 5 5 3 8 10 7 4 1 1 12 17 15 13 2 2 12 17 19 18 10 4 31 29 25 24 12 12 4 1 2 0 4 37 0 4 0 1 8 25 1 3 3 7 8 15 4 4 8 17 12 21 3 12 8 11 20 20 11 8 12 11 12 33 ] *** STATISTICS *** Number of iterations = 8 Converge criterion = 0.0000007701 X-squared = 98.8974 (0.1438) L-squared = 108.1864 (0.0457) Cressie-Read = 99.3489 (0.1369) Dissimilarity index = 0.0894 Degrees of freedom = 85 Log-likelihood = -6395.07676 Number of parameters = 58 (+1) Sample size = 1399.0 BIC(L-squared) = -507.5122 AIC(L-squared) = -61.8136 BIC(log-likelihood) = 13210.2773 AIC(log-likelihood) = 12906.1535 Eigenvalues information matrix 3494.6794 3047.0249 2465.4976 1487.9665 1423.4166 1181.6992 1053.1398 691.7399 556.0756 407.4707 392.3357 336.2621 327.1280 308.9386 302.4187 297.9285 288.8009 280.7291 273.5078 267.5752 248.3926 239.5201 227.6456 222.8331 216.3940 204.0767 186.0671 175.6963 171.0997 167.0254 159.8356 144.0320 142.9599 137.2334 134.3028 130.8009 128.7430 123.9319 110.6050 86.6838 75.3012 68.9159 49.6844 43.3039 41.3859 39.9368 37.8094 36.9394 35.7518 33.4848 30.6370 28.6774 24.7808 23.3667 21.7152 20.6495 16.5095 13.2570 *** FREQUENCIES *** L H X Y observed estimated std. res. 1 1 1 1 44.000 39.031 0.795 1 1 1 2 4.000 6.167 -0.873 1 1 1 3 9.000 9.488 -0.158 1 1 1 4 7.000 5.854 0.473 1 1 1 5 6.000 5.567 0.184 1 1 1 6 7.000 10.893 -1.179 1 1 2 1 13.000 13.597 -0.162 1 1 2 2 30.000 24.751 1.055 1 1 2 3 20.000 19.287 0.162 1 1 2 4 8.000 12.671 -1.312 1 1 2 5 14.000 14.952 -0.246 1 1 2 6 7.000 6.743 0.099 1 1 3 1 9.000 12.466 -0.982 1 1 3 2 23.000 20.183 0.627 1 1 3 3 17.000 15.433 0.399 1 1 3 4 17.000 14.289 0.717 1 1 3 5 3.000 4.239 -0.602 1 1 3 6 0.000 2.390 -1.546 1 1 4 1 16.000 14.028 0.527 1 1 4 2 17.000 17.549 -0.131 1 1 4 3 10.000 12.507 -0.709 1 1 4 4 20.000 17.609 0.570 1 1 4 5 2.000 3.114 -0.631 1 1 4 6 2.000 2.193 -0.130 1 1 5 1 5.000 5.053 -0.023 1 1 5 2 4.000 7.624 -1.312 1 1 5 3 9.000 7.447 0.569 1 1 5 4 10.000 7.054 1.109 1 1 5 5 4.000 3.624 0.198 1 1 5 6 0.000 1.200 -1.095 1 1 6 1 3.000 3.150 -0.085 1 1 6 2 3.000 2.626 0.231 1 1 6 3 0.000 2.395 -1.547 1 1 6 4 1.000 2.093 -0.755 1 1 6 5 4.000 0.987 3.033 1 1 6 6 1.000 0.750 0.289 1 2 1 1 7.000 10.503 -1.081 1 2 1 2 4.000 2.375 1.054 1 2 1 3 5.000 4.390 0.291 1 2 1 4 5.000 4.005 0.497 1 2 1 5 14.000 11.771 0.650 1 2 1 6 69.000 70.956 -0.232 1 2 2 1 5.000 3.958 0.524 1 2 2 2 7.000 10.310 -1.031 1 2 2 3 13.000 9.653 1.077 1 2 2 4 15.000 9.375 1.837 1 2 2 5 38.000 34.197 0.650 1 2 2 6 37.000 47.507 -1.524 1 2 3 1 6.000 3.573 1.284 1 2 3 2 7.000 8.279 -0.445 1 2 3 3 8.000 7.607 0.143 1 2 3 4 10.000 10.412 -0.128 1 2 3 5 10.000 9.548 0.146 1 2 3 6 15.000 16.581 -0.388 1 2 4 1 4.000 4.140 -0.069 1 2 4 2 12.000 7.413 1.685 1 2 4 3 5.000 6.348 -0.535 1 2 4 4 13.000 13.212 -0.058 1 2 4 5 6.000 7.223 -0.455 1 2 4 6 14.000 15.665 -0.421 1 2 5 1 2.000 0.727 1.492 1 2 5 2 3.000 1.570 1.141 1 2 5 3 1.000 1.843 -0.621 1 2 5 4 1.000 2.581 -0.984 1 2 5 5 3.000 4.098 -0.543 1 2 5 6 5.000 4.180 0.401 1 2 6 1 0.000 0.447 -0.669 1 2 6 2 0.000 0.534 -0.731 1 2 6 3 1.000 0.585 0.543 1 2 6 4 1.000 0.756 0.281 1 2 6 5 1.000 1.102 -0.097 1 2 6 6 3.000 2.577 0.264 2 1 1 1 8.000 8.617 -0.210 2 1 1 2 2.000 1.362 0.547 2 1 1 3 2.000 2.095 -0.065 2 1 1 4 1.000 1.293 -0.257 2 1 1 5 0.000 1.229 -1.109 2 1 1 6 4.000 2.405 1.029 2 1 2 1 5.000 4.138 0.424 2 1 2 2 5.000 7.533 -0.923 2 1 2 3 5.000 5.870 -0.359 2 1 2 4 5.000 3.856 0.582 2 1 2 5 5.000 4.551 0.211 2 1 2 6 3.000 2.052 0.662 2 1 3 1 8.000 5.601 1.014 2 1 3 2 10.000 9.068 0.310 2 1 3 3 7.000 6.934 0.025 2 1 3 4 4.000 6.420 -0.955 2 1 3 5 1.000 1.905 -0.655 2 1 3 6 1.000 1.074 -0.071 2 1 4 1 12.000 12.772 -0.216 2 1 4 2 17.000 15.978 0.256 2 1 4 3 15.000 11.387 1.071 2 1 4 4 13.000 16.032 -0.757 2 1 4 5 2.000 2.835 -0.496 2 1 4 6 2.000 1.996 0.003 2 1 5 1 12.000 12.632 -0.178 2 1 5 2 17.000 19.059 -0.472 2 1 5 3 19.000 18.617 0.089 2 1 5 4 18.000 17.634 0.087 2 1 5 5 10.000 9.059 0.313 2 1 5 6 4.000 2.999 0.578 2 1 6 1 31.000 34.915 -0.663 2 1 6 2 29.000 29.102 -0.019 2 1 6 3 25.000 26.541 -0.299 2 1 6 4 24.000 23.196 0.167 2 1 6 5 12.000 10.939 0.321 2 1 6 6 12.000 8.307 1.281 2 2 1 1 4.000 4.847 -0.385 2 2 1 2 1.000 1.096 -0.092 2 2 1 3 2.000 2.026 -0.019 2 2 1 4 0.000 1.848 -1.360 2 2 1 5 4.000 5.433 -0.615 2 2 1 6 37.000 32.749 0.743 2 2 2 1 0.000 1.308 -1.144 2 2 2 2 4.000 3.407 0.321 2 2 2 3 0.000 3.190 -1.786 2 2 2 4 1.000 3.098 -1.192 2 2 2 5 8.000 11.300 -0.982 2 2 2 6 25.000 15.698 2.348 2 2 3 1 1.000 2.361 -0.886 2 2 3 2 3.000 5.470 -1.056 2 2 3 3 3.000 5.026 -0.904 2 2 3 4 7.000 6.879 0.046 2 2 3 5 8.000 6.308 0.673 2 2 3 6 15.000 10.955 1.222 2 2 4 1 4.000 5.060 -0.471 2 2 4 2 4.000 9.060 -1.681 2 2 4 3 8.000 7.758 0.087 2 2 4 4 17.000 16.147 0.212 2 2 4 5 12.000 8.828 1.068 2 2 4 6 21.000 19.146 0.424 2 2 5 1 3.000 3.588 -0.310 2 2 5 2 12.000 7.747 1.528 2 2 5 3 8.000 9.093 -0.362 2 2 5 4 11.000 12.732 -0.485 2 2 5 5 20.000 20.219 -0.049 2 2 5 6 20.000 20.621 -0.137 2 2 6 1 11.000 6.488 1.772 2 2 6 2 8.000 7.739 0.094 2 2 6 3 12.000 8.480 1.209 2 2 6 4 11.000 10.956 0.013 2 2 6 5 12.000 15.972 -0.994 2 2 6 6 33.000 37.366 -0.714 *** LOG-LINEAR PARAMETERS *** * TABLE LHXY [or P(LHXY)] * effect beta std err z-value exp(beta) Wald df prob main 1.8190 6.1654 L 1 -0.0793 0.0357 -2.223 0.9238 2 0.0793 1.0825 4.94 1 0.026 H 1 0.0501 0.0390 1.283 1.0514 2 -0.0501 0.9511 1.65 1 0.199 X 1 -0.1885 0.0904 -2.084 0.8282 2 0.2350 0.0723 3.249 1.2649 3 0.0493 0.0738 0.669 1.0506 4 0.3247 0.0665 4.881 1.3836 5 -0.0563 0.0821 -0.685 0.9453 6 -0.3643 0.6947 38.65 5 0.000 Y 1 -0.1143 0.0793 -1.442 0.8920 2 0.0131 0.0752 0.174 1.0132 3 0.0147 0.0714 0.206 1.0148 4 0.0723 0.0708 1.022 1.0750 5 -0.0717 0.0745 -0.963 0.9308 6 0.0859 1.0896 4.35 5 0.500 LH 1 1 0.1020 0.0357 2.859 1.1074 1 2 -0.1020 0.9030 2 1 -0.1020 0.9030 2 2 0.1020 1.1074 8.17 1 0.004 LX 1 1 0.6503 0.0744 8.741 1.9160 1 2 0.6535 0.0684 9.560 1.9223 1 3 0.3830 0.0714 5.367 1.4666 1 4 0.0526 0.0631 0.834 1.0540 1 5 -0.5488 0.0802 -6.840 0.5777 1 6 -1.1906 0.3040 2 1 -0.6503 0.5219 2 2 -0.6535 0.5202 2 3 -0.3830 0.6818 2 4 -0.0526 0.9488 2 5 0.5488 1.7311 2 6 1.1906 3.2890 270.01 5 0.000 HX 1 1 -0.1682 0.0887 -1.896 0.8452 1 2 -0.0437 0.0764 -0.572 0.9573 1 3 -0.1119 0.0791 -1.414 0.8942 1 4 -0.1037 0.0705 -1.470 0.9015 1 5 0.1590 0.0872 1.823 1.1724 1 6 0.2685 1.3080 2 1 0.1682 1.1832 2 2 0.0437 1.0446 2 3 0.1119 1.1184 2 4 0.1037 1.1093 2 5 -0.1590 0.8530 2 6 -0.2685 0.7646 12.86 5 0.025 HY 1 1 0.5901 0.0763 7.737 1.8042 1 2 0.4109 0.0689 5.967 1.5082 1 3 0.3191 0.0693 4.606 1.3759 1 4 0.1237 0.0656 1.885 1.1316 1 5 -0.4406 0.0710 -6.203 0.6436 1 6 -1.0032 0.3667 2 1 -0.5901 0.5543 2 2 -0.4109 0.6630 2 3 -0.3191 0.7268 2 4 -0.1237 0.8837 2 5 0.4406 1.5537 2 6 1.0032 2.7270 275.10 5 0.000 XY 1 1 0.9209 0.1360 6.770 2.5115 1 2 -0.8724 0.2264 -3.853 0.4179 1 3 -0.3514 0.1869 -1.881 0.7037 1 4 -0.6964 0.2105 -3.309 0.4984 1 5 -0.0384 0.1727 -0.223 0.9623 1 6 1.0378 2.8230 2 1 -0.5211 0.1667 -3.126 0.5938 2 2 0.1297 0.1336 0.971 1.1385 2 3 -0.0295 0.1390 -0.212 0.9709 2 4 -0.3118 0.1514 -2.059 0.7321 2 5 0.5621 0.1238 4.542 1.7543 2 6 0.1707 1.1861 3 1 -0.1594 0.1675 -0.952 0.8526 3 2 0.3743 0.1395 2.683 1.4539 3 3 0.1962 0.1456 1.347 1.2167 3 4 0.2570 0.1418 1.812 1.2930 3 5 -0.2499 0.1724 -1.450 0.7789 3 6 -0.4181 0.6583 4 1 0.0282 0.1449 0.195 1.0286 4 2 0.3040 0.1313 2.316 1.3553 4 3 0.0555 0.1395 0.398 1.0571 4 4 0.5354 0.1217 4.400 1.7082 4 5 -0.4886 0.1700 -2.875 0.6135 4 6 -0.4346 0.6475 5 1 -0.3696 0.1717 -2.152 0.6910 5 2 0.0936 0.1458 0.642 1.0981 5 3 0.1603 0.1426 1.124 1.1739 5 4 0.2439 0.1391 1.753 1.2762 5 5 0.2861 0.1458 1.963 1.3313 5 6 -0.4143 0.6608 6 1 0.1011 1.1064 6 2 -0.0292 0.9712 6 3 -0.0311 0.9694 6 4 -0.0280 0.9724 6 5 -0.0713 0.9312 6 6 0.0585 1.0602 178.48 25 0.000 LHX 1 1 1 0.0823 0.0744 1.107 1.0858 1 1 2 -0.0815 0.0684 -1.191 0.9218 1 1 3 -0.0056 0.0714 -0.078 0.9944 1 1 4 -0.0284 0.0631 -0.450 0.9720 1 1 5 0.0679 0.0802 0.847 1.0703 1 1 6 -0.0348 0.9658 1 2 1 -0.0823 0.9210 1 2 2 0.0815 1.0849 1 2 3 0.0056 1.0056 1 2 4 0.0284 1.0288 1 2 5 -0.0679 0.9343 1 2 6 0.0348 1.0354 2 1 1 -0.0823 0.9210 2 1 2 0.0815 1.0849 2 1 3 0.0056 1.0056 2 1 4 0.0284 1.0288 2 1 5 -0.0679 0.9343 2 1 6 0.0348 1.0354 2 2 1 0.0823 1.0858 2 2 2 -0.0815 0.9218 2 2 3 -0.0056 0.9944 2 2 4 -0.0284 0.9720 2 2 5 0.0679 1.0703 2 2 6 -0.0348 0.9658 3.21 5 0.667