Title: | Biodiversity Index Calculation and Bootstrap Confidence Interval Estimation |
---|---|
Description: | Provides tools for the calculation of common biodiversity indices from count data. Additionally, it incorporates bootstrapping techniques to generate multiple samples, facilitating the estimation of confidence intervals around these indices. Furthermore, the package allows for the exploration of how variation in these indices changes with differing numbers of sites, making it a useful tool with which to begin an ecological analysis. Methods are based on the following references: Chao et al. (2014) <doi:10.1890/13-0133.1>, Chao and Colwell (2022) <doi:10.1002/9781119902911.ch2>, Hsieh, Ma,` and Chao (2016) <doi:10.1111/2041-210X.12613>. |
Authors: | Craig Eric Simpkins [aut, cre]
|
Maintainer: | Craig Eric Simpkins <[email protected]> |
License: | MIT + file LICENSE |
Version: | 1.0.4 |
Built: | 2025-03-10 03:18:14 UTC |
Source: | https://github.com/csim063/biosampler |
This dataset contains tree counts from Barro-Colorado Island. It has 50 rows each representing the counts taken from a separate one hectare plot for each of the 225 species (columns)
BCI
BCI
A data frame with 50 rows and 225 columns:
Count for Abarema.macradenia
Count for Vachellia.melanoceras
Count for Acalypha.diversifolia
Count for Acalypha.macrostachya
Count for Adelia.triloba
Count for Aegiphila.panamensis
Count for Alchornea.costaricensis
Count for Alchornea.latifolia
Count for Alibertia.edulis
Count for Allophylus.psilospermus
Count for Alseis.blackiana
Count for Amaioua.corymbosa
Count for Anacardium.excelsum
Count for Andira.inermis
Count for Annona.spraguei
Count for Apeiba.glabra
Count for Apeiba.tibourbou
Count for Aspidosperma.desmanthum
Count for Astrocaryum.standleyanum
Count for Astronium.graveolens
Count for Attalea.butyracea
Count for Banara.guianensis
Count for Beilschmiedia.pendula
Count for Brosimum.alicastrum
Count for Brosimum.guianense
Count for Calophyllum.longifolium
Count for Casearia.aculeata
Count for Casearia.arborea
Count for Casearia.commersoniana
Count for Casearia.guianensis
Count for Casearia.sylvestris
Count for Cassipourea.guianensis
Count for Cavanillesia.platanifolia
Count for Cecropia.insignis
Count for Cecropia.obtusifolia
Count for Cedrela.odorata
Count for Ceiba.pentandra
Count for Celtis.schippii
Count for Cespedesia.spathulata
Count for Chamguava.schippii
Count for Chimarrhis.parviflora
Count for Maclura.tinctoria
Count for Chrysochlamys.eclipes
Count for Chrysophyllum.argenteum
Count for Chrysophyllum.cainito
Count for Coccoloba.coronata
Count for Coccoloba.manzinellensis
Count for Colubrina.glandulosa
Count for Cordia.alliodora
Count for Cordia.bicolor
Count for Cordia.lasiocalyx
Count for Coussarea.curvigemma
Count for Croton.billbergianus
Count for Cupania.cinerea
Count for Cupania.latifolia
Count for Cupania.rufescens
Count for Cupania.seemannii
Count for Dendropanax.arboreus
Count for Desmopsis.panamensis
Count for Diospyros.artanthifolia
Count for Dipteryx.oleifera
Count for Drypetes.standleyi
Count for Elaeis.oleifera
Count for Enterolobium.schomburgkii
Count for Erythrina.costaricensis
Count for Erythroxylum.macrophyllum
Count for Eugenia.florida
Count for Eugenia.galalonensis
Count for Eugenia.nesiotica
Count for Eugenia.oerstediana
Count for Faramea.occidentalis
Count for Ficus.colubrinae
Count for Ficus.costaricana
Count for Ficus.insipida
Count for Ficus.maxima
Count for Ficus.obtusifolia
Count for Ficus.popenoei
Count for Ficus.tonduzii
Count for Ficus.trigonata
Count for Ficus.yoponensis
Count for Garcinia.intermedia
Count for Garcinia.madruno
Count for Genipa.americana
Count for Guapira.myrtiflora
Count for Guarea.fuzzy
Count for Guarea.grandifolia
Count for Guarea.guidonia
Count for Guatteria.dumetorum
Count for Guazuma.ulmifolia
Count for Guettarda.foliacea
Count for Gustavia.superba
Count for Hampea.appendiculata
Count for Hasseltia.floribunda
Count for Heisteria.acuminata
Count for Heisteria.concinna
Count for Hirtella.americana
Count for Hirtella.triandra
Count for Hura.crepitans
Count for Hieronyma.alchorneoides
Count for Inga.acuminata
Count for Inga.cocleensis
Count for Inga.goldmanii
Count for Inga.laurina
Count for Inga.semialata
Count for Inga.nobilis
Count for Inga.oerstediana
Count for Inga.pezizifera
Count for Inga.punctata
Count for Inga.ruiziana
Count for Inga.sapindoides
Count for Inga.spectabilis
Count for Inga.umbellifera
Count for Jacaranda.copaia
Count for Lacistema.aggregatum
Count for Lacmellea.panamensis
Count for Laetia.procera
Count for Laetia.thamnia
Count for Lafoensia.punicifolia
Count for Licania.hypoleuca
Count for Licania.platypus
Count for Lindackeria.laurina
Count for Lonchocarpus.heptaphyllus
Count for Luehea.seemannii
Count for Macrocnemum.roseum
Count for Maquira.guianensis.costaricana
Count for Margaritaria.nobilis
Count for Marila.laxiflora
Count for Maytenus.schippii
Count for Miconia.affinis
Count for Miconia.argentea
Count for Miconia.elata
Count for Miconia.hondurensis
Count for Mosannona.garwoodii
Count for Myrcia.gatunensis
Count for Myrospermum.frutescens
Count for Nectandra.cissiflora
Count for Nectandra.lineata
Count for Nectandra.purpurea
Count for Ochroma.pyramidale
Count for Ocotea.cernua
Count for Ocotea.oblonga
Count for Ocotea.puberula
Count for Ocotea.whitei
Count for Oenocarpus.mapora
Count for Ormosia.amazonica
Count for Ormosia.coccinea
Count for Ormosia.macrocalyx
Count for Pachira.quinata
Count for Pachira.sessilis
Count for Perebea.xanthochyma
Count for Cinnamomum.triplinerve
Count for Picramnia.latifolia
Count for Piper.reticulatum
Count for Platymiscium.pinnatum
Count for Platypodium.elegans
Count for Posoqueria.latifolia
Count for Poulsenia.armata
Count for Pourouma.bicolor
Count for Pouteria.fossicola
Count for Pouteria.reticulata
Count for Pouteria.stipitata
Count for Prioria.copaifera
Count for Protium.costaricense
Count for Protium.panamense
Count for Protium.tenuifolium
Count for Pseudobombax.septenatum
Count for Psidium.friedrichsthalianum
Count for Psychotria.grandis
Count for Pterocarpus.rohrii
Count for Quararibea.asterolepis
Count for Quassia.amara
Count for Randia.armata
Count for Sapium.broadleaf
Count for Sapium.glandulosum
Count for Schizolobium.parahyba
Count for Senna.dariensis
Count for Simarouba.amara
Count for Siparuna.guianensis
Count for Siparuna.pauciflora
Count for Sloanea.terniflora
Count for Socratea.exorrhiza
Count for Solanum.hayesii
Count for Sorocea.affinis
Count for Spachea.membranacea
Count for Spondias.mombin
Count for Spondias.radlkoferi
Count for Sterculia.apetala
Count for Swartzia.simplex.var.grandiflora
Count for Swartzia.simplex.continentalis
Count for Symphonia.globulifera
Count for Handroanthus.guayacan
Count for Tabebuia.rosea
Count for Tabernaemontana.arborea
Count for Tachigali.versicolor
Count for Talisia.nervosa
Count for Talisia.princeps
Count for Terminalia.amazonia
Count for Terminalia.oblonga
Count for Tetragastris.panamensis
Count for Tetrathylacium.johansenii
Count for Theobroma.cacao
Count for Thevetia.ahouai
Count for Tocoyena.pittieri
Count for Trattinnickia.aspera
Count for Trema.micrantha
Count for Trichanthera.gigantea
Count for Trichilia.pallida
Count for Trichilia.tuberculata
Count for Trichospermum.galeottii
Count for Triplaris.cumingiana
Count for Trophis.caucana
Count for Trophis.racemosa
Count for Turpinia.occidentalis
Count for Unonopsis.pittieri
Count for Virola.multiflora
Count for Virola.sebifera
Count for Virola.surinamensis
Count for Vismia.baccifera
Count for Vochysia.ferruginea
Count for Xylopia.macrantha
Count for Zanthoxylum.ekmanii
Count for Zanthoxylum.juniperinum
Count for Zanthoxylum.panamense
Count for Zanthoxylum.setulosum
Count for Zuelania.guidonia
https://www.science.org/doi/10.1126/science.1066854
Calculate the change in variance with increasing number of sites
calc_delta_var( data, col_name, site_name = "num_sites", rep_name = "rep", visualize = FALSE )
calc_delta_var( data, col_name, site_name = "num_sites", rep_name = "rep", visualize = FALSE )
data |
A data frame containing the biodiversity indices to analyze, for a different number of sites over multiple repetitions. |
col_name |
The name of the column containing the biodiversity index to analyze. |
site_name |
The name of the column containing the number of sites. |
rep_name |
The name of the column containing the repetition number. |
visualize |
A logical indicating whether to visualize the results. |
A data frame with the number of sites and the variance and standard deviation of the mean of the biodiversity index for each number of sites.
ss <- generate_subsamples(BCI, min_sites = 1, max_sites = 5, step = 1, reps = 5) data <- ss data <- unlist(data, recursive = FALSE) data <- do.call(rbind, data) calc_delta_var(data, col_name = "richness", site_name = "num_sites", rep_name = "rep", visualize = TRUE)
ss <- generate_subsamples(BCI, min_sites = 1, max_sites = 5, step = 1, reps = 5) data <- ss data <- unlist(data, recursive = FALSE) data <- do.call(rbind, data) calc_delta_var(data, col_name = "richness", site_name = "num_sites", rep_name = "rep", visualize = TRUE)
Calculate biodiversity summary indices from count data
calc_diversity_indices(data)
calc_diversity_indices(data)
data |
A data frame of count data, with sites as rows and species as columns. |
A data frame with sites as rows and diversity indices as columns. The columns are: abundance, species richness, Shannon diveristy index, Simpson diversity index, Chao1, Difference between Choa1 and species richness.
ind <- calc_diversity_indices(BCI)
ind <- calc_diversity_indices(BCI)
Create multiple resamples of a data set.
create_resample(data, reps = 100, summary = TRUE, seed = sample(0:9999, 1))
create_resample(data, reps = 100, summary = TRUE, seed = sample(0:9999, 1))
data |
A data frame of count data, with sites as rows and species as columns. |
reps |
The number of resamples to create. |
summary |
A logical indicating whether to calculate summary indices
using |
seed |
A random seed to use for reproducibility. |
A list of data frames, if summary = FALSE
, each data frame
is a resample of the original data set. If summary = TRUE
,
each data frame is a resample of the original data set with
diversity indices calculated using
calc_diversity_indices
.
rs <- create_resample(BCI, reps = 10, summary = TRUE)
rs <- create_resample(BCI, reps = 10, summary = TRUE)
Generate subsamples of a data frame with a number of sites between a minimum and maximum value.
generate_subsamples( data, min_sites = 1, max_sites = nrow(data), step = 1, reps = 100, summary = TRUE, seed = sample(0:9999, 1) )
generate_subsamples( data, min_sites = 1, max_sites = nrow(data), step = 1, reps = 100, summary = TRUE, seed = sample(0:9999, 1) )
data |
A data frame of count data, with sites as rows and species as columns. |
min_sites |
The minimum number of sites to include in a subsample. |
max_sites |
The maximum number of sites to include in a subsample. Defaults to the number of sites in the original data set. |
step |
The number of sites to increase by at each iteration. |
reps |
The number of subsamples with a given number of sites to create. |
summary |
A logical indicating whether to calculate summary indices
using |
seed |
A random seed to use for reproducibility. |
A list of lists of data frames, if summary = FALSE
, each data
frame is a subsample of the original data set. If summary =
TRUE
, each data frame is a subsample of the original data set with
diversity indices calculated using
calc_diversity_indices
.
ss <- generate_subsamples(BCI, min_sites = 1, max_sites = 5, step = 1, reps = 2)
ss <- generate_subsamples(BCI, min_sites = 1, max_sites = 5, step = 1, reps = 2)
Calculate biodiversity measures and summary statistics for a data set using repeated sampling
get_sample_stats(data, sites_col = 1, reps = 100, indices = "all")
get_sample_stats(data, sites_col = 1, reps = 100, indices = "all")
data |
A data frame of count data, with sites as rows and species as columns. |
sites_col |
The column number of column containing site IDs. |
reps |
The number of resamples to create. |
indices |
A vector of indices to calculate. Use "all" to calculate all indices. Available indices are: abundance, richness, shannon, simpson, chao1, and chao_diff. |
A list of two data frames. The first data frame contains site specific data with sites as rows and summary statistics as columns. The second contains an overall summary of the data.
stats <- get_sample_stats(BCI, reps = 5)
stats <- get_sample_stats(BCI, reps = 5)