spmprojDet: spmprojDet conducts forward projections from known conditions

In MQMF: Modelling and Quantitative Methods in Fisheries

spmprojDet conducts forward projections from known conditions

Description

spmprojDet conducts deterministic forward projections of the dynamics of a fitted surplus production model. The parameters and original data need to be put through the function spm to form the spmobj, i.e. the list generated by the function spm. This contains all required information except for details of the projection. The application of spm is where the dynamics are defined as either Schaefer or Fox. If no plot is generated then the projected dynamics are output invisibly, where the biomass and predCE are matrices of years vs projcatch.

Usage

spmprojDet(spmobj, projcatch, projyr = 10, plotout = FALSE, useft = 7)

Arguments

spmobj	the list generated by the function spm
projcatch	the projected constant catch levels as a vector
projyr	the number of years of projection. default = 10
plotout	should the projection be plotted? default=FALSE
useft	which font used in a plot? default=7 = bold times

Value

the projected biomass, CPUE, and the projected years

Examples

 data(abdat)
 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=FALSE)
 out <- spm(bestmod$estimate,indat=abdat,schaefer=FALSE)
 catches <- seq(700,1000,50)
 spmprojDet(spmobj = out,projcatch=catches,projyr=10,plotout=TRUE)

MQMF documentation built on Sept. 8, 2023, 5:14 p.m.