Skip to content
Prev 3040 / 20628 Next

Convergence diagnostics for MCMCglmm

Dear R users,

I've managed to get to grips with using MCMCglmm for phylogenetic  
meta-analysis, thanks to help on here from Jarrod Hadfield. However, I  
haven't yet worked out how to use convergence diagnostics on MCMCglmm  
output.

I am currently running models for 50,000 iterations, with a long  
burnin of 25000, and a thinning interval of 10. Whilst my trace plots  
look good, I would like to use the Gelman-Rubin or Raftery-Lewis  
diagnostics in the coda package to check that my burnin/no. of  
iterations are adequate. I've read relevant sections in Ch. 7 of Ben  
Bolkers ecological models book, and the coda/MCMCglmm package pdfs,  
but I'm still not sure what mcmc output from MCMCglmm I am supposed to  
pass to the coda diagnostic functions. Apologies for the simple  
question, but hopefully there is a simple answer!

Any suggestions appreciated as always,

Thanks,

Wayne