getBM (biomaRt) sometimes returns NULL
1
0
Entering edit mode
@giovanni-coppola-893
Last seen 10.2 years ago
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"
Annotation limma marray convert biomaRt Annotation limma marray convert biomaRt • 1.7k views
ADD COMMENT
0
Entering edit mode
@wolfgang-huber-3550
Last seen 3 months ago
EMBL European Molecular Biology Laborat…
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"
ADD COMMENT
0
Entering edit mode
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"
ADD REPLY
0
Entering edit mode
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="", &gt; 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
ADD REPLY
0
Entering edit mode
An embedded and charset-unspecified text was scrubbed... Name: not available Url: https://stat.ethz.ch/pipermail/bioconductor/attachments/20060330/ 1befd947/attachment.pl
ADD REPLY
0
Entering edit mode
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
ADD REPLY

Login before adding your answer.

Traffic: 617 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6