Question: ReactomePA : error with gsePathway
gravatar for samuel collombet
3.2 years ago by
samuel collombet120 wrote:


I would like to performed a continuous gsea with the ReactomePA, but I get the following result:

> y <- gsePathway(geneList_entrez, organism="mouse", nPerm = 100, minGSSize = 20,pvalueCutoff = 1, pAdjustMethod = "BH", exponent=1, verbose = T)
[1] "calculating observed enrichment scores..."
[1] "calculating permutation scores..."
  |============================================== |  99%
[1] "calculating p values..."
Error in .checkKeys(value, Lkeys(x), x@ifnotfound) : 
  value for "5604596" not found

my geneList seems well formated:

> head(geneList_entrez)
    18223     52428     67080     13056    225908    380684 
0.9999980 0.9999966 0.9999965 0.9999948 0.9999936 0.9999923 

I can reproduce the example in the vignette...
I thought it might be because my score are <1, but it doesn't work if I change them into the rank of the value.

> sessionInfo()
R version 3.1.0 RC (2014-04-05 r65382)
Platform: x86_64-pc-linux-gnu (64-bit)

 [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C               LC_TIME=en_US.UTF-8        LC_COLLATE=en_US.UTF-8     LC_MONETARY=en_US.UTF-8   
 [6] LC_MESSAGES=en_US.UTF-8    LC_PAPER=en_US.UTF-8       LC_NAME=C                  LC_ADDRESS=C               LC_TELEPHONE=C            

attached base packages:
 [1] tcltk     grid      parallel  stats4    stats     graphics  grDevices utils     datasets  methods   base     

other attached packages:
 [1] reactome.db_1.50.0        topGO_2.18.0              SparseM_1.6               GO.db_3.0.0              
 [6] graph_1.44.1              biomaRt_2.22.0            Mfuzz_2.26.0              DynDoc_1.44.0             widgetTools_1.44.0       
[11] e1071_1.6-4               wq_0.4-1                  zoo_1.7-12                reshape2_1.4.1            ggplot2_1.0.1            
[16] RColorBrewer_1.1-2        DESeq2_1.6.3              RcppArmadillo_0.4.650.1.1 Rcpp_0.11.5               GenomicRanges_1.18.4     
[21]        ReactomePA_1.10.1         DOSE_2.4.0                AnnotationDbi_1.28.2      GenomeInfoDb_1.2.4       
[26] IRanges_2.0.1             S4Vectors_0.4.0           Biobase_2.26.0            BiocGenerics_0.12.1       RSQLite_1.0.0            
[31] DBI_0.3.1                

loaded via a namespace (and not attached):
 [1] acepack_1.3-3.3     annotate_1.44.0     base64enc_0.1-2     BatchJobs_1.6       BBmisc_1.9          BiocParallel_1.0.3  bitops_1.0-6       
 [8] brew_1.0-6          checkmate_1.5.2     class_7.3-12        cluster_2.0.1       codetools_0.2-11    colorspace_1.2-6    digest_0.6.8       
[15] DO.db_2.8.0         fail_1.2            foreach_1.4.2       foreign_0.8-63      Formula_1.2-0       genefilter_1.48.1   geneplotter_1.44.0 
[22] GOSemSim_1.24.1     graphite_1.12.0     gtable_0.1.2        Hmisc_3.15-0        igraph_0.7.1        iterators_1.0.7     lattice_0.20-30    
[29] latticeExtra_0.6-26 locfit_1.5-9.1      MASS_7.3-40         munsell_0.4.2       nnet_7.3-9          plyr_1.8.1          proto_0.3-10       
[36] qvalue_1.43.0       RCurl_1.95-4.5      rpart_4.1-9         scales_0.2.4        sendmailR_1.2-1     splines_3.1.0       stringr_0.6.2      
[43] survival_2.38-1     tkWidgets_1.44.0    tools_3.1.0         XML_3.98-1.1        xtable_1.7-4        XVector_0.6.0     
ADD COMMENTlink modified 3.2 years ago by Guangchuang Yu990 • written 3.2 years ago by samuel collombet120
gravatar for Guangchuang Yu
3.2 years ago by
Hong Kong
Guangchuang Yu990 wrote:

Hi Samuel,

This is an known issue, please refer to reactome.db is not updated? .

This issue is not caused by ReactomePA, but due to reactome.db. I don't know why reactome are not responsive this time. The current release of reactome.db is not correct and hopefully they will offer right data in next release.

Currently, you can manual install the old version of reactome.db,, and everything will be OK.




ADD COMMENTlink written 3.2 years ago by Guangchuang Yu990

Now it works fine :) Thanks!

ADD REPLYlink modified 3.2 years ago • written 3.2 years ago by samuel collombet120
gravatar for Guangchuang Yu
3.2 years ago by
Hong Kong
Guangchuang Yu990 wrote:

BTW: you may interesting in GSEA analysis with latest online KEGG data using clusterProfiler, .

ADD COMMENTlink written 3.2 years ago by Guangchuang Yu990

Great Thanks! I haven't tried clusterProfiler yet, will do :)

ADD REPLYlink written 3.2 years ago by samuel collombet120
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: 133 users visited in the last hour