[BioC] GOstats After Running Lumi

Shawn Westaway westaway at ohsu.edu
Fri May 4 03:29:52 CEST 2007


Using Lumi on Illumina Mouse-Ref6 chip microarray data. I tried to use
GOstats to generate a nice table like the one shown in lumi.pdf using
the following script

> if (require(GOstats) & require(lumiMouseV1)) {
sigLL <- unique(unlist(mget(sigGene, env=lumiMouseV1ENTREZID,
ifnotfound=NA)))
sigLL <- as.character(sigLL[!is.na(sigLL)])
params <- new("GOHyperGParams",
geneIds= sigLL,
annotation="lumiMouseV1",
ontology="BP",
pvalueCutoff= 0.01,
conditional=FALSE,
testDirection="over")
hgOver <- hyperGTest(params)
gGhyp.pv <- pvalues(hgOver)
sigGO.ID <- names(gGhyp.pv[gGhyp.pv < 0.001])
sigGO.Term <- getGOTerm(sigGO.ID)[["BP"]]
}

The script does not seem to output anything. Is there a command missing
to generate the data?
Thanks,
Shawn



More information about the Bioconductor mailing list