is.broadpeak FALSE read C2C12-myocyte-DNAse.pooled.sorted.pseudoRep1.MACS-2.1.0.p1e-1_peaks.narrowPeak.top100K-p-value-sort : 100000 peaks 100000 peaks are left after cleaning C2C12-myocyte-DNAse.pooled.sorted.pseudoRep2.MACS-2.1.0.p1e-1_peaks.narrowPeak.top100K-p-value-sort : 100000 peaks 100000 peaks are left after cleaninghalf.width= reported significant measure= p.value URI is done URI is saved at: IDR-C2C12-myocyte-DNAse.pooled.sorted.pseudoRep-uri.sav EM is done EM is saved at: IDR-C2C12-myocyte-DNAse.pooled.sorted.pseudoRep-em.sav EM estimation for the following files C2C12-myocyte-DNAse.pooled.sorted.pseudoRep1.MACS-2.1.0.p1e-1_peaks.narrowPeak.top100K-p-value-sort C2C12-myocyte-DNAse.pooled.sorted.pseudoRep2.MACS-2.1.0.p1e-1_peaks.narrowPeak.top100K-p-value-sort $p [1] 0.7082079 $rho1 [1] 0.9877519 $rho2 [1] 0 Write overlapped peaks and local idr to: IDR-C2C12-myocyte-DNAse.pooled.sorted.pseudoRep-overlapped-peaks.txt Write number of peaks above IDR cutoff [0.01, 0.25]: npeaks-aboveIDR.txt Marginal mean of two components: $x.mean1 [1] 33724.6 $x.mean2 [1] 72435.25 $y.mean1 [1] 33734.47 $y.mean2 [1] 72411.31 $x.sd1 [1] 21436.05 $x.sd2 [1] 11779.01 $y.sd1 [1] 21487.13 $y.sd2 [1] 11631.11