Hello everybody,
I am using biomaRt to retrieve worm annotation.
Sometimes (apparently unpredictably) I get a NULL output and a
warning:
> getBM (attributes=c
("wormbase_gene","go_description"),filter="wormbase_gene",values=wormi
ds
,mart=mart)->gobiomart
Warning message:
getBM returns NULL. in: getBM(attributes = c("wormbase_gene",
"go_description"), filter = "wormbase_gene",
The exact same line works later. Has anyone experienced this? Is
there a timeout issue or something like that?
Thanks
Giovanni
> sessionInfo()
R version 2.2.1, 2005-12-20, powerpc-apple-darwin7.9.0
attached base packages:
[1] "tools" "methods" "stats" "graphics" "grDevices"
"utils"
[7] "datasets" "base"
other attached packages:
biomaRt RCurl XML convert marray limma Biobase
"1.5.12" "0.6-0" "0.99-6" "1.4.0" "1.8.0" "2.4.7" "1.8.0"
Ciao Giovanni,
can you provide a script that allows us to reproduce this error?
In any case, I think getBM should produce a useful error message
rather
than just quietly returning NULL and we will fix this.
To be even more helpful, you could also say "debug(getBM)" before
calling getBM, and report to us the content of "xmlQuery" and
"postRes"
after the call to the "postForm" function (this is the main action of
the function getBM, when it talks to the Wormbase webservice).
Best wishes
Wolfgang
-------------------------------------
Wolfgang Huber
European Bioinformatics Institute
European Molecular Biology Laboratory
Cambridge CB10 1SD
England
Phone: +44 1223 494642
Fax: +44 1223 494486
Http: www.ebi.ac.uk/huber
-------------------------------------
Giovanni Coppola wrote:
> Hello everybody,
> I am using biomaRt to retrieve worm annotation.
> Sometimes (apparently unpredictably) I get a NULL output and a
warning:
>
> > getBM (attributes=c
> ("wormbase_gene","go_description"),filter="wormbase_gene",values=wor
mids
> ,mart=mart)->gobiomart
> Warning message:
> getBM returns NULL. in: getBM(attributes = c("wormbase_gene",
> "go_description"), filter = "wormbase_gene",
>
> The exact same line works later. Has anyone experienced this? Is
> there a timeout issue or something like that?
>
> Thanks
> Giovanni
>
>
> > sessionInfo()
> R version 2.2.1, 2005-12-20, powerpc-apple-darwin7.9.0
>
> attached base packages:
> [1] "tools" "methods" "stats" "graphics" "grDevices"
"utils"
> [7] "datasets" "base"
>
> other attached packages:
> biomaRt RCurl XML convert marray limma Biobase
> "1.5.12" "0.6-0" "0.99-6" "1.4.0" "1.8.0" "2.4.7" "1.8.0"
Hello Wolfgang,
I am actually getting a different error now... the wormids list is
usually much longer.
> wormids<-c("ZK666.6","Y41C4A.
11","T02G5.7.1","K08F4.7","A_12_P11025","ZC410.5b","C32H11.2","F54D11.
1"
)
> getBM (attributes=c
("wormbase_gene","go_description"),filter="wormbase_gene",values=wormi
ds
,mart=mart)->gobiomart
Error: ncol(result) == length(attributes) is not TRUE
>
Thank you very much
Giovanni
debug(getBM)
[...]
Browse[1]>
debug: xmlQuery = paste("<query virtualschemaname="default" count="0">
<dataset name="",
mart at dataset, "">", sep = "")
Browse[1]>
debug: attributeXML = paste("<attribute name="", attributes, ""/>",
collapse = "", sep = "")
> sessionInfo()
R version 2.1.0, 2005-05-12, powerpc-apple-darwin7.9.0
attached base packages:
[1] "tools" "methods" "stats" "graphics" "grDevices"
"utils"
[7] "datasets" "base"
other attached packages:
biomaRt RCurl XML convert marray limma Biobase
"1.5.12" "0.6-1" "0.99-6" "1.1.7" "1.6.3" "2.0.8" "1.5.12"
>
On Mar 24, 2006, at 10:26 AM, Wolfgang Huber wrote:
> Ciao Giovanni,
>
> can you provide a script that allows us to reproduce this error?
>
> In any case, I think getBM should produce a useful error message
> rather than just quietly returning NULL and we will fix this.
>
> To be even more helpful, you could also say "debug(getBM)" before
> calling getBM, and report to us the content of "xmlQuery" and
> "postRes" after the call to the "postForm" function (this is the
> main action of the function getBM, when it talks to the Wormbase
> webservice).
>
> Best wishes
> Wolfgang
>
>
> -------------------------------------
> Wolfgang Huber
> European Bioinformatics Institute
> European Molecular Biology Laboratory
> Cambridge CB10 1SD
> England
> Phone: +44 1223 494642
> Fax: +44 1223 494486
> Http: www.ebi.ac.uk/huber
> -------------------------------------
>
>
> Giovanni Coppola wrote:
>> Hello everybody,
>> I am using biomaRt to retrieve worm annotation.
>> Sometimes (apparently unpredictably) I get a NULL output and a
>> warning:
>> > getBM (attributes=c
>>
("wormbase_gene","go_description"),filter="wormbase_gene",values=worm
>> ids ,mart=mart)->gobiomart
>> Warning message:
>> getBM returns NULL. in: getBM(attributes = c("wormbase_gene",
>> "go_description"), filter = "wormbase_gene",
>> The exact same line works later. Has anyone experienced this? Is
>> there a timeout issue or something like that?
>> Thanks
>> Giovanni
>> > sessionInfo()
>> R version 2.2.1, 2005-12-20, powerpc-apple-darwin7.9.0
>> attached base packages:
>> [1] "tools" "methods" "stats" "graphics" "grDevices"
>> "utils"
>> [7] "datasets" "base"
>> other attached packages:
>> biomaRt RCurl XML convert marray limma Biobase
>> "1.5.12" "0.6-0" "0.99-6" "1.4.0" "1.8.0" "2.4.7" "1.8.0"
Dear Giovanni,
can you send us a sufficient code example to reproduce this error
please? Otherwise nobody can really do anything about your problem.
I.e. how did you obtain your "mart" object?
Best
Wolfgang
Giovanni Coppola wrote:
> Hello Wolfgang,
> I am actually getting a different error now... the wormids list is
> usually much longer.
>
> > wormids<-c("ZK666.6","Y41C4A.
> 11","T02G5.7.1","K08F4.7","A_12_P11025","ZC410.5b","C32H11.2","F54D1
1.1" )
> > getBM (attributes=c
> ("wormbase_gene","go_description"),filter="wormbase_gene",values=wor
mids
> ,mart=mart)->gobiomart
> Error: ncol(result) == length(attributes) is not TRUE
> >
>
> Thank you very much
> Giovanni
>
>
> debug(getBM)
> [...]
> Browse[1]>
> debug: xmlQuery = paste(" DOCTYPE Query><query virtualschemaname="default" count="0">
> <dataset name="",
> mart at dataset, "">", sep = "")
> Browse[1]>
> debug: attributeXML = paste("<attribute name="", attributes,
""/>",
> collapse = "", sep = "")
>
> > sessionInfo()
> R version 2.1.0, 2005-05-12, powerpc-apple-darwin7.9.0
>
> attached base packages:
> [1] "tools" "methods" "stats" "graphics" "grDevices"
"utils"
> [7] "datasets" "base"
>
> other attached packages:
> biomaRt RCurl XML convert marray limma Biobase
> "1.5.12" "0.6-1" "0.99-6" "1.1.7" "1.6.3" "2.0.8" "1.5.12"
> >
>
>
> On Mar 24, 2006, at 10:26 AM, Wolfgang Huber wrote:
>
>> Ciao Giovanni,
>>
>> can you provide a script that allows us to reproduce this error?
>>
>> In any case, I think getBM should produce a useful error message
>> rather than just quietly returning NULL and we will fix this.
>>
>> To be even more helpful, you could also say "debug(getBM)" before
>> calling getBM, and report to us the content of "xmlQuery" and
>> "postRes" after the call to the "postForm" function (this is the
main
>> action of the function getBM, when it talks to the Wormbase
webservice).
>>
>> Best wishes
>> Wolfgang
>>
>>
>> -------------------------------------
>> Wolfgang Huber
>> European Bioinformatics Institute
>> European Molecular Biology Laboratory
>> Cambridge CB10 1SD
>> England
>> Phone: +44 1223 494642
>> Fax: +44 1223 494486
>> Http: www.ebi.ac.uk/huber
>> -------------------------------------
>>
>>
>> Giovanni Coppola wrote:
>>
>>> Hello everybody,
>>> I am using biomaRt to retrieve worm annotation.
>>> Sometimes (apparently unpredictably) I get a NULL output and a
warning:
>>> > getBM (attributes=c
>>>
("wormbase_gene","go_description"),filter="wormbase_gene",values=worm
>>> ids ,mart=mart)->gobiomart
>>> Warning message:
>>> getBM returns NULL. in: getBM(attributes = c("wormbase_gene",
>>> "go_description"), filter = "wormbase_gene",
>>> The exact same line works later. Has anyone experienced this? Is
>>> there a timeout issue or something like that?
>>> Thanks
>>> Giovanni
>>> > sessionInfo()
>>> R version 2.2.1, 2005-12-20, powerpc-apple-darwin7.9.0
>>> attached base packages:
>>> [1] "tools" "methods" "stats" "graphics" "grDevices"
"utils"
>>> [7] "datasets" "base"
>>> other attached packages:
>>> biomaRt RCurl XML convert marray limma Biobase
>>> "1.5.12" "0.6-0" "0.99-6" "1.4.0" "1.8.0" "2.4.7" "1.8.0"
>
>
--
Best regards
Wolfgang
-------------------------------------
Wolfgang Huber
European Bioinformatics Institute
European Molecular Biology Laboratory
Cambridge CB10 1SD
England
Phone: +44 1223 494642
Fax: +44 1223 494486
Http: www.ebi.ac.uk/huber
Hi Giovanni,
> sorry Wolfgang, here it is.
> BTW, I just tried and it works today.
library(biomaRt)
mart=useMart("ensembl")
mart=useDataset(dataset = "celegans_gene_ensembl", mart = mart)
wormids=c("ZK666.6","Y41C4A.11","T02G5.7.1","K08F4.7","A_12_P11025","Z
C410.5b","C32H11.2","F54D11.1")
gbm=getBM(attributes=c("wormbase_gene","go_description"),filter="wormb
ase_gene",values=wormids,mart=mart)
thanks. I also just tried this just now and it worked fine, but I
understand there may have been errors previously. It appears that
there
are sporadic problems in the interaction between the biomaRt R package
and the biomart webservices at Ensembl and/or Wormbase, it is not
clear
at this point whether the error is on the side of the biomaRt package
or
of the webservice provider, or whether it is a miscommunication
between
them. Any hints from users that experience such a problem (including
enough info to reproduce it) is welcome!
Cheers
Wolfgang