Cannot download org.Hs.egPATH
3
0
Entering edit mode
ladypurrsia ▴ 10
@ladypurrsia-11377
Last seen 3.1 years ago
United States

Hello:

I have a list of Kegg Orthology (K) numbers that I would like to covert from K to entrezID to use clusterProfiler package. Currently, I am trying to install the package org.Hs.egPATH to use with the bitr command from clusterProfiler. My inquiry is this:

I am getting the following error after I run the following command: source("https://bioconductor.org/biocLite.R")
biocLite("org.Hs.egPATH”)

1: package ‘org.Hs.egPATH’ is not available (for R version 3.3.1)

2: In install.packages(update[instlib == l, "Package"], l, repos = repos,  :
  installation of package ‘Matrix’ had non-zero exit status
3: In install.packages(update[instlib == l, "Package"], l, repos = repos,  :
  installation of package ‘rgl’ had non-zero exit status

 

I was able to download the org.Hs.eg.db package without a problem. Any recommendations as to when it will be available? Is there a work-around on how to obtain the org.Hs.egPATH package? 

Thank you very much!

Joany 

org.Hs.egPATH • 1.8k views
ADD COMMENT
0
Entering edit mode

Thank you so very much, Martin!

Joany

ADD REPLY
0
Entering edit mode

pls tag the post with clusterprofiler.
 

ADD REPLY
1
Entering edit mode
@martin-morgan-1513
Last seen 4 months ago
United States

'org.Hs.egPATH' is not a package, but a variable in the org.Hs.eg.db package, so

biocLite("org.Hs.eg.db")
library(org.Hs.eg.db)
org.Hs.egPATH

but unfortunately the information in this object is now several years old (KEGG no longer freely distributes the aggregate information. You might try the KEGGREST package instead.

ADD COMMENT
0
Entering edit mode

Is the KEGGREST package needed? It seems to me that one can read the info directly from the KEGG website without needing an intermediary package:

EG.KEGG <- read.table("http://rest.kegg.jp/link/pathway/hsa",stringsAsFactors=FALSE)
ADD REPLY
1
Entering edit mode
@gordon-smyth
Last seen 2 hours ago
WEHI, Melbourne, Australia

You can read the latest Entrez ID to KEGG links from the following URL at any time:

  http://rest.kegg.jp/link/pathway/hsa

The getGeneKEGGLinks() function in the limma package automates this for you:

> library(limma)
> EG.KEGG <- getGeneKEGGLinks("hsa")
> head(EG.KEGG)
  GeneID     PathwayID
1  10327 path:hsa00010
2    124 path:hsa00010
3    125 path:hsa00010
4    126 path:hsa00010
5    127 path:hsa00010
6    128 path:hsa00010

The first column in Entrez ID and the second is KEGG pathway ID.

However I rather doubt that you need to do any ID conversion to use the clusterProfiler package. I'm pretty sure that the clusterProfiler package knows how to access the KEGG website and that it works fine with Entrez IDs. The bitr() function is intended to convert between gene ID systems, not between gene IDs and pathway IDs. If you have Entrez IDs then you're already ready to go, because KEGG is Entrez Id based.

ADD COMMENT
0
Entering edit mode

He doesn't have entrezgene ID and only have a list of K number (very curious and I also have no idea).

ADD REPLY
0
Entering edit mode
Guangchuang Yu ★ 1.2k
@guangchuang-yu-5419
Last seen 26 days ago
China/Guangzhou/Southern Medical Univer…

1. bitr doesn't support converting K number to gene ID.

2. You can use enrichKEGG with a list of K number as input gene, see also #clusterProfiler# KEGG enrichment with non-reference species

 

ADD COMMENT

Login before adding your answer.

Traffic: 650 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6