[R-sig-ME] convert into ORs

raphaëlle métras metras.raphaelle at choisyclub.org
Fri Nov 13 14:23:49 CET 2009


Dear all,

I run a GLMM binomial, using glmer function. Please see the command and
the ouputs at the end of my message.

I wanted to know if there is a command that would convert my estimates
into odds ratios, e.g. similar to the command logistic.display() that is
working for 'glm'.

I've been looking for a while but I really cant find it.

Thanks very much for your help,

Regards,

Raphaelle


> k8b <- glmer(status  ~ no_tradcat + dayscat + (1|fm_id) ,
family=binomial, data=basket)
> summary(k8b)
Generalized linear mixed model fit by the Laplace approximation
Formula: status ~ no_tradcat + dayscat + (1 | fm_id)
   Data: basket
 AIC BIC logLik deviance
 120 134  -54.8      110
Random effects:
 Groups Name        Variance Std.Dev.
 fm_id  (Intercept)  1        1
Number of obs: 117, groups: fm_id, 48

Fixed effects:
                  Estimate Std. Error z value Pr(>|z|)
(Intercept)         -4.501      1.166   -3.86  0.00011 ***
no_tradcat1-2trad    2.686      0.752    3.57  0.00036 ***
dayscat1-2days       2.896      1.072    2.70  0.00689 **
dayscatsame day      3.563      1.015    3.51  0.00045 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
            (Intr) n_t1-2 dys1-2
n_trdct1-2t -0.617
dysct1-2dys -0.878  0.339
dayscatsmdy -0.745  0.067  0.793




More information about the R-sig-mixed-models mailing list