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)
## Not run:
# # 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)
# ## End(Not run)
Run the code above in your browser using DataLab