dlproj: This function performs projections for dbsra and catchmsy...

In fishmethods: Fishery Science Methods and Models

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

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 [e-mail]

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

## 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)

fishmethods documentation built on April 27, 2023, 9:10 a.m.