But I can use the code to make a barplot, like this
SVs<-read_qza("table-dada2-3.qza")$d
metadata<-read_q2metadata("qiime2/metadata.tsv")
head(metadata)
taxonomy<-read_qza("taxonomy.qza")$data %>% parse_taxonomy()submeta<-subset(metadata, Source=="Plastic")
subSVs<-SVs[,submeta$SampleID]
tsum<-summarize_taxa(subSVs,taxonomy)$Phylum
taxa_barplot(tsum, submeta, "Plastic-type", ntoplot=30)
How can I change my subgroup? such as this:
because I have a subgroup"Fertilizer-type" under the group of "Plastic-type".
I know how to change it in the original code,but i dont know in taxa_barplot.
In addition, how can I change the x axis to display the group and y axis display the log10?