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

Adriaan De Jong Adr|@@n@de@Jong @end|ng |rom @|u@@e
Fri Feb 1 14:49:05 CET 2019


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

-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: Arter 2002-2015 Visits Impact Short.txt
URL: <https://stat.ethz.ch/pipermail/r-sig-mixed-models/attachments/20190201/3c92833f/attachment-0001.txt>


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