model_height: Fit a height-diameter model.
In forestgeo/fgeo.biomass: Calculate Biomass with the Best Available Allometries

Description Usage Arguments Details Value See Also Examples

This functions wraps BIOMASS::modelHD(), which you can see for more options and details.

1	model_height(data, method = c("log1", "log2", "weibull", "michaelis"))

data

A ForestGEO-like census-dataframe.

method

A character string. One of these:

log 1 (equivalent to a power model): log(H) = a + b * log(D)
log 2: log(H) = a + b * log(D) + c * log(D)^2
weibull: H = a * (1 - exp(-(D / b) ^ c))
michaelis: H = (A * D) / (B + D)

All the back transformations for log-log models are done using the Baskerville correction (0.5 * RSE^2, where RSE is the Residual Standard Error).

A list with these elements:

input: list of the data used to construct the model (list(H, D)).
model: outputs of the model (same outputs as given by stats::lm(), stats::nls()).
RSE: Residual Standard Error of the model.
RSElog: Residual Standard Error of the log model (NULL if other model).
residuals: Residuals of the model.
coefficients: Coefficients of the model.
R.squared: R^2 of the model.
formula: Formula of the model.
method: Name of the method used to construct the model.
predicted: Predicted height values.

BIOMASS::modelHD().

data <- fgeo.biomass::scbi_tree1[1:30, ]

str(model_height(data))

out <- model_height(data, method = "log2")
# Avoid too much output
str(out[-2])
# Glimpse what wasn't shown above
lapply(out[2], names)