Skip to content
Prev 14498 / 20628 Next

Help with MCMC fitting in R

Janelle Sylvester <jsylves92 at ...> writes:
snip
snip
[snip]

  It's possible that you're just having problems with an out-of-date
binary: some people reported difficulty like this, and solved it
using instructions at:

https://stat.ethz.ch/pipermail/r-sig-mixed-models/2016q1/024490.html

  Data and figures get stripped by the mailing list software, so we
haven't seen that (you could send it to me, but I can't guarantee
I'll have time to take a look at it).

  I'm not sure what you're up to with defining mcmc.control there ...

  As far as other solutions go: do you absolutely need the post-hoc
MCMC?  It is nice, but I would generally say that if you do standard
model diagnostics (examine residuals and model predictions, ideally
graphically) it's not an ironclad requirement. (Among other things,
most other non-Bayesian model-fitting methods don't offer this 
feature ...)  Other ways to go to cross-check your model would be:

- fit a zero-inflated Poisson-logNormal with MCMCglmm (a bit of a nuisance,
but doable: search for "owls NCEAS bolker" to find an example)
- use the relatively new/experimental glmmTMB package (install
via devtools::install_github("glmmTMB/glmmTMB",sub="glmmTMB") , 
then library("glmmTMB"); ?Owls for an example)

  Ben Bolker