[BioC] GOstats and org.EcK12.eg.db

Robert Castelo robert.castelo at upf.edu
Wed Oct 22 15:48:22 CEST 2008


dear list,

I cannot get to work GOstats with the annotation for E. coli in
org.EcK12.eg.db. Please find below the code that reproduces the problem
including the error message, and my sessionInfo() at the end of this
email. I have included the same exercise with the human annotation
package org.Hs.eg.db which runs fine in my system. Any help with this
will be very much appreciated.

thanks!!
robert.
==========CODE STARTS HERE===========

library(org.Hs.eg.db)
library(org.EcK12.eg.db)
library(GOstats)

geneuniverse <- mappedkeys(org.EcK12.egSYMBOL)
set.seed(12345)
geneset <- sample(geneuniverse, size=100, replace=FALSE)


goHypGparams <- new("GOHyperGParams",
                    geneIds=geneset,
                    universeGeneIds=geneuniverse,
                    annotation="org.EcK12.eg.db", ontology="BP",
                    pvalueCutoff=1.0, conditional=TRUE,
                    testDirection="over")

goHypGcond <- hyperGTest(goHypGparams)

Error in get(mapName, envir = pkgEnv, inherits = FALSE) : 
  variable "org.EcK12.egENTREZID" was not found
Error in mget(probes, ID2EntrezID(datPkg)) : 
  error in evaluating the argument 'envir' in selecting a method for
function 'mget'

geneuniverse <- mappedkeys(org.Hs.egSYMBOL)
set.seed(12345)
geneset <- sample(geneuniverse, size=100, replace=FALSE)


goHypGparams <- new("GOHyperGParams",
                    geneIds=geneset,
                    universeGeneIds=geneuniverse,
                    annotation="org.Hs.eg.db", ontology="BP",
                    pvalueCutoff=1.0, conditional=TRUE,
                    testDirection="over")

goHypGcond <- hyperGTest(goHypGparams)


sessionInfo()
R version 2.8.0 beta (2008-10-05 r46601) 
x86_64-unknown-linux-gnu 

locale:
LC_CTYPE=en_US.UTF-8;LC_NUMERIC=C;LC_TIME=en_US.UTF-8;LC_COLLATE=en_US.UTF-8;LC_MONETARY=C;LC_MESSAGES=en_US.UTF-8;LC_PAPER=en_US.UTF-8;LC_NAME=C;LC_ADDRESS=C;LC_TELEPHONE=C;LC_MEASUREMENT=en_US.UTF-8;LC_IDENTIFICATION=C

attached base packages:
[1] splines   tools     stats     graphics  grDevices utils
datasets 
[8] methods   base     

other attached packages:
 [1] org.Hs.eg.db_2.2.6    GOstats_2.7.0         Category_2.7.6       
 [4] genefilter_1.21.5     survival_2.34-1       RBGL_1.17.2          
 [7] annotate_1.19.2       xtable_1.5-4          GO.db_2.2.5          
[10] graph_1.19.6          org.EcK12.eg.db_2.2.6 AnnotationDbi_1.3.12 
[13] RSQLite_0.7-0         DBI_0.2-4             Biobase_2.1.7        

loaded via a namespace (and not attached):
[1] cluster_1.11.11 GSEABase_1.3.6  XML_1.98-1



More information about the Bioconductor mailing list