起因是笔者选取了otu table中的部分otu,另外有全部otu的注释表(taxonomy),有些分析需要保证otu table中的otu数目和taxonomy中otu数目一致,不然就会出现以下报错
接下来是我原始的数据格式
otu table
taxonomy表格
所以目的就是筛选taxonomy中的ASV,让它跟otu表格是能对应上的
用了如下代码
#导入两个表格
species <- read.csv("otutable.txt",header=T,sep="",row.names=1)
taxa <-read.csv("taxonomy.txt",header=T,sep="",row.names=1)
#筛选
taxa1=taxa[c(rownames(species)),]
略做学习记录,欢迎交流指正