production_estimates: Estimate production for each node in a food web

In jdyen/trophic: Trophic Models of Energetics to Estimate Carrying Capacity

Description

Estimation updates replicate production values based on a trophic_dynamics object

Usage

estimate_production(trophic_dynamics, primary_producers,
  stochastic = c("food_web", "efficiency", "primary_production"),
  nsim = 100)

is.production_estimates(x)

## S3 method for class 'production_estimates'
print(x, ...)

## S3 method for class 'production_estimates'
plot(x, nodes = NULL, settings = list(), ...)

Arguments

trophic_dynamics	a dynamics object - food web, trophic efficiency, trophic dominance, and P:B ratios
primary_producers	a primary_producers object containing estimates of primary production
stochastic	a character vector naming the stochastic elements of the trophic projection; can be any combination of "food_web", "efficiency", and "primary_production"
nsim	number of estimates to generate
x	a production_estimates object
...	further arguments passed to or from other methods
nodes	integer or integer vector of nodes to plot (indexed by food_web row)
settings	plot settings passed directly to plot

Value

An object of class production_estimates

Examples

library(trophic)
library(future)
plan(multiprocess)

# Construct the component objects
test_fw <- build_food_web(interaction_matrix = food_web)
test_efficiency_matrix <- build_efficiency_matrix(efficiency_mean = efficiency_mean,
                                                  efficiency_sd = 0.01)
test_dominance <- build_dominance_matrix(dominance = dominance_matrix)
test_primary_producers <- build_primary_producers(production_mean = c(1, 2),
                                                  production_sd = c(0.5, 0.5))

# Construct the trophic_dynamics object
test_trophic_dynamics <- build_trophic_dynamics(food_web = test_fw,
                                                efficiency_matrix = test_efficiency_matrix,
                                                dominance_matrix = test_dominance)

# Estimate production values from constructed trophic_dynamics object
production_estimates <- estimate_production(test_trophic_dynamics, test_primary_producers)

# Test if object is of the type 'production_estimates'
## Not run: 
is.production_estimates(x)

## End(Not run)

# Print information about the 'production_estimates' object
## Not run: 
print(x)

## End(Not run)

# Plot a 'production_estimates' object

## Not run: 
plot(x)

## End(Not run)

jdyen/trophic documentation built on May 15, 2019, 3:19 p.m.