Map gene IDs to gene names
1
0
Entering edit mode
6.9 years ago
bisansamara ▴ 20

Hi, I have a huge list of gene names, and I want to assign the ID for each gene. I tried the following code, and it kind of worked, but the problem was that it created extra IDs (more than the available genes) --> causing misalignment between names and IDs

CODE

library("org.Hs.eg.db")

a <- read.csv("filename.csv",TRUE,",")

y=as.character(a$Gene.refGene) #column name is Gene.refGene

gene=y

output = unlist(mget(x=gene,envir=org.Hs.egALIAS2EG,ifnotfound=NA))

write.csv(output, file = "output.csv") #write the IDs which is the output in a csv file

Is there any other way to get the gene IDs? or any suggestion on how to modify the code? your help is highly appreciated!

gene R Ensembl GeneID • 1.8k views
ADD COMMENT
1
Entering edit mode

Some gene names have more IDs and vice versa. There is no way around it. One way to deal with it is to keep only one.

ADD REPLY
0
Entering edit mode
ADD COMMENT

Login before adding your answer.

Traffic: 2637 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6