Estimate Steady State Biomass

Usage

steady_state(x = NULL, pars = NULL, model = NULL)

Arguments

x: A litfit object. If provided, pars and model parameters are extracted from this object.
pars: A numeric vector of parameters for the model. Only needed if x is not provided.
model: A character string specifying the decomposition model. Must be one of the following: "neg.exp", "weibull", "discrete.parallel", "discrete.series", or "cont.quality2". Only needed if x is not provided.

Value

A named numeric value representing the estimated steady state biomass from the specified model.

Details

Computes the steady state biomass, as a proportion of the annual input, based on a given model fit or parameters.

Currently, the function supports a subset of decomposition models. New model support is planned for future updates.

Author

Will Cornwell

Examples

# Example with litfit object
fit <- fit_litter(
  time = c(0, 1, 2, 3, 4, 5, 6),
  mass.remaining = c(1, 0.9, 1.01, 0.4, 0.6, 0.2, 0.01),
  model = "neg.exp",
  iters = 250
)
#> Number of successful fits:  248  out of 250 
steady_state(fit)
#> neg.exp 
#> 4.60503 

# Example with specific model and parameter values
steady_state(pars = c(6, 2), model = "weibull")
#>  weibull 
#> 5.317362