[BioC] biomaRt error in getGene

Georg Otto georg.otto at tuebingen.mpg.de
Thu Jul 13 16:36:19 CEST 2006


Besides "affy_zebrafish" there is also an attribute
"affy_zebrafish_primary_id". Could this have something to do with the
problem, or could I use "affy_zebrafish_primary_id" instead of
"affy_zebrafish"? Sorry for the naive question, but I am not very
familiar with ensembl/biomart...

Cheers,

Georg

Steffen Durinck <durincks at mail.nih.gov> writes:

> Hi Georg,
>
> I've repeated your query and unfortunately get the same error. 
> When using biomaRt with Ensembl,  it should return the same output as 
> the Ensembl BioMart web application (http://www.ensembl.org/Multi/martview).
> I often use this to verify that biomaRt is properly functioning or to 
> find out if the reported error is due to a database error or due to a 
> bug in the biomaRt package.
>
> After doing a query using zebrafish affy ids, the MartView web 
> application also returns no values for affy ids (e.g. querying for all 
> affy ids on zebrafish chromosome 1).
> It looks thus like there is an error in the zebrafish dataset of 
> Ensembl.  Any problems of  this type should be reported at the 
> corresponding database so they can be fixed in the next database 
> release. For Ensembl this is helpdesk at ensembl.org. 
> <mailto:helpdesk at ensembl.org>
> I'll document this in the biomaRt vignette.
>
> Best,
> Steffen
>
>
>
> Georg Otto wrote:
>> Hi Steffen,
>>
>> thanks for your help. I updated to biomaRt 1.7.3 and followed your
>> advice. What I got is this:
>>
>>   
>>> mart <- useMart("ensembl")
>>> mart<-useDataset("drerio_gene_ensembl", mart)
>>>     
>>
>>   
>>> getGene(id="Dr.10336.1.S1_at", array="affy_zebrafish", mart=mart)
>>>     
>> NULL
>>
>> then I tried with getBM
>>
>>   
>>> getBM(attributes=c("ensembl_transcript_id", "affy_zebrafish", "refseq_dna"), filters="affy_zebrafish", values="Dr.10336.1.S1_at", mart=mart)
>>>     
>> NULL
>>
>> I tried with several different affymetrix probes that should be
>> annotated in ensembl, but all of them gave me a "NULL"
>>
>>   
>>> sessionInfo()
>>>     
>> Version 2.3.0 (2006-04-24) 
>> x86_64-redhat-linux-gnu 
>>
>> attached base packages:
>> [1] "stats"     "graphics"  "grDevices" "utils"     "datasets"  "methods"  
>> [7] "base"     
>>
>> other attached packages:
>>  biomaRt    RCurl      XML 
>>  "1.7.3"  "0.6-2" "0.99-7" 
>>
>> Any idea what could be wrong here?
>>
>> Best,
>>
>> Georg



More information about the Bioconductor mailing list