# 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) %>%
filter(day == 7) ## Only including the first week
## Bray-Curtis beta diversity
bray <- set %>% beta_div(table = "Family")
## Morisita-Horn beta diversity
horn <- set %>% beta_div(table = "Family", method = "horn")
# }
Run the code above in your browser using DataLab