stability.value {fluxweb} | R Documentation |
Estimates network stability
Description
Computes resilience of the system through Jacobian matrix eigenvalues.
Usage
stability.value(
val.mat,
biomasses,
efficiencies,
metabolic.types,
ef.level = "prey",
full.output = FALSE
)
Arguments
val.mat |
A matrix describing fluxes between species (usually a result of |
biomasses |
A vector of species biomasses. |
efficiencies |
A vector or an array of conversion efficiencies of species in the adjacency matrix. These values describe the proportion of consumed energy that is converted to biomass of the consumer. |
metabolic.types |
A vector containing information on species type ( |
ef.level |
Set to |
full.output |
Logical, if |
Details
efficiencies
: Determines how efficient species are to convert energy (seeef.level
for more details). Providing an array will assume values depending on both prey and predator identity.ef.level
: If"prey"
(resp"pred"
), the total amount of energy that can be metabolized from a trophic link will be determined by prey (resp pred) identity."link.specific"
assumes that efficiencies are defined for each trophic interaction and impliesefficiencies
parameter to be a matrixfull.output
: IfTRUE
, function result is a list of eigenvalues and eigenvectors of the Jacobian matrix.
Value
Maximum eigenvalue of the Jacobian matrix of a Lotka Voltera like system of equations. If full.output, Jacobian eigenvalues and eigenvectors are returned.
Author(s)
Benoit Gauzens, benoit.gauzens@gmail.com
Examples
losses = 0.15 * groups.level$bodymasses^(-0.25)
# define metbolic types:
met.types = rep('animal', length(losses))
met.types[groups.level$efficiencies == 0.545] = 'plant'
val.mat = fluxing(groups.level$mat,
groups.level$biomasses,
losses,
groups.level$efficiencies,
bioms.pref = TRUE,
ef.level = "pred")
stability.value(val.mat,
groups.level$biomasses,
groups.level$efficiencies,
metabolic.types = met.types,
ef.level = "pred")