fillell2: fillell2 runs the criteria of success on each SRA

In haddonm/datalowSA: A Package to Faciliate Application of Data poor Stock Assessments

fillell2 runs the criteria of success on each SRA

Description

fillell2 runs the acceptance criteria over the biomass trajectories obtained from the Stock Reduction Analysis for each combination of r, K, and starting depletion conducted in makebiomC. The criteria agsint which each trajectory are tested include that the final depletion be less than the maximum of the expected range of final depletions and greater than the minimum of final depletion levels. In addition I have included that the initial biomass level be greater than the final biomass level. This differs from the extra criteria used by Martell and Froese (2013), who selected the smallest K still able to generate the mean predicted MSY.

Usage

fillell2(biot, intheta, mult, ct, Hmax = 1, Fyear)

Arguments

biot	the matrix of biomass trajectories obtained from makebiomC, containing the SRA trajectories for the aprticular combination of r and K across the range of starting biomass trajectories.
intheta	the array of parameters associated with the biomass trajectories. Includes r, K, the range of depletions and the process error added to productivity - sigR.
mult	a multiplier for K to allow for stock biomasses to rise above K
ct	the vector of catches per year
Hmax	upper limit of harvest rate included in the constraints; defaults to 1.0, which implies no upper limit.
Fyear	is the index to the year in which a rnage of harvest rates is to be used to constrina the acceptable trajectories.

Value

a vector of 1 and 0 relating to an acceptable trajectory (one that met the criteria) and unacceptable (one that failed the criteria)

Examples

## Not run: 
print("This is the function that imposes constraints on the biomass")
print("trajectories. Add new ones to the end of the long if statement")
print("They need to define what is required to keep a trajectory")

## End(Not run)    #   biot=biot; intheta=itheta[tick,]; mult=mult; ct=ct; Hmax=maxH;Fyear=NA 

haddonm/datalowSA documentation built on Nov. 5, 2023, 6:40 p.m.