getPredRate: 'getPredRate' method for the size based model

In pwoodworth-jefcoats/Size-Based-Modeling: Mizer With The Effects Of Temperature Included

getPredRate method for the size based model

Description

Calculates the predation rate of each predator species at size on prey size. In formulas

\phi_i(w_p/w) (1-f_i(w)) \gamma_i w^q N_i(w) dw

This method is used by the project_therMizer method for performing simulations. In the simulations, it is combined with the interaction matrix (see therMizerParams) to calculate the realised predation mortality (see getM2).

getPredRate method with feeding_level argument.

Usage

getPredRate(object, n, n_pp, feeding_level)

## S4 method for signature 'therMizerParams,matrix,numeric,matrix'
getPredRate(object, n, n_pp,
  feeding_level)

Arguments

object	A therMizerParams object.
n	A matrix of species abundance (species x size).
n_pp	A vector of the background abundance by size.
feeding_level	The current feeding level (optional). A matrix of size no. species x no. size bins. If not supplied, is calculated internally using the getFeedingLevel() method.
ocean_temp	A numeric vectory of the ocean temperature by realm or a single numeric ocean_temp value which is used for all realms.

Value

A three dimensional array (predator species x predator size x prey size), where the predator size runs over the community size range only and prey size runs over community plus background spectrum.

See Also

project_therMizer, getM2, getFeedingLevel and therMizerParams

Examples

## Not run: 
data(NS_species_params_gears)
data(inter)
params <- therMizerParams(NS_species_params_gears, inter)
# With constant fishing effort for all gears for 20 time steps
sim <- project_therMizer(params, t_max = 20, effort = 0.5)
# Get the feeding level at one time step
n <- sim@n[21,,]
n_pp <- sim@n_pp[21,]
getPredRate(params,n,n_pp)

## End(Not run)

pwoodworth-jefcoats/Size-Based-Modeling documentation built on Sept. 15, 2023, 8:13 a.m.