[BioC] Error in summary(hyperGtest)
Fraser Sim
fjsim at buffalo.edu
Mon Jul 13 21:28:48 CEST 2009
I am using bioconductor to perform a hyperGtest on GO BP annotations. This
code has worked in the past but I am getting a strange error message from
the summary method of a hyperGtest-results object.
Here's the code:
> params = new("GOHyperGParams",
+ geneIds = as.numeric(geneIds),
+ universeGeneIds = as.numeric(GeneIDUniverse), annotation =
"hgu133plus2",
+ ontology = "BP", pvalueCutoff = 0.001, conditional = FALSE,
+ testDirection = "over")
> x <- hyperGTest(params)
> x
Gene to GO BP test for over-representation
652 GO BP ids tested (22 have p < 0.001)
Selected gene set size: 74
Gene universe size: 11847
Annotation package: hgu133plus2
> summary(x)
Error in if (nrow(df) == 0) { : argument is of length zero
Why does summary have this error when it looks like there are 22 categories
which should be returned?
Here's my sessionInfo()
R version 2.9.1 (2009-06-26)
i386-pc-mingw32
locale:
LC_COLLATE=English_United States.1252;LC_CTYPE=English_United
States.1252;LC_MONETARY=English_United
States.1252;LC_NUMERIC=C;LC_TIME=English_United States.1252
attached base packages:
[1] stats graphics grDevices datasets utils methods base
other attached packages:
[1] GOstats_2.10.0 graph_1.22.2 Category_2.10.1
[4] hgu133plus2.db_2.2.11 annaffy_1.16.0 KEGG.db_2.2.11
[7] GO.db_2.2.11 RSQLite_0.7-1 DBI_0.2-4
[10] AnnotationDbi_1.6.1 limma_2.18.2 Biobase_2.4.1
loaded via a namespace (and not attached):
[1] annotate_1.22.0 genefilter_1.24.2 GSEABase_1.6.0 RBGL_1.20.0
[5] splines_2.9.1 survival_2.35-4 tools_2.9.1 XML_2.3-0
[9] xtable_1.5-5
Thanks,
Fraser
More information about the Bioconductor
mailing list