Description Usage Arguments Value Examples
Estimation updates replicate production values based on a trophic_dynamics object
1 2 3 4 5 6 7 8 9 10 11 | 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(), ...)
|
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 |
An object of class production_estimates
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 | 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)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.