is.broadpeak FALSE [1] "Fix 6 stops\n" [1] "Fix 11 stops\n" read /broad/hptmp/anshul/idr_1580616803/dPeak3-P300-K562-Rep1.relaxed.narrowPeak_1580616803 : 75220 peaks 75220 peaks are left after cleaning /broad/hptmp/anshul/idr_1580616803/dPeak3-P300-K562-Rep2.relaxed.narrowPeak_1580616803 : 95247 peaks 95247 peaks are left after cleaninghalf.width= reported significant measure= signal.value URI is done URI is saved at: /broad/compbio/anshul/projects/encode3/evaluation/tfs/phase1/idrResults/origWidth/dPeak3/idr/dPeak3-P300-K562-Rep1_VS_dPeak3-P300-K562-Rep2-uri.sav EM is done EM is saved at: /broad/compbio/anshul/projects/encode3/evaluation/tfs/phase1/idrResults/origWidth/dPeak3/idr/dPeak3-P300-K562-Rep1_VS_dPeak3-P300-K562-Rep2-em.sav EM estimation for the following files /broad/hptmp/anshul/idr_1580616803/dPeak3-P300-K562-Rep1.relaxed.narrowPeak_1580616803 /broad/hptmp/anshul/idr_1580616803/dPeak3-P300-K562-Rep2.relaxed.narrowPeak_1580616803 $p [1] 0.6604739 $rho1 [1] 0.7790013 $rho2 [1] 0 Write overlapped peaks and local idr to: /broad/compbio/anshul/projects/encode3/evaluation/tfs/phase1/idrResults/origWidth/dPeak3/idr/dPeak3-P300-K562-Rep1_VS_dPeak3-P300-K562-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] 16246.6 $x.mean2 [1] 28774.38 $y.mean1 [1] 16586.84 $y.mean2 [1] 28113.67 $x.sd1 [1] 10811.96 $x.sd2 [1] 9032.037 $y.sd1 [1] 11078.2 $y.sd2 [1] 9281.883