MQMF (version 0.1.1)

spmprojDet: 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

# NOT RUN {
 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)
# }