Question: Retrieving KEGG pathway diagrams via KEGGSOAP
0
gravatar for Daren Tan
10.1 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
keggsoap • 565 views
ADD COMMENTlink modified 10.1 years ago by Paolo Sonego140 • written 10.1 years ago by Daren Tan120
Answer: Retrieving KEGG pathway diagrams via KEGGSOAP
0
gravatar for Saroj Mohapatra
10.1 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 10.1 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 10.1 years ago by Daren Tan120
Answer: Retrieving KEGG pathway diagrams via KEGGSOAP
0
gravatar for Paolo Sonego
10.1 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 10.1 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 16.09
Traffic: 211 users visited in the last hour