[BioC] hyperGTest, KEGG

Nianhua Li nli at fhcrc.org
Sat Apr 14 02:27:16 CEST 2007


Hi, Ivan,

library(KEGG, lib="~/annoprep/bin-13")
x <- c("04360", "05214", "01030")
## Assume Hu19K8Build17102006 is the annotation for human genes,
## you prefix the KEGG ID with human KEGG pathway code "hsa"
x <- paste("hsa", x, sep="")
g <- mget(x, KEGGPATHID2EXTID)
## then maybe you want to intersect g with selectedEntrezIds ...

hope this helps

nianhua



More information about the Bioconductor mailing list