Skip to content

Monte Carlo simulations for lmer with binomial link

3 messages · Marco van de Ven, Ben Bolker, Jarrod Hadfield

#
Marco van de Ven wrote:
Not extremely easily.  See if anything at
<http://glmm.wikidot.com/faq> helps.
#
Hi Marco,

The newest version of MCMCglmm calculates p-values when you call  
summary. These should be reasonably close to what should be obtained  
under mcmcsamp/pvals.func unless there is a lot of over-dispersion. In  
this case the p-values will be larger  - as they should be.

For testing groups of effects as in anova, I recently recommended to  
someone that they calculate the posterior covariance matrix and pass  
it to Wald.test in the aod package. I am not advocating this, but I  
would be interested in other people's thoughts. If the posterior  
distribution of the effects is close to multivariate normal (its  
multivariate-t for Gaussian models) I can't see a problem, but perhaps  
there are some issues....

Cheers,

Jarrod



Cheers,

Jarrod


Quoting Marco van de Ven <m.a.m.vande.ven at gmail.com>: