create.jacob {fluxweb} | R Documentation |
Compute the Jacobian matrix
Description
Compute the Jacobian matrix
Usage
create.jacob(
val.mat,
biomasses,
efficiencies,
metabolic.types,
ef.level = "prey"
)
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. |
metabolic.types |
A vector containing information on species type ( |
ef.level |
Set to |
Value
A matrix representing the Jacobian matrix of the dynamical system associated with the fluxes from the fluxing
function.
Author(s)
Benoit Gauzens, benoit.gauzens@gmail.com
Examples
# First compute species per unit biomass metabolic rates:
losses = 0.15 * groups.level$bodymasses^(-0.25)
val.mat = fluxing(groups.level$mat,
groups.level$biomasses,
losses,
groups.level$efficiencies,
bioms.pref = TRUE,
ef.level = "prey")
# define metabolic types
met.types = rep("animal", nrow(val.mat))
met.types[groups.level$efficiencies == 0.545] = "plant"
met.types[groups.level$efficiencies == 0.158] = "detritus"
create.jacob(val.mat,
groups.level$biomasses,
groups.level$efficiencies,
met.types,
ef.level = "prey")
[Package fluxweb version 2.0.1 Index]