[ADMB Users] Boostrap
SteveM at iphc.int
Tue Nov 26 10:51:19 PST 2013
You can use the mceval phase to do such at thing; something along the lines of
ofs = ofstream("ssb.mcmc",ios::app);
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:
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 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,
> Users mailing list
> Users at admb-project.org
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