[ADMB Users] Boostrap

Steve Martell SteveM at iphc.int
Tue Nov 26 10:51:19 PST 2013

Hi Luis,

You can use the mceval phase to do such at thing;  something along the lines of


                ofs = ofstream("ssb.mcmc",ios::app);
                ofs <<ssb<<endl;

Then run your model:
        mymodel -mcmc 10000 -mcsave 100

        note tha thte -mcsave command line option will save every nth sample from the posterior
THen run it again using:

        mymodel -mceval
This will then create a file named ssb.mcmc, where each row is a vector of your ssb estimates.  Then read into R and make your CI.

There are other ways as well, this is just off the top of my head.

On Nov 26, 2013, at 10:43 AM, Luis Ridao <luridao at gmail.com> wrote:

> ADMB-comunity,
> ADMB has full capabilities to address parameter uncertainty (e.g. MCMC and likelihood profiles)
> but i was wondering how to implement a simple bootstrap on a fisheries catch-at-age model
> My main interest is to display the uncertainty in ssb and recruitment in a time series fashion with confidence bands and so on.
> is anyone tried this? any simple script on how to perform  this in admb?
> thanks in advance,
> Luis
> _______________________________________________
> Users mailing list
> Users at admb-project.org
> http://lists.admb-project.org/mailman/listinfo/users


This internet e-mail message, and any files transmitted with it, contains confidential, privileged information that is intended only for the addressee. If you have received this e-mail message in error, please call us at (206) 634-1838 collect if necessary) and ask to speak to the message sender. Nothing in this e-mail or the act of transmitting it, is to be construed as a waiver of any rights or privileges enjoyed by the sender or the International Pacific Halibut Commission pursuant to the International Organizations Immunities Act, 22 U.S.C. Sec. 288 et seq.

More information about the Users mailing list