___ ____ ____ ____ ____ tm /__ / ____/ / ____/ ___/ / /___/ / /___/ 8.2 Copyright 1984-2004 Statistics/Data Analysis StataCorp 4905 Lakeway Drive Special Edition College Station, Texas 77845 USA 800-STATA-PC http://www.stata.com 979-696-4600 stata@stata.com 979-696-4601 (fax) 46-user Stata for Windows (network) perpetual license: Serial number: 8198044710 Licensed to: SRC - ISR University of Michigan Notes: 1. (/m# option or -set memory-) 97.66 MB allocated to data 2. (/v# option or -set maxvar-) 5000 maximum variables 3. Stata running in batch mode running C:\Program Files\stata8SE\profile.do ... . do W:\yuxie\web\files\assoc\clogg-rc2.do . *log using clogg-rc2.log, replace . . input f1-f4 f1 f2 f3 f4 1. 44 11 38 62 2. 59 41 147 293 3. 23 11 13 27 4. 27 8 16 27 5. 258 57 105 110 6. end . . input r1-r4 r1 r2 r3 r4 1. 1 1 1 1 2. 2 2 2 2 3. 3 3 3 3 4. 4 4 4 4 5. 5 5 5 5 . . input c1-c4 c1 c2 c3 c4 1. 1 2 3 4 2. 1 2 3 4 3. 1 2 3 4 4. 1 2 3 4 5. 1 2 3 4 . . stack f1 r1 c1 f2 r2 c2 f3 r3 c3 f4 r4 c4, into(freq row column) clear . . tab row column, sum(freq) nost nofreq Means of freq | column row | 1 2 3 4 | Total -----------+--------------------------------------------+---------- 1 | 44 11 38 62 | 38.75 2 | 59 41 147 293 | 135 3 | 23 11 13 27 | 18.5 4 | 27 8 16 27 | 19.5 5 | 258 57 105 110 | 132.5 -----------+--------------------------------------------+---------- Total | 82.2 25.6 63.8 103.8 | 68.85 . . xi:glm freq i.row i.column , family(poisson) i.row _Irow_1-5 (naturally coded; _Irow_1 omitted) i.column _Icolumn_1-4 (naturally coded; _Icolumn_1 omitted) Iteration 0: log likelihood = -191.97426 Iteration 1: log likelihood = -173.76424 Iteration 2: log likelihood = -173.70926 Iteration 3: log likelihood = -173.70926 Generalized linear models No. of obs = 20 Optimization : ML: Newton-Raphson Residual df = 12 Scale parameter = 1 Deviance = 236.3378892 (1/df) Deviance = 19.69482 Pearson = 222.6175106 (1/df) Pearson = 18.55146 Variance function: V(u) = u [Poisson] Link function : g(u) = ln(u) [Log] Standard errors : OIM Log likelihood = -173.7092614 AIC = 18.17093 BIC = 200.3891019 ------------------------------------------------------------------------------ freq | Coef. Std. Err. z P>|z| [95% Conf. Interval] -------------+---------------------------------------------------------------- _Irow_2 | 1.248144 .0911233 13.70 0.000 1.069546 1.426743 _Irow_3 | -.73936 .141298 -5.23 0.000 -1.016299 -.462421 _Irow_4 | -.6867163 .1388241 -4.95 0.000 -.9588065 -.4146261 _Irow_5 | 1.229452 .0913149 13.46 0.000 1.050478 1.408426 _Icolumn_2 | -1.166563 .1012205 -11.52 0.000 -1.364951 -.9681744 _Icolumn_3 | -.2534021 .0746183 -3.40 0.001 -.3996512 -.107153 _Icolumn_4 | .2333107 .0660293 3.53 0.000 .1038956 .3627258 _cons | 3.834356 .0903243 42.45 0.000 3.657323 4.011388 ------------------------------------------------------------------------------ . . rc2, row(row) col(column) Estimating rc2 effects between row and column ---------------------------------------------------------------------- iteration log likelihood sub-changes main changes ---------------------------------------------------------------------- 1.1 -59.8190 113.8903 -59.8190 1.2 -58.4032 1.4158 -58.4032 2.1 -58.3234 0.0798 1.4956 2.2 -58.3193 0.0041 0.0839 3.1 -58.3190 0.0002 0.0043 3.2 -58.3190 0.0000 0.0002 4.1 -58.3190 0.0000 0.0000 4.2 -58.3190 0.0000 0.0000 ---------------------------------------------------------------------- Convergence criterion .0001 reached in 4 iterations Sigma scale for row: ---------------------- row | sigma ----------+----------- 1 | 0.0000 2 | -1.0407 3 | 0.2577 4 | 0.3419 5 | 1.0000 ---------------------- Phi scale for column: ---------------------- column | phi ----------+----------- 1 | 0.0000 2 | 0.4597 3 | 0.7562 4 | 1.0000 ---------------------- RC2 association parameters: ------------------------------------------------------------------------------- | Coef. Std. Err. z P>|z| -------------------------------------+----------------------------------------- mu: RC2 assoc column & row | -1.1823 0.0850 -13.9097 0.0000 ------------------------------------------------------------------------------- Full parameter listing: Poisson regression Number of obs = 20 LR chi2(8) = 1352.39 Prob > chi2 = 0.0000 Log likelihood = -58.319034 Pseudo R2 = 0.9206 ------------------------------------------------------------------------------ freq | Coef. Std. Err. z P>|z| [95% Conf. Interval] -------------+---------------------------------------------------------------- _Irow_2 | .3779616 .1167082 3.24 0.001 .1492177 .6067054 _Irow_3 | -.5588933 .1419795 -3.94 0.000 -.837168 -.2806187 _Irow_4 | -.4507749 .1399833 -3.22 0.001 -.725137 -.1764127 _Irow_5 | 1.837073 .1004113 18.30 0.000 1.64027 2.033876 _Icolumn_2 | -.9784989 .1050099 -9.32 0.000 -1.184315 -.7726832 _Icolumn_3 | -.059287 .0819176 -0.72 0.469 -.2198425 .1012686 _Icolumn_4 | .3596105 .0738297 4.87 0.000 .2149068 .5043141 __mu | -1.182253 .0849947 -13.91 0.000 -1.34884 -1.015667 _cons | 3.721381 .0941497 39.53 0.000 3.536851 3.905911 ------------------------------------------------------------------------------ Normalized Sigma scale for row: ---------------------- row | sigma ----------+----------- 1 | 0.0753 2 | 0.7763 3 | -0.0983 4 | -0.1550 5 | -0.5983 ---------------------- Normalized Phi scale for column: ---------------------- column | phi ----------+----------- 1 | -0.7432 2 | -0.1265 3 | 0.2713 4 | 0.5984 ---------------------- RC2 association parameters: ------------------------------------------------------------------------------- | Coef. Std. Err. z P>|z| -------------------------------------+----------------------------------------- mu: RC2 assoc column & row | 1.3083 0.0941 13.9097 0.0000 ------------------------------------------------------------------------------- Full parameter listing: Poisson regression Number of obs = 20 LR chi2(8) = 1352.39 Prob > chi2 = 0.0000 Log likelihood = -58.319034 Pseudo R2 = 0.9206 ------------------------------------------------------------------------------ freq | Coef. Std. Err. z P>|z| [95% Conf. Interval] -------------+---------------------------------------------------------------- _Irow_2 | 1.059565 .0947077 11.19 0.000 .8739414 1.245189 _Irow_3 | -.7276807 .1413508 -5.15 0.000 -1.004723 -.4506381 _Irow_4 | -.6747102 .1389152 -4.86 0.000 -.9469791 -.4024414 _Irow_5 | 1.182134 .092873 12.73 0.000 1.000107 1.364162 _Icolumn_2 | -1.039248 .1040451 -9.99 0.000 -1.243173 -.8353236 _Icolumn_3 | -.159227 .0798743 -1.99 0.046 -.3157778 -.0026762 _Icolumn_4 | .2274536 .072154 3.15 0.002 .0860344 .3688728 __mu | 1.308277 .0940548 13.91 0.000 1.123933 1.492621 _cons | 3.794592 .0934024 40.63 0.000 3.611527 3.977658 ------------------------------------------------------------------------------ ------------------------------------------------------------------------------------------ Model summary ------------------------------------------------------------------------------------------ Baseline model log likelihood -173.7093 Baseline model model df 7 rc2 model log likelihood -58.3190 rc2 model df 13 Chi-square improvement 230.7805 df change 6 significance 0.0000 rc2 model deviance 5.5570 rc2 model residual df 6 rc2 model prob 0.4746 rc2 model bic -37.8090 rc2 model aic -6.4430 Number of cases 1377 ------------------------------------------------------------------------------------------ . . . . . end of do-file