Jarrod Hadfield wrote: > Hi Wayne, > > You should pass the mcmc objects from >2 models as a list. For > example, the fixed effects: > > diag.gelman(mcmc.list(m1$Sol, m2$Sol)) PS I think that's gelman.diag()