注释表中提取特定的OTU注释行

起因是笔者选取了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)),] 

略做学习记录,欢迎交流指正

评论 2
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值