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]]

