why DiffBind, diffferentially binding site is the for different samples
1
0
Entering edit mode
24 months ago
iridha • 0

I have data for 3 histone marks (2 for silencing and 1 for activation) each mark has three replicates. when I run the diffBind package I have three contrast:

Silence1-Activation Silence1-Silence2 Silence2-Activation When I run dba.plotBox for any of contrasts I will have the same distribution of reads.like below enter image description here

I even tried to start with 1000 peaks of each one that showed higher fold change but again got the same pattern for plotbox. Even the outliers are the same, How it's possible?

here is the code:

peaks10kthreemarksobj <- dba(sampleSheet="samplesheets10kPeaks.csv") peaks10kthreemarksobj <- dba.count(peaks10kthreemarksobj, bUseSummarizeOverlaps=TRUE,bParallel=TRUE) peaks10kthreemarksobj <- dba.contrast(peaks10kthreemarksobj, categories=DBAFACTOR) peaks10kthreemarksobj <- dba.analyze(peaks10kthreemarksobj, method=DBAALLMETHODS, bParallel=T) dba.plotBox(peaks10kthreemarksobj, contrast=3)

ChIP-Seq alignment diffBind • 692 views
ADD COMMENT
0
Entering edit mode

Cross-posted to Bioconductor https://support.bioconductor.org/p/132502/

ADD REPLY
0
Entering edit mode
ADD COMMENT

Login before adding your answer.

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