Issue with wavClusteR
0
0
Entering edit mode
btharada • 0
@btharada-22272
Last seen 4.5 years ago

I am having an issue with the getAllSub function in the wavClusteR package. I have read my bam file using the readSortedBam function,

library(wavClusteR)
bam <- readSortedBam("~/parclip/bam/replicate_1.bam")
countTable <- getAllSub(bam, minCov=10)

but when I attempt to run the getAllSub function on the data, I get the following error:

Loading required package: doMC
Loading required package: foreach
Loading required package: iterators
Considering substitutions, n = 28204, processing in 29 chunks
   chunk #: 1
   chunk #: 2
   chunk #: 3
   chunk #: 4
   chunk #: 5
   chunk #: 6
   chunk #: 7
   chunk #: 8
   chunk #: 9
   chunk #: 10
   chunk #: 11
   chunk #: 12
   chunk #: 13
   chunk #: 14
   chunk #: 15
   chunk #: 16
   chunk #: 17
   chunk #: 18
   chunk #: 19
   chunk #: 20
   chunk #: 21
   chunk #: 22
   chunk #: 23
   chunk #: 24
   chunk #: 25
   chunk #: 26
   chunk #: 27
   chunk #: 28
   chunk #: 29
Error in { : 
  task 1 failed - "solving row 7017: range cannot be determined from the supplied arguments (too many NAs)"

My bam file has been sorted and indexed via samtools as directed in the instructions. The error seems to occur during the "at <- IRanges(unlist(offset), width = 1)" call within the processChunk function within the processMD function called by getAllSub:

>   at <- IRanges(unlist(offset), width = 1)
Error in .Call2("solve_user_SEW0", start, end, width, PACKAGE = "IRanges") : 
  solving row 7017: range cannot be determined from the supplied arguments (too many NAs)

Please advise as to how to get around this error message.

Here is the my session info for reference

R version 3.6.0 (2019-04-26)
Platform: x86_64-pc-linux-gnu (64-bit)
Running under: Ubuntu 18.04.2 LTS

Matrix products: default
BLAS:   /usr/lib/x86_64-linux-gnu/openblas/libblas.so.3
LAPACK: /usr/lib/x86_64-linux-gnu/libopenblasp-r0.2.20.so

locale:
 [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C               LC_TIME=en_US.UTF-8       
 [4] LC_COLLATE=en_US.UTF-8     LC_MONETARY=en_US.UTF-8    LC_MESSAGES=en_US.UTF-8   
 [7] LC_PAPER=en_US.UTF-8       LC_NAME=C                  LC_ADDRESS=C              
[10] LC_TELEPHONE=C             LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C       

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

other attached packages:
 [1] doMC_1.3.6           iterators_1.0.12     foreach_1.4.7        wavClusteR_2.18.0   
 [5] Rsamtools_2.0.0      Biostrings_2.52.0    XVector_0.24.0       GenomicRanges_1.36.0
 [9] GenomeInfoDb_1.20.0  IRanges_2.18.2       S4Vectors_0.22.0     BiocGenerics_0.30.0 

loaded via a namespace (and not attached):
 [1] Biobase_2.44.0              httr_1.4.1                  bit64_0.9-7                
 [4] splines_3.6.0               Formula_1.2-3               assertthat_0.2.1           
 [7] latticeExtra_0.6-28         blob_1.2.0                  GenomeInfoDbData_1.2.1     
[10] progress_1.2.2              pillar_1.4.2                RSQLite_2.1.2              
[13] backports_1.1.4             lattice_0.20-38             glue_1.3.1                 
[16] digest_0.6.20               checkmate_1.9.4             RColorBrewer_1.1-2         
[19] colorspace_1.4-1            htmltools_0.4.0             Matrix_1.2-17              
[22] XML_3.98-1.20               pkgconfig_2.0.2             biomaRt_2.40.4             
[25] zlibbioc_1.30.0             purrr_0.3.2                 scales_1.0.0               
[28] ifultools_2.0-5             BiocParallel_1.18.1         htmlTable_1.13.1           
[31] tibble_2.1.3                ggplot2_3.2.1               SummarizedExperiment_1.14.1
[34] GenomicFeatures_1.36.4      nnet_7.3-12                 lazyeval_0.2.2             
[37] survival_2.44-1.1           magrittr_1.5                crayon_1.3.4               
[40] mclust_5.4.5                memoise_1.1.0               MASS_7.3-51.4              
[43] foreign_0.8-72              data.table_1.12.2           tools_3.6.0                
[46] prettyunits_1.0.2           hms_0.5.1                   matrixStats_0.54.0         
[49] stringr_1.4.0               munsell_0.5.0               cluster_2.1.0              
[52] DelayedArray_0.10.0         AnnotationDbi_1.46.1        wmtsa_2.0-3                
[55] ade4_1.7-13                 compiler_3.6.0              rlang_0.4.0                
[58] grid_3.6.0                  RCurl_1.95-4.12             rstudioapi_0.10            
[61] htmlwidgets_1.5.1           bitops_1.0-6                base64enc_0.1-3            
[64] gtable_0.3.0                codetools_0.2-16            DBI_1.0.0                  
[67] R6_2.4.0                    gridExtra_2.3               GenomicAlignments_1.20.1   
[70] knitr_1.24                  dplyr_0.8.3                 rtracklayer_1.44.3         
[73] seqinr_3.6-1                bit_1.1-14                  zeallot_0.1.0              
[76] Hmisc_4.2-0                 stringi_1.4.3               Rcpp_1.0.2                 
[79] rpart_4.1-15                vctrs_0.2.0                 acepack_1.4.1              
[82] splus2R_1.2-2               xfun_0.9                    tidyselect_0.2.5
wavClusteR • 470 views
ADD COMMENT

Login before adding your answer.

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