is.broadpeak FALSE [1] "Fix 1 stops\n" read H3K27ac_mH2A1i_MT_Rep1-GSM1972405.mm10.36mers.MACS-2.1.0.p1e-1_peaks.narrowPeak.top100K-p-value-sort : 100000 peaks 100000 peaks are left after cleaning H3K27ac_mH2A1i_MT_Rep2-GSM1972406.mm10.36mers.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-H3K27ac_mH2A1i_MT.MACS-2.1.0.p1e-1.narrowPeak-uri.sav EM is done EM is saved at: IDR-H3K27ac_mH2A1i_MT.MACS-2.1.0.p1e-1.narrowPeak-em.sav EM estimation for the following files H3K27ac_mH2A1i_MT_Rep1-GSM1972405.mm10.36mers.MACS-2.1.0.p1e-1_peaks.narrowPeak.top100K-p-value-sort H3K27ac_mH2A1i_MT_Rep2-GSM1972406.mm10.36mers.MACS-2.1.0.p1e-1_peaks.narrowPeak.top100K-p-value-sort $p [1] 0.6065876 $rho1 [1] 0.6767851 $rho2 [1] 0 Write overlapped peaks and local idr to: IDR-H3K27ac_mH2A1i_MT.MACS-2.1.0.p1e-1.narrowPeak-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] 21796.77 $x.mean2 [1] 25649.88 $y.mean1 [1] 23729.36 $y.mean2 [1] 22672.94 $x.sd1 [1] 14556.67 $x.sd2 [1] 11167.05 $y.sd1 [1] 14307.01 $y.sd2 [1] 12005.36