hyperGTest on KEGG and PFAM with org.XX.eg annotations
0
0
Entering edit mode
Marc Carlson ★ 7.2k
@marc-carlson-2264
Last seen 5.5 years ago
United States
Hi James, The current release and devel versions of Category have been patched. Barring any incident, these should now work more reliably and should also represent things more consistently internally. The latest versions of Category to look for are: 2.9.8 for devel and 2.8.2 for release. Thank you for reporting the bug. Marc James F. Reid wrote: > Dear list, > > hyperGTest behaves differently when using org.XX.eg.db packages > compared to microarray based ones, like hgu95av2.db for example, for > doing a KEGG analysis. hyperGTest complains if the annotation string > does not end with the suffix ".db", it works if you add it but then > you can't run a summary on the result. A quick fix is to re-assign the > ".db"-less string to the annotation slot of the hyperGTest result. > So I am wondering if I am doing something wrong of if it is a bug. > > For the PFAM analysis everything works fine except that in the summary > output the Term (Description) is just the PFAMID which is not very > useful for interpretation. I think this could easily be fixed by using > the same approach as for the KEGG output in the PFAMHyperGResult > summary method: > ## implicit require("PFAM.db") > pfamEnv <- getAnnMap("DE", "PFAM", load=TRUE) > pfamTerms <- unlist(mget(pfamIds, pfamEnv, ifnotfound=NA)) > > > Many thanks, > James. > > Here is a session reporting the problem: > > library("Category") > library("org.Hs.eg.db") > > set.seed(123) > geneBackground <- Lkeys(org.Hs.egPATH) > geneList <- sample(geneBackground, 500) > > params <- new("KEGGHyperGParams", > geneIds = geneList, > universeGeneIds = geneBackground, > annotation = "org.Hs.eg") > hgKEGG <- hyperGTest(params) > # Error in get(paste(lib, name, sep = "")) : > # variable "org.Hs.egPATH2PROBE" was not found > > params at annotation <- "org.Hs.eg.db" > hgKEGG <- hyperGTest(params) > summary(hgKEGG) > # Error in get(paste(annotation(object), "ORGANISM", sep = "")) : > # variable "org.Hs.eg.dbORGANISM" was not found > > hgKEGG at annotation <- "org.Hs.eg" > summary(hgKEGG) > # KEGGID Pvalue OddsRatio ExpCount Count Size > #1 05130 0.003282103 7.314332 0.6239536 4 51 > #2 05131 0.003282103 7.314332 0.6239536 4 51 > # Term > #1 Pathogenic Escherichia coli infection - EHEC > #2 Pathogenic Escherichia coli infection - EPEC > > > > sessionInfo() > R version 2.8.0 (2008-10-20) > i486-pc-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_IDEN TIFICATION=C > > > attached base packages: > [1] splines tools stats graphics grDevices utils datasets > [8] methods base > > other attached packages: > [1] KEGG.db_2.2.5 org.Hs.eg.db_2.2.6 RSQLite_0.7-1 > [4] DBI_0.2-4 Category_2.8.1 genefilter_1.22.0 > [7] survival_2.34-1 annotate_1.20.1 xtable_1.5-4 > [10] AnnotationDbi_1.4.1 graph_1.20.0 Biobase_2.2.1 > > loaded via a namespace (and not attached): > [1] cluster_1.11.11 GSEABase_1.4.0 RBGL_1.18.0 XML_1.98-1 > > _______________________________________________ > Bioconductor mailing list > Bioconductor at stat.math.ethz.ch > https://stat.ethz.ch/mailman/listinfo/bioconductor > Search the archives: > http://news.gmane.org/gmane.science.biology.informatics.conductor >
Microarray Annotation Escherichia coli hgu95av2 Category Microarray Annotation hgu95av2 • 489 views
ADD COMMENT

Login before adding your answer.

Traffic: 455 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6