[R-sig-ME] Poisson GLMM: where is the deviance?

Luciano La Sala lucianolasala at yahoo.com.ar
Sat Mar 17 20:25:50 CET 2012


Dear R experts:

I need to fit the following GLMM with Poisson distribution: glmer(yobs~year
+ (year|stake), data=mibase,family=poisson)


Where "stake" (plots where the number of birds is measured) is as random
effect, and "yobs" is the abundance of birds. I need an estimation of slope
an intercept for the random term.

I specified my model as follows:

glmer(yobs~year + (year|stake), data=mibase,family=poisson)

and I get the following (it does not calculate the deviance)

  Generalized linear mixed model fit by the Laplace approximation

  Formula: yobs ~ year + (year | stake)
     Data: mibase
   AIC BIC logLik deviance
   NaN NaN    NaN      NaN

  Random effects:
   Groups Name        Variance   Std.Dev.   Corr
   stake  (Intercept) 1.4165e-01 0.37636766
          year        3.5469e-08 0.00018833 0.000
  Number of obs: 1186, groups: stake, 63
  Fixed effects:
                Estimate Std. Error z value Pr(>|z|)
  (Intercept)  3.213e+01  6.445e+08       0        1
  year        -1.952e-01  3.225e+05       0        1

- Ignored:

  Correlation of Fixed Effects:
       (Intr)
  year -1.000
  Mensajes de aviso perdidos
  In mer_finalize(ans) : singular convergence (7)


Q. What does this error mean and how can I solve it?

By the way, if I run the same model but only estimating the intercept of
random effects as
glmer(yobs~year + (1|stake), data=mibase,family=poisson) it seems to work
properly.


I look forward for your response and thanks in advance.

Luciana




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