Update Practical 3 Processing 16S rRNA amplicon data authored by Taylor Priest's avatar Taylor Priest
...@@ -278,18 +278,18 @@ This looks much better, but it is hard to distinguish between some of the colour ...@@ -278,18 +278,18 @@ This looks much better, but it is hard to distinguish between some of the colour
> marmic_filtered_family_barplot_refined <- ggplot(marmic_filtered_family_barplot_data, > marmic_filtered_family_barplot_refined <- ggplot(marmic_filtered_family_barplot_data,
aes(x=Sample, y=Abundance)) + aes(x=Sample, y=Abundance)) +
geom_bar(aes(fill=Family), stat="identity", position="stack") + geom_bar(aes(fill=Family), stat="identity", position="stack") +
facet_grid(~Sample_group, scales="free_x") + facet_grid(~Sample_group, scales="free_x") +
guides(fill=guide_legend(ncol=1)) + guides(fill=guide_legend(ncol=1)) +
labs(y="Relative abundance", x="Sample") + labs(y="Relative abundance", x="Sample") +
scale_fill_manual(values=colorRampPalette(brewer.pal(12,"Paired"))(colourCount)) + scale_fill_manual(values=colorRampPalette(brewer.pal(12,"Paired"))(colourCount)) +
theme_bw() + theme_bw() +
theme(axis.title.x = element_text(size=12,face="bold",colour="black"), theme(axis.title.x = element_text(size=12,face="bold",colour="black"),
axis.title.y = element_text(size=12,face="bold",colour="black"), axis.title.y = element_text(size=12,face="bold",colour="black"),
axis.text.x = element_text(size=10, colour="black", angle=45, hjust=1), axis.text.x = element_text(size=10, colour="black", angle=45, hjust=1),
legend.title = element_text(size=12,face="bold",colour="black"), legend.title = element_text(size=12,face="bold",colour="black"),
strip.background.x = element_rect(fill="white"), strip.background.x = element_rect(fill="white"),
strip.text.x = element_text(size=11,face="bold",colour="black")) strip.text.x = element_text(size=11,face="bold",colour="black"))
> marmic_filtered_family_barplot_refined > marmic_filtered_family_barplot_refined
> pdf(file="Marmic_NGS2022_amplicons_family_barplot.pdf", height=10, width=14) > pdf(file="Marmic_NGS2022_amplicons_family_barplot.pdf", height=10, width=14)
... ...
......