[BioC] R: Pathway Information

John Zhang jzhang at jimmy.harvard.edu
Thu Jun 16 17:55:42 CEST 2005

Sorry for the typo. steriods should be steroids.

>X-Original-To: jzhang at jimmy.harvard.edu
>Delivered-To: jzhang at jimmy.harvard.edu
>From: "Alberto Goldoni" <alberto.goldoni at eurogene.org>
>To: "'John Zhang'" <jzhang at jimmy.harvard.edu>, <michael.auer at meduniwien.ac.at>, 
<rhaddad at genetics.wayne.edu>
>Cc: <bioconductor at stat.math.ethz.ch>
>Subject: R: [BioC] Pathway Information
>Date: Thu, 16 Jun 2005 17:43:49 +0200
>MIME-Version: 1.0
>Content-Transfer-Encoding: quoted-printable
>X-Priority: 3 (Normal)
>X-MSMail-Priority: Normal
>X-MimeOLE: Produced By Microsoft MimeOLE V6.00.2800.1506
>Importance: Normal
>X-Virus-Scanned: Cineca AppOs 0.86 at poster.unibo.it
>X-Spam-Checker-Version: SpamAssassin 3.0.1 (2004-10-22) on  
>X-Spam-Status: No, score=-1.3 required=3.0 tests=AWL,BAYES_00 autolearn=ham  
>Sorry but it doesn't work:
>> get(paste("hsa", get("Biosynthesis of steriods", KEGGPATHNAME2ID),
>Error in get(x, envir, mode, inherits) : variable "Biosynthesis of
>steriods" not find
>-----Messaggio originale-----
>Da: bioconductor-bounces at stat.math.ethz.ch
>[mailto:bioconductor-bounces at stat.math.ethz.ch] Per conto di John Zhang
>Inviato: giovedì 16 giugno 2005 16.55
>A: michael.auer at meduniwien.ac.at; rhaddad at genetics.wayne.edu
>Cc: bioconductor at stat.math.ethz.ch
>Oggetto: Re: [BioC] Pathway Information
>Suppose you are looking for gene that are involved in human
>"Biosynthesis of 
>steriods" pathway, the following will give you a vector of Entrez Gene
>>get(paste("hsa", get("Biosynthesis of steriods", KEGGPATHNAME2ID), sep
>= ""), 
>>X-Original-To: jzhang at jimmy.harvard.edu
>>Delivered-To: jzhang at jimmy.harvard.edu
>>From: Ramsi Haddad <rhaddad at genetics.wayne.edu>
>>To: michael.auer at meduniwien.ac.at
>>Mime-Version: 1.0
>>Date: 16 Jun 2005 09:58:13 -0400
>>Received-SPF: pass (hypatia.math.ethz.ch: is
>authenticated by a 
>trusted mechanism)
>>Received-SPF: none (hypatia: domain of rhaddad at genetics.wayne.edu does
>designate permitted sender hosts)
>>X-Virus-Scanned: by amavisd-new at stat.math.ethz.ch
>>Cc: bioconductor mailing list <bioconductor at stat.math.ethz.ch>
>>Subject: Re: [BioC] Pathway Information
>>X-BeenThere: bioconductor at stat.math.ethz.ch
>>X-Mailman-Version: 2.1.6
>>List-Id: The Bioconductor Project Mailing List
>>List-Unsubscribe: <https://stat.ethz.ch/mailman/listinfo/bioconductor>,
><mailto:bioconductor-request at stat.math.ethz.ch?subject=unsubscribe>
>>List-Archive: <https://stat.ethz.ch/pipermail/bioconductor>
>>List-Post: <mailto:bioconductor at stat.math.ethz.ch>
>>List-Help: <mailto:bioconductor-request at stat.math.ethz.ch?subject=help>
>>List-Subscribe: <https://stat.ethz.ch/mailman/listinfo/bioconductor>, 
><mailto:bioconductor-request at stat.math.ethz.ch?subject=subscribe>
>>Content-Transfer-Encoding: 7bit
>>X-Spam-Checker-Version: SpamAssassin 3.0.1 (2004-10-22) on  
>>X-Spam-Status: No, score=-2.2 required=3.0 tests=AWL,BAYES_00
>>Dear Michael,
>>	There is a pretty good description in the globaltest vignette.
>Most of
>>the stuff shown below is based on it:
>>### this next line generates a list of kegg pathways and all
>>### the probe sets that belong in a given pathway.
>>mouse.kegg.sets <- as.list(mgu74av2PATH2PROBE)
>>### the next 3 lines give the names of the pathways from their
>>### kegg ID numbers
>>all.kegg.ids <- ls(KEGGPATHID2NAME)
>>all.kegg.names <- mget(all.kegg.ids, KEGGPATHID2NAME)
>>table.of.keggs <- as.matrix(unlist(all.kegg.names))
>>write.table(table.of.keggs, "keggs.csv", sep=",",col.names=NA,
>>### here's how to get the name of a given kegg pathway:
>>what.is.00010 <- get("00010", KEGGPATHID2NAME) ## this is glycolysis.
>>### here's how to get the probesets in the Glycolysis Path from the
>>mouse array:
>>glycolysis.affyIDs <- get("00010", mgu74av2PATH2PROBE) ## there are 68
>>such probe sets.
>>another.glycolysis.affyIDs <- mouse.kegg.sets[["00010"]] ## more than
>>one way to skin a cat
>>### here's how to get the kegg pathway numbers for anything with
>>colera.kegg.paths.index <- grep("Cholera", table.of.keggs) ## this
>>provides the index
>>colera.keggs <- table.of.keggs[colera.kegg.paths]
>>### here it is in one step:
>>### table.of.keggs[grep("Cholera", table.of.keggs)]
>>hope that gets you started.
>>> Hi
>>> Can anybody tell me the following
>>> I want to identify genes which are located on a certain pathway NF
>>> How can I obtain such an information. Which package is appropriate
>>> which commands are needed. If any body ever encounterd the same
>>> please let me know.
>>> Michael
>>Ramsi Haddad, Ph.D.
>>Center for Molecular Medicine and Genetics.
>>Functional Genomics Laboratory,
>>Perinatology Research Branch, NICHD, NIH.
>>259 Mack Avenue,
>>Room 3146 Applebaum Bldg.
>>Detroit, MI 48201, USA.
>>phone:(313) 577-2569 / fax:(313) 577-7736
>>Bioconductor mailing list
>>Bioconductor at stat.math.ethz.ch
>Jianhua Zhang
>Department of Medical Oncology
>Dana-Farber Cancer Institute
>44 Binney Street
>Boston, MA 02115-6084
>Bioconductor mailing list
>Bioconductor at stat.math.ethz.ch

Jianhua Zhang
Department of Medical Oncology
Dana-Farber Cancer Institute
44 Binney Street
Boston, MA 02115-6084

More information about the Bioconductor mailing list