[BioC] KEGG: gene ids for nodes in a pathway

Marc Carlson mcarlson at fhcrc.org
Tue May 5 00:05:52 CEST 2009


Hi Tim,

I think that the mapping you are using below already maps the entrez
gene IDs associated with a particular pathway.  All you need to do is
use mget() instead of toTable().

So for pathway "04630", you can just get the associated entrez gene IDs
like this:

library(org.Hs.eg.db)
mget("04630", revmap(org.Hs.egPATH), ifnotfound=NA)


  Marc





Tim Smith wrote:
> Hi,
>
> I wanted a list of genes for a particular pathway arranged nodewise. For example, if I select the Jak-stat pathway ("http://www.genome.jp/kegg/pathway/hsa/hsa04630.html"), how do I get the entrez ids of genes associated with the node 'STAT' ? Currently, I use the following code:
>
> x <- toTable(org.Hs.egPATH)
>
> and then select genes associated with a particular pathway (e.g. for Jak-stat: "04630") . But this gives the entire set of genes associated with the pathway. Is there a way to get the entrez ids of the genes associated with each of the nodes ('JAK', 'STAT', 'STAM','PIAS' etc.) in the pathway?
>
> thanks!
>
>
>
>       
> 	[[alternative HTML version deleted]]
>
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/bioconductor
> Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor
>
>



More information about the Bioconductor mailing list