spmproj {MQMF}

spmproj calculates biomass trajectories for replicate parameters

Description

spmproj uses a matrix of parameter vectors to project surplus production dynamics forward including future projection years under constant catches. This is used to conduct risk assessments for different constant catches allowing a search for optimum future catch levels.

Usage

spmproj(
  parmat,
  indat,
  constC,
  projyr = 10,
  year = "year",
  cats = "catch",
  index = "cpue"
)

Arguments

parmat

a matrix of N parameter vectors obtained from either asymptotic errors (parasympt), bootstraps (bootpar), or from a Bayesian analysis (parsbayes).

indat

the fisheries data used during model fitting

constC

the constant catch level imposed in the projection years

projyr

the number of years of projection, default = 10

year

name of the year variable within indat, default=year

cats

name of the catch variable within indat, default=catch

index

name of the cpue variable within indat, default=cpue

Value

an N x years matrix of biomass trajectories, one for each parameter vector

Examples

  data(abdat)
  schf <- FALSE
  param <- log(c(r=0.3,K=11500,Binit=3300,sigma=0.05))
  bestmod <- nlm(f=negLL1,p=param,funk=simpspm,logobs=log(abdat$cpue),
                 indat=abdat,typsize=magnitude(param),iterlim=1000,
                 schaefer=schf,hessian = TRUE)
  out <- spm(bestmod$estimate,indat=abdat,schaefer=schf)
  matpar <- parasympt(bestmod,10) # normally use 1000 or more
  projs <- spmproj(matpar,abdat,projyr=10,constC=900)
  plotproj(projs,out)