Estimate PD (Faith 1992) or FD (Petchey and Gaston 2002) from abundance or probability data and a phylogenetic or functional dendrogram.
Usage
div_pd(x, tree, ...)
# S3 method for class 'numeric'
div_pd(x, tree, prune = FALSE, as_numeric = FALSE, ..., check_arguments = TRUE)
# S3 method for class 'species_distribution'
div_pd(
x,
tree,
prune = FALSE,
gamma = FALSE,
as_numeric = FALSE,
...,
check_arguments = TRUE
)
Arguments
- x
An object, that may be a named numeric vector (names are species names) containing abundances or probabilities, or an object of class abundances or probabilities.
- tree
an ultrametric, phylogenetic tree. May be an object of class phylo_divent, ape::phylo, ade4::phylog or stats::hclust.
- ...
Unused.
- prune
What to do when some species are in the tree but not in
x
? IfTRUE
, the tree is pruned to keep species ofx
only. The height of the tree may be changed if a pruned branch is related to the root. IfFALSE
(default), the length of branches of missing species is not summed but the height of the tree is never changed.- as_numeric
if
TRUE
, a number or a numeric vector is returned rather than a tibble.- check_arguments
if
TRUE
, the function arguments are verified. Should be set toFALSE
to save time when the arguments have been checked elsewhere.- gamma
if
TRUE
, \(\gamma\) diversity, i.e. diversity of the metacommunity, is computed.
Details
Estimators to deal with incomplete sampling are not implemented.
Use function div_hill with argument q = 0
if they are needed.
PD and FD are defined as the total length of the branches of the tree.
All species of the species_distribution
must be found in the tips of the
tree
.
References
Faith DP (1992).
“Conservation Evaluation and Phylogenetic Diversity.”
Biological Conservation, 61(1), 1–10.
doi:10.1016/0006-3207(92)91201-3
.
Petchey OL, Gaston KJ (2002).
“Functional Diversity (FD), Species Richness and Community Composition.”
Ecology Letters, 5, 402–411.
doi:10.1046/j.1461-0248.2002.00339.x
.
Examples
# diversity of each community
div_pd(paracou_6_abd, tree = paracou_6_taxo)
#> # A tibble: 4 × 5
#> site weight estimator order diversity
#> <chr> <dbl> <chr> <dbl> <dbl>
#> 1 subplot_1 1.56 naive 0 337
#> 2 subplot_2 1.56 naive 0 373
#> 3 subplot_3 1.56 naive 0 350
#> 4 subplot_4 1.56 naive 0 339
# gamma diversity
div_pd(paracou_6_abd, tree = paracou_6_taxo, gamma = TRUE)
#> # A tibble: 1 × 4
#> site estimator order diversity
#> <chr> <chr> <dbl> <dbl>
#> 1 Metacommunity naive 0 555