[R-sig-ME] Is there an option to not include the correlation matrix in output of the glmmadmb function?

Ben Bolker bbolker at gmail.com
Fri Jul 6 16:46:18 CEST 2012


  I'm a little confused.  Can you specify which version of glmmADMB
you're using (i.e. results of sessionInfo()) ?  Can you give a small
reproducible example?

  In the current version it doesn't seem that the correlation matrix of
the fixed-effect parameters is printed at all?

For example, on my system:

  data(bacteria,package="MASS")
     bacteria$present <- as.numeric(bacteria$y)-1
     (bfit <-  glmmadmb(present ~ trt + I(week > 2), random = ~ 1 | ID,
                          family = "binomial", data = bacteria))
summary(bfit)

Call:
glmmadmb(formula = present ~ trt + I(week > 2), data = bacteria,
    family = "binomial", random = ~1 | ID)


Coefficients:
                Estimate Std. Error z value Pr(>|z|)
(Intercept)        3.548      0.696    5.10  3.5e-07 ***
trtdrug           -1.367      0.677   -2.02  0.04355 *
trtdrug+          -0.783      0.683   -1.15  0.25197
I(week > 2)TRUE   -1.599      0.476   -3.36  0.00078 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Number of observations: total=220, ID=50
Random effect variance(s):
Group=ID
            Variance StdDev
(Intercept)    1.554  1.246

Log-likelihood: -96.1307



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