is.broadpeak FALSE read HepG2-CTCF-Rep1-HA.hg20-male.36mers.unique.MACS-2.1.0.p1e-1_peaks.narrowPeak.top100K-p-value-sort : 100000 peaks 100000 peaks are left after cleaning HepG2-CTCF-Rep2-HA.hg20-male.36mers.unique.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-HepG2-CTCF-HA-Rep1-vs-Rep2-uri.sav EM is done EM is saved at: IDR-HepG2-CTCF-HA-Rep1-vs-Rep2-em.sav EM estimation for the following files HepG2-CTCF-Rep1-HA.hg20-male.36mers.unique.MACS-2.1.0.p1e-1_peaks.narrowPeak.top100K-p-value-sort HepG2-CTCF-Rep2-HA.hg20-male.36mers.unique.MACS-2.1.0.p1e-1_peaks.narrowPeak.top100K-p-value-sort $p [1] 0.8264497 $rho1 [1] 0.8336794 $rho2 [1] 0 Write overlapped peaks and local idr to: IDR-HepG2-CTCF-HA-Rep1-vs-Rep2-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] 24821.56 $x.mean2 [1] 47119.33 $y.mean1 [1] 24630.45 $y.mean2 [1] 48025.79 $x.sd1 [1] 15277.72 $x.sd2 [1] 7686.694 $y.sd1 [1] 15085.14 $y.sd2 [1] 6738.035