is.broadpeak TRUE read H3K27ac_MB_Rep1-GSM1972399.mm10.36mers.MACS-2.1.0.p1e-1.broad_peaks.broadPeak.top100K-p-value-sort : 78287 peaks 78287 peaks are left after cleaning H3K27ac_MB_Rep2-GSM1972400.mm10.36mers.MACS-2.1.0.p1e-1.broad_peaks.broadPeak.top100K-p-value-sort : 47586 peaks 47586 peaks are left after cleaninghalf.width= reported significant measure= p.value URI is done URI is saved at: IDR-H3K27ac_MB.MACS-2.1.0.p1e-1.broad_peaks-uri.sav EM is done EM is saved at: IDR-H3K27ac_MB.MACS-2.1.0.p1e-1.broad_peaks-em.sav EM estimation for the following files H3K27ac_MB_Rep1-GSM1972399.mm10.36mers.MACS-2.1.0.p1e-1.broad_peaks.broadPeak.top100K-p-value-sort H3K27ac_MB_Rep2-GSM1972400.mm10.36mers.MACS-2.1.0.p1e-1.broad_peaks.broadPeak.top100K-p-value-sort $p [1] 0.4853467 $rho1 [1] 0.06886322 $rho2 [1] 0 Write overlapped peaks and local idr to: IDR-H3K27ac_MB.MACS-2.1.0.p1e-1.broad_peaks-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] 13222.8 $x.mean2 [1] 13603.12 $y.mean1 [1] 13744.49 $y.mean2 [1] 13111.03 $x.sd1 [1] 8422.583 $x.sd2 [1] 7043.796 $y.sd1 [1] 8405.247 $y.sd2 [1] 7054.478