美文网首页学习
Gene id 的转换方式

Gene id 的转换方式

作者: 超级无敌大蜗牛 | 来源:发表于2020-02-25 22:07 被阅读0次

先把基因集对应的注释包安装了,并且加载出来

library(org.Hs.eg.db) #基因注释包。

基因集准备

gene <- write.csv("gene.csv")

第一种方法

安装clusterProfiler包,bitr转换ID函数属于这个包
BiocManager::install(clusterProfiler)
library("clusterProfiler")

gene.df <- bitr(gene, fromType = "ENTREZID", #fromType是指你的数据ID类型是属于哪一类的
                toType = c("ENSEMBL", "SYMBOL"), #toType是指你要转换成哪种ID类型,可以写多种,也可以只写一种
                OrgDb = org.Hs.eg.db)

第二种方法

library(AnnotationDbi)
mySymbols <- mget(gene,
                  org.Hs.egSYMBOL, #这个是可以选择的,选择不同,转换的ID类型也不一样
                  ifnotfound=NA)
转换成Symbol ID
head(mySymbols)
class(mySymbols)

第三种方法

geneIDselect <-select(org.Hs.eg.db, #.db是这个芯片数据对应的注释包
                      keys=gene,
                      columns=c("SYMBOL","ENSEMBL","GENENAME"), #clolumns参数是你要转换的ID类型是什么,这里选择三个。
                      keytype="ENTREZID" )#函数里面的keytype与keys参数是对应的,keys是你输入的那些数据,keytype是指这些数据是属于什么类型的数据。 
head(geneIDselect)

第四种方法

gene <-mapIds(org.Hs.eg.db, #.db是这个芯片数据对应的注释包
                      keys=gene,
                      column="SYMBOL", #clolumns参数是你要转换的ID类型是什么,只能选择一个。
                      keytype="ENTREZID" )#函数里面的keytype与keys参数是对应的,keys是你输入的那些数据,keytype是指这些数据是属于什么类型的数据。 
head(geneIDselect)

原文链接

相关文章

网友评论

    本文标题:Gene id 的转换方式

    本文链接:https://www.haomeiwen.com/subject/vywqchtx.html