2021-12-14 gene ID 转换

genes<-data$gene
head(genes)
图片.png
# clusterProfile包的bitr转化
library(clusterProfiler)
gene_an<-bitr(data$gene,
              fromType = "ENSEMBL",
               toType = c("ENTREZID","SYMBOL"),
              OrgDb = "org.Mm.eg.db")
##转化的gene_an为data.frame
图片.png
## 
library(org.Mm.eg.db)
library(AnnotationDbi)
gene_ann2<-mapIds(org.Mm.eg.db,
                  keys=genes, ##输入的ID
                  column ="SYMBOL",  #要转化为的类型,只能选1种
                  keytype ="ENSEMBL", ##输入ID的类型
                  multiVals = "first" # 多种转换结果选择第一个)
##  gene_ann2    为字符向量
图片.png
library(org.Mm.eg.db)
library(AnnotationDbi)
k<-keys(org.Mm.eg.db,keytype = "ENSEMBL")
##得到所有的ensembl_ID以及其对应的SYMBOL和ENTRZID
gene_list<-AnnotationDbi::select(org.Mm.eg.db,
                  keys=k, ##所有的ensembl_ID
                  columns = c("ENTREZID","SYMBOL"),#要转化为的类型
                   keytype="ENSEMBL")
## 在list的ENSEMBL中匹配我们要转化的ensemble_ID(其他类型类似)
gene_ann3<-gene_list[match(genes,gene_list[,"ENSEMBL"]),]
head(gene_ann3)    
temp2<-na.omit(gene_ann3) # 去除未注释到的genes
图片.png
# 方法同上,但是更为简单
gene_ann4<-AnnotationDbi::select(org.Mm.eg.db,
                                 keys=genes,
                                 columns = c("ENTREZID","SYMBOL"),
                                 keytype="ENSEMBL")
最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

推荐阅读更多精彩内容