A range of stock-recruitment (SR) models commonly used in fisheries science are provided in FLCore.
ricker() bevholt() segreg() geomean() shepherd() cushing() rickerSV() bevholtSV() shepherdSV() bevholtAR1() rickerAR1() segregAR1() rickerCa()
spr0
.Beverton, R.J.H. and Holt, S.J. (1957) On the dynamics of exploited fish populations. MAFF Fish. Invest., Ser: II 19, 533.
Needle, C.L. Recruitment models: diagnosis and prognosis. Reviews in Fish Biology and Fisheries 11: 95-111, 2002.
Ricker, W.E. (1954) Stock and recruitment. J. Fish. Res. Bd Can. 11, 559-623.
Shepherd, J.G. (1982) A versatile new stock-recruitment relationship for fisheries and the construction of sustainable yield curves. J. Cons. Int. Explor. Mer 40, 67-75.
# inspect the output of one of the model functions bevholt()#> $logl #> function (a, b, rec, ssb) #> loglAR1(log(rec), log(a * ssb/(b + ssb))) #> <bytecode: 0x1046d068> #> <environment: 0xb9fd090> #> #> $model #> rec ~ a * ssb/(b + ssb) #> <environment: 0xb9fd090> #> #> $initial #> function (rec, ssb) #> { #> a <- max(quantile(c(rec), 0.75, na.rm = TRUE)) #> b <- max(quantile(c(rec)/c(ssb), 0.9, na.rm = TRUE)) #> return(FLPar(a = a, b = a/b)) #> } #> <bytecode: 0xf13cb28> #> <environment: 0xb9fd090> #> attr(,"lower") #> [1] -Inf -Inf #> attr(,"upper") #> [1] Inf Inf #>names(bevholt())#> [1] "logl" "model" "initial"bevholt()$logl#> function(a, b, rec, ssb) #> loglAR1(log(rec), log(a*ssb/(b+ssb))) #> <bytecode: 0x1046d068> #> <environment: 0x9fe7b38># once an FLSR model is in the workspace ... data(nsher) # the three model-definition slots can be modified # at once by calling 'model<-' with # (1) a list model(nsher) <- bevholt() # (2) the name of the function returning this list model(nsher) <- 'bevholt' # or (3) the function itself that returns this list model(nsher) <- bevholt