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]]

