数据a
数据b
即提取a数据前两列中的基因都属于b的数据及其所在的行
x <- matrix(NA,ncol = 5,nrow = 6793)
x <- as.data.frame(x)
j<- 1
for (i in 1:dim(a)[1]){
if((a$gene[i]%in%b$gene_name)&&(a$gene.1[i]%in%b$gene_name))
{x[j,] <- a[i,]
j <- j+1
}
}
最后删除掉多余的行
或者
x <- a[((a$gene)%in%(b$gene_name))&((a$gene.1)%in%(b$gene_name)),]