hill_tuckey_pq: Calculate hill number and compute Tuckey post-hoc test

In MiscMetabar: Miscellaneous Functions for Metabarcoding Analysis

Calculate hill number and compute Tuckey post-hoc test

Description

Note that, by default, this function use a sqrt of the read numbers in the linear model in order to correct for uneven sampling depth.

Usage

hill_tuckey_pq(
  physeq,
  modality,
  q = c(0, 1, 2),
  hill_scales = lifecycle::deprecated(),
  silent = TRUE,
  correction_for_sample_size = TRUE,
  ...
)

Arguments

physeq	(required) a phyloseq-class object obtained using the phyloseq package.
modality	(required) the variable to test
q	(numeric vector) Hill diversity orders to compute (q values). Default computes Hill number 0 (species richness), Hill number 1 (exponential of Shannon index) and Hill number 2 (inverse of Simpson index). Formerly hill_scales. Hill numbers are more appropriate in DNA metabarcoding studies when q > 0 (Alberdi & Gilbert, 2019; Calderón-Sanou et al., 2019).
hill_scales	Use q instead.
silent	(logical) If TRUE, no message are printing.
correction_for_sample_size	(logical, default TRUE) This function use a sqrt of the read numbers in the linear model in order to correct for uneven sampling depth.
...	Additional arguments passed to divent_hill_matrix_pq() and hence to divent::div_hill() (e.g. estimator = "naive" to match vegan-style results).

Value

A ggplot2 object

Author(s)

Adrien Taudière

References

Alberdi, A., & Gilbert, M. T. P. (2019). A guide to the application of Hill numbers to DNA-based diversity analyses. Molecular Ecology Resources. \Sexpr[results=rd]{tools:::Rd_expr_doi("10.1111/1755-0998.13014")}

Calderón-Sanou, I., Münkemüller, T., Boyer, F., Zinger, L., & Thuiller, W. (2019). From environmental DNA sequences to ecological conclusions: How strong is the influence of methodological choices? Journal of Biogeography, 47. \Sexpr[results=rd]{tools:::Rd_expr_doi("10.1111/jbi.13681")}

Examples

data_f <- prune_samples(
  sample_names(data_fungi_mini)[1:20],
  data_fungi_mini
)
hill_tuckey_pq(data_f, "Height")

MiscMetabar documentation built on June 8, 2026, 5:07 p.m.