DESeq2_Error in .call2, what does it mean?
2
0
Entering edit mode
tg369 ▴ 10
@tg369-9320
Last seen 4.7 years ago
United Kingdom

I was running DESeq2. However after running DESEq() I got the following:

> dds<-DESeq(dds);
using pre-existing size factors
estimating dispersions
gene-wise dispersion estimates
Error in .Call2("solve_user_SEW0", start, end, width, PACKAGE = "IRanges") : 
  solving row 463827: negative widths are not allowed

The sessionInfo is given below:

> sessionInfo()
R version 3.2.3 (2015-12-10)
Platform: x86_64-apple-darwin13.4.0 (64-bit)
Running under: OS X 10.9.2 (Mavericks)

locale:
[1] en_GB.UTF-8/en_GB.UTF-8/en_GB.UTF-8/C/en_GB.UTF-8/en_GB.UTF-8

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

other attached packages:
 [1] DESeq2_1.10.0              RcppArmadillo_0.6.300.2.0  Rcpp_0.12.2                SummarizedExperiment_1.0.1
 [5] Biobase_2.30.0             GenomicRanges_1.22.2       GenomeInfoDb_1.6.1         IRanges_2.4.6             
 [9] S4Vectors_0.8.5            BiocGenerics_0.16.1        BiocInstaller_1.20.1      

loaded via a namespace (and not attached):
 [1] RColorBrewer_1.1-2   futile.logger_1.4.1  plyr_1.8.3           XVector_0.10.0       futile.options_1.0.0 tools_3.2.3         
 [7] zlibbioc_1.16.0      rpart_4.1-10         digest_0.6.8         RSQLite_1.0.0        annotate_1.48.0      gtable_0.1.2        
[13] lattice_0.20-33      DBI_0.3.1            proto_0.3-10         gridExtra_2.0.0      genefilter_1.52.0    cluster_2.0.3       
[19] stringr_1.0.0        locfit_1.5-9.1       nnet_7.3-11          grid_3.2.3           AnnotationDbi_1.32.2 XML_3.98-1.3        
[25] survival_2.38-3      BiocParallel_1.4.0   foreign_0.8-66       latticeExtra_0.6-26  Formula_1.2-1        geneplotter_1.48.0  
[31] ggplot2_1.0.1        reshape2_1.4.1       lambda.r_1.1.7       magrittr_1.5         scales_0.3.0         Hmisc_3.17-0        
[37] MASS_7.3-45          splines_3.2.3        xtable_1.8-0         colorspace_1.2-6     stringi_1.0-1        acepack_1.3-3.3     
[43] munsell_0.4.2       

 

Please advice me.  Thank you in advance.

 

deseq2 • 1.0k views
ADD COMMENT
0
Entering edit mode
@mikelove
Last seen 23 hours ago
United States
Was dds built in a previous version of Bioconductor? The base class, SummarizedExperiment underwent some refactoring, you may need to call updateObject on dds first.
ADD COMMENT
0
Entering edit mode
tg369 ▴ 10
@tg369-9320
Last seen 4.7 years ago
United Kingdom

Thank you Mike. That was the case. doing updateObject (dds) fixed the problem.

ADD COMMENT

Login before adding your answer.

Traffic: 988 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