[R-sig-ME] Combining MCMCglmm estimates

Hans Ekbrand hans at sociologi.cjb.net
Thu Oct 11 10:05:44 CEST 2012

On Wed, Oct 10, 2012 at 11:56:52AM +0100, Paul Johnson wrote:
> > Why is different starting values important?
> > Shouldn't burning make the 10 chains independent enough?
> ...
> > The idea that different starting points are needed would, if I
> > understand the rationale correctly, imply that the chains are
> > better in the end than in the beginning. Is that the point?
> I think it's a precaution. Assuming that you have ended up with homogenous looking samples from different runs, you'll be more confident (but never certain of course) that they've converged if they started from different points in parameter space. E.g. there might be local optima where chains could get stuck, and this problem would be much more likely to be discovered starting from different values. I don't know how likely local optima are in practice with a typical MCMCglmm model.
> However, I don't see that 10 samples of (effective) size 1000 from the same starting values, with sufficient burnin, are any worse than 10,000 samples from a single run. So my feeling is that using different starting values is always worthwhile (given how easy this is), but not strictly essential.

Thanks for your answer, Paul.

10 different runs with the same starting point would then be just as
good (or bad) as 1 run, while 10 different runs from different
starting points would be better. Is that a correct conclusion?

More information about the R-sig-mixed-models mailing list