# NOT RUN {
data(phy); data(cla); data(ord); data(fam); data(clin)
otu_tabs = list(Phylum = phy, Class = cla, Order = ord, Family = fam)
set <- tidy_micro(otu_tabs = otu_tabs, clinical = clin) %>%
filter(day == 7) ## Only including the first week
## Bray-Curtis beta diversity
bray <- set %>% beta_div(table = "Family")
bray %>% beta_heatmap(micro_set = set, bpd1)
# }
Run the code above in your browser using DataLab