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.
hilldiv3 measures and compares the diversity of
biological communities — OTU, ASV or MAG count tables — using
Hill numbers. Hill numbers are a single, intuitive
family of diversity metrics: each one is an effective number of
taxa (“how many equally-abundant taxa would give this much
diversity”), and one parameter, the diversity order
q, slides smoothly between counting all taxa equally
(richness), weighting them by abundance (Shannon) and focusing on the
common ones (Simpson). Because everything is expressed in the same
currency, results are directly comparable across samples, studies and
methods.
From that one foundation, hilldiv3 provides a unified
toolkit for neutral diversity (abundances only),
phylogenetic diversity (accounting for how related taxa
are) and functional diversity (accounting for how
different their traits are), covering measurement, partitioning,
(dis)similarity, profiles, evenness and redundancy. You call the same
functions for all three — the diversity type is chosen by whether you
supply a tree or a distance matrix.
# install.packages("devtools")
devtools::install_github("alberdilab/hilldiv3")library(hilldiv3)
# Bundled simulated gut-microbiome MAG data.
hilldiv(gut_counts) # neutral Hill numbers q = 0, 1, 2
hilldiv(gut_counts, tree = gut_tree) # neutral + phylogenetic
dist <- traits2dist(gut_traits)
hilldiv(gut_counts, dist = dist) # neutral + functional
hilldiv(gut_counts, tree = gut_tree, dist = dist) # all three types at once
# Results are tidy by default and plot directly.
plot(hillprof(gut_counts)) # diversity profile
hilldiv(gut_counts, out = "matrix") # matrix: samples x q ordersFull documentation lives on the package website: https://alberdilab.github.io/hilldiv3/
vignette("hilldiv3"), a
gentle introduction for anyone using Hill numbers for the first
time.If you have used hilldiv2, v3 is a
full redesign that keeps the familiar function names
(hilldiv(), hillpart(),
hilldiss(), hillsim(),
hillpair(), hillred(), tss(),
traits2dist()) while changing how they work underneath:
match_data() helper.phyloseq and TreeSummarizedExperiment
objects.ape post-order
traversal in place of geiger::tips();
hillpair() computes the shared structure once and reuses it
across all sample pairs.hill* function
returns a long-format data.frame with
print()/plot()/autoplot()
methods; pass out = "matrix" for a plain matrix.type = c("auto", "neutral", "phylogenetic", "functional")
argument that asserts and validates the diversity type (auto-detected by
default).hillprof() (diversity profiles) and
hilleven() (evenness), hierarchical multi-scale
partitioning in hillpart(), plus bundled example data
(gut_counts, gut_tree,
gut_traits).See NEWS.md for the full
changelog.
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.