I have a list of about 800 KEGG IDs for my genes of interest. I want to plot those genes and get a nice KEGG mapped plot for each of my 800 genes. How can I do this in R package? Thanks
You may use the `pathview` package:

<s>1) use `kegg.gsets()` to download the KEGG pathways of your species.

2) map your gene ids to the pathways - I use `ids2indices()` from limma.</s>

3) use `pathview()` to download KEGG pathways figures and color your genes on them.

Example:

    library(pathview)
    test.uid <- c("5494737" , "5495078", "5495093", "5494418", "5495039")
    pv.out <- pathview(gene.data = test.uid, pathway.id = "ssl00190", species = "ssl", 
                       out.suffix = "kegg.get.all", kegg.native = T)