基因ID转换

2022-08-09  本文已影响0人  Shift_shift

#BiocManager::install("clusterProfiler")

#BiocManager::install("org.Hs.eg.db")

library(org.Hs.eg.db)

library("clusterProfiler")

data <- read.table("C:/Users/cendy/Desktop/counts.txt", header=TRUE, quote="\t",skip = 1)

geneList <- data$Geneid

source("http://bioconductor.org/biocLite.R")

options(BioC_mirror="http://mirrors.ustc.edu.cn/bioc/")

gene.df <- bitr(geneList, fromType = "ENSEMBL", #fromType是指你的数据ID类型是属于哪一类的

                toType = c("ENSEMBL", "SYMBOL"), #toType是指你要转换成哪种ID类型,可以写多种,也可以只写一种

                OrgDb = org.Hs.eg.db)#Orgdb是指对应的注释包是哪个

head(gene.df)

上一篇下一篇

猜你喜欢

热点阅读