Entering edit mode
Ha
•
0
@Ha-24678
Last seen 3.8 years ago
Hello,
I'm trying to run DiffBind with my ChIP-seq data. I do not want to subtract control reads from ChIP reads, which should be controlled by bSubControl
in dba.count()
. However, the parameter does not seem to be available to me. I also did not run grey list analysis.
> ?dba.count
Usage:
dba.count(DBA, peaks, minOverlap=2, score=DBA_SCORE_NORMALIZED,
fragmentSize=DBA$config$fragmentSize,
summits=200, filter=5, bRemoveDuplicates=FALSE, bScaleControl=TRUE,
mapQCth=DBA$config$mapQCth, filterFun=max, minCount=0,
bLog=FALSE, bUseSummarizeOverlaps=TRUE,
readFormat=DBA_READS_DEFAULT, bParallel=DBA$config$RunParallel)
> sessionInfo()
R version 4.0.2 (2020-06-22)
Platform: x86_64-conda_cos6-linux-gnu (64-bit)
Running under: Ubuntu 16.04 LTS
Matrix products: default
BLAS/LAPACK: /opt/miniconda3/envs/diffbind/lib/libopenblasp-r0.3.10.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] parallel stats4 stats graphics grDevices utils datasets
[8] methods base
other attached packages:
[1] forcats_0.5.0 stringr_1.4.0
[3] purrr_0.3.4 readr_1.4.0
[5] tidyr_1.1.2 tibble_3.0.4
[7] ggplot2_3.3.3 tidyverse_1.3.0
[9] dplyr_1.0.2 DiffBind_3.0.3
[11] SummarizedExperiment_1.20.0 Biobase_2.50.0
[13] MatrixGenerics_1.2.0 matrixStats_0.57.0
[15] GenomicRanges_1.42.0 GenomeInfoDb_1.26.0
[17] IRanges_2.24.0 S4Vectors_0.28.0
[19] BiocGenerics_0.36.0
loaded via a namespace (and not attached):
[1] readxl_1.3.1 backports_1.2.1 GOstats_2.56.0
[4] BiocFileCache_1.14.0 plyr_1.8.6 GSEABase_1.52.0
[7] splines_4.0.2 BiocParallel_1.24.0 amap_0.8-18
[10] digest_0.6.27 invgamma_1.1 GO.db_3.12.1
[13] fansi_0.4.1 SQUAREM_2021.1 magrittr_2.0.1
[16] checkmate_2.0.0 memoise_1.1.0 BSgenome_1.58.0
[19] base64url_1.4 limma_3.46.0 Biostrings_2.58.0
[22] annotate_1.68.0 modelr_0.1.8 systemPipeR_1.24.2
[25] askpass_1.1 bdsmatrix_1.3-4 prettyunits_1.1.1
[28] jpeg_0.1-8.1 colorspace_2.0-0 rvest_0.3.6
[31] blob_1.2.1 rappdirs_0.3.1 apeglm_1.12.0
[34] ggrepel_0.9.0 haven_2.3.1 crayon_1.3.4
[37] RCurl_1.98-1.2 jsonlite_1.7.2 graph_1.68.0
[40] genefilter_1.72.0 brew_1.0-6 survival_3.2-7
[43] VariantAnnotation_1.36.0 glue_1.4.2 gtable_0.3.0
[46] zlibbioc_1.36.0 XVector_0.30.0 DelayedArray_0.16.0
[49] V8_3.4.0 Rgraphviz_2.34.0 scales_1.1.1
[52] pheatmap_1.0.12 mvtnorm_1.1-1 DBI_1.1.0
[55] edgeR_3.32.0 Rcpp_1.0.5 xtable_1.8-4
[58] progress_1.2.2 emdbook_1.3.12 bit_4.0.4
[61] rsvg_2.1 AnnotationForge_1.32.0 truncnorm_1.0-8
[64] httr_1.4.2 gplots_3.1.1 RColorBrewer_1.1-2
[67] ellipsis_0.3.1 pkgconfig_2.0.3 XML_3.99-0.5
[70] dbplyr_2.0.0 locfit_1.5-9.4 tidyselect_1.1.0
[73] rlang_0.4.10 AnnotationDbi_1.52.0 cellranger_1.1.0
[76] munsell_0.5.0 tools_4.0.2 cli_2.2.0
[79] generics_0.1.0 RSQLite_2.2.2 broom_0.7.3
[82] yaml_2.2.1 fs_1.5.0 bit64_4.0.5
[85] caTools_1.18.1 RBGL_1.66.0 xml2_1.3.2
[88] biomaRt_2.46.0 rstudioapi_0.13 compiler_4.0.2
[91] curl_4.3 png_0.1-7 reprex_0.3.0
[94] stringi_1.5.3 ps_1.5.0 GenomicFeatures_1.42.0
[97] lattice_0.20-41 Matrix_1.3-2 vctrs_0.3.6
[100] pillar_1.4.7 lifecycle_0.2.0 data.table_1.13.6
[103] bitops_1.0-6 irlba_2.3.3 rtracklayer_1.50.0
[106] R6_2.5.0 latticeExtra_0.6-29 hwriter_1.3.2
[109] ShortRead_1.48.0 KernSmooth_2.23-18 MASS_7.3-53
[112] gtools_3.8.2 assertthat_0.2.1 openssl_1.4.3
[115] Category_2.56.0 rjson_0.2.20 withr_2.3.0
[118] GenomicAlignments_1.26.0 batchtools_0.9.15 Rsamtools_2.6.0
[121] GenomeInfoDbData_1.2.4 hms_0.5.3 grid_4.0.2
[124] DOT_0.1 coda_0.19-4 GreyListChIP_1.22.0
[127] ashr_2.2-47 mixsqp_0.3-43 bbmle_1.0.23.1
[130] lubridate_1.7.9.2 numDeriv_2016.8-1.1
I really appreciate any help. Thank you!