is.broadpeak FALSE [1] "Fix 1 stops\n" [1] "Fix 2 stops\n" read SPP-300K-G1E-ER4-CTCF-Rep1--vs-Rep2-PSU-PooledPseudoRep1/G1E-ER4-CTCF-Rep1--vs-Rep2-PSU.ChIP.pooled.sorted.pseudoRep1_VS_G1E-ER4-CTCF-Rep1--vs-Rep2-PSU.Control.pooled.sorted.pseudoRep1.regionPeak : 291742 peaks 291742 peaks are left after cleaning SPP-300K-G1E-ER4-CTCF-Rep1--vs-Rep2-PSU-PooledPseudoRep2/G1E-ER4-CTCF-Rep1--vs-Rep2-PSU.ChIP.pooled.sorted.pseudoRep2_VS_G1E-ER4-CTCF-Rep1--vs-Rep2-PSU.Control.pooled.sorted.pseudoRep2.regionPeak : 292171 peaks 292171 peaks are left after cleaninghalf.width= reported significant measure= signal.value URI is done URI is saved at: IDR-SPP-G1E-ER4-CTCF-Rep1--vs-Rep2-PSU-PooledPseudoReps-uri.sav EM is done EM is saved at: IDR-SPP-G1E-ER4-CTCF-Rep1--vs-Rep2-PSU-PooledPseudoReps-em.sav EM estimation for the following files SPP-300K-G1E-ER4-CTCF-Rep1--vs-Rep2-PSU-PooledPseudoRep1/G1E-ER4-CTCF-Rep1--vs-Rep2-PSU.ChIP.pooled.sorted.pseudoRep1_VS_G1E-ER4-CTCF-Rep1--vs-Rep2-PSU.Control.pooled.sorted.pseudoRep1.regionPeak SPP-300K-G1E-ER4-CTCF-Rep1--vs-Rep2-PSU-PooledPseudoRep2/G1E-ER4-CTCF-Rep1--vs-Rep2-PSU.ChIP.pooled.sorted.pseudoRep2_VS_G1E-ER4-CTCF-Rep1--vs-Rep2-PSU.Control.pooled.sorted.pseudoRep2.regionPeak $p [1] 0.480705 $rho1 [1] 0.9513139 $rho2 [1] 0 Write overlapped peaks and local idr to: IDR-SPP-G1E-ER4-CTCF-Rep1--vs-Rep2-PSU-PooledPseudoReps-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] 17861.11 $x.mean2 [1] 47478.69 $y.mean1 [1] 17846.91 $y.mean2 [1] 47491.83 $x.sd1 [1] 12862.42 $x.sd2 [1] 11584.22 $y.sd1 [1] 12780.19 $y.sd2 [1] 11634.75