想要实现的目标如题。
善用属于符号:%in%
,以及不要忘记XX$V1
。
我想要提取数据框cpm2中行名在vector interest中的行,使用如下指令解决。
interest<-read.table("result_coding_gene_id_2.txt")
data<-cpm2[which(rownames(cpm2)%in%interest$V1),]
- interest
V1
1 ENSG00000000460
2 ENSG00000001084
3 ENSG00000001497
4 ENSG00000001626
5 ENSG00000002746
6 ENSG00000002822
- rownames(cpm2)
head(row.names(cpm2))
[1] “ENSG00000000003” “ENSG00000000005” “ENSG00000000419” “ENSG00000000457”
[5] “ENSG00000000460” “ENSG00000000938”
如上,实现需求。