ZarnackGroup / BindingSiteFinder

Package for the definition of biniding sites for iCLIP data
https://www.bioconductor.org/packages/release/bioc/html/BindingSiteFinder.html
6 stars 1 forks source link

! Can't recycle `quantiles` (size 101) to size 0. #11

Closed dongyangzheng closed 2 months ago

dongyangzheng commented 2 months ago

Hi profess when i repeat the workflow from this repo https://github.com/malewins/Plant-iCLIPseq , come across some problem Estimation at: 0% Error in dplyr_col_modify(): ! Can't recycle quantiles (size 101) to size 0. Backtrace: ▆

  1. ├─BindingSiteFinder::estimateBsWidth(...)
  2. │ └─base::lapply(...)
  3. │ └─BindingSiteFinder (local) FUN(X[[i]], ...)
  4. │ └─BindingSiteFinder::pureClipGeneWiseFilter(...)
  5. │ └─... %>% ...
  6. ├─dplyr::filter(., as.character(cutoffs) == as.character(cutoff))
  7. ├─dplyr::mutate(., cutoffs = potentialCutoffs)
  8. ├─dplyr::group_by(., queryHits)
  9. ├─dplyr::reframe(., quantiles = quantile(score, probs = potentialCutoffs))
    1. ├─dplyr:::reframe.data.frame(., quantiles = quantile(score, probs = potentialCutoffs))
    2. │ └─dplyr:::summarise_build(by, cols)
    3. │ ├─dplyr::dplyr_col_modify(out, cols$new)
    4. │ └─dplyr:::dplyr_col_modify.data.frame(out, cols$new)
    5. │ └─vctrs::vec_recycle_common(!!!cols, .size = nrow(data))
    6. └─vctrs:::stop_recycle_incompatible_size(...)
    7. └─vctrs:::stop_vctrs(...)
    8. └─rlang::abort(message, class = c(class, "vctrs_error"), ..., call = call) Warning messages: 1: In .merge_two_Seqinfo_objects(x, y) : The 2 combined objects have no sequence levels in common. (Use suppressWarnings() to suppress this warning.) 2: In .merge_two_Seqinfo_objects(x, y) : The 2 combined objects have no sequence levels in common. (Use suppressWarnings() to suppress this warning.) Execution halted R version 4.2.2 (2022-10-31) Platform: x86_64-conda-linux-gnu (64-bit) Running under: CentOS Linux 7 (Core)

Matrix products: default BLAS/LAPACK: /home/zhengdy/miniconda3/envs/diffscan/lib/libopenblasp-r0.3.21.so

locale: [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C
[3] 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
[7] LC_PAPER=en_US.UTF-8 LC_NAME=C
[9] LC_ADDRESS=C LC_TELEPHONE=C
[11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C

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

other attached packages: [1] BindingSiteFinder_2.0.0 GenomicRanges_1.50.2 GenomeInfoDb_1.34.9
[4] IRanges_2.32.0 S4Vectors_0.36.2 BiocGenerics_0.44.0

loaded via a namespace (and not attached): [1] bitops_1.0-7 matrixStats_1.3.0
[3] bit64_4.0.5 filelock_1.0.3
[5] doParallel_1.0.17 RColorBrewer_1.1-3
[7] progress_1.2.3 httr_1.4.7
[9] tools_4.2.2 utf8_1.2.4
[11] R6_2.5.1 DBI_1.2.2
[13] colorspace_2.1-0 ggdist_3.3.2
[15] GetoptLong_1.0.5 withr_3.0.0
[17] tidyselect_1.2.1 prettyunits_1.2.0
[19] bit_4.0.5 curl_5.2.1
[21] compiler_4.2.2 cli_3.6.2
[23] Biobase_2.58.0 xml2_1.3.6
[25] DelayedArray_0.24.0 rtracklayer_1.58.0
[27] scales_1.3.0 rappdirs_0.3.3
[29] systemfonts_1.0.6 stringr_1.5.1
[31] digest_0.6.35 Rsamtools_2.14.0
[33] svglite_2.1.3 rmarkdown_2.26
[35] XVector_0.38.0 htmltools_0.5.8.1
[37] pkgconfig_2.0.3 MatrixGenerics_1.10.0
[39] dbplyr_2.5.0 fastmap_1.1.1
[41] rlang_1.1.3 GlobalOptions_0.1.2
[43] rstudioapi_0.16.0 RSQLite_2.3.6
[45] farver_2.1.1 shape_1.4.6.1
[47] BiocIO_1.8.0 generics_0.1.3
[49] BiocParallel_1.32.6 dplyr_1.1.4
[51] distributional_0.4.0 RCurl_1.98-1.14
[53] magrittr_2.0.3 kableExtra_1.4.0
[55] GenomeInfoDbData_1.2.9 Matrix_1.5-3
[57] Rcpp_1.0.12 munsell_0.5.1
[59] fansi_1.0.6 lifecycle_1.0.4
[61] stringi_1.8.3 yaml_2.3.8
[63] MASS_7.3-58.2 SummarizedExperiment_1.28.0 [65] zlibbioc_1.44.0 plyr_1.8.9
[67] BiocFileCache_2.6.1 grid_4.2.2
[69] blob_1.2.4 parallel_4.2.2
[71] crayon_1.5.2 lattice_0.22-6
[73] Biostrings_2.66.0 GenomicFeatures_1.50.4
[75] circlize_0.4.16 hms_1.1.3
[77] KEGGREST_1.38.0 knitr_1.46
[79] ComplexHeatmap_2.15.4 pillar_1.9.0
[81] rjson_0.2.21 codetools_0.2-20
[83] biomaRt_2.54.1 XML_3.99-0.13
[85] glue_1.7.0 evaluate_0.23
[87] tweenr_2.0.3 png_0.1-8
[89] vctrs_0.6.5 foreach_1.5.2
[91] polyclip_1.10-6 gtable_0.3.5
[93] purrr_1.0.2 tidyr_1.3.1
[95] clue_0.3-65 cachem_1.0.8
[97] ggplot2_3.5.1 xfun_0.43
[99] ggforce_0.4.2 restfulr_0.0.15
[101] viridisLite_0.4.2 tibble_3.2.1
[103] iterators_1.0.14 GenomicAlignments_1.34.1
[105] AnnotationDbi_1.60.2 memoise_2.0.1
[107] cluster_2.1.6

dongyangzheng commented 2 months ago

i have solved this thank u

MelinaKlostermann commented 2 months ago

Ok perfect!