Skip to contents

Get diffusion rate from predation

Source: R/diffusion.R
getDiffusion.Rd

Calculates the diffusion rate \(D_i(w)\) (grams^2/year) for each species. This is the rate at which the abundance density is diffused along the size axis due to the variability in prey sizes. This is the diffusion term from the jump-growth equation.

Usage

getDiffusion(
  params,
  n = initialN(params),
  n_pp = initialNResource(params),
  n_other = initialNOther(params),
  t = 0,
  ...
)

Arguments

params

A MizerParams object

n

A matrix of species abundances (species x size). Defaults to the initial abundances in params.

n_pp

A vector of the resource abundance by size. Defaults to the initial resource abundances in params.

n_other

A list of abundances for other dynamical components.

t

The time for which to do the calculation.

...

Unused

Value

An array of dimensions species x size holding the diffusion rates.

References

Datta, S., Delius, G. W. and Law, R. (2010). A jump-growth model for predator-prey dynamics: derivation and application to marine ecosystems. Bulletin of Mathematical Biology, 72(6):1361–1382

See also

Other rate functions: getEGrowth(), getERepro(), getEReproAndGrowth(), getEncounter(), getFMort(), getFMortGear(), getFeedingLevel(), getFlux(), getMort(), getPredMort(), getPredRate(), getRDD(), getRDI(), getRates(), getResourceMort()