[BioC] Retrieving KEGG pathway diagrams via KEGGSOAP
Paolo Sonego
paolo.sonego at gmail.com
Fri Apr 17 16:52:20 CEST 2009
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
More information about the Bioconductor
mailing list