I have 2 Chipseq, including IP and input samples.
for example, the the number of reads in the chr1:1-500bp is 1000 and 50 in IP and input samples, respectively.
I want to compare the statistic differences for that. Does anyone know how to work on that?
thanks in advance.
1 answer
The simplest answer would be to use the Poisson distribution. Assuming your data is normalized for read depth, you have a background estimate for your window of interest (lambda = 50), and an observed count for that window in your IP (1000), so in R you could estimate this with:
ppois(100, lambda=50, lower.tail=FALSE)
However, this is only one window, and there maybe many other factors involved. The MACS paper is an easy read which resembles this approach (Zhang et al. 2008, Genome Biology). However, opinions differ on the best distribution to use, as well as other factors involved. See for example: Kuan et al., (2011). A Statistical Framework for the Analysis of ChIP-Seq Data.
Log in to answer this question.
-
Normalization and Genome-Wide Comparisons with Chipseq
written by zbeddingfield3 •
-
Create bed file where each column is the number of intersections with another bed file
written by rxs1018 •
-
statistic test
written by BenHu •
-
csaw filterWindows unused argument error
written by jennifer.d •
-
Finding the cause for corrupt ChIP-Seq data
written by Aspire
-
Dealing with spike-in reads that overlap between two genomes
written by juanpjlozano •
-
Does MetDiff support paired-end reads input?
written by camelbbs
-
Precision and recall of binned contigs
written by kimsumin94 •
-
Where, When and How to use the "breadth of coverage" values reported per bin in BEDtools coverageBE…
written by Zee_S
-
actual coverage vs. macs2
written by blur
