Hi Edgar,<div>You should be able to set a different random number seed for MCMC using the -mcseed input. You can also theoretically change the starting point for the MCMC using -mcpin. However, multiple short MCMC chains are not a replacement for a single long chain, and if the convergence diagnostics suggest 6 000 000 samples, I think you are very unlikely to get a good approximation to the posterior with 12 runs of 500 000.</div>

<div><br></div><div>I would say that a more fruitful path would be to look at which parameters were the most correlated with each other or had the most autocorrelation in the shorter MCMC and look at reparameterizing the model or reducing the number of parameters.</div>

<div>-Ian<br><br><div class="gmail_quote">On Thu, Jul 5, 2012 at 9:00 AM, Edgar Gonzalez <span dir="ltr"><<a href="mailto:edgarjgonzalez@ymail.com" target="_blank">edgarjgonzalez@ymail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

<div><div style="font-size:10pt;font-family:times new roman,new york,times,serif"><div><span>Hi everyone:</span></div><div><span><br></span></div><div><span>I'm trying to run mcmc on ADMB. I ran a pre-run to estimate the Total N required to obtain the posterior distribution of my parameters. I got a rounded N of 6 000 000. As this would take weeks (months?) to do, I'd like to know if there is a way to partition this into 12 runs with N = 500 000 and to assemble them into a single run. As far as I know, the mcmc involves random walks, but when I ran two mcmc's I got the same walk (the traces look the same).</span></div>

<div><span><br></span></div><div><span>Thanks in advance,</span></div><div> </div><div>Edgar J. González<br>Ecology and Natural Resources Dept.<br>Science Faculty, UNAM</div></div></div><br>_______________________________________________<br>


Users mailing list<br>
<a href="mailto:Users@admb-project.org" target="_blank">Users@admb-project.org</a><br>
<a href="http://lists.admb-project.org/mailman/listinfo/users" target="_blank">http://lists.admb-project.org/mailman/listinfo/users</a><br>
<br></blockquote></div><br></div>