Skip to content

MCMCglmm: priors for ordinal regression

2 messages · m.fenati at libero.it, Jarrod Hadfield

3 days later
#
Dear Massimo,

Do you mean the chain did not converge or the chain did not mix?  
Generally the former is rare, and is usually only seen with  
ordinal/categorical data with complete (or near complete) separation.   
Sometimes a prior that constrains the linear predictor away from  
extreme values on the logit/probit scale can fix this with a  
relatively minor prior influence on inferences made on the data scale.  
Sometimes not. Its not clear to me what the motivation is behind your  
prior - is it that the sum of your variance components is close to  
100? If so I would be careful. Use pl=TRUE in your call to MCMCglmm  
and make sure your latent variables are in the range -7 to 7.

Cheers,

Jarrod





Quoting "m.fenati at libero.it" <m.fenati at libero.it> on Wed, 4 Jul 2012  
16:48:18 +0200 (CEST):