proof of principle for revision code

analysis/16S/16S_Analysis_pipeline_v4-3.Rmd

@@ -853,6 +853,13 @@ set$topntax_tbl <- psmelt(set$ps.topnTax) %>%
                     arrange(ndvi) %>%
                     mutate(ndvi_rank = c(1:length(ndvi))) 
                     
+set$other <- set$topntax_tbl %>% 
+  group_by(Alias) %>% 
+  summarise(sumabu = sum(Abundance)) %>% 
+  mutate(Genus = 'Other', Abundance = 1-sumabu) %>% 
+  select(-sumabu)
+
+set$topntax_tbl <- bind_rows(set$topntax_tbl, set$other)
 
 plots$topn_tax_custom <- ggplot(set$topntax_tbl, aes(x = fct_reorder(Alias, ndvi),
                            y = Abundance,