Get consumption rate for each species

For each species, this function returns the rate at which food is consumed by all individuals of that species. This corresponds to the consumption rate \(Q\) used by Ecopath models.

Usage

getConsumption(params, min_w_pred = 0, max_w_pred = Inf)

Arguments

params

A MizerParams object

min_w_pred

The minimum weight of predators to include in the diet matrix. Default is 0.

max_w_pred

The maximum weight of predators to include in the diet matrix. Default is Inf.

Value

A named vector of consumption rate for each species

Details

The consumption is calculated as $$Q_i = \int_{w_0}^{w_{max}} E_{e.i}(w) (1 - f_i(w)) N_i(w) dw$$ where \(E_{e.i}(w)\) is the encounter rate of an individual of species \(i\) and weight \(w\) (calculated with getEncounter), \(f_i(w)\) is the feeding level of an individual of species \(i\) and weight \(w\) (calculated with getFeedingLevel) and \(N_i(w)\) is the number density of species \(i\) at weight \(w\).

See also

Other rate functions: getDietMatrix(), getEcotrophicEfficiency(), getGonadicProduction(), getM0B(), getM2B(), getMetabolicRespiration(), getOffspringProduction(), getProduction(), getReproductiveLoss(), getRespiration(), getSomaticProduction(), getTotalProduction(), getUnassimilated(), getZB()

Examples

getConsumption(NS_params)
#>        Sprat      Sandeel       N.pout      Herring          Dab      Whiting 
#> 2.395035e+12 3.648862e+13 7.082083e+12 1.004841e+13 7.802687e+10 1.338066e+12 
#>         Sole      Gurnard       Plaice      Haddock          Cod       Saithe 
#> 5.798710e+11 5.472728e+11 9.504280e+12 3.801462e+12 2.031014e+12 2.467478e+12