DivNet: an R package to estimate diversity when taxa in the community cooccur via a ecological network.
Willis, A.D. and Martin, B.D. (2020) DivNet: Estimating diversity in networked communities. Biostatistics. doi.org/10.1093/biostatistics/kxaa015
If you knew the exact composition of a community (such as all of the microbes living on your skin), you could calculate the Shannon diversity of the community, or the Bray-Curtis distance between your skin’s microbes and your cat’s fur microbes. However, you will only ever observe a small fraction of the microbes on your skin from your experiment, so you have to estimate these diversity indices based on the data, and ideally you will get a confidence interval for those diversity indices.
Almost all ecologists use “plug-in” diversity estimates, obtained by taking the observed relative abundances of each taxon and plugging them into the formula for the true diversity. This is problematic for a number of reasons, including that there may be unobserved taxa, and the observed relative abundances probably don’t exactly equal the true relative abundances due to random sampling. Furthermore, typical approaches to obtaining a confidence interval assume that all taxa behave independently: if taxon A is there, this tells you nothing about whether taxon B is likely to be there too. DivNet takes care of all of these issues, and gives you a confidence interval that incorporates all of these challenges.
Want a confidence interval for a different diversity index? Let us know by posting an issue!
DivNet is currently available via GitHub, and can be installed as follows. Please note that you will need to install the latest version of
breakaway to install
DivNet; this is included in the following commands:
# install.packages("remotes") ## run this line if you do not already have remotes installed remotes::install_github("adw96/breakaway") remotes::install_github("adw96/DivNet") library(DivNet)
Let’s analyse the Lee et al. dataset using DivNet. We could analyse at the ASV level, but to keep things simple let’s analyse the phylum level diversity.
library(phyloseq) data(Lee) divnet_phylum <- divnet(tax_glom(Lee, taxrank="Phylum"), X = "char", ncores = 4) divnet_phylum
To test if alpha-diversity (by default, Shannon) is equal across the values of the covariate X, we can run
and to plot the DivNet estimates of diversity we can run
Advanced Usage (for large datasets)
One limitation of DivNet is that it is slow when you have many taxa (thousands and upwards). Fortunately for us, the amazing Ryan Moore implemented a faster version of DivNet in Rust called divnet-rs. Check out the documentation here.
Unfortunately the StatDivLab do not have expertise with
divnet-rs, so please direct your questions about
divnet-rs to Ryan and the
divnet-rs team (
divnet-rs Issues page here).
DivNet integrates with