...
| Code Block |
|---|
asvtable <- read.table("../illumina/results/dada2/ASV_table.tsv", check.names = FALSE, sep = "\t", skip = 1, stringsAsFactors = FALSE, comment.char = "", header = T)
colnames(asvtable)[1] <- "ASV" |
...
| Code Block |
|---|
asv_table <- cbind(asvtable, mytaxsplitmytax) # Also remove the ASV_ID, Confidence and sequence columns, as they are not needed asv_table <- subset(asv_table, select=-c(ASV_ID, Confidenceconfidence, sequence)) |
| Code Block |
|---|
head(asv_table) |
...