joey711/joey711.github.com

Official GitHub central page for joey711

joey711.github.com's Issues

Top 10 phyla with 10th categorised as 'Other' for stacked relative abundance bar plot

Hi there,

I managed to create a 100% stacked proportional bar chart representing relative abundance of phyla across 3 forest age categories using this code. However, what I want is this relative abundance figure except with only 10 phyla represented, including the 9 most abundant phyla and a 10th category merging all the remaining taxa.

Here's my code so far, where Bacteria_physeq is my phyloseq object and age.category is the variable of interest:

remove OTUs with no assigned phylum

cleaned_phyla = subset_taxa(Bacteria_physeq, Phylum!="NA" & Phylum!="")
cleaned_phyla
otu_table() OTU Table: [ 17566 taxa and 19 samples ]
sample_data() Sample Data: [ 19 samples by 6 sample variables ]
tax_table() Taxonomy Table: [ 17566 taxa by 7 taxonomic ranks ]

melt phyloseq object to use with ggplot2

Bacteriaphyseq <- psmelt(cleaned_phyla)

plot stacked bar plot for % relative abundance

barstackprop <- ggplot(Bacteriaphyseq, aes(x = age.category, y = Abundance, fill = Phylum)) +
geom_bar(position = "fill",stat = "identity") +
scale_y_continuous(labels = scales::percent_format()) +
labs(x = "Age category", y = "Relative abundance")
plot(barstackprop)

As you can see, there are far too many colours and phyla represented.

I've tried a couple of tutorials and other forum recommendations but can't seem to get it to work for what I want to achieve here. Any help would be really useful.