[ADMB Users] How to use the mcr command?

[Jim Ianelli]

Hi,

as a minimalist test I made the model:
DATA_SECTION

PARAMETER_SECTION
  init_number a
  sdreport_number b
  objective_function_value f
PROCEDURE_SECTION
  b=a;
  f = square(a);

and invoked mcmc as:

program -mcmc 1000 -mcsave 10

and using readbin:
readbin program.psv

I found there was a chain of 100 elements.

restarting as:

program -mcr -mcmc 1000 -mcsave

I got 200 elements in the chain...

So it should work just to add -mcr with the previous arguments.

8 days to get 150000 MCMCs...wow.

Cheers,
Jim


On Tue, Aug 7, 2012 at 4:58 PM, Edgar Gonzalez <edgarjgonzalez at ymail.com>wrote:

> Hi everyone,
>
> Does anyone know how the -mcr command works? I've searched the web and
> couldn't find its use, apart from the vague:
>
> -mcr resume previous mcmc
>
> I ran my code (which is not a random effects model) as:
>
> program -mcmc 500000 -mcsave -1500
>
> which took about 8 days to run. My problem is that I can only run codes
> for 10 days and then they're automatically stopped. So I'll have to
> partition my mcmc run as I require 16 million steps (calculated from a
> previous diagnostics run).
>
> Could anyone help me?
>
> Thanks in advance,
>
> Edgar J. González
> Ecology Department
> Science Faculty, UNAM
>
> _______________________________________________
> Users mailing list
> Users at admb-project.org
> http://lists.admb-project.org/mailman/listinfo/users
>
>


-- 

James Ianelli
REFM Division
Alaska Fisheries Science Center
NMFS/NOAA Building 4
7600 Sand Pt Way NE
Seattle WA 98115

 206 526 6510
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.admb-project.org/pipermail/users/attachments/20120807/021cbbee/attachment.html>