[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