I've several refseq id that I want to convert to gene symbol. So in R with biomaRt :
library(biomaRt) mart<- useDataset("hsapiens_gene_ensembl", useMart("ensembl")) refseq <- c("NM_001271003") getBM(filters="refseq_mrna", attributes="external_gene_id", values=refseq, mart=mart)
gives me logical(0)
But the gene is well defined in pubmed (TFPI2) : http://www.ncbi.nlm.nih.gov/nuccore/NM_001271003
The gene is also in ensembl : http://www.ensembl.org/Homo_sapiens/Gene/Matches?db=core;g=ENSG00000105825;r=7:93514709-93520303 but in external references, the refseq id : NM_001271003 is not reported..
So is there an another way to get the gene symbol from refseq id using R ?