DiffBind—错误:没有站点的活动大于minMaxval
1
0
进入编辑模式
@095e334e
最后一次出现是3小时前

在这里输入正文

运行时dba.count ()在我的dba对象上:

samples <- dba(sampleSheet = "../ sampleSheet .csv",peakFormat = "bed",scoreCol = 5) ep_293T_siNeg_R1 ep_293T_siNeg_R1 ep_293t_sinegg 2床ep_293T_siFUS_R1 ep_293T_siFUS_R2 ep_293T_siFUS 1床ep_293T_siFUS_R2 ep_293T_siFUS 2床

我得到以下错误:

<- dba.count(samples)计算峰会…示例:切分0 /执行/ 10249 _sample1。bam125 Sample:分片-1/execution/10249_sample2。bam125 Sample:分片-2/execution/10249_sample3。bam125 Sample:分片-3/execution/10249_sample4。bam125 Sample:分片-4/execution/10249_sample5。bam125 Sample:分片-5/execution/10249_sample6。bam125 Sample:分片-6/execution/10249_sample7。bam125 Sample:分片-7/execution/10249_sample8。重新定心山峰…示例:切分0 /执行/ 10249 _sample1。bam125 read将被算作成对端。示例:切分1 /执行/ 10249 _sample2。bam125 read将被算作成对端。示例:切分2 /执行/ 10249 _sample3。bam125 read将被算作成对端。示例:shard-3 /执行/ 10249 _sample4。bam125 Reads will be counted as Paired-end. Sample: shard-4/execution/10249_sample5.bam125 Reads will be counted as Paired-end. Sample: shard-5/execution/10249_sample6.bam125 Reads will be counted as Paired-end. Sample: shard-6/execution/10249_sample7.bam125 Reads will be counted as Paired-end. Sample: shard-7/execution/10249_sample8.bam125 Reads will be counted as Paired-end. Error: No sites have activity greater than minMaxval In addition: Warning message: Parallel execution unavailable: executing serially.

我哪里搞错了?谢谢。我注意到有一个奇怪的数字连在我的本演讲文件即125

附:从github错误源于这一节:https://github.com/aeron15/DiffBind/blob/master/R/counts.R#L579但我不太熟悉这句话的意思。

R版本4.1.0(2021-05-18)平台:x86_64-w64-mingw32/x64(64位)在Windows 10 x64 (build 19042)下运行矩阵产品:default locale: [1] LC_COLLATE=English_United States.1252 LC_CTYPE=English_United States.1252 LC_MONETARY=English_United States.1252 LC_NUMERIC=C [5] LC_TIME=English_United States.1252附带的基本包:[1]parallel stats4 stats graphics grDevices utils datasets methods基于其他附加包:[1] DiffBind_3.2.1 SummarizedExperiment_1.22.0 Biobase_2.52.0 MatrixGenerics_1.4.0 matrixStats_0.58.0 [6] GenomicRanges_1.44.0 GenomeInfoDb_1.28.0 IRanges_2.26.0 S4Vectors_0.30.0 BiocGenerics_0.38.0通过命名空间加载(未附加):[1] backports_1.2.1 GOstats_2.58.0 BiocFileCache_2.0.0 plyr_1.8.6 gseabase_1.54 4.0 splines_4.1.0 [7] BiocParallel_1.26.0 ggplot2_3.3.3 amap_0.8-18 digest_0.6.27 invgamma_1. 1.1 GO.db_3.13.0 [13] SQUAREM_2021.1 fansi_0.4.2 magrittr_2.0.1 checkmate_2.0.0 memoise_2.0.0 BSgenome_1.60.0 [19] base64url_1. 1.4 limma_3.48.0 Biostrings_2.60.0 annotate_1.70.0 systemPipeR_1.26.2 bdsmatrix_1.3-4 [25] prettyunits_1.1.1 jpeg_0.1-8.1 colorspace_2.0-1 blob_1.2.1 rappdirs_0.3.3 apeglm_1.14.0 [31] ggrepel_0.9.1dplyr_1.0.6 crayon_1.4.1 RCurl_1.98-1.3 jsonlite_1.7.2 graph_1.70.0 [37] genefilter_1.74.0 brew_1. 38.0 glue_1.4.2 gtable_0.3.0 [43] zlibbioc_1.38.0 XVector_0.32.0 DelayedArray_0.18.0 v83.4.2 Rgraphviz_2.36.0 scales_1.1.1 [49] pheatmap_1.0.12 mvtnorm_1. 1.1-1 DBI_1.1.1 edgeR_3.34.0 Rcpp_1.0.6 xtable_1.8-4 [55] progress_1.2.2 emdbook_1.3.12 bit_4.0.4 rsvg_2.1.2 AnnotationForge_1.34.0 truncnorm_1.0-8 [61] httr_1.4.2 gplots_3.1.1 RColorBrewer_1.1-2ellipsis_0.3.2 pkgconfig_2.0.3 XML_3.99-0.6 [67] dbplyr_2.1.1 locfit_1.5-9.4 utf8_1.2.1 tidyselect_1.1.1 rlang_0.4.11 annotationdbi_1.54 4.0 [73] munsell_0.5.0 tools_4.1.0 cachem_1.0.5 generics_0.1.0 RSQLite_2.2.7 string_1 .4.0 [79] fastmap_1.1.0 yaml_2.2.1 bit64_4.0.5 caTools_1.18.2 purrr_0.3.4 KEGGREST_1.32.0 [85] RBGL_1.68.0 biomaRt_2.48.0 compiler_4.1.0 rstudioapi_0.13 filelock_1.0.2 curl_4.3.1 [91] png_0.1-7 tibble_3.1.2 stringi_1.6.1 GenomicFeatures_1.44.0 lattice_0.20-44 Matrix_1.3-3 [97]vctrs_0.3.8 pillar_1.6.1 lifecycle_1.0.0 data.table_1.14.0 bitops_1.0-7 irlba_2.3.3 [103] rtracklayer_1.2.0 R6_2.5.0 BiocIO_1.2.0 latticeExtra_0.6-29 hwriter_1.3.2 ShortRead_1.50.0 [109] KernSmooth_2.23-20 MASS_7.3-54 gtools_3.8.2 assertthat_0.2.1 Category_2.58.0 rjson_0.2.20 [115] withrs_4.2 GenomicAlignments_1.28.0 batchtools_0.9.15 Rsamtools_2.8.0 GenomeInfoDbData_1.2.6 hms_1.1.0 [121] grid_4.1.0 coda_0.19-4 GreyListChIP_1.24.0 ashr2.2 -47 mixsqp_0.3-43
DiffBind•19个视图
添加评论
0
进入编辑模式
@rory -斯塔克- 5741
最后一次出现是15分钟前
克鲁克,剑桥,英国

错误消息需要修复,因为它引用了一个内部值(minMaxval)而不是参数过滤器

这表示所有计算的RPKM值都小于的默认筛选器值1.0

你可能想在阅读实验之后和计数之前检查一下矩阵是什么样子的。如下图所示:

dba.peakset(样本)

这种方法可以检查您是否有一组间隔,您希望在其中找到唯一重叠的读取。

添加评论

登录在回答之前。

流量:最近一小时访问了502个用户

使用本网站即表示接受我们的用户协议和隐私政策

2.3.6版本