# NOT RUN {
data(bpd_phy); data(bpd_cla); data(bpd_ord); data(bpd_fam); data(bpd_clin)
otu_tabs = list(Phylum = bpd_phy, Class = bpd_cla,
Order = bpd_ord, Family = bpd_fam)
set <- tidy_micro(otu_tabs = otu_tabs, clinical = bpd_clin) %>%
mutate(bpd1 = factor(bpd1))
## Summarize each taxa by Table
set %>% taxa_summary
## Summarize each taxa by a categorical variable of interest
set %>% taxa_summary(bpd1)
## Summarize each taxa by a categorical variable of interest within a Table
set %>% taxa_summary(bpd1, table = "Phylum")
## Summarize within group or table only
set %>% taxa_summary(taxa = FALSE)
# }
Run the code above in your browser using DataLab