Search
Question: Retrieving KEGG pathway diagrams via KEGGSOAP
0
gravatar for Daren Tan
9.6 years ago by
Daren Tan120
Daren Tan120 wrote:
There was an earlier example about KEGGSOAP, which I am looking for all this while. My work centers around hsa, so I tried one of the pathway annotated in human and sound genes. However, NULL is always returned. What is wrong ? > get.genes.by.pathway("path:hsa03050") [1] "hsa:100132108" "hsa:10197" "hsa:10213" "hsa:121906" "hsa:122706" "hsa:143471" "hsa:23198" "hsa:345645" "hsa:3458" [10] "hsa:51371" "hsa:5682" "hsa:5683" "hsa:5684" "hsa:5685" "hsa:5686" "hsa:5687" "hsa:5688" "hsa:5689" [19] "hsa:5690" "hsa:5691" "hsa:5692" "hsa:5693" "hsa:5694" "hsa:5695" "hsa:5696" "hsa:5698" "hsa:5699" [28] "hsa:5700" "hsa:5701" "hsa:5702" "hsa:5704" "hsa:5705" "hsa:5706" "hsa:5707" "hsa:5708" "hsa:5709" [37] "hsa:5710" "hsa:5713" "hsa:5714" "hsa:5717" "hsa:5718" "hsa:5719" "hsa:5720" "hsa:5721" "hsa:652826" [46] "hsa:7979" "hsa:9491" "hsa:9861" > url <- mark.pathway.by.objects("path:hsa03050", "hsa:10197") > url NULL > url <- mark.pathway.by.objects("path:hsa03050", "hsa:5684") > url NULL > sessionInfo() R version 2.8.1 (2008-12-22) i386-pc-mingw32 locale: LC_COLLATE=English_United States.1252;LC_CTYPE=English_United States.1252;LC_MONETARY=English_United States.1252;LC_NUMERIC=C;LC_TIME=English_United States.1252 attached base packages: [1] stats graphics grDevices utils datasets methods base other attached packages: [1] KEGGSOAP_1.16.0 SSOAP_0.4-6 RCurl_0.94-1 XML_2.1-0 limma_2.16.4 RSQLite_0.7-1 DBI_0.2-4 loaded via a namespace (and not attached): [1] gdata_2.4.2 gplots_2.6.0 gtools_2.5.0-1 tools_2.8.1
ADD COMMENTlink modified 9.6 years ago by Paolo Sonego140 • written 9.6 years ago by Daren Tan120
0
gravatar for Saroj Mohapatra
9.6 years ago by
Saroj Mohapatra450 wrote:
Hi Daren, I do not know why but if you try with more than one node, the code works. > url <- mark.pathway.by.objects("path:hsa03050", c("hsa:5709", "hsa:10197")) > url [1] "http://soap.genome.jp/tmp/mark_pathway_www2_api15820/hsa03050.gif" > if(interactive()){ browseURL(url) } Best, Saroj ----- Original Message ----- From: "Daren Tan" <darentan76@gmail.com> To: bioconductor at stat.math.ethz.ch Sent: Thursday, April 16, 2009 8:46:30 PM GMT -05:00 US/Canada Eastern Subject: [BioC] Retrieving KEGG pathway diagrams via KEGGSOAP There was an earlier example about KEGGSOAP, which I am looking for all this while. My work centers around hsa, so I tried one of the pathway annotated in human and sound genes. However, NULL is always returned. What is wrong ? > get.genes.by.pathway("path:hsa03050") [1] "hsa:100132108" "hsa:10197" "hsa:10213" "hsa:121906" "hsa:122706" "hsa:143471" "hsa:23198" "hsa:345645" "hsa:3458" [10] "hsa:51371" "hsa:5682" "hsa:5683" "hsa:5684" "hsa:5685" "hsa:5686" "hsa:5687" "hsa:5688" "hsa:5689" [19] "hsa:5690" "hsa:5691" "hsa:5692" "hsa:5693" "hsa:5694" "hsa:5695" "hsa:5696" "hsa:5698" "hsa:5699" [28] "hsa:5700" "hsa:5701" "hsa:5702" "hsa:5704" "hsa:5705" "hsa:5706" "hsa:5707" "hsa:5708" "hsa:5709" [37] "hsa:5710" "hsa:5713" "hsa:5714" "hsa:5717" "hsa:5718" "hsa:5719" "hsa:5720" "hsa:5721" "hsa:652826" [46] "hsa:7979" "hsa:9491" "hsa:9861" > url <- mark.pathway.by.objects("path:hsa03050", "hsa:10197") > url NULL > url <- mark.pathway.by.objects("path:hsa03050", "hsa:5684") > url NULL > sessionInfo() R version 2.8.1 (2008-12-22) i386-pc-mingw32 locale: LC_COLLATE=English_United States.1252;LC_CTYPE=English_United States.1252;LC_MONETARY=English_United States.1252;LC_NUMERIC=C;LC_TIME=English_United States.1252 attached base packages: [1] stats graphics grDevices utils datasets methods base other attached packages: [1] KEGGSOAP_1.16.0 SSOAP_0.4-6 RCurl_0.94-1 XML_2.1-0 limma_2.16.4 RSQLite_0.7-1 DBI_0.2-4 loaded via a namespace (and not attached): [1] gdata_2.4.2 gplots_2.6.0 gtools_2.5.0-1 tools_2.8.1 _______________________________________________ Bioconductor mailing list Bioconductor at stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/bioconductor Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor -- Computational Biology VBI @ Virginia Tech
ADD COMMENTlink written 9.6 years ago by Saroj Mohapatra450
Grateful for the help. Now I can get some nice KEGG pictures. On hand, I have about 20 pathways for visualization, how can I download the 20 gifs or possibly higher resolution pictures and collate them into a pdf ? Possible to do all these steps via R ? On Fri, Apr 17, 2009 at 10:02 AM, Saroj K Mohapatra <smohapat at="" vbi.vt.edu=""> wrote: > Hi Daren, > > I do not know why but if you try with more than one node, the code works. > >> url <- mark.pathway.by.objects("path:hsa03050", > ?c("hsa:5709", "hsa:10197")) > >> url > [1] "http://soap.genome.jp/tmp/mark_pathway_www2_api15820/hsa03050.gif" > >> if(interactive()){ browseURL(url) } > > Best, > > Saroj > > > ----- Original Message ----- > From: "Daren Tan" <darentan76 at="" gmail.com=""> > To: bioconductor at stat.math.ethz.ch > Sent: Thursday, April 16, 2009 8:46:30 PM GMT -05:00 US/Canada Eastern > Subject: [BioC] Retrieving KEGG pathway diagrams via KEGGSOAP > > There was an earlier example about KEGGSOAP, which I am looking for > all this while. My work centers around hsa, so I tried one of the > pathway annotated in human and sound genes. However, NULL is always > returned. What is wrong ? > >> get.genes.by.pathway("path:hsa03050") > ?[1] "hsa:100132108" "hsa:10197" ? ? "hsa:10213" ? ? "hsa:121906" > "hsa:122706" ? ?"hsa:143471" ? ?"hsa:23198" ? ? "hsa:345645" > "hsa:3458" > [10] "hsa:51371" ? ? "hsa:5682" ? ? ?"hsa:5683" ? ? ?"hsa:5684" > "hsa:5685" ? ? ?"hsa:5686" ? ? ?"hsa:5687" ? ? ?"hsa:5688" > "hsa:5689" > [19] "hsa:5690" ? ? ?"hsa:5691" ? ? ?"hsa:5692" ? ? ?"hsa:5693" > "hsa:5694" ? ? ?"hsa:5695" ? ? ?"hsa:5696" ? ? ?"hsa:5698" > "hsa:5699" > [28] "hsa:5700" ? ? ?"hsa:5701" ? ? ?"hsa:5702" ? ? ?"hsa:5704" > "hsa:5705" ? ? ?"hsa:5706" ? ? ?"hsa:5707" ? ? ?"hsa:5708" > "hsa:5709" > [37] "hsa:5710" ? ? ?"hsa:5713" ? ? ?"hsa:5714" ? ? ?"hsa:5717" > "hsa:5718" ? ? ?"hsa:5719" ? ? ?"hsa:5720" ? ? ?"hsa:5721" > "hsa:652826" > [46] "hsa:7979" ? ? ?"hsa:9491" ? ? ?"hsa:9861" >> url <- mark.pathway.by.objects("path:hsa03050", "hsa:10197") >> url > NULL >> url <- mark.pathway.by.objects("path:hsa03050", "hsa:5684") >> url > NULL > >> sessionInfo() > R version 2.8.1 (2008-12-22) > i386-pc-mingw32 > > locale: > LC_COLLATE=English_United States.1252;LC_CTYPE=English_United > States.1252;LC_MONETARY=English_United > States.1252;LC_NUMERIC=C;LC_TIME=English_United States.1252 > > attached base packages: > [1] stats ? ? graphics ?grDevices utils ? ? datasets ?methods ? base > > other attached packages: > [1] KEGGSOAP_1.16.0 SSOAP_0.4-6 ? ? RCurl_0.94-1 ? ?XML_2.1-0 > limma_2.16.4 ? ?RSQLite_0.7-1 ? DBI_0.2-4 > > loaded via a namespace (and not attached): > [1] gdata_2.4.2 ? ?gplots_2.6.0 ? gtools_2.5.0-1 tools_2.8.1 > > _______________________________________________ > Bioconductor mailing list > Bioconductor at stat.math.ethz.ch > https://stat.ethz.ch/mailman/listinfo/bioconductor > Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor > > -- > Computational Biology > VBI @ Virginia Tech >
ADD REPLYlink written 9.6 years ago by Daren Tan120
0
gravatar for Paolo Sonego
9.6 years ago by
Paolo Sonego140
Italy
Paolo Sonego140 wrote:
Dear Daren, If you are in a Unix environment (Linux/Mac OS X), below there is a trivial example taking advantage of EBImage package and ImageMagick toolkit (sorry I can figure out how to use EBImage to combine the different images in a single pdf). library(KEGGSOAP) url <- mark.pathway.by.objects("path:hsa03050", c("hsa:5709", "hsa:10197")) urls <- c(url, url, url) ## set different urls here library(EBImage) im.gif <- NULL for (i in 1:length(urls)){ im.gif <- readImage(urls[i], TrueColor) writeImage(im.gif, paste("pathway_", i, ".pdf", sep="")) } to combine the files in a single pdf: system("convert pathway_* pathways.pdf") HIH. Paolo
ADD COMMENTlink written 9.6 years ago by Paolo Sonego140
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.2.0
Traffic: 405 users visited in the last hour