Diversity Estimation of a metacommunity
DivEst.Rd
Estimates diversity of a metacommunity.
Usage
DivEst(q = 0, MC, Biased = TRUE, Correction = "Best", Tree = NULL,
Normalize = TRUE, Z = NULL, Simulations = 100,
ShowProgressBar = TRUE, CheckArguments = TRUE)
is.DivEst(x)
# S3 method for DivEst
plot(x, ..., main = NULL, Which = "All",
Quantiles = c(0.025, 0.975), colValue = "red", lwdValue = 2, ltyValue = 2,
colQuantiles = "black", lwdQuantiles = 1, ltyQuantiles = 2)
# S3 method for DivEst
autoplot(object, ..., main = NULL, Which = "All",
labels = NULL, font.label = list(size=11, face="plain"),
Quantiles = c(0.025, 0.975), colValue = "red",
colQuantiles = "black", ltyQuantiles = 2)
# S3 method for DivEst
summary(object, ...)
Arguments
- q
A number: the order of diversity.
- MC
A
MetaCommunity
object.- Biased
Logical; if
FALSE
, a bias correction is appplied.- Correction
A string containing one of the possible corrections. The correction must be accepted by
DivPart
."Best"
is the default value.- Tree
An object of class
hclust
,phylo
,phylog
orPPtree
. The tree must be ultrametric.- Normalize
If
TRUE
(default), diversity is not affected by the height of the tree..
IfFALSE
, diversity is proportional to the height of the tree.- Z
A relatedness matrix, i.e. a square matrix whose terms are all positive, strictly positive on the diagonal. Generally, the matrix is a similarity matrix, i.e. the diagonal terms equal 1 and other terms are between 0 and 1.
- Simulations
The number of simulations to build confidence intervals.
- ShowProgressBar
If
TRUE
(default), a progress bar is shown.- CheckArguments
Logical; if
TRUE
, the function arguments are verified. Should be set toFALSE
to save time when the arguments have been checked elsewhere.- x
An object to be tested or plotted.
- main
The title of the plot.
- Which
May be
"Alpha"
,"Beta"
or"Gamma"
to respectively plot the metacommunity's alpha, beta or gamma diversity. If"All"
(default), all three plots are shown.- labels
Vector of labels to be added to multiple plots.
"auto"
is the same asc("a", "b", "c", "d)"
.- font.label
A list of arguments to customize labels. See
ggarrange
.- object
A
MCdiversity
object to be summarized or plotted.- Quantiles
A vector containing the quantiles of interest.
- colValue
The color of the line representing the real value on the plot.
- lwdValue
The width of the line representing the real value on the plot.
- ltyValue
The line type of the line representing the real value on the plot.
- colQuantiles
The color of the lines representing the quantiles on the plot.
- lwdQuantiles
The width of the lines representing the quantiles on the plot.
- ltyQuantiles
The line type of the lines representing the quantiles on the plot.
- ...
Additional arguments to be passed to the generic methods.
Details
Divest
estimates the diversity of the metacommunity and partitions it into alpha and beta components.
If Tree
is provided, the phylogenetic diversity is calculated else if Z
is not NULL
, then similarity-based entropy is calculated.
Bootstrap confidence intervals are calculated by drawing simulated communities from a multinomial distribution following the observed frequencies (Marcon et al, 2012; 2014).
Value
A Divest
object which is a DivPart
object with an additional item in its list:
- SimulatedDiversity
A matrix containing the simulated values of alpha, beta and gamma diversity.
Divest
objects can be summarized and plotted.
References
Marcon, E., Herault, B., Baraloto, C. and Lang, G. (2012). The Decomposition of Shannon's Entropy and a Confidence Interval for Beta Diversity. Oikos 121(4): 516-522.
Marcon, E., Scotti, I., Herault, B., Rossi, V. and Lang, G. (2014). Generalization of the partitioning of Shannon diversity. PLOS One 9(3): e90289.
Marcon, E., Herault, B. (2015). Decomposing Phylodiversity. Methods in Ecology and Evolution 6(3): 333-339.
Examples
# Load Paracou data (number of trees per species in two 1-ha plot of a tropical forest)
data(Paracou618)
# Estimate Shannon diversity.
Estimation <- DivEst(q = 1, Paracou618.MC, Biased = FALSE, Correction = "UnveilJ",
Simulations = 20)
#>
plot(Estimation)
#> Warning: NAs introduced by coercion
#> Warning: NAs introduced by coercion
#> Warning: NAs introduced by coercion
#> Warning: NAs introduced by coercion
#> Warning: NAs introduced by coercion
#> Warning: NAs introduced by coercion
summary(Estimation)
#> Diversity partitioning of order 1 of MetaCommunity Paracou618.MC
#> with correction: UnveilJ
#> Alpha diversity of communities:
#> P006 P018
#> 83.7268 118.2713
#> Total alpha diversity of the communities:
#> [1] 97.06467
#> Beta diversity of the communities:
#> UnveilJ
#> 1.422843
#> Gamma diversity of the metacommunity:
#> UnveilJ
#> 138.1078
#> Quantiles of simulations (alpha, beta and gamma diversity):
#> 0% 1% 2.5% 5% 10% 25% 50% 75%
#> 89.64103 90.10350 90.79720 91.95337 92.86043 93.85074 96.91485 98.49928
#> 90% 95% 97.5% 99% 100%
#> 104.11466 104.28565 105.25133 105.83074 106.21701
#> 0% 1% 2.5% 5% 10% 25% 50% 75%
#> 1.347289 1.354156 1.364457 1.381624 1.384839 1.407170 1.433999 1.442474
#> 90% 95% 97.5% 99% 100%
#> 1.450413 1.459604 1.466567 1.470745 1.473531
#> 0% 1% 2.5% 5% 10% 25% 50% 75%
#> 129.3010 129.7963 130.5392 131.7774 133.6159 134.6326 138.3291 141.3829
#> 90% 95% 97.5% 99% 100%
#> 143.4091 146.1764 146.4703 146.6466 146.7641