Entering edit mode
Tim Smith
★
1.1k
@tim-smith-1532
Last seen 10.2 years ago
Hi,
I was trying to map gene symbols to gene names using the org.Hs.eg.db
package. I first convert the gene symbol to an entrez id, and then
convert that to a gene name (example code below). However, during this
process I can't get the gene names for some of the genes:
--------------------------
library(org.Hs.eg.db)
### First two genes are ok...
symbols <- c('JAK1','AATK','A2BP1','A2LD1')
for(sym in symbols){
print('***')
entrezid <-
as.character(mget(sym,org.Hs.egSYMBOL2EG,ifnotfound=NA))
print(sym)
print(entrezid)
genename <-
as.character(mget(entrezid,org.Hs.egGENENAME,ifnotfound=NA))
print(genename)
}
-------------------------
I don't get any gene names for the last two genes in this example. Is
there a better way of doing this? Perhaps with some other
function/package?
thanks.
[[alternative HTML version deleted]]