fillell2 | R Documentation |
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.
fillell2(biot, intheta, mult, ct, Hmax = 1, Fyear)
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. |
a vector of 1 and 0 relating to an acceptable trajectory (one that met the criteria) and unacceptable (one that failed the criteria)
## 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
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.