soilfoodwebs (version 1.0.2)

corrstoich: Correct stoichiometry

Description

Correct stoichiometry

Usage

corrstoich(
  usin,
  forceProd = FALSE,
  dietlimits = c(NA),
  Immobilizationlimit = Inf
)

Value

The modified community with new production efficiencies and diets.

Arguments

usin

# Community in which to correct stoichiometry.

forceProd

Boolean: Should we force organisms to only change their production efficiency and not their diet?

dietlimits

A matrix the same size as imat that gives the diet limits as a proportion of the total diet. All values must be between 0 and 1. Leaving it as NA sets the limits of all diet items to 1.

Immobilizationlimit

Set a limit for the maximum rate of immobilization of inorganic nitrogen for the food web per time step.

Details

A function that corrects trophic species production efficiency and/or diet composition to balance carbon and nitrogen demands.

Examples

corrstoich(intro_comm)
# To force the correction to modify production efficiency only.

corrstoich(intro_comm, forceProd = TRUE)

# Set limits on the composition of the animal diet
DL = intro_comm$imat # copy over the feeding matrix to use for the diet limits.
# Oribatid 1 can only have up to 20% of its diet from each fungal species
DL["Orib1",] = c(0,0,0,0.2,0.2,1,1)
# Oribatid 2 can only have up to 10% of its diet from fungi 1
DL["Orib2",] = c(0,0,0,0.1,1,1,1)

# Run them with the limits:
corrstoich(intro_comm, dietlimits = DL)