[R] odds ratio: how to create reference
Eik Vettorazzi
E.Vettorazzi at uke.uni-hamburg.de
Tue Sep 23 11:36:30 CEST 2008
Hi,
You can state a probability p as odds p/(1-p) and vice versa. To get an
odds ratio you need actually two odds. Then you can get the odds ration
of being/having "a" instead of "b" by odds(a)/odds(b), where "b" is the
reference level.
If you fit a logistic regression model (which means that your outcome is
dichotomous) then the estimated coefficients are actually
log(oddsratios) - which you can transform to odds by exp() .
You can use a factor-variable with three levels for race and
treatment-contrasts to get odds ratios for not being white against being
white - make sure, that either your factor has "white" as first level or
specify the contrast with the "base" argument.
If you create 3 dummy variables and involve an intercept in your model
your model will be perfectly collinear - the so called "dummy variable
trap" - you can use an intercept and create two dummies for the
covariate levels you are actually interested in and put this in your
logistic model - the result will be the same as with the treatment
contrasts.
hth.
Bunny, lautloscrew.com schrieb:
> HI there,
>
> i know this is a basic question, though i need some help because this
> is somewhat away from my current issue, but nevertheless interesting
> to me... Lets assume i have some estimated probabilities, say
> estimated by a logit model. i know i can also state them as an odds
> ratio.
>
> Now i´d like to state these odds ratios as a reference to a specific
> outcome of my investigated variable.
>
> for example, if my covariate of interest is race and possible outcomes
> are white, black and hispanic, whereas the latter are minorities in my
> case - how can i state the odds ratio in such a way that white is the
> reference (always 1) and other races' odds ratio are relative to the
> reference. e.g. hispanics are 1.5 times more likely to ...
>
> Is creating 3 binary dummies for race the right way ? And if so how
> can i go on.
> As i said, i know this is rather basic, i am thankful for any links /
> references...
>
> thanks in advance !
>
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide
> http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
--
Eik Vettorazzi
Institut für Medizinische Biometrie und Epidemiologie
Universitätsklinikum Hamburg-Eppendorf
Martinistr. 52
20246 Hamburg
T ++49/40/42803-8243
F ++49/40/42803-7790
More information about the R-help
mailing list