[R-sig-ME] MCMCglmm, priorR and binomial distribution

Camille Madec camille.madec at ebc.uu.se
Thu Jan 31 15:21:12 CET 2013


Dear everyone,

I have a model with 2 fixed factors, 1 random factor and a binary  
response variable. I ran a MCMCglmm with family=”categorical” and the  
prior for the residual being R=list(V=1, nu=0.002). In the summary of  
the model I got high post.mean values (around 50 for fixed effects and  
 >1000 for random effects and sometimes up to 14000).
I ran the same model with R=list(V=1, fix=1) which means that the  
variance of the residual is fixed to 1, so the residual becomes a  
fixed factor (if I understand correctly). In that case my post.mean  
values are smaller (between zero and 24).

My questions are:
1) Are the large values in the first case normal?
2) How do I know which prior is the more appropriate for the residual?

Bests,
Camille

-----
Camille Madec
PhD student
Plant Ecology and Evolution
Uppsala University
Sweden



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