Dear All
I have fitted the following glmm:
cmai ~ time.f * intrv.f + (1 | nhome.f/Res_Code.f)
with poisson distribution, using both glmer and glmmadmb.
But the estimation for the fixed and random effects were different, i.e.
> summary(lmer.AGGREG.cmai.out3)
Call:
glmmadmb(formula = cmai ~ time.f * intrv.f + (1 | nhome.f/Res_Code.f),
data = beam.AGGREG.cmai.long, family = "poisson", link = "log",
zeroInflation = F, admb.opts = admbControl(impSamp = 0, run = F),
save.dir = "tmp")
AIC: 1032.2
Coefficients:
Estimate Std. Error z value Pr(>|z|)
(Intercept) 0.542 3.105 0.17 0.86
time.f2 0.104 5.177 0.02 0.98
time.f3 -0.526 3.230 -0.16 0.87
intrv.f1 0.929 2.712 0.34 0.73
time.f2:intrv.f1 -0.416 5.302 -0.08 0.94
time.f3:intrv.f1 0.177 3.261 0.05 0.96
Number of observations: total=1032, nhome.f=35, nhome.f:Res_Code.f=344
Random effect variance(s):
Group=nhome.f
Variance StdDev
(Intercept) 0.7118 0.8437
Group=nhome.f:Res_Code.f
Variance StdDev
(Intercept) 1.454 1.206
Log-likelihood: -508.108
> summary(lmer.AGGREG.cmai.out2)
Generalized linear mixed model fit by the Laplace approximation
Formula: cmai ~ time.f * intrv.f + (1 | nhome.f/Res_Code.f)
Data: beam.AGGREG.cmai.long
AIC BIC logLik deviance
1835 1874 -909.5 1819
Random effects:
Groups Name Variance Std.Dev.
Res_Code.f:nhome.f (Intercept) 0.040125 0.20031
nhome.f (Intercept) 0.033702 0.18358
Number of obs: 1032, groups: Res_Code.f:nhome.f, 344; nhome.f, 35
Fixed effects:
Estimate Std. Error z value Pr(>|z|)
(Intercept) 3.62040 0.04749 76.23 <2e-16 ***
time.f2 -0.01964 0.01706 -1.15 0.2496
time.f3 0.01643 0.01691 0.97 0.3310
intrv.f1 0.07540 0.06819 1.11 0.2689
time.f2:intrv.f1 0.02148 0.02395 0.90 0.3698
time.f3:intrv.f1 -0.04835 0.02394 -2.02 0.0435 *
---
Signif. codes: 0 *** 0.001 ** 0.01 * 0.05 . 0.1 1
Correlation of Fixed Effects:
(Intr) tim.f2 tim.f3 intr.1 t.2:.1
time.f2 -0.178
time.f3 -0.179 0.500
intrv.f1 -0.696 0.124 0.125
tm.f2:ntr.1 0.127 -0.712 -0.356 -0.175
tm.f3:ntr.1 0.127 -0.353 -0.706 -0.175 0.498
Does anyone have similar experience, which can be explained for the reason
behind it?
Regards,
Kenny
[[alternative HTML version deleted]]
______________________________________________
R-help@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.
