[BioC] biomaRt query
Sebastien Gerega
seb at gerega.net
Wed Dec 10 03:55:42 CET 2008
I thought that might be the case, I was just surprised by how many genes
didn't show up in Ensembl....
Sean Davis wrote:
> On Tue, Dec 9, 2008 at 7:21 PM, Sebastien Gerega <seb at gerega.net> wrote:
>
>> Hi,
>> there is probably a very simple answer to this question but why does the
>> following code produce the output NULL?
>>
>> ensembl = useMart("ensembl", dataset="hsapiens_gene_ensembl")
>> getBM(attributes = c("entrezgene", "hgnc_symbol"), filters = "entrezgene",
>> values = "117581", mart = ensembl, na.value = NA)
>>
>> The NCBI site has a hgnc symbol listed for the gene:
>> http://www.ncbi.nlm.nih.gov/sites/entrez?db=gene&cmd=Retrieve&dopt=Graphics&list_uids=117581
>>
>
> If you do a similar search at Ensembl, it appears that Ensembl does
> not have an entry for TWIST2. (I did the search quickly, so I could
> be mistaken.) Remember that these genome databases do not always
> agree about what a gene represents (or even which ones exist).
>
> Sean
>
>
More information about the Bioconductor
mailing list