modelHD | R Documentation |
This function fits and compares (optional) height-diameter models.
modelHD(D, H, method = NULL, useWeight = FALSE, drawGraph = FALSE, plot = NULL)
D |
Vector with diameter measurements (in cm). NA values are accepted but a minimum of 10 valid entries (i.e. having a corresponding height in H) is required. |
H |
Vector with total height measurements (in m). NA values are accepted but a minimum of 10 valid entries (i.e. having a corresponding diameter in D) is required. |
method |
Method used to fit the relationship. To be chosen between:
If |
useWeight |
If weight is |
drawGraph |
If |
plot |
(optional) Plot ID, must be either one value, or a vector of the same length as D. This argument is used to build stand-specific HD models. |
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).
If plot
is NULL or has a single value, a single list is returned. If there is more than one plot,
multiple embedded lists are returned with plots as the list names.
If model
is not null (model comparison), returns a list :
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
fitPlot
: a ggplot object containing the model fitting plot
If the parameter model is null, the function return a plot with all the methods for comparison, the function also returns a data.frame with:
method
: The method that had been used to construct the plot
RSE
: Residual Standard Error of the model
RSElog
: Residual Standard Error of the log model (NULL
if other model)
Average_bias
: The average bias for the model
Maxime REJOU-MECHAIN, Arthur PERE, Ariane TANGUY, Arthur Bailly
retrieveH()
# Load a data set
data(NouraguesHD)
# Fit H-D models for the Nouragues dataset
HDmodel <- modelHD(D = NouraguesHD$D, H = NouraguesHD$H, drawGraph = TRUE)
# For a chosen model
HDmodel <- modelHD(D = NouraguesHD$D, H = NouraguesHD$H,
method = "log2", drawGraph = TRUE)
# Using weights
HDmodel <- modelHD(
D = NouraguesHD$D, H = NouraguesHD$H,
method = "log2", useWeight = TRUE,
drawGraph = TRUE)
# With multiple stands (plots)
HDmodel <- modelHD(
D = NouraguesHD$D, H = NouraguesHD$H,
method = "log2", useWeight = TRUE,
plot = NouraguesHD$plotId, drawGraph = TRUE)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.