| importProj {scapeMCMC} | R Documentation |
Import Coleraine MCMC traces for spawning biomass and catch, projected into the near future.
importProj(dir, info="", coda=FALSE, quiet=TRUE)
dir |
directory containing the files ‘strategy.out’, ‘projspbm.pst’ and ‘procatch.pst’. |
info |
optional string containing information to store with MCMC projections. |
coda |
whether data frames should be coerced to class mcmc
using the coda package. |
quiet |
whether to report progress while parsing files in directory. |
A list containing:
B |
biomass by catch policy and year |
Y |
catch by catch policy and year |
as lists of data frames, or mcmc objects if coda=TRUE.
MCMC projections can be used to evaluate the short-term outcome of harvest policies (constant catch or constant harvest rate), given the uncertainty about parameter values and random future recruitment.
The example dataset xproj was generated using
importProj.
The functions ll (package gdata) and head
are recommended for browsing MCMC projections, e.g.
ll(xproj); ll(xproj$B); ll(xproj$B$"0.2"); head(xproj$B$"0.2").
Arni Magnusson arnima@u.washington.edu.
Hilborn, R., M. Maunder, A. Parma, B. Ernst, J. Payne, and P. Starr. 2003. Coleraine: A generalized age-structured stock assessment model. User's manual version 2.0. University of Washington Report SAFS-UW-0116. Available at http://fish.washington.edu/research/coleraine/coleraine.pdf.
importMCMC, read.table,
readLines, scan, xproj.
scapeMCMC-package gives an overview of the package.
## Not run:
path <- paste(.find.package("scapeMCMC"), "/example", sep="")
xproj <- importProj(path)
## End(Not run)