Question: chipenrich mm10 error
gravatar for jamie.gearing
8 months ago by
jamie.gearing0 wrote:

I have encountered an error while trying to use the chipenrich package with the mm10 genome. Here is a small example:

peaks <- data(peaks_E2F4, package = '')
results <- chipenrich(head(peaks), out_name=NULL, genome = "mm10", genesets = c("GOBP", "GOCC", "GOMF"))

# Error in chipenrich(head(peaks), out_name = NULL, genome = "mm10") : 
#   Invalid organism / geneset(s) combination requested: mmu TRUE
# In addition: Warning message:
#   In if (file.exists(genesets)) { :
#       the condition has length > 1 and only the first element will be used​

I think the problem might lie on line 33 of the chipenrich() function. For "mmu" and these gene sets:

any(supported_genesets()$organism == "mmu" & c("GOBP", "GOCC", "GOMF") == supported_genesets()$geneset)
# [1] FALSE

Perhaps this should be:

any(supported_genesets()$organism == "mmu" & supported_genesets()$geneset %in% c("GOBP", "GOCC", "GOMF"))
# [1] TRUE




# R version 3.3.2 (2016-10-31)
# Platform: x86_64-apple-darwin13.4.0 (64-bit)
# Running under: macOS Sierra 10.12.3
# locale:
#   [1] en_AU.UTF-8/en_AU.UTF-8/en_AU.UTF-8/C/en_AU.UTF-8/en_AU.UTF-8
# attached base packages:
#   [1] stats     graphics  grDevices utils     datasets  methods   base     
# other attached packages:
#   [1] chipenrich_1.12.0
# loaded via a namespace (and not attached):
#   [1] httr_1.2.1             tidyr_0.6.1            jsonlite_1.3           viridisLite_0.1.3      splines_3.3.2          Formula_1.2-1         
# [7] shiny_1.0.0            assertthat_0.1         stats4_3.3.2           latticeExtra_0.6-28    backports_1.0.5        lattice_0.20-34       
# [13] quantreg_5.29          limma_3.30.12          digest_0.6.12          GenomicRanges_1.26.3   RColorBrewer_1.1-2     XVector_0.14.0        
# [19] checkmate_1.8.2        colorspace_1.3-2       sandwich_2.3-4         rms_5.1-0              htmltools_0.3.5        httpuv_1.3.3          
# [25] Matrix_1.2-8           plyr_1.8.4             SparseM_1.74           chipenrich.data_1.10.0 zlibbioc_1.20.0        purrr_0.2.2           
# [31] xtable_1.8-2           mvtnorm_1.0-6          scales_0.4.1           MatrixModels_0.4-1     htmlTable_1.9          tibble_1.2            
# [37] mgcv_1.8-17            IRanges_2.8.1          ggplot2_2.2.1          TH.data_1.0-8          nnet_7.3-12            BiocGenerics_0.20.0   
# [43] lazyeval_0.2.0         survival_2.40-1        magrittr_1.5           mime_0.5               polspline_1.1.12       MASS_7.3-45           
# [49] nlme_3.1-131           foreign_0.8-67         colourpicker_0.3       tools_3.3.2            data.table_1.10.4      multcomp_1.4-6        
# [55] stringr_1.2.0          plotly_4.5.6.9000      S4Vectors_0.12.1       munsell_0.4.3          cluster_2.0.5          GenomeInfoDb_1.10.3   
# [61] grid_3.3.2             RCurl_1.95-4.8         htmlwidgets_0.8        miniUI_0.1.1           bitops_1.0-6           base64enc_0.1-3       
# [67] codetools_0.2-15       gtable_0.2.0           DBI_0.5-1              R6_2.2.0               gridExtra_2.2.1        zoo_1.7-14            
# [73] knitr_1.15.1           dplyr_0.5.0            Hmisc_4.0-2            stringi_1.1.2          parallel_3.3.2         Rcpp_0.12.9           
# [79] rpart_4.1-10           acepack_1.4.1         


ADD COMMENTlink modified 8 months ago by rcavalca40 • written 8 months ago by jamie.gearing0
gravatar for rcavalca
8 months ago by
United States
rcavalca40 wrote:

Yikes! Thank you for pointing this out!

I just pushed a patch (v1.12.1) to the release branch of Bioconductor. You can wait for it to work through the build system, or you can install the version with devtools::install_github('Bioconductor-mirror/chipenrich@release-3.4').

Thanks for using chipenrich.

ADD COMMENTlink written 8 months ago by rcavalca40

Thanks very much for such a fast response! This package is just what I need, so I am very grateful.

ADD REPLYlink written 8 months ago by jamie.gearing0
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.2.0
Traffic: 249 users visited in the last hour