getProductionC: getProductionC estimates MSY from output of ASPM

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

getProductionC estimates MSY from output of ASPM

Description

getProductionC takes the optimum estimate of R0 from ASPM and estimates the production curve, from which it gains the MSY, Bmsy, Hmsy, and Dmsy (depletion at MSY). It generate the production curve by stepping through the dynamics of the fishery over a series of constant harvest rates for however many iterations of nyr required. It is identical to getProduction but uses an Rcpp routine to speed up the iterative steps by about 26 times. That is very useful when conducting bootstraps

Usage

getProductionC(
  inR0,
  infish,
  inglb,
  inprops,
  Hrg = c(0.025, 0.4, 0.025),
  nyr = 50,
  maxiter = 3
)

Arguments

inR0	the optimum estimate of R0 from ASPM
infish	the fish object from a data set
inglb	the glb object from a data set
inprops	the props object from a data set
Hrg	the desired sequence of harvest rates to be used to define the production curve. The default is c(0.025,0.4,0.025), but it is a good idea to plot teh harvest rate against yield and manually adjust the sequence values to focus attention and detail on where the MSY falls.
nyr	The number of years making up a single iteration while searching for the equilibrium yield. default = 50.
maxiter	defaults to 3. Determines how many runs through the nyr steps are conducted to reach equilibrium. One can trial different values but 3 is a reasonable compromise. It is best to test different maxiter values to ensure equilibria are reached in each yield estimate.

Value

a data.frame containing the sequence of harvest rates, the related spawning biomass, the exploitable biomass, the yield, and depletion.

Examples

## Not run: 
data(fishdat)
fish <- fishdat$fish
glb <- fishdat$glb
props <- fishdat$props
pop1 <- unfished(glb,props,exp(glb$R0))
B0 <-pop1$B0
M <- glb$M
nyrs <- 30
years <- 0:nyrs
pars <- c(13.6,0.3)
bestL <- optim(par,aspmLL,method="Nelder-Mead",infish=fish,inglb=glb,
               inprops=props,
               control=list(maxit = 1000, parscale = c(10,0.1)))
prod <- getProduction(exp(bestL$par[1])infish,inglb,inprops,
                      Hrg=c(0.0005,0.07,0.0005),nyr=50)
print(prod)

## End(Not run)

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