summaryByPlot {BIOMASS} | R Documentation |
Summarise by plot the posterior distribution of AGB values
Description
This function summarizes the matrix AGB_val
given by the function AGBmonteCarlo()
by plot.
Usage
summaryByPlot(AGB_val, plot, drawPlot = FALSE)
Arguments
AGB_val |
Matrix resulting from the |
plot |
Vector corresponding to the plots code (plots ID) |
drawPlot |
A logic indicating whether the graphic should be displayed or not |
Details
If some trees belong to an unknown plot (i.e. NA value in the plot arguments), their AGB values are randomly assigned to a plot at each iteration of the AGB monte Carlo approach.
Value
a data frame where:
-
plot
: the code of the plot -
AGB
: AGB value at the plot level -
Cred_2.5
: the 2.5\ -
Cred_97.5
: the 97.5\
Examples
# Load a database
data(NouraguesHD)
data(NouraguesTrees)
# Modelling height-diameter relationship
HDmodel <- modelHD(D = NouraguesHD$D, H = NouraguesHD$H, method = "log2")
# Retrieving wood density values
NouraguesWD <- getWoodDensity(NouraguesTrees$Genus, NouraguesTrees$Species,
stand = NouraguesTrees$plotId)
# Propagating errors
resultMC <- AGBmonteCarlo(
D = NouraguesTrees$D, WD = NouraguesWD$meanWD,
errWD = NouraguesWD$sdWD, HDmodel = HDmodel )
# The summary by plot
summaryByPlot(AGB_val = resultMC$AGB_simu, plot = NouraguesTrees$Plot)
[Package BIOMASS version 2.2.4 Index]