[R-sig-ME] NaN output from mcse on a glmm model

Ben Bolker bbo|ker @end|ng |rom gm@||@com
Fri Feb 1 18:26:05 CET 2019

  I can confirm this on my system. I looked and nothing obviously
fishy pops out about the model or the data. Inside the mcvcov()
function, after calling the C-code guts of the computation, stuff[[4]]
(the "numsum" component of the list passed to the C code) is full of
Inf and NaN values.  If I were you I'd (1) try a couple of very simple
examples, one with bernoulli and one with poisson data, to support
your idea that there's something wrong with the Poisson case; (2) look
around for examples of people using the package, e.g.
and see if there are Poisson examples; (3) contact the maintainer ...

On Fri, Feb 1, 2019 at 8:49 AM Adriaan De Jong <Adriaan.de.Jong using slu.se> wrote:
> Dear list members,
> I ran the following glmm on the attached data file, guided by Christina Knudson's "An introduction to Model-Fitting with the R package glmm, 11th of December 2018. (Fresh R download and all packages recently updated. I'm aware of the fact that the cluster part of the script I redundant)
> Everything appears to work fine, but when I try to extract the Monte Carlo standard errors with mcse I only receive NaN's for each of the parameters in the model.
> The example in Christina Knudson's text uses Bernoulli data, while I use count data (Poisson). Is that the cause of the NaN's?
> Grateful for any explanation or suggestion. Thanks in advance.
> Have a nice weekend,
> Adjan
> Adriaan "Adjan" de Jong
> Senior researcher
> Dept. of Wildlife, Fish, and Environmental Studies
> Swedish University of Agricultural Sciences
> set.seed(2019)
> clust<-makeCluster(2)
> mm01<-glmm(Arter ~ Status, random=list(~0+Site), varcomps.names=c("Site"), data=ArterVis, family.glmm=poisson.glmm, m=10000, cluster=clust)
> stopCluster(clust)
> summary(mm01)
> coef(mm01)
> confint(mm01)
> mcse(mm01)
> se(mm01)
> ---
> N?r du skickar e-post till SLU s? inneb?r detta att SLU behandlar dina personuppgifter. F?r att l?sa mer om hur detta g?r till, klicka h?r <https://www.slu.se/om-slu/kontakta-slu/personuppgifter/>
> E-mailing SLU will result in SLU processing your personal data. For more information on how this is done, click here <https://www.slu.se/en/about-slu/contact-slu/personal-data/>
> _______________________________________________
> R-sig-mixed-models using r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models

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