[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,
Adriaan "Adjan" de Jong
Dept. of Wildlife, Fish, and Environmental Studies
Swedish University of Agricultural Sciences
mm01<-glmm(Arter ~ Status, random=list(~0+Site), varcomps.names=c("Site"), data=ArterVis, family.glmm=poisson.glmm, m=10000, cluster=clust)
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
More information about the R-sig-mixed-models