Practical-3-Processing-16S-rRNA-amplicon-data.md

@@ -157,8 +157,7 @@ plotTaxon <- function(seqtab, taxtable, taxon, min_abund, taxonomicLevel) {
   seqtab.m$Order <- sapply(strsplit(seqtab.m$variable, "-"), "[[", 4)
   seqtab.m$Family <- sapply(strsplit(seqtab.m$variable, "-"), "[[", 5)
   seqtab.m$Genus <- sapply(strsplit(seqtab.m$variable, "-"), "[[", 6)
-  seqtab.m <- cbind(seqtab.m, replicate(1,seqtab.m[colnames(seqtab.m) == 
-taxonomicLevel]))
+  seqtab.m <- cbind(seqtab.m, replicate(1,seqtab.m[colnames(seqtab.m) == taxonomicLevel]))
   colnames(seqtab.m) <- c(colnames(seqtab.m)[1:length(colnames(seqtab.m)) -1], "displayLevel")
   nb.cols <- length(seqtab) - 1
   mycolours <- rep(brewer.pal(12, "Paired"), ceiling(nb.cols/12))