The suitability function for predation used in the R model is: $$S_{pred,prey}(L,l) = \frac{\delta}{1+e^{-\alpha-\beta l-\gamma L}}$$ With one predator, one prey and otherfood the equation becomes: $$C_{L,l}=N_{L}M_{L}\Psi_{L}\frac{F_{L,l}}{\sum_lF_{L,l}+OA}$$ $$=N_{L}M_{L}\frac{F_{L,l}}{\sum_lF_{L,l}+OA+HA}$$ where $O$ is the density of otherfood.

suitability(
  params,
  l,
  L = c(0),
  type = "exponential",
  to.data.frame = FALSE,
  normalize = FALSE
)

Arguments

params

suitability paramters

l

prey length group(s)

L

predator length group(s)

type

suitability function

to.data.frame

Logical, return a data.frame instead of an array

normalize

Logical, should the output be normalized

Value

matrix of suitabilities, columns prey length, lines predator length