DiffBind: same numbers in dba.overlap
0
0
Entering edit mode
Udi Landau ▴ 30
@udi-landau-9726
Last seen 4.2 years ago

Hello,

I have two conditions with 2 and 3 biological replicates, I made a dba object like this:

k4 <- dba.count(k4, minOverlap = 2)

k4 <- dba.contrast(k4, categories = DBA_CONDITION, minMembers = 2)

k4 <- dba.analyze(k4, bSubControl = FALSE)

olap.rate <- dba.overlap(k4, mode = DBA_OLAP_RATE)

 olap.rate
[1] 17534 17534 17534 17534 17534

for some reason I get the same number of peaks for each replicate/treatment...

when I use plotVenn 

  dba.plotVenn(k4 ,k4$masks$Consensus & k4$masks$WT)

I get different numbers for each replicate.

Don't know if it matters but I have duplicated peaks in my BED file (with different scores).

Hope that Dr. Stark could help me with that...

> sessionInfo()
R version 3.2.2 (2015-08-14)
Platform: x86_64-apple-darwin13.4.0 (64-bit)
Running under: OS X 10.11.3 (El Capitan)

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] stats4    parallel  stats     graphics  grDevices utils     datasets  methods  
[9] base     

other attached packages:
 [1] DiffBind_1.16.3            RSQLite_1.0.0              DBI_0.3.1                 
 [4] locfit_1.5-9.1             GenomicAlignments_1.6.3    Rsamtools_1.22.0          
 [7] Biostrings_2.38.3          XVector_0.10.0             limma_3.26.7              
[10] SummarizedExperiment_1.0.2 Biobase_2.30.0             GenomicRanges_1.22.4      
[13] GenomeInfoDb_1.6.3         IRanges_2.4.6              S4Vectors_0.8.11          
[16] BiocGenerics_0.16.1        ggplot2_2.0.0             

loaded via a namespace (and not attached):
 [1] Rcpp_0.12.3             lattice_0.20-33         GO.db_3.2.2            
 [4] gtools_3.5.0            digest_0.6.9            plyr_1.8.3             
 [7] futile.options_1.0.0    BatchJobs_1.6           ShortRead_1.28.0       
[10] gplots_2.17.0           zlibbioc_1.16.0         GenomicFeatures_1.22.12
[13] annotate_1.48.0         gdata_2.17.0            Matrix_1.2-3           
[16] checkmate_1.7.0         systemPipeR_1.4.8       labeling_0.3           
[19] GOstats_2.36.0          splines_3.2.2           BiocParallel_1.4.3     
[22] stringr_1.0.0           pheatmap_1.0.8          RCurl_1.95-4.7         
[25] biomaRt_2.26.1          munsell_0.4.2           sendmailR_1.2-1        
[28] rtracklayer_1.30.1      base64enc_0.1-3         BBmisc_1.9             
[31] fail_1.3                edgeR_3.12.0            XML_3.98-1.3           
[34] AnnotationForge_1.12.2  bitops_1.0-6            grid_3.2.2             
[37] RBGL_1.46.0             xtable_1.8-0            GSEABase_1.32.0        
[40] gtable_0.1.2            magrittr_1.5            scales_0.3.0           
[43] graph_1.48.0            KernSmooth_2.23-15      amap_0.8-14            
[46] stringi_1.0-1           hwriter_1.3.2           genefilter_1.52.1      
[49] latticeExtra_0.6-26     futile.logger_1.4.1     brew_1.0-6             
[52] rjson_0.2.15            lambda.r_1.1.7          RColorBrewer_1.1-2     
[55] tools_3.2.2             Category_2.36.0         survival_2.38-3        
[58] AnnotationDbi_1.32.3    colorspace_1.2-6        caTools_1.17.1         

 

diffbind dba.overlap • 1.2k views
ADD COMMENT
0
Entering edit mode

You should change your tag to diffbind so the author sees your question. Also show the output of sessionInfo().

Valerie

ADD REPLY
0
Entering edit mode

Thanks Valerie. 

ADD REPLY

Login before adding your answer.

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