diff --git a/analysis/16S/16S_Analysis_pipeline_v4-3.Rmd b/analysis/16S/16S_Analysis_pipeline_v4-3.Rmd index 02f27bebb6d2b9e1f99a9bd791978d5ea3c99e09..566e0f464a62eeb50a6f354df5a8faeb5b2cce0e 100644 --- a/analysis/16S/16S_Analysis_pipeline_v4-3.Rmd +++ b/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,