KEGG pathway analysis with ReportingTools
Entering edit mode
Last seen 7.1 years ago
United States
Dear All, I am successfully using the ReportingTools package to detect over- expressed categories using hypergeometric Test for Gene Ontologies and PFAM follwoing the package vignette. When I tried to do the same using KEGG pathways, the function publish does not work. I am surprise because class PFAMHyperGResult and KEGGHyperGResult have the same structure. Can anyone shed light onto this? Mayte ### CODE > keggParams<-new("KEGGHyperGParams", geneIds = DEGi_entrez, universeGeneIds=UnivEntrez, + annotation='hgu133plus2.db', pvalueCutoff = 0.1, testDirection = "over") Warning in makeValidParams(.Object) : removing duplicate IDs in universeGeneIds > keggResults<-hyperGTest(keggParams) > keggtab <-HTMLReport(shortName = paste("KEGG Analysis",cni_original,sep=''), + title='KEGG Analysis', reportDirectory="./reports") > publish(keggResults, keggtab, selectedIDs=DEGi_entrez, annotation.db="",categorySize=10, makePlot=FALSE) Error in, "data.frame") : cannot coerce class "structure("KEGGHyperGResult", package = "Category")" to a data.frame > sessionInfo() R version 3.0.1 (2013-05-16) Platform: x86_64-apple-darwin10.8.0 (64-bit) locale: [1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8 attached base packages: [1] splines grid parallel stats graphics grDevices utils datasets methods base other attached packages: [1] EMA_1.4.4 PFAM.db_2.9.0 proto_0.3-10 Category_2.26.0 hwriter_1.3 [6] lattice_0.20-24 affycoretools_1.32.1 KEGG.db_2.9.1 GO.db_2.9.0 ReportingTools_2.0.1 [11] statmod_1.4.18 MASS_7.3-29 hgu133plus2.db_2.9.0 RSQLite_0.11.4 [16] DBI_0.2-7 genefilter_1.42.0 car_2.0-19 Rcmdr_2.0-4 RcmdrPlugin.KMggplot2_0.2-0 [21] ggplot2_0.9.3.1 mgcv_1.7-28 nlme_3.1-113 rgl_0.93.996 hgu133plus2probe_2.12.0 [26] hgu133plus2cdf_2.12.0 AnnotationDbi_1.22.6 plyr_1.8 gcrma_2.32.0 gdata_2.13.2 [31] limma_3.16.8 affy_1.38.1 Biobase_2.20.1 BiocGenerics_0.6.0 loaded via a namespace (and not attached): [1] affyio_1.28.0 annaffy_1.32.0 annotate_1.38.0 AnnotationForge_1.2.2 BiocInstaller_1.10.4 biomaRt_2.16.0 [7] Biostrings_2.28.0 biovizBase_1.8.1 bit_1.1-11 bitops_1.0-6 BSgenome_1.28.0 caTools_1.16 [13] cluster_1.14.4 codetools_0.2-8 colorspace_1.2-4 dichromat_2.0-0 digest_0.6.4 edgeR_3.2.4 [19] FactoMineR_1.25 ff_2.2-12 foreach_1.4.1 Formula_1.1-1 GenomicFeatures_1.12.4 GenomicRanges_1.12.5 [25] ggbio_1.8.8 ggthemes_1.6.0 GOstats_2.26.0 gplots_2.12.1 graph_1.38.3 gridExtra_0.9.1 [31] GSA_1.03 GSEABase_1.22.0 gtable_0.1.2 gtools_3.3.0 heatmap.plus_1.3 Hmisc_3.14-0 [37] IRanges_1.18.4 iterators_1.0.6 KernSmooth_2.23-10 labeling_0.2 latticeExtra_0.6-26 Matrix_1.1-2 [43] multtest_2.16.0 munsell_0.4.2 nnet_7.3-7 oligoClasses_1.22.0 preprocessCore_1.22.0 R.methodsS3_1.6.1 [49] R.oo_1.17.0 R.utils_1.29.8 R2HTML_2.2.1 RBGL_1.36.2 RColorBrewer_1.0-5 RCurl_1.95-4.1 [55] reshape2_1.2.2 Rsamtools_1.12.4 rtracklayer_1.20.4 scales_0.2.3 siggenes_1.34.0 stats4_3.0.1 [61] stringr_0.6.2 survival_2.37-7 tcltk_3.0.1 tcltk2_1.2-10 tools_3.0.1 VariantAnnotation_1.6.8 [67] XML_3.95-0.2 xtable_1.7-1 zlibbioc_1.6.0 Mayte Suarez-Farinas, PhD Research Assistant Professor Laboratory of Investigative Dermatology Biostatistician, Center for Clinical and Translational Science The Rockefeller University 1230 York Ave, Box 178 New York, NY 10065 Phone: +1(212) 327-8213 Fax: +1(212) 327-8232
Pathways GO ReportingTools Pathways GO ReportingTools • 901 views

Login before adding your answer.

Traffic: 623 users visited in the last hour
Help About
Access RSS

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

Powered by the version 2.3.6