Calculates the energy rate (grams/year) available for reproduction after growth and metabolism have been accounted for. You would not usually call this function directly but instead use getERepro(), which then calls this function unless an alternative function has been registered, see below.

mizerERepro(params, n, n_pp, n_other, t, e, ...)

Arguments

params

A MizerParams object

n

A matrix of species abundances (species x size).

n_pp

A vector of the resource abundance by size

n_other

A list of abundances for other dynamical components of the ecosystem

t

The time for which to do the calculation (Not used by standard mizer rate functions but useful for extensions with time-dependent parameters.)

e

A two dimensional array (species x size) holding the energy available for reproduction and growth as calculated by mizerEReproAndGrowth().

...

Unused

Value

A two dimensional array (species x size) holding $$\psi_i(w)E_{r.i}(w)$$ where \(E_{r.i}(w)\) is the rate at which energy becomes available for growth and reproduction, calculated with mizerEReproAndGrowth(), and \(\psi_i(w)\) is the proportion of this energy that is used for reproduction. This proportion is taken from the params object and is set with setReproduction().

Your own reproduction rate function

By default getERepro() calls mizerERepro(). However you can replace this with your own alternative reproduction rate function. If your function is called "myERepro" then you register it in a MizerParams object params with

params <- setRateFunction(params, "ERepro", "myERepro")

Your function will then be called instead of mizerERepro(), with the same arguments.

See also