欢迎光临散文网 会员登陆 & 注册

【clusterProfiler包】1行代码将SYMBOL转换为entrezID_比DAVID高效多了_20211027

2021-10-27 17:08 作者:天马行空的坦克兵  | 我要投稿

#有gene symbol怎么拿到entrezid ?

rm(list=ls())

library(org.Rn.eg.db)   #大鼠数据包

library(clusterProfiler)  

#看一下待转换数据

图1  ID数据框中  id列有219个基因,idid列有22个基因



#使用id列进行转换,没有成功。后面解释为什么不成功

,AA=bitr(ID$id,fromType = "SYMBOL",

              toType =c('ENTREZID') ,

              OrgDb = org.Rn.eg.db)



#使用idid列转换,成功了,但有些没有entrezid

(1)BB=bitr(ID$idid,fromType = "SYMBOL",

                toType =c('ENTREZID') ,

                OrgDb = org.Rn.eg.db)

图2  使用idid列转换,成功了,但有些没有entrezid

(2)CC=bitr(ID$idid,fromType = "SYMBOL",

        toType =c('ENSEMBL') ,

        OrgDb = org.Rn.eg.db)

图3  将idid列转换为ENSEMBL  ID

#-----------------

BBCC=bitr(ID$idid,fromType = "SYMBOL",

          toType =c('ENTREZID','ENSEMBL') ,

          OrgDb = org.Rn.eg.db)

图4  一次性获得图2和图3

#为什么只有idid列转换成功了?

#因为id列基因名称不对,id列基因全是大写字母,R语言系统识别不了

------------------------------下期继续---------------------------------------

2021年10月27日

up主:天马行空的坦克兵





【clusterProfiler包】1行代码将SYMBOL转换为entrezID_比DAVID高效多了_20211027的评论 (共 条)

分享到微博请遵守国家法律