R/Bioconductor与转录组数据分析(5)
2018-12-3 15:30
阅读:3268
标签:转录组数据分析
从转录组中提取感兴趣的基因子集
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
收藏
当前推荐数:0
推荐到博客首页
网友评论0 条评论