data(ITS1)
# filter.OTU() returns a list
otu <- filter.OTU(list(ITS1=ITS1), percent=0.001)[[1]]
# split and keep taxonomy column
otu.split <- col.splitup(otu, col="taxonomy", sep="; ",
drop=FALSE)
# give new column names
tax.classes <- c("kingdom", "phylum", "class",
"order", "family", "genus")
otu.split <- col.splitup(otu, col="taxonomy", sep="; ",
drop=TRUE, names=tax.classes)Run the code above in your browser using DataLab