___ ____ ____ ____ ____ 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.do . . *log using clogg.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 . . mclgen column (60 observations created) Your response factor is column with 4 categories. Its main effects form the intercept of a multinomial logistic model, interactions with independent variables form their effects. . . xi: mclest i.column [fweight=freq], rc2(row) i.column _Icolumn_1-4 (naturally coded; _Icolumn_1 omitted) Estimating rc2 effects for row ---------------------------------------------------------------------- iteration log likelihood sub-changes main changes ---------------------------------------------------------------------- 1.1 -1660.0652 -1660.0652 -1660.0652 1.2 -1658.6494 1.4158 -1658.6494 2.1 -1658.5696 0.0798 1.4956 2.2 -1658.5655 0.0041 0.0839 3.1 -1658.5653 0.0002 0.0043 3.2 -1658.5653 0.0000 0.0002 4.1 -1658.5653 0.0000 0.0000 4.2 -1658.5653 -0.0000 0.0000 ---------------------------------------------------------------------- Convergence criterion .0001 reached in 4 iterations Phi scale for column: ---------------------- column | phi ----------+----------- 1 | 0.0000 2 | 0.4597 3 | 0.7562 4 | 1.0000 ---------------------- Sigma scale for row: ---------------------- row | sigma ----------+----------- 1 | 0.0000 2 | -1.0407 3 | 0.2577 4 | 0.3419 5 | 1.0000 ---------------------- RC2 association parameters: ------------------------------------------------------------------------------- | Coef. Std. Err. z P>|z| -------------------------------------+----------------------------------------- mu: RC2 assoc column and row | -1.1823 0.0850 -13.9097 0.0000 ------------------------------------------------------------------------------- Full parameter listing: Conditional (fixed-effects) logistic regression Number of obs = 5508 LR chi2(4) = 500.72 Prob > chi2 = 0.0000 Log likelihood = -1658.5653 Pseudo R2 = 0.1312 ------------------------------------------------------------------------------ __didep | Coef. Std. Err. z P>|z| [95% Conf. Interval] -------------+---------------------------------------------------------------- _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 ------------------------------------------------------------------------------ Normalized Solution: Normalized Phi scale for column: ---------------------- column | phi ----------+----------- 1 | -0.7432 2 | -0.1265 3 | 0.2713 4 | 0.5984 ---------------------- Normalized Sigma scale for row: ---------------------- row | sigma ----------+----------- 1 | 0.0753 2 | 0.7763 3 | -0.0983 4 | -0.1550 5 | -0.5983 ---------------------- RC2 association parameters: ------------------------------------------------------------------------------- | Coef. Std. Err. z P>|z| -------------------------------------+----------------------------------------- mu: RC2 assoc column and row | 1.3083 0.0941 13.9097 0.0000 ------------------------------------------------------------------------------- Full parameter listing: Conditional (fixed-effects) logistic regression Number of obs = 5508 LR chi2(4) = 500.72 Prob > chi2 = 0.0000 Log likelihood = -1658.5653 Pseudo R2 = 0.1312 ------------------------------------------------------------------------------ __didep | Coef. Std. Err. z P>|z| [95% Conf. Interval] -------------+---------------------------------------------------------------- _Icolumn_2 | -1.039248 .1040429 -9.99 0.000 -1.243169 -.8353281 _Icolumn_3 | -.159227 .0798745 -1.99 0.046 -.3157783 -.0026758 _Icolumn_4 | .2274536 .0721542 3.15 0.002 .0860341 .3688731 __mu | 1.308277 .094055 13.91 0.000 1.123932 1.492621 ------------------------------------------------------------------------------ ------------------------------------------------------------------------------------------ Model fit information ------------------------------------------------------------------------------------------ log likelihood -1658.5653 model df 9 ------------------------------------------------------------------------------------------ . . end of do-file