Calculates the predation mortality rate \(\mu_p(w)\) on the resource spectrum by resource size (in units 1/year). You would not usually call this function directly but instead use getResourceMort(), which then calls this function unless an alternative function has been registered, see below.

mizerResourceMort(params, n, n_pp, n_other, t, pred_rate, ...)



A MizerParams object


A matrix of species abundances (species x size).


A vector of the resource abundance by size


A list of abundances for other dynamical components of the ecosystem


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


A two dimensional array (predator species x prey size) with the predation rate, where the prey size runs over fish community plus resource spectrum.




A vector of mortality rate by resource size.

Your own resource mortality function

By default getResourceMort() calls mizerResourceMort(). However you can replace this with your own alternative resource mortality function. If your function is called "myResourceMort" then you register it in a MizerParams object params with

params <- setRateFunction(params, "ResourceMort", "myResourceMort")

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

See also