xvtyzn / nyankomicro Goto Github PK

microbiome_barplot(physeq, "Genus", c("sampling_date","colony"),plot_percent = 10)

taxotu_table <- otu_table(agg_phylo_rel) %>% as.data.frame() %>% rownames_to_column("ASV") %>%
        left_join(level_tax, by = ("ASV" = "ASV")) %>% select(-ASV) %>% group_by(across(all_of(level))) %>%
        summarise_all(sum) %>% remove_rownames() %>% column_to_rownames(level)
    taxotu_table$mean <- apply(taxotu_table, 1, mean)

physeq <- prune_samples(sample_sums(physeq)>=1000, physeq)

══ Failed tests ════════════════════════════════════════════════════════════════
─ Error (test_microbiome_barplot.R:7:1): (code run outside of `test_that()`) ─
Error: `.data` must be a data frame without row names.
Backtrace:
    █
 1. ├─nyankomicro::microbiome_barplot(...) test_microbiome_barplot.R:7:0
 2. │ └─`%>%`(...)
 3. └─tibble::column_to_rownames(., level)