fudc: fudc()

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

fudc()

Description

generates an up-down-constant F-pattern

Usage

fudc(
  object,
  fref = 0.2,
  fhi = 2.5,
  flo = 0.8,
  sigmaF = 0.2,
  breaks = c(0.5, 0.75)
)

Arguments

object	An *FLStock*
fref	reference denominator for fbar
fhi	factor for high F as fhi = fbar/fref
flo	factor for low F as flo = fbar/fref
sigmaF	variation on fbar
breaks	relative location of directional change

Value

FLQuant

Examples

data(ple4)
sr <- srrTMB(as.FLSR(ple4,model=bevholtSV),spr0=mean(spr0y(ple4)))
brp = computeFbrp(ple4,sr,proxy="msy")
fmsy = Fbrp(brp)["Fmsy"]
stki = propagate(ple4,100)
fy = fudc(ple4,fhi=2,flo=0.9,fref=fmsy,sigmaF=0)
fyi = fudc(stki,fhi=2,flo=0.9,fref=fmsy,sigmaF=0.2)
plot(fy,fyi)+ylab("F")
#Forcasting
om <- FLStockR(ffwd(stki,sr,fbar=fyi))
om@refpts = Fbrp(brp)
plotAdvice(window(om,start=1960))

Henning-Winker/FLRef documentation built on July 14, 2024, 7:28 a.m.