从转录组中提取感兴趣的基因子集

1、读入转录组数据

expdata <- read.table("GSE119382expr.txt",header = T)

expdata <- expdata[-c(1:5),]   #去除前5行不需要的数据

2、读入自己的基因列表

genelist <- read.table("gene_id_name.txt",header = T,sep = "\t")

3、提取子集

方法一：用merge，可保留两个表的信息

colnames(genelist)[1] <- "ID"   #把合并的依据列名称统一

terp_expdata1 <- merge(expdata,genelist,by = "ID")

【%in%，网页上有乱码，只好贴图】

转载本文请联系原作者获取授权，同时请注明本文来自陈新科学网博客。

链接地址：https://wap.sciencenet.cn/blog-62701-1149618.html?mobile=1

收藏

分享