[BioC] Error in package goseq
Joydeep Mitra
joy.mitra at yahoo.com
Tue Aug 17 20:52:04 CEST 2010
Hi,
I am a grad student at Institute of Bioinformatics, University of Georgia, USA.
I have been trying to use the bioconductor package goseq for go-term/pathway
enrichment of my expression data.
Since my organism of interest isn't included in the goseq database, I am
manually entering the data for gene lengths and categories following the
directions in the vignette provided with the package. I have tried representing
the gene-category mappings in both data frame and list formats as specified.
While the 'nullp' function worked fine for my data, when I run goseq() with the
following command, I get an error message:
> eco1.gs<-goseq(DEgenes=eco1.de,pwf=eco1.pwf,gene2cat=ids.go.list)
Using manually entered categories.
Error in summary(map)[, 1] : incorrect number of dimensions
Here's my session info:
> sessionInfo()
R version 2.11.1 (2010-05-31)
i386-apple-darwin9.8.0
locale:
[1] en_US.UTF-8/en_US.UTF-8/C/C/en_US.UTF-8/en_US.UTF-8
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
[1] goseq_1.0.3 geneLenDataBase_0.99.4 BiasedUrn_1.03
loaded via a namespace (and not attached):
[1] AnnotationDbi_1.10.2 Biobase_2.8.0 biomaRt_2.4.0
Biostrings_2.16.9 BSgenome_1.16.5 DBI_0.2-5
[7] GenomicFeatures_1.0.6 GenomicRanges_1.0.6 grid_2.11.1
IRanges_1.6.8 lattice_0.18-8 Matrix_0.999375-39
[13] mgcv_1.6-2 nlme_3.1-96 RCurl_1.4-2
RSQLite_0.9-1 rtracklayer_1.8.1 tools_2.11.1
[19] XML_3.1-0
Has anyone else faced this problem in goseq? Any help is appreciated.
Thanks,
Joydeep
More information about the Bioconductor
mailing list