Skip to content
Prev 18513 / 20628 Next

Specifying starting point for MCMCglmm()

Hi everyone~

I am running a series of large models in R using MCMCglmm. After running the models for 14 days (complicated, multivariate models), I found that the models have not converged yet and I need to run more iterations. I know you can use the start= to specify a starting function, but I?m having difficulties getting the model to run.

model_EC <- parLapply(cl=cl,listEC, function(i) {
  MCMCglmm(i[[1]],
           random = ~us(trait):phylo, family=rep("gaussian", 6), rcov=~us(trait):units,
           ginverse=list(phylo=inv.phyloEC$Ainv),prior=prior.1,data=i[[2]], start=1820000,
           nitt=2600000,thin=1400,burnin=0)})

As the code is written, I have the start=1820000, which is where the last iteration left off. Can anyone explain what I?m doing wrong in specifying the start function?

Any help or suggestions would be greatly appreciated!

Best,
Emily Cornelius Ruhs
Postdoctoral Scholar
University of South Florida