dmrseq icon indicating copy to clipboard operation
dmrseq copied to clipboard

Rbind error during permutations

Open plbngl opened this issue 8 months ago • 1 comments

Hello! Thanks for the tool! I am encountering an error that seems to be recurrent but should have been fixed according to the forums.

I am using dmrseq_1.26.0 , installed from BiocManager and having the following error (testCovariate is multi-factor):

blocks <- dmrseq(bs=bsseq_obj,
                 cutoff = 0.05,
                 testCovariate="Genotype",
                 adjustCovariate = 'Replicate',
                 block = TRUE,
                 minInSpan = 500,
                 bpSpan = 5e4,
                 maxGapSmooth = 1e6,
                 maxGap = 5e3)
 Beginning permutation 5
...Chromosome chr1: Smoothed (0.04 min). 430 regions scored (0.47 min). 
...Chromosome chr10: Smoothed (0.05 min). 272 regions scored (0.25 min). 
...Chromosome chr11: Smoothed (0.02 min). Error in rbind(deparse.level, ...) : 
  numbers of columns of arguments do not match

R version 4.4.0 (2024-04-24)
Platform: x86_64-pc-linux-gnu
Running under: Red Hat Enterprise Linux 8.8 (Ootpa)

Matrix products: default
BLAS/LAPACK: /usr/lib64/libopenblasp-r0.3.12.so;  LAPACK version 3.9.0

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

time zone: Europe/Rome
tzcode source: system (glibc)

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

other attached packages:
 [1] BiocParallel_1.40.2         data.table_1.17.0           RColorBrewer_1.1-3          dmrseq_1.26.0              
 [5] bsseq_1.42.0                SummarizedExperiment_1.36.0 Biobase_2.66.0              MatrixGenerics_1.18.1      
 [9] matrixStats_1.5.0           GenomicRanges_1.58.0        GenomeInfoDb_1.42.3         IRanges_2.40.1             
[13] S4Vectors_0.44.0            BiocGenerics_0.52.0        

loaded via a namespace (and not attached):
  [1] DBI_1.2.3                 bitops_1.0-9              permute_0.9-7             rlang_1.1.6              
  [5] magrittr_2.0.3            compiler_4.4.0            RSQLite_2.3.9             GenomicFeatures_1.58.0   
  [9] DelayedMatrixStats_1.28.1 png_0.1-8                 vctrs_0.6.5               reshape2_1.4.4           
 [13] stringr_1.5.1             pkgconfig_2.0.3           crayon_1.5.3              fastmap_1.2.0            
 [17] dbplyr_2.5.0              XVector_0.46.0            Rsamtools_2.22.0          tzdb_0.5.0               
 [21] UCSC.utils_1.2.0          bit_4.6.0                 zlibbioc_1.52.0           cachem_1.1.0             
 [25] jsonlite_2.0.0            blob_1.2.4                rhdf5filters_1.18.1       DelayedArray_0.32.0      
 [29] Rhdf5lib_1.28.0           parallel_4.4.0            R6_2.6.1                  stringi_1.8.7            
 [33] limma_3.62.2              rtracklayer_1.66.0        iterators_1.0.14          Rcpp_1.0.14              
 [37] R.utils_2.13.0            readr_2.1.5               splines_4.4.0             Matrix_1.7-0             
 [41] tidyselect_1.2.1          abind_1.4-8               yaml_2.3.10               codetools_0.2-20         
 [45] curl_6.2.2                doRNG_1.8.6.2             lattice_0.22-6            tibble_3.2.1             
 [49] regioneR_1.38.0           plyr_1.8.9                KEGGREST_1.46.0           BiocFileCache_2.14.0     
 [53] Biostrings_2.74.1         pillar_1.10.2             BiocManager_1.30.25       filelock_1.0.3           
 [57] rngtools_1.5.2            foreach_1.5.2             generics_0.1.3            RCurl_1.98-1.17          
 [61] BiocVersion_3.20.0        hms_1.1.3                 ggplot2_3.5.2             sparseMatrixStats_1.18.0 
 [65] munsell_0.5.1             scales_1.3.0              bumphunter_1.48.0         gtools_3.9.5             
 [69] glue_1.8.0                tools_4.4.0               AnnotationHub_3.14.0      BiocIO_1.16.0            
 [73] BSgenome_1.74.0           locfit_1.5-9.12           GenomicAlignments_1.42.0  XML_3.99-0.18            
 [77] rhdf5_2.50.2              grid_4.4.0                AnnotationDbi_1.68.0      colorspace_2.1-1         
 [81] nlme_3.1-164              GenomeInfoDbData_1.2.13   HDF5Array_1.34.0          restfulr_0.0.15          
 [85] annotatr_1.32.0           cli_3.6.4                 rappdirs_0.3.3            S4Arrays_1.6.0           
 [89] dplyr_1.1.4               gtable_0.3.6              outliers_0.15             R.methodsS3_1.8.2        
 [93] digest_0.6.37             SparseArray_1.6.2         rjson_0.2.23              memoise_2.0.1            
 [97] R.oo_1.27.0               lifecycle_1.0.4           httr_1.4.7                statmod_1.5.0            
[101] bit64_4.6.0-1            

This error did not occur when using just a subset of the regions to test. Can you please advise on how to fix it? Thank you!!

plbngl avatar Apr 22 '25 12:04 plbngl