data("Three_group")
# Check meta file
print(Three_group$Groupfile)
# Subsetting tax summary objects
# Select BF and OF groups
sub_testtax_summary <- sub_tax_summary(Three_group, Group %in% c("BF", "OF"))
print(sub_testtax_summary$Groupfile)
# Subsetting according to taxonomy
Proteo <- sub_tax_summary(
Three_group,
taxnum = which(Three_group$Base_taxonomy$Phylum == "p__Proteobacteria")
)
print(Proteo$Phylum_percent) # Check phylum table
print(Proteo$Genus_percent) # Check genus table
Run the code above in your browser using DataLab