This function performs projections for dbsra and catchmsy objects

Description

Make biomass projections by using inputted catch and results of dbsra or catchmsy functions

Usage

1
2
3
4
dlproj(dlobj = NULL, projyears = NULL, projtype = 1, projcatch = NULL, 
grout = 1, grargs = list(lwd = 1, unit = "MT", mains = " ", cex.main = 1, 
cex.axis = 1, cex.lab = 1), grtif = list(zoom = 4, width = 11, height = 13, 
pointsize = 10))

Arguments

dlobj

function dbsra or catchmsy output object

projyears

the number of years for projection. The first year will be the last year of catch data plus one in the original dbsra or catchmsy run.

projtype

the type of catch input. 0 = use median MSY from dbsra or catchmsy object, 1 = use mean MSY from dbsra or catchmsy object, 2 = user-inputted catch

projcatch

if projtype = 2, a single catch value used over all projection years or a vector of catch values (length is equal to projyears).

grout

numeric argument specifying whether projection graph should be shown on the console only (grout==1) or shown on the console and exported to a TIF graph file (grout==2). No graph (grout== 0). If plotted, the median (solid line), mean (dashed line), and 2.5th and 97.5 percentiles(dotted lines) are displayed. Use setwd before running function to direct .tif file to a specific directory. The name of .tif file is automatically determined.

grargs

list control arguments for plotting functions. lwd is the line width, unit is the biomass unit for the y-axis label,mains and cex.main are the title and character expansion value for the graph, cex.axis is the character expansion value for the x and y-axis tick labels and cex.lab is the character expansion value(s) for the x and y-axis labels.

grtif

list control arguments for the .TIF graph file. See tiff help file in R.

Details

The biomass estimate of the last year+1 is used as the starting biomass (year 1 in projections) and leading parameters from each plausible (accepted) run are used to project biomass ahead projyears years using either the MSY estimate (median or mean) from all plausible runs or inputted catch values. The biomass estimates are loaded from either the "Biotraj-dbsra.csv" or "Biotroj-cmsy.csv" files that were automatically saved in functions "dbsra" and "catchmsy".

Use setwd() before running the function to change the directory where .csv files are stored.

Value

type

object projection type

ProjBio

dataframe of biomass projections for each plausible run

Author(s)

Gary A. Nelson, Massachusetts Division of Marine Fisheries gary.nelson@state.ma.us

References

Martell, S. and R. Froese. 2012. A simple method for estimating MSY from catch and resilience. Fish and Fisheries 14:504-514.

Dick, E. J. and A. D. MacCall. 2011. Depletion-based stock reduction analysis: a catch-based method for determining sustainable yield for data-poor fish stocks. Fisheries Research 110: 331-341.

See Also

catchmsy dbsra

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
## Not run: 
  data(lingcod)
   outs<-catchmsy(year=lingcod$year,
    catch=lingcod$catch,catchCV=NULL,
    catargs=list(dist="none",low=0,up=Inf,unit="MT"),
    l0=list(low=0.8,up=0.8,step=0),
    lt=list(low=0.01,up=0.25,refyr=2002),sigv=0,
    k=list(dist="unif",low=4333,up=433300,mean=0,sd=0),
    r=list(dist="unif",low=0.015,up=0.5,mean=0,sd=0),
    bk=list(dist="unif",low=0.5,up=0.5,mean=0,sd=0),
    M=list(dist="unif",low=0.24,up=0.24,mean=0.00,sd=0.00),
    nsims=30000)
   outbio<-dlproj(dlobj = outs, projyears = 20, projtype = 0, grout = 1)
 
## End(Not run)

Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.