rffwd: rffwd() Project forward an FLStock with evolutionary Fbar

In Henning-Winker/FLRef: Reference point computation for advice rules

rffwd() Project forward an FLStock with evolutionary Fbar

Description

rffwd() Project forward an FLStock with evolutionary Fbar

Usage

rffwd(object, sr, fbar = control, control = fbar, deviances = "missing")

Arguments

object	An *FLStock*
sr	A stock-recruit relationship, *FLSR* or *predictModel*.
fbar	Yearly target for average fishing mortality, *FLQuant*.
control	Yearly target for average fishing mortality, *FLPar*.
deviances	Deviances for the strock-recruit relationsip, *FLQuant*.

Value

The projected *FLStock* object.

Examples

data(ple4)
sr <- srrTMB(as.FLSR(ple4,model=bevholtSV),spr0=mean(spr0y(ple4)))
brp = computeFbrp(ple4,sr,proxy="msy") 
fbar(brp) = FLQuant(rep(0.01,70))
stk = as(brp,"FLStock")
units(stk) = standardUnits(stk)
its = 100
stk <- FLStockR(propagate(stk, its))
stk@refpts= Fbrp(brp)
b0=an(Fbrp(brp)["B0"])
control = FLPar(Feq=0.15,Frate=0.1,Fsigma=0.15,SB0=b0,minyear=2,maxyear=70,its=its)
run <- rffwd(stk, sr=sr,control=control,deviances=ar1rlnorm(0.3, 1:70, its, 0, 0.6))
plotAdvice(run)

Henning-Winker/FLRef documentation built on Dec. 16, 2024, 7:23 p.m.