Problem reading Bed file using dba.count in DiffBind
2
1
Entering edit mode
urjaswita ▴ 40
@urjaswita-13128
Last seen 19 months ago

Hi folks,

I am having a trouble with dba.count function in diffbind. It fails with warnings but without any specific information. Please see below:

> obj <- dba(sampleSheet='analysis_H3K4me3_samples.csv')
S1 Eye  A None 1 Bed
S2 Eye  A None 2 Bed
S3 Eye  B None 1 Bed
S4 Eye  B None 2 Bed
S5 Eye  C None 1 Bed
S6 Eye  C None 2 Bed


> obj
6 Samples, 19333 sites in matrix (23582 total):
  ID  Tissue Condition Treatment Replicate Caller Intervals
1 S1     Eye         A      None         1    Bed     17211
2 S2     Eye         A      None         2    Bed     16219
3 S3     Eye         B      None         1    Bed     18846
4 S4     Eye         B      None         2    Bed     19704
5 S5     Eye         C      None         1    Bed     17983
6 S6     Eye         C      None         2    Bed     18471


> obj <- dba.count(obj)
Error: Error processing one or more read files. Check warnings().
In addition: Warning messages:
1:  
2:  
3:  
4:  
5:  
6:  
7:  
8:  
9:  


> warnings()
Warning messages:
1: 
2: 
3: 
4: 
5: 
6: 
7: 
8: 
9: 

I checked and the path of the files in sample sheet is fine. I am using Bed files as input. Thanks a lot in advance for your help!!

Best,
U


Here is the output of sessioninfo if you need:

> sessionInfo()
R version 3.4.1 (2017-06-30)
Platform: x86_64-apple-darwin15.6.0 (64-bit)
Running under: macOS Sierra 10.12.4

Matrix products: default
BLAS: /System/Library/Frameworks/Accelerate.framework/Versions/A/Frameworks/vecLib.framework/Versions/A/libBLAS.dylib
LAPACK: /Library/Frameworks/R.framework/Versions/3.4/Resources/lib/libRlapack.dylib

