MCMCglmm does not produce more than one chain. So the way to check convergence is to run the model twice, each with different starting values, then conduct the diagnostic with:
gelman.diag(m1$Sol, m2$Sol)
2 comments:
Anonymous
said...
This did not work - it came up with
Error in gelman.diag(m1$Sol, m2$Sol) : You need at least two chains
2 comments:
This did not work - it came up with
Error in gelman.diag(m1$Sol, m2$Sol) : You need at least two chains
I got it
you have to do
gelman.diag(list(m1$Sol,m2$Sol))
Post a Comment