The hardware and bandwidth for this mirror is donated by METANET, the Webhosting and Full Service-Cloud Provider.
If you wish to report a bug, or if you are interested in having us mirror your free-software or open-source project, please feel free to contact us at mirror[@]metanet.ch.
This introduction shows the list of biodiversity patterns that can be
derived using mobsim
. In this vignette the patterns are
evaluated on a simulated data set, but of course the same patterns can
be also derived from real data.
First, we generate community data by simulating 1,000 individuals
from 30 species with a cluster size of sigma = 0.01
and one
cluster per species (mother_points = 1
). For more
information on the simulation of communities see the vignette
Simulating communities with mobsim
.
library(mobsim)
sim_dat1 <- sim_thomas_community(s_pool = 30, n_sim = 1000, sad_type = "lnorm",
sad_coef = list("meanlog" = 2, "sdlog" = 1),
sigma = 0.1, mother_points = 1)
Then we explore the generated community object. In the plot each dot represents one individual and the colour indicates the species identity.
## No. of individuals: 1000
## No. of species: 30
## x-extent: 0 1
## y-extent: 0 1
##
## x y species
## Min. :0.0017 Min. :0.00079 species_01:217
## 1st Qu.:0.5826 1st Qu.:0.38876 species_02:184
## Median :0.7224 Median :0.60217 species_03:137
## Mean :0.6797 Mean :0.58890 species_04: 57
## 3rd Qu.:0.8326 3rd Qu.:0.84126 species_05: 55
## Max. :1.0000 Max. :0.99942 species_06: 46
## (Other) :304
## List of 3
## $ census :'data.frame': 1000 obs. of 3 variables:
## ..$ x : num [1:1000] 0.914 0.822 0.869 0.824 0.799 ...
## ..$ y : num [1:1000] 0.859 0.861 0.98 0.933 0.794 ...
## ..$ species: Factor w/ 30 levels "species_01","species_02",..: 1 1 1 1 1 1 1 1 1 1 ...
## $ x_min_max: num [1:2] 0 1
## $ y_min_max: num [1:2] 0 1
## - attr(*, "class")= chr "community"
A fundamental non-spatial pattern of a community is the abundance
distribution, i.e. the distribution of commonness and rarity in a
community. The abundance distribution can be extracted from a community
object using the function communtiy_to_sad()
A standard plot in community ecology is the rank-abundance plot, where the abundance of each species is plotted vs. its rank from highest to lowest abundance. As a standard log-scaling is used for the abundance axis in this plot.
Of course the abundance distribution can be also visualized as a histogram. By tradition for the binning of species abundances logarithms with base 2 are used, following the suggestion of Preston (1948). This means the first abundance class includes species with just one individual, the second class with two individuals, the third class with 3-4, the fourth with 5-8 etc. These abundance classes are called “octaves”.
Another important biodiversity pattern is the rarefaction curve, which estimates how the number of observed species increases with sample size of individuals. The rarefaction curve assumes that individuals are sampled randomly and independently (Gotelli & Colwell 2001). The rarefaction curve only depends on the abundance distribution of species.
## Classes 'spec_sample_curve' and 'data.frame': 1000 obs. of 2 variables:
## $ n : int 1 2 3 4 5 6 7 8 9 10 ...
## $ spec_rarefied: num 1 1.89 2.68 3.39 4.04 ...
The main strengths of mobsim
are the simulation and
analysis of spatial biodiversity patterns. In the following all
functions for spatial pattern evaluation included in mobsim
are introduced.
Closely related to the rarefaction curve is the spatial species accumulation curve. In contrast to the rarefaction curve for the derivation of the species accumulation curve individuals are not sampled randomly, but starting from a focal individual always the closest neighbour individual is sampled and the number of encountered species is counted. The resulting curve is derived as the average of the curves for each focal individual. The species accumulation curve is influenced by the abundance distribution, but also by the spatial distribution of individuals and species. Therefore the species accumulation curves requires individual’s positions in addition to species abundances.
Due to the close relationship of the accumulation and rarefaction curves they are calculated and plotted with the same function.
It is comprehensive to plot the rarefaction and the species accumulation curve together. The difference between the two curves indicates aggregation or overdispersion of conspecific individuals.
spec_curves1 <- spec_sample_curve(sim_dat1, method = c("accumulation", "rarefaction"))
plot(spec_curves1)
The most well-known spatial biodiversity pattern is the
species-area relationship (SAR). In mobsim
the
function divar
(diversity-area relationships) calculates
the species richness in randomly located subplots (quadrats or
rectangles) of different sizes. However, the function divar
calculates additional indices, including the Shannon and Simpson
diversity indices for each subplot, as well as the number of endemic
species, which are the species that only occur within, but not
outside the subplot. See ?div_rect
for detailed information
on the diversity indices.
For the Shannon and Simpson diversity indices also the corresponding Effective Number of Species (ENS) is calculated (Jost 2006). This measure corresponds to the species number in a community of equally abundant species, which results in the same Shannon and Simpson indices as the observed community (with unequal abundances).
The endemics-area relationship (EAR) has been suggested as important tool to investigate the consequences of habitat loss for biodiversity (He & Hubbell 2011, Keil et al. 2015). For the evaluation of the diversity-area relationships a vector with subplot sizes measured as proportion of the total community size has to be defined.
subplot_size <- c(0.01, 0.05, 0.1, 0.2, 0.3, 0.4, 0.5,
0.6, 0.7, 0.8, 0.9, 0.95, 0.99, 1)
divar1 <- divar(sim_dat1, prop_area = subplot_size)
head(divar1)
## prop_area m_species sd_species m_endemics sd_endemics m_shannon sd_shannon
## 1 0.01 3.236559 1.527984 0.01075269 0.1036952 0.8476667 0.4793897
## 2 0.05 6.270000 2.960532 0.13000000 0.3379977 1.1884362 0.4647952
## 3 0.10 9.630000 3.807635 0.24000000 0.4292347 1.4896940 0.3931594
## 4 0.20 15.190000 2.901393 0.84000000 0.9610600 1.8984476 0.1649661
## 5 0.30 18.810000 2.236949 1.72000000 1.5832456 2.1193808 0.1468489
## 6 0.40 22.990000 1.817396 3.77000000 1.9114488 2.2311739 0.0557845
## m_ens_shannon sd_ens_shannon m_simpson sd_simpson m_ens_simpson
## 1 2.587596 1.0892155 0.5097275 0.21666918 2.396748
## 2 3.591544 1.3352842 0.5690234 0.21363593 2.814330
## 3 4.737762 1.5380408 0.6569855 0.15664511 3.454018
## 4 6.768003 1.1572716 0.7672961 0.05792136 4.592310
## 5 8.417319 1.2792297 0.8048639 0.04562801 5.460863
## 6 9.325188 0.5238804 0.8180267 0.02089588 5.570515
## sd_ens_simpson
## 1 0.8927158
## 2 1.1167927
## 3 1.3461093
## 4 1.2639255
## 5 1.5076464
## 6 0.6749775
The last spatial pattern provided by mobsim
is the
distance decay of community similarity, which quantifies how
quickly the similarity in species composition between two subplots
decreases with the distance between two subplots. The function
dist_decay
distributes non-overlapping subplots with
user-defined size and number in the community and calculates all
pairwise similarities.
The function dist_decay
makes use of the function
vegdist
from the package vegan. Therefore you
can consult ?vegdist
for a list of available similarity
indices.
Here is a demonstration how the distance decay can be estimated and visualized
## distance similarity
## 114 0.1418358 0.7000000
## 16 0.1554387 0.5000000
## 45 0.1585154 0.0000000
## 99 0.1641721 0.0000000
## 76 0.1654742 0.2222222
## 17 0.1688944 0.5405405
F. W. Preston 1948. The Commonness, and Rarity, of Species. Ecology 29:254-283.
Gotelli & Colwell 2001. Quantifying biodiversity: procedures and pitfalls in the measurement and comparison of species richness. Ecology Letters 4: 379–391.
He & Hubbell 2001. Species-area relationships always overestimate extinction rates from habitat loss. Nature 473:368–371
Jost 2006. Entropy and diversity. Oikos, 113:363–375.
Keil et al. 2015. On the decline of biodiversity due to area loss. Nature communications 6.
Morlon et al. 2008. A general framework for the distance–decay of similarity in ecological communities. Ecology Letters 9:904–917.
These binaries (installable software) and packages are in development.
They may not be fully stable and should be used with caution. We make no claims about them.