Question: Map Gene names/IDs on Pathway with KEGG
0
gravatar for Shelle
7 months ago by
Shelle0
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. https://www.bioconductor.org/packages/devel/bioc/vignettes/KEGGgraph/inst/doc/KEGGgraph.pdf

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 • 290 views
ADD COMMENTlink modified 7 months ago by h.mon24k • written 7 months ago by Shelle0
0
gravatar for h.mon
7 months ago by
h.mon24k
Brazil
h.mon24k 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 7 months ago • written 7 months ago by h.mon24k
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.3.0
Traffic: 1110 users visited in the last hour