INFCvpeSUM {ForIT} | R Documentation |
Estimate tree bole volume or phytomass for stems groups, with associated accuracy info
Description
Cumulative estimation of the volume or phytomass of groups of trees is just
the summation of the values computed with INFCvpe()
, but the computation of
accuracy estimates is improved using these summation functions.
Two approaches are available.
Via
INFCvpe_summarise()
that processes and returns a data frameor by following aggregation functions within a standard
summarise()
:-
INFCvpe_sum()
-
INFCvpe_ConfInt()
-
INFCvpe_OutOfDomain()
-
Usage
INFCvpe_summarise(
in.data,
EPPOcode_C,
dbh_C,
h_tot_C,
quantity = "vol",
p = 0.95
)
INFCvpe_sum(EPPOcode, dbh, h_tot, quantity = "vol")
INFCvpe_ConfInt(EPPOcode, dbh, h_tot, quantity = "vol", p = 0.95)
INFCvpe_OutOfDomain(EPPOcode, dbh, h_tot)
Arguments
in.data |
A dataframe (or tibble) containing tally data to be matched with "EPPOcode_C", "dbh_C" and "htot_C" arguments |
EPPOcode_C |
A string, the name of the column in |
dbh_C |
A string, the name of the column in |
h_tot_C |
A string, the name of the column in |
quantity |
(default =
|
p |
(default OPZIONE 2 |
EPPOcode |
A character vector with the species EPPO codes (with length = 1 or length = length(dbh)) |
dbh |
A numeric vector with the brest height diameter values |
h_tot |
A numeric vector with the tree total height values (with length = 1 or length = length(dbh)) |
Details
Functions developed following Tabacchi et al. (2011), pages 23-26.
Value
-
INFCvpe_summarise()
returns a dataframe (tibble) with the grouping columns defined withgroup_by()
, and the following columns:-
quantity
: as additional grouping column, -
n
: number of trees in the group, -
n_out
: the number of(dbh, htot)
pairs that are 'out of the domain', -
est
: the estimated value, -
cihw
: confidence interval half width -
p
: probability used computingcihw
-
-
INFCvpe_SUM
- the functions of this family return a numeric vector, aggregating rows within the same group,-
INFCvpe_sum()
returns the sum of the estimated quantities, -
INFCvpe_ConfInt()
returns 'confidence interval half width', -
INFCvpe_OutOfDomain()
returns the number of 'out of domain'(dhb, h_tot)
pairs included in the summation
-
See Also
INFCvpe()
to compute individual estimates, with detailed accuracy evaluation
Examples
## Not run:
Sezione <- function(EPPOcodes){
# retrive 'Sezione' name, decoding EPPO codes
INFCspecies %>%
dplyr::filter(EPPOcode %in% EPPOcodes) %>%
dplyr::left_join(INFCcatalog,by = "pag")%>%
dplyr::select(section) %>%
purrr::pluck(1)
}
tst <- ForIT_test_data %>%
dplyr::filter(UC != "U0")
# select Tabachi et al. example data
tst %>%
dplyr::group_by(specie) %>%
INFCvpe_summarise("specie", "d130", "h_dendro") %>%
dplyr::ungroup() %>%
dplyr::mutate(specie = Sezione(specie),
dplyr::across(c("est", "cihw"), ~round(.x, 1))
) %>%
dplyr::select(specie, est, cihw) %>%
dplyr::arrange(specie) %>%
dplyr::slice(2, 1, 3) %>%
t() %>%
provideDimnames(base = list(dimnames(.)[[1]], ""), unique=FALSE)
# Compare ForIT (ver 2) output
## specie "Aceri" "Abete bianco" "Faggio"
## est "4623.0" "4044.2" "1079.4"
## cihw "567.5" "661.2" "275.4"
# with 'Tabella 2' in Tabacchi et al. (2011, pag. 27)
## specie "aceri" "abete bianco" "faggio"
## est "4623.0" "4044.2" "1079.4"
## cihw "567.4" "662.4" "279.2"
# Using 'INFCvpe_summarise()'
## Overall totals
tst %>%
INFCvpe_summarise("specie", "d130", "h_dendro", quantity = c("vol", "dw4"))
## Group by dbh class ('cld')
tst %>%
dplyr::mutate(cld = ceiling(d130/5)*5) %>%
dplyr::group_by(UC, specie, cld) %>%
INFCvpe_summarise("specie", "d130", "h_dendro")
## Group by sampling unit ('UC')
tst %>%
dplyr::group_by(UC) %>%
INFCvpe_summarise("specie", "d130", "h_dendro", quantity = "dw4")
# Using 'INFCvpeSUM' aggregation functions
## Esitmate 'dw4' phytomass, by sampling unit ('UC')
tst %>%
dplyr::group_by(UC) %>%
dplyr::summarise(
n_stems = dplyr::n(),
OoD = INFCvpe_OutOfDomain(specie, d130, h_dendro),
dw4 = INFCvpe_sum(specie, d130, h_dendro, quantity = "dw4"),
dw4_ConfInt = INFCvpe_ConfInt(specie, d130, h_dendro, quantity = "dw4")
)
## Esitmate volume, by sampling unit ('UC')
tst %>%
dplyr::group_by(UC) %>%
dplyr::summarise(
n_stems = dplyr::n(),
OoD = INFCvpe_OutOfDomain(specie, d130, h_dendro),
vol = INFCvpe_sum(specie, d130, h_dendro),
vol_ConfInt = INFCvpe_ConfInt(specie, d130, h_dendro)
)
rm(tst, Sezione)
## End(Not run)