Microbiome-oriented wrapper around agglomerate_taxa() that summarizes
abundance at a chosen taxonomy rank.
Usage
summarise_taxa(
object,
taxonomic_rank = c("Kingdom", "Phylum", "Class", "Order", "Family", "Genus", "Species"),
what = c("sum_intensity", "mean_intensity", "median_intensity"),
na_remove = TRUE
)
Arguments
- object
A microbiome_dataset object.
- taxonomic_rank
A taxonomy rank in variable_info.
- what
Aggregation method.
- na_remove
Should taxa with missing rank be removed before agglomeration?
Value
A microbiome_dataset object.
Examples
data("global_patterns", package = "microbiomedataset")
x0 <- prune_taxa(global_patterns, variable_id = global_patterns@variable_info$variable_id[1:200])
x <- summarise_taxa(
object = x0,
taxonomic_rank = "Family",
what = "sum_intensity"
)
dim(x@expression_data)
#> [1] 16 26