tax_dat <- read_tax(system.file("extdata", "test.taxonomy",
package = "schtools"
))
shared_dat <- readr::read_tsv(system.file("extdata", "test.shared",
package = "schtools"
))
pool_taxon_counts(shared_dat, tax_dat, "genus")
pool_taxon_counts(shared_dat, tax_dat, "family")
pool_taxon_counts(shared_dat, tax_dat, "phylum")
Run the code above in your browser using DataLab