Notice: On March 31, it was **announced** that Statalist is moving from an email list to a **forum**. The old list will shut down on April 23, and its replacement, **statalist.org** is already up and running.

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

From |
maarten buis <maartenlbuis@googlemail.com> |

To |
statalist@hsphsun2.harvard.edu |

Subject |
Re: st: computing odds ratios for models with interraction terms |

Date |
Fri, 6 May 2011 09:39:34 +0200 |

On Thu, May 5, 2011 at 11:38 PM, Fitzmaurice, Ann E. wrote: > I use the following for example > > Logistic outcome i.varA i.varB b i.varA#i.varB other variables > > And I obtain the odds ratios for var a var b and the interaction terms > > > Var a has four levels and var b has 3 levels > > What I would like to do , is to generate the odds ratios for the 12 cells in the table var a by var b > > Is there a way of doing this in stata, I would also like the confidence intervals You can quite easily do so with the new factor variable notation, the trick is to use # instead of ##. However you won't be able to fill the entire 4 by 3 table with odds ratios as you need to define one of these cells as your reference category. An odds ratio is a comparison of groups, so it needs to have a group to compare it with. this means that one of your cells will be fixed at 1. This is what I have done in the first -logit- model. The reference category is divorced women with less than highschool. The odds ratio for that category is reported as (base). Alternatively you can replace it with 1 and leave the standard error empty, as the odds ratio for that group is the odds of attaining a high occupation within that group divided by the odds of attaining a high occupation within that group, which is trivially equal to 1. The coefficient of baseline gives you the baseline odds, that is the odds of attaining a high job for divorced women with less than high school education. Within that group we expect to find .12 women with a high job for every women who does not have a high job. All odds ratios tell you by what factor the group that belongs to that odds ratio differs from this reference category. For example the odds increases by a factor 1.78 [i.e. (1.78 - 1)*100% = 78%] when a divorced woman gets a high school diploma. You can get coefficients for all your cells, but than you will get odds not odds ratios. The trick is to precede the categorical variables with ibn. instead of i., meaning that you don't want Stata to leave out the reference category, and you must make sure there is no constant estimated, i.e. specify the -nocons- option and leave out the variable baseline. This is what I did in the second -logit- model. So for divorced women without high school we find that there are .12 women with a high occupation per woman without a high occupation, which is exactly the same as we found above. For divorced women with high school there are .22 women with a high occupation per woman without a high occupation. Notice that if we divide those two odds [ -di .2193878 /.1235955-] we get exactly the odds ratio from our first -logit- model. This should be true, as the two models are just two different ways of displaying the same model. *------------------------- begin example ----------------------- // data preparation sysuse nlsw88, clear gen byte marst = never_married + 2*married label define marst 1 "divorced/widowed" /// 2 "never married" /// 3 "married" label value marst marst label variable marst "marital status" gen byte edcat = cond(grade < 12, 1, /// cond(grade == 12, 2, 3)) /// if grade < . label define edcat 1 "less than high school" /// 2 "high school" /// 3 "more than high school" label value edcat edcat label variable edcat "education in categories" gen byte high_occ = occupation < 3 /// if occupation < . label define high_occ 1 "proffesionals and managers" /// 2 "other" label value high_occ high_occ gen byte baseline = 1 // get all odds ratios, // reference = divorced and less than highschool logit high_occ i.marst#i.edcat baseline, /// or nocons baselevels // get all odds logit high_occ ibn.marst#ibn.edcat, nocons or *------------------------ end example ---------------------------- (For more on examples I sent to the Statalist see: http://www.maartenbuis.nl/example_faq ) Hope this helps, Maarten -------------------------- Maarten L. Buis Institut fuer Soziologie Universitaet Tuebingen Wilhelmstrasse 36 72074 Tuebingen Germany http://www.maartenbuis.nl -------------------------- * * For searches and help try: * http://www.stata.com/help.cgi?search * http://www.stata.com/support/statalist/faq * http://www.ats.ucla.edu/stat/stata/

**Follow-Ups**:**RE: st: computing odds ratios for models with interraction terms***From:*"Fitzmaurice, Ann E." <a.e.fitzmaurice@abdn.ac.uk>

**References**:**st: computing odds ratios for models with interraction terms***From:*"Fitzmaurice, Ann E." <a.e.fitzmaurice@abdn.ac.uk>

- Prev by Date:
**st: dummy variables with option cluster()** - Next by Date:
**Re: st: how to generate mean for a secquence of five year of a variable** - Previous by thread:
**st: computing odds ratios for models with interraction terms** - Next by thread:
**RE: st: computing odds ratios for models with interraction terms** - Index(es):