Reference points based on life histories

lhRef calculates a variety of reference points i.e. population growth rate at small population sizs (r), and at B~MSY (rc), ratio of virgin biomass to B~MSY (sk), life time reproductive output (srp0) and reproductive output at B~MSY (sprmsy)

lhRef(
  params,
  m = function(length, params) exp(0.55) * (length^-1.61) %*% (params["linf"]^1.44)
    %*% params["k"],
  sr = "bevholt",
  range = c(min = 0, max = 40, minfbar = 1, maxfbar = 40, plusgroup = 40),
  what = c("r", "rc", "msy", "lopt", "sk", "spr0", "sprmsy"),
  msy = "msy"
)

Arguments

params

FLPar

m

function for natural mortality

sr

character with stock recruitment relationship, e.g. "bevholt","ricker",...

range

ages used i.e. c(min=0,max=40,minfbar=1,maxfbar=40,plusgroup=40)

what

quantities to calculate "r","lopt","rc","sk","spr0","sprmsy"

msy,

character with "msy", "f0.1", ...

Value

object of class FLPar with reference points, i.e r, rc, sk, lopt,

See also

lhPar, lhEql

Examples

if (FALSE) {
library(FLBRP)
params=FLPar(linf=100,t0=0,k=.4)
params=lhPar(params)
lhRef(params)
}