Question: Map Gene names/IDs on Pathway with KEGG
gravatar for Shelle
17 months ago by
Shelle0 wrote:

I want to map the list of gene IDs that i have to pathway and see if there is a pathway between the genes. I decided to go with KEGG using the code below to see what can be the name of nodes given the pathway id but for some reason in one line it gives me error so i am not able to run the code in the documentation website.

Does anyone know how this problem can be solved?

> library(KEGG.db)
> tmp <- tempfile()
> pName <- "Metabolic pathways"
> pId <- mget(pName, KEGGPATHNAME2ID)[[1]]
> retrieveKGML(pId, organism="cel", destfile=tmp, method="wget", quiet=TRUE)
> mapkKGML <- system.file("extdata/hsa01100.xml", package="KEGGgraph")

# Error was in the line above: 
Error in if (fileSize < 100L) msg <- paste(msg, "[WARNING] File size (",  :
  missing value where TRUE/FALSE needed

> mapkG <- parseKGML2Graph(mapkKGML,expandGenes=TRUE)
> mapkNodes <- nodes(mapkG)
> nodes(mapkG)[1:3]

The only change here comparing to the documentation is the pathway identifier i have also tried their exact example and same issue.

kegg gene pathway • 562 views
ADD COMMENTlink modified 17 months ago by h.mon29k • written 17 months ago by Shelle0
gravatar for h.mon
17 months ago by
h.mon29k wrote:

A quick careful reading of the KEGGgraph manual will tell you the cause of the failure, which can be confirmed with:

list.files( system.file( "extdata/", package = "KEGGgraph" ) )
ADD COMMENTlink modified 17 months ago • written 17 months ago by h.mon29k
Please log in to add an answer.


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