propagate_errors: Propagate errors in above-ground biomass (AGB).
In forestgeo/fgeo.biomass: Calculate Biomass with the Best Available Allometries

Description Usage Arguments Details Value References Examples

This functions is warps BIOMASS::AGBmonteCarlo(), which you may use directly for more options.

1	propagate_errors(data, n = 1000, dbh_sd = NULL, height_model = NULL)

`data`	The output of `add_tropical_biomass()`.
`n`	Number of iterations. Cannot be smaller than 50 or larger than 1000. By default `n = 1000`
`dbh_sd`	This variable can take three kind of values, indicating how to propagate the errors on diameter measurements: a single numerical value or a vector of the same size as D, both representing the standard deviation associated with the diameter measurements or "chave2004" (an important error on 5 percent of the measures, a smaller error on 95 percent of the trees).
`height_model`	Model used to estimate tree height from tree diameter (output from `model_height()`, see example).

See Rejou-Mechain et al. (2017) for all details on the error propagation procedure.

A list with the following elements:

meanAGB: Mean stand AGB value following the error propagation.
medAGB: Median stand AGB value following the error propagation.
sdAGB: Standard deviation of the stand AGB value following the error propagation.
credibilityAGB: Credibility interval at 95 following the error propagation.
AGB_simu: Matrix with the AGB of the trees (rows) times the n iterations (columns).

Chave, J. et al. (2004). Error propagation and scaling for tropical forest biomass estimates. Philosophical Transactions of the Royal Society B: Biological Sciences, 359(1443), 409-420.

Rejou-Mechain et al. (2017). BIOMASS: An R Package for estimating above-ground biomass and its uncertainty in tropical forests. Methods in Ecology and Evolution, 8 (9), 1163-1167.

library(dplyr)

data <- fgeo.biomass::scbi_tree1 %>%
  slice(1:100)
species <- fgeo.biomass::scbi_species

# Using `region` (default)
biomass <- add_tropical_biomass(data, species)
model <- model_height(biomass)
str(
  propagate_errors(biomass, n = 50, height_model = model)
)

# Using `latitude` and `longitude`
biomass <- add_tropical_biomass(
  data = data,
  species = species,
  latitude = 4,
  longitude = -52
)

model <- model_height(biomass)

# Asks to confirm using the model instead of coordinates
if (interactive()) {
  str(
    propagate_errors(biomass, n = 50, height_model = model)
  )
}