[ADMB Users] Boostrap

[Steve Martell]

Hi Luis,

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


PROCEDURE_SECTION

        if(mceval_phase())
        {
                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.

Steve
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.