Thank you. I used getURL function of RCurl to retrieve the information from the site.
Hi.
Like other identifiers, entrez ids are also change as time goes by.
I used two library, 'org.Hs.eg.db' and 'annotate', to convert entrez ids into gene symbols.
Some of entrez ids were not changed because of updates.
Take entrez id 164022 as an example,
library(org.Hs.eg.db);library(annotate)
getSYMBOL('164022',data = 'org.Hs.eg')
164022
NA
If you search the entrez id at NCBI, it says that 164022 was replaced with 653505.
(See. http://www.ncbi.nlm.nih.gov/gene/?term=164022)
Therefore, you should use the newest entrez id to get its gene symbol.
getSYMBOL('653505',data = 'org.Hs.eg')
653505
"PPIAL4A"
There are about two hundreds entrez ids whose symbol-matching failed.
As manual searching requires huge time, I need a solution.
How to update old entrez ids into newest entrez id?
Is there a function or library for this?
1 answer
This is a very good question.
You may try with the mygene.info service:
$: curl mygene.info/v2/gene/164022?fields=entrezgene
{
_id: "653505",
entrezgene: 653505
}
Remove the "fields" parameter to get more information.
For more documentation on mygene.info, check http://mygene.info/v2/api#MyGene.info-gene-annotation-services-GET-Gene-annotation-service
Log in to answer this question.
Hello agicict!
It appears that your post has been cross-posted to another site: http://stackoverflow.com/questions/38366725
This is typically not recommended as it runs the risk of annoying people in both communities.
Thank you for your comment. I made a mistake.