Skip to content
Prev 8366 / 20628 Next

mcmcglmm and parallel chains

I am learning mcmcglmm in order to use it on a beowulf cluster.

In https://stat.ethz.ch/pipermail/r-sig-mixed-models/2011q3/006558.html

Jarrod Hadfield writes:

"You can merge MCMC chains from multiple runs, although you should make
sure you start them from different initial values"

Is it sufficient to provide differents random seeds for each run,
or does this refer to the start parameter of mcmcglmm()?

   start: optional list having 4 possible elements: ?R? (R-structure)
          ?G? (G-structure) and ?liab? (latent variables or
          liabilities) should contain the starting values where ?G?
          itself is also a list with as many elements as random effect
          components. The fourth element ?QUASI? should be logical: if
          ?TRUE? starting latent variables are obtained heuristically,
          if ?FALSE? then they are sampled from a Z-distribution