数据框筛选特定的子集

a=read.table(file = "H_V75_ENSG,T,P_type-mart_export.txt",header = TRUE,sep = "\t")#以数据框格式读取原文件
pro_cod=a[a$Gene.Biotype=="protein_coding",]#筛选出数据框中第四列为protein_coding的子数据框
write.table(pro_cod,file = "pro_cod_selected.txt")#重新写入筛选后文件