plotconstC | R Documentation |
plotconstC plots and summarizes the outcome of constant catch projections made on each of the successful trajectories from the catch-MSY analysis. The catch-MSY analysis provides an uncertain estimate of MSY and of final depletion. But by conducting constant catch projectios the implications in terms of what proportion of trajectories increase and what proportion decrease can aid decisions and could form he basis for the development of formal harvest control rules. It will be noticed that a large proportion of trajectories can be below 20 notice also that at the same time a high proportion can be above 48 assuming that this Commonwealth default target is used. It is the case that a target of 0.4 or 40 species, as a means of preventing the management of such species from prevent the primary economic drivers of a fishery from being caught.
plotconstC(
deplete,
endyear,
constC = 0,
limit = 0.2,
target = 0.48,
console = TRUE,
intensity = NA,
contours = TRUE
)
deplete |
the output of 'gettraject', a matrix of biomass trajectories |
endyear |
the final year of the known catches and biomass trajectories |
constC |
the constant catch applied to generate the proections |
limit |
biomass depletion level used as a limit for the species |
target |
the depletion level used as a target for the species concerned. |
console |
a boolean determinng whether the projection years results are printed to the console |
intensity |
defaults to NA but if contains a value this is the density of trajectories that lead to full colour |
contours |
default to TRUE, should the 80 be plotted? |
a matrix of all years with the proportion < 20 (the input target), the mean and median depletion, and the proportion of trajectories that were increasing relative to the endyear of data.
## Not run:
data(invert)
fish <- invert$fish
glb <- invert$glb
reps <- 10000 # one would run at least 20000, preferably more
answer <- run_cMSY(fish,glb,n=reps,sigpR=1e-8)
traject <- gettraject(answer$R1,projn=5)
newtraj <- doproject(traject,constC=300)
trajdepl <- makedeplet(newtraj)
plotconstC(trajdepl,endyear=2017,constC=300,limit=0.2,target=0.40)
## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.