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
Dr. Wayne Dawson Institute of Plant Sciences University of Bern Altenbergrain 21 3013 Bern Switzerland +41 (0)31 631 49 25