locale:
[1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8

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

other attached packages:
 [1] DiffBind_2.4.8             SummarizedExperiment_1.6.4
 [3] DelayedArray_0.2.7         matrixStats_0.52.2        
 [5] Biobase_2.36.2             GenomicRanges_1.28.5      
 [7] GenomeInfoDb_1.12.2        IRanges_2.10.3            
 [9] S4Vectors_0.14.4           BiocGenerics_0.22.0       

loaded via a namespace (and not attached):
 [1] edgeR_3.18.1             bit64_0.9-7              splines_3.4.1           
 [4] gtools_3.5.0             assertthat_0.2.0         latticeExtra_0.6-28     
 [7] amap_0.8-14              RBGL_1.52.0              blob_1.1.0              
[10] GenomeInfoDbData_0.99.0  Rsamtools_1.28.0         ggrepel_0.6.5           
[13] Category_2.42.1          RSQLite_2.0              backports_1.1.0         
[16] lattice_0.20-35          glue_1.1.1               limma_3.32.6            
[19] digest_0.6.12            RColorBrewer_1.1-2       XVector_0.16.0          
[22] checkmate_1.8.3          colorspace_1.3-2         Matrix_1.2-11           
[25] plyr_1.8.4               GSEABase_1.38.1          XML_3.98-1.9            
[28] pkgconfig_2.0.1          pheatmap_1.0.8           ShortRead_1.34.1        
[31] biomaRt_2.32.1           genefilter_1.58.1        zlibbioc_1.22.0         
[34] xtable_1.8-2             GO.db_3.4.1              scales_0.5.0            
[37] brew_1.0-6               gdata_2.18.0             BiocParallel_1.10.1     
[40] tibble_1.3.4             annotate_1.54.0          ggplot2_2.2.1           
[43] GenomicFeatures_1.28.4   lazyeval_0.2.0           magrittr_1.5            
[46] survival_2.41-3          memoise_1.1.0            systemPipeR_1.10.2      
[49] fail_1.3                 gplots_3.0.1             hwriter_1.3.2           
[52] GOstats_2.42.0           graph_1.54.0             tools_3.4.1             
[55] BBmisc_1.11              stringr_1.2.0            sendmailR_1.2-1         
[58] munsell_0.4.3            locfit_1.5-9.1           bindrcpp_0.2            
[61] AnnotationDbi_1.38.2     Biostrings_2.44.2        compiler_3.4.1          
[64] caTools_1.17.1           rlang_0.1.2              grid_3.4.1              
[67] RCurl_1.95-4.8           rjson_0.2.15             AnnotationForge_1.18.2  
[70] bitops_1.0-6             base64enc_0.1-3          gtable_0.2.0            
[73] DBI_0.7                  R6_2.2.2                 GenomicAlignments_1.12.2
[76] dplyr_0.7.3              rtracklayer_1.36.4       bit_1.1-12              
[79] bindr_0.1                KernSmooth_2.23-15       stringi_1.1.5           
[82] BatchJobs_1.6            Rcpp_0.12.12     


 

DiffBind dba.count Chip-Seq Differential analysis • 1.1k views
ADD COMMENT
0
Entering edit mode
Rory Stark ★ 4.1k
@rory-stark-5741
Last seen 2 days ago
CRUK, Cambridge, UK

So the warning messages are all blank? I've never seen that one before!

Could you run the following script and report the output?

> obj <- dba(sampleSheet='analysis_H3K4me3_samples.csv')
> obj$class[10:11,]
> file.exists(unique(obj$class[10,]))
> file.exists(unique(obj$class[11,]))

Regards-

Rory

ADD COMMENT
0
Entering edit mode

Thanks so much Rory! Please find the output below. Everything seems okay though. I am really stuck at this step and your help is very much appreciated.

I can share the data etc. with you if you wish to have a look. Thanks again for your help!

> obj <- dba(sampleSheet='liver_H3K27Ac_samples.csv')

> obj$class[10:11,]

           S1                                                                                                                                                 
bamRead    "/Volumes/ChIP-seq/Bed_resample/S1_H3K27ac.bed"
bamControl "/Volumes/ChIP-seq/Bed_resample/S1_H3.bed"     
           S2                                                                                                                                                 
bamRead    "/Volumes/ChIP-seq/Bed_resample/S2_H3K27ac.bed"
bamControl "/Volumes/ChIP-seq/Bed_resample/S2_H3.bed"     
           S3                                                                                                                                                 
bamRead    "/Volumes/ChIP-seq/Bed_resample/S3_H3K27ac.bed"
bamControl "/Volumes/ChIP-seq/Bed_resample/S3_H3.bed"     
           S4                                                                                                                                                 
bamRead    "/Volumes/ChIP-seq/Bed_resample/S4_H3K27ac.bed"
bamControl "/Volumes/ChIP-seq/Bed_resample/S4_H3.bed"     
           S5                                                                                                                                                 
bamRead    "/Volumes/ChIP-seq/Bed_resample/S5_H3K27ac.bed"
bamControl "/Volumes/ChIP-seq/Bed_resample/S5_H3.bed"     
           S6                                                                                                                                                 
bamRead    "/Volumes/ChIP-seq/Bed_resample/S6_H3K27ac.bed"
bamControl "/Volumes/ChIP-seq/Bed_resample/S6_H3.bed"     
> file.exists(unique(obj$class[10,]))
[1] TRUE TRUE TRUE TRUE TRUE TRUE
> file.exists(unique(obj$class[11,]))
[1] TRUE TRUE TRUE TRUE TRUE TRUE

ADD REPLY
0
Entering edit mode

The next step is to check the bed files to make sure they are well-formed. You are using 12 files, and there is some issue with 9 of them. I expect the warnings to state the name of the file which makes it easy to see which file shave a problem, but in this case we are seeing completely blank warning messages (!).

We could start with seeing the head (first six lines or so) of each bed file to make sure it is well-formed. If the issue is somewhere deep in the bed files, we'll need the whole file. If you could give me access to at least four of your bed files we can be sure at least one of them has an issue. Could you put these up somewhere I can access them, as well as the DBA object obj.

It is rare these days to see sequencing alignments in a bed file, indeed we have been considering removing support for bed alignments in an upcoming release.

-Rory

ADD REPLY
0
Entering edit mode
urjaswita ▴ 40
@urjaswita-13128
Last seen 19 months ago

Thanks Rory! Okay, so it seems to be the issue with Bed alignment file. I am not sure it is the file format or those files in particular, but I converted them in Bam and everything works now. Thanks again for your help!

ADD COMMENT
0
Entering edit mode

Hi

I am getting same error messages. Can you please help. Thanks.

dbObj <- dba.count(dbObj)

Computing summits...

Error: Error processing one or more read files. Check warnings().

In addition: There were 50 or more warnings (use warnings() to see the first 50)

dbObj$class[10:11,]

XX_ChipH3Ia_DA_Rep1 bamRead "/XX_ChipH3Ia_DA_Rep1.bam" bamControl "/XX_input_DA_Rep1.bam" XX_ChipH3Ia_DA_Rep2 bamRead "/XX_ChipH3Ia_DA_Rep2.bam" bamControl "/XX_input_DA_Rep2.bam" XX_ChipH3Ia_DB_Rep1 bamRead "/XX_ChipH3Ia_DB_Rep1.bam" bamControl "/XX_input_DB_Rep1.bam" XX_ChipH3Ia_DB_Rep2 bamRead "/XX_ChipH3Ia_DB_Rep2.bam" bamControl "/XX_input_DB_Rep2.bam" XX_ChipH3Ia_DC_Rep1 bamRead "/XX_ChipH3Ia_DC_Rep1.bam" bamControl "/XX_input_DC_Rep1.bam" XX_ChipH3Ia_DC_Rep2 bamRead "/XX_ChipH3Ia_DC_Rep2.bam" bamControl "/XX_input_DC_Rep2.bam" XX_ChipH3Ib_DA_Rep1 bamRead "/XX_ChipH3Ib_DA_Rep1.bam" bamControl "/XX_input_DA_Rep1.bam" XX_ChipH3Ib_DA_Rep2 bamRead "/XX_ChipH3Ib_DA_Rep2.bam" bamControl "/XX_input_DA_Rep2.bam" XX_ChipH3Ib_DB_Rep1 bamRead "/XX_ChipH3Ib_DB_Rep1.bam" bamControl "/XX_input_DB_Rep1.bam" XX_ChipH3Ib_DB_Rep2 bamRead "/XX_ChipH3Ib_DB_Rep2.bam" bamControl "/XX_input_DB_Rep2.bam" XX_ChipH3Ib_DC_Rep1 bamRead "/XX_ChipH3Ib_DC_Rep1.bam" bamControl "/XX_input_DC_Rep1.bam" XX_ChipH3Ib_DC_Rep2 bamRead "/XX_ChipH3Ib_DC_Rep2.bam" bamControl "/XX_input_DC_Rep2.bam" XX_ChipH3Ie_DA_Rep1 bamRead "/XX_ChipH3Ie_DA_Rep1.bam" bamControl "/XX_input_DA_Rep1.bam" XX_ChipH3Ie_DA_Rep2 bamRead "/XX_ChipH3Ie_DA_Rep2.bam" bamControl "/XX_input_DA_Rep2.bam" XX_ChipH3Ie_DB_Rep1 bamRead "/XX_ChipH3Ie_DB_Rep1.bam" bamControl "/XX_input_DB_Rep1.bam" XX_ChipH3Ie_DB_Rep2 bamRead "/XX_ChipH3Ie_DB_Rep2.bam" bamControl "/XX_input_DB_Rep2.bam" XX_ChipH3Ie_DC_Rep1 bamRead "/XX_ChipH3Ie_DC_Rep1.bam" bamControl "/XX_input_DC_Rep1.bam" XX_ChipH3Ie_DC_Rep2 bamRead "/XX_ChipH3Ie_DC_Rep2.bam" bamControl "/XX_input_DC_Rep2.bam" XX_ChipH3Ic_DA_Rep1 bamRead "/XX_ChipH3Ic_DA_Rep1.bam" bamControl "/XX_input_DA_Rep1.bam" XX_ChipH3Ic_DA_Rep2 bamRead "/XX_ChipH3Ic_DA_Rep2.bam" bamControl "/XX_input_DA_Rep2.bam" XX_ChipH3Ic_DB_Rep1 bamRead "/XX_ChipH3Ic_DB_Rep1.bam" bamControl "/XX_input_DB_Rep1.bam" XX_ChipH3Ic_DB_Rep2 bamRead "/XX_ChipH3Ic_DB_Rep2.bam" bamControl "/XX_input_DB_Rep2.bam" XX_ChipH3Ic_DC_Rep1 bamRead "/XX_ChipH3Ic_DC_Rep1.bam" bamControl "/XX_input_DC_Rep1.bam" XX_ChipH3Ic_DC_Rep2 bamRead "/XX_ChipH3Ic_DC_Rep2.bam" bamControl "/XX_input_DC_Rep2.bam" XX_ChipH3Id_DA_Rep1 bamRead "/XX_ChipH3Id_DA_Rep1.bam" bamControl "/XX_input_DA_Rep1.bam" XX_ChipH3Id_DA_Rep2 bamRead "/XX_ChipH3Id_DA_Rep2.bam" bamControl "/XX_input_DA_Rep2.bam" XX_ChipH3Id_DB_Rep1 bamRead "/XX_ChipH3Id_DB_Rep1.bam" bamControl "/XX_input_DB_Rep1.bam" XX_ChipH3Id_DB_Rep2 bamRead "/XX_ChipH3Id_DB_Rep2.bam" bamControl "/XX_input_DB_Rep2.bam" XX_ChipH3Id_DC_Rep1 bamRead "/XX_ChipH3Id_DC_Rep1.bam" bamControl "/XX_input_DC_Rep1.bam" XX_ChipH3Id_DC_Rep2 bamRead "/XX_ChipH3Id_DC_Rep2.bam" bamControl "/XX_input_DC_Rep2.bam" YY_ChipH3Ia_DA_Rep1 bamRead "/YY_ChipH3Ia_DA_Rep1.bam" bamControl "/YY_input_DA_Rep1.bam" YY_ChipH3Ia_DA_Rep2 bamRead "/YY_ChipH3Ia_DA_Rep2.bam" bamControl "/YY_input_DA_Rep2.bam" YY_ChipH3Ia_DB_Rep1 bamRead "/YY_ChipH3Ia_DB_Rep1.bam" bamControl "/YY_input_DB_Rep1.bam" YY_ChipH3Ia_DB_Rep2 bamRead "/YY_ChipH3Ia_DB_Rep2.bam" bamControl "/YY_input_DB_Rep2.bam" YY_ChipH3Ia_DC_Rep1 bamRead "/YY_ChipH3Ia_DC_Rep1.bam" bamControl "/YY_input_DC_Rep1.bam" YY_ChipH3Ia_DC_Rep2 bamRead "/YY_ChipH3Ia_DC_Rep2.bam" bamControl "/YY_input_DC_Rep2.bam" YY_ChipH3Ib_DA_Rep1 bamRead "/YY_ChipH3Ib_DA_Rep1.bam" bamControl "/YY_input_DA_Rep1.bam" YY_ChipH3Ib_DA_Rep2 bamRead "/YY_ChipH3Ib_DA_Rep2.bam" bamControl "/YY_input_DA_Rep2.bam" YY_ChipH3Ib_DB_Rep1 bamRead "/YY_ChipH3Ib_DB_Rep1.bam" bamControl "/YY_input_DB_Rep1.bam" YY_ChipH3Ib_DB_Rep2 bamRead "/YY_ChipH3Ib_DB_Rep2.bam" bamControl "/YY_input_DB_Rep2.bam" YY_ChipH3Ib_DC_Rep1 bamRead "/YY_ChipH3Ib_DC_Rep1.bam" bamControl "/YY_input_DC_Rep1.bam" YY_ChipH3Ib_DC_Rep2 bamRead "/YY_ChipH3Ib_DC_Rep2.bam" bamControl "/YY_input_DC_Rep2.bam" YY_ChipH3Ie_DA_Rep1 bamRead "/YY_ChipH3Ie_DA_Rep1.bam" bamControl "/YY_input_DA_Rep1.bam" YY_ChipH3Ie_DA_Rep2 bamRead "/YY_ChipH3Ie_DA_Rep2.bam" bamControl "/YY_input_DA_Rep2.bam" YY_ChipH3Ie_DB_Rep1 bamRead "/YY_ChipH3Ie_DB_Rep1.bam" bamControl "/YY_input_DB_Rep1.bam" YY_ChipH3Ie_DB_Rep2 bamRead "/YY_ChipH3Ie_DB_Rep2.bam" bamControl "/YY_input_DB_Rep2.bam" YY_ChipH3Ie_DC_Rep1 bamRead "/YY_ChipH3Ie_DC_Rep1.bam" bamControl "/YY_input_DC_Rep1.bam" YY_ChipH3Ie_DC_Rep2 bamRead "/YY_ChipH3Ie_DC_Rep2.bam" bamControl "/YY_input_DC_Rep2.bam" YY_ChipH3Ic_DA_Rep1 bamRead "/YY_ChipH3Ic_DA_Rep1.bam" bamControl "/YY_input_DA_Rep1.bam" YY_ChipH3Ic_DA_Rep2 bamRead "/YY_ChipH3Ic_DA_Rep2.bam" bamControl "/YY_input_DA_Rep2.bam" YY_ChipH3Ic_DB_Rep1 bamRead "/YY_ChipH3Ic_DB_Rep1.bam" bamControl "/YY_input_DB_Rep1.bam" YY_ChipH3Ic_DB_Rep2 bamRead "/YY_ChipH3Ic_DB_Rep2.bam" bamControl "/YY_input_DB_Rep2.bam" YY_ChipH3Ic_DC_Rep1 bamRead "/YY_ChipH3Ic_DC_Rep1.bam" bamControl "/YY_input_DC_Rep1.bam" YY_ChipH3Ic_DC_Rep2 bamRead "/YY_ChipH3Ic_DC_Rep2.bam" bamControl "/YY_input_DC_Rep2.bam" YY_ChipH3Id_DA_Rep1 bamRead "/YY_ChipH3Id_DA_Rep1.bam" bamControl "/YY_input_DA_Rep1.bam" YY_ChipH3Id_DA_Rep2 bamRead "/YY_ChipH3Id_DA_Rep2.bam" bamControl "/YY_input_DA_Rep2.bam" YY_ChipH3Id_DB_Rep1 bamRead "/YY_ChipH3Id_DB_Rep1.bam" bamControl "/YY_input_DB_Rep1.bam" YY_ChipH3Id_DB_Rep2 bamRead "/YY_ChipH3Id_DB_Rep2.bam" bamControl "/YY_input_DB_Rep2.bam" YY_ChipH3Id_DC_Rep1 bamRead "/YY_ChipH3Id_DC_Rep1.bam" bamControl "/YY_input_DC_Rep1.bam" YY_ChipH3Id_DC_Rep2 bamRead "/YY_ChipH3Id_DC_Rep2.bam" bamControl "/YY_input_DC_Rep2.bam"

file.exists(unique(dbObj$class[10,]))

[1] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE [16] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE [31] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE [46] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE

file.exists(unique(dbObj$class[11,]))

[1] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE

I checked quality of bam files using

$ samtools quickcheck -v *.bam $

but did not get any error.

sessionInfo()

R version 4.0.3 (2020-10-10) Platform: x86_64-pc-linux-gnu (64-bit) Running under: CentOS Linux 7 (Core)

Matrix products: default BLAS/LAPACK: /spack/apps/linux-centos7-x86_64/gcc-8.3.0/openblas-0.3.8-2no6mfziiclwxb7lstxoos335gnhjpes/lib/libopenblasp-r0.3.8.so

locale: [1] C

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

other attached packages: [1] forcats_0.5.1 stringr_1.4.0 [3] dplyr_1.0.2 purrr_0.3.4 [5] readr_1.4.0 tidyr_1.1.3 [7] tibble_3.0.4 ggplot2_3.3.3 [9] tidyverse_1.3.0 DiffBind_3.0.9 [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.2 [17] IRanges_2.24.1 S4Vectors_0.28.1 [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.1 [7] splines_4.0.3 BiocParallel_1.24.1 amap_0.8-18 [10] digest_0.6.27 invgamma_1.1 GO.db_3.12.1 [13] fansi_0.4.1 SQUAREM_2020.5 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_1.0.0 [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.2.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_1.3 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.3 cli_2.2.0 [79] generics_0.1.0 RSQLite_2.2.2 broom_0.7.5 [82] yaml_2.2.1 fs_1.5.0 bit64_4.0.5 [85] caTools_1.18.0 RBGL_1.66.0 xml2_1.3.2 [88] biomaRt_2.46.0 rstudioapi_0.13 compiler_4.0.3 [91] curl_4.3 png_0.1-7 reprex_1.0.0 [94] stringi_1.5.3 ps_1.5.0 GenomicFeatures_1.42.1 [97] lattice_0.20-41 Matrix_1.3-2 vctrs_0.3.6 [100] pillar_1.4.7 lifecycle_1.0.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.14 Rsamtools_2.6.0 [121] GenomeInfoDbData_1.2.4 hms_0.5.3 grid_4.0.3 [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.10 numDeriv_2016.8-1.1 >

ADD REPLY

Login before adding your answer.

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