multiSA (version 0.1.0)

calc_growth: Calculate von Bertalanffy length-at-age

Description

Returns an array of length-at-age with seasonal dimension. Useful for MSAdata inputs.

Usage

calc_growth(
  Linf_s,
  K_s,
  t0_s,
  ns = length(Linf_s),
  nm = 4,
  ny = 20,
  a = seq(1, 10) - 1
)

Value

Array [y, m, a, s]

Arguments

Linf_s

Vector by stock s of asymptotic length

K_s

Vector by stock s of the growth coefficient

t0_s

Vector by s of the age at length zero.

ns

Integer, number of stocks

nm

Integer, number of seasons

ny

Integer, number of years

a

Integer vector of ages

Examples

len_ymas <- calc_growth(c(30, 40), c(0.4, 0.2), c(-1, -1))

# Calculate stock weight at age
a_s <- rep(1e-6, 2)
b_s <- c(3, 3.1)

ns <- length(a_s)
swt_ymas <- sapply(1:ns, function(s) {
  a_s[s]*len_ymas[, , , s]^b_s[s]
}, simplify = "array")