getRoot5D: Get the Root of the Full Ecosystem Model in 5D

In cboettig/bs-tipping: Bifurcation or State Tipping in a Trophic Cascade Experiment

Description

Numerically finds the roots of the ecosystem model given starting values (adults, [sun]fish, juveniles, herbivores, and phytoplankton) and parameter (harvest of adults, qE)

Usage

getRoot5D(x, pars, maxiter = 1000, ...)

Arguments

x	a vector of length 5 with names A0 (adult bass), F0 (planktivorous fish), J0 (juvenile bass), H0 (herbivorous zooplankton), and P0 (phytoplankton)
pars	vector of parameters to pass to model
maxiter	maximum number of iterations to use when finding root
...	additional arguments to pass to ecoStep

Value

numeric vector with roots

Examples

getRoot5D(c(A0=15, F0=50, J0=10, H0=100, P0=100), pars=c(qE=-0.02)) # stocking, high starting A
getRoot5D(c(A0=2, F0=50, J0=5, H0=100, P0=100), pars=c(qE=-0.02)) # stocking, but start low A
getRoot5D(c(A0=15, F0=1, J0=15, H0=100, P0=100), pars=c(qE=0.02)) # harvesting, start high A

cboettig/bs-tipping documentation built on May 5, 2019, 7:08 a.m.