getx: Estimate an initial vector of depletion values

In cttedwards/bdm: Bayesian biomass dynamics model

Estimate an initial vector of depletion values

Description

This function can be applied to a list or bdmData object to provide a rough estimate of the depletion time series x.

Usage

getx(object, ...)

## S3 method for class 'bdm'
getx(object, ...)

## S3 method for class 'list'
getx(object, r, logK, ...)

Arguments

object	either a bdm class object or a list object containing a harvest vector
...	additional arguments to generic function
r	an assumed value for the intrinsic growth rate r
logK	an assumed value for the carrying capacity ln(K)

Details

This function assumes a logistic production function and uses input values for the intrinsic growth and carrying capacity. These would usually be obtained through calls to rCalc in the lhm package and getlogK.

Examples

# get some data
data(albio)
dat <- bdmData(harvest = albio$catch, index = albio$cpue, time = rownames(albio))

# default model
mdl <- bdm()

# extract depletion from
# bdm object
getx(mdl)

# extract with data
mdl@data <- dat
getx(mdl)

# calculate depletion from
# catches and assumed parameters
getx(dat, logK = 6, r = exp(-1))

cttedwards/bdm documentation built on Oct. 11, 2022, 7:52 p.m.