def __overlap_fe (self, gold_peaks, sample_peaks, fe_low, fe_up, top, bottom, step):
ret = ["%d-%d" % (fe_low,fe_up)]
gp = PeakIO()
gp.init_from_dict(gold_peaks)
gp.filter_fc(fe_low,fe_up)
gptotal = gp.total()
if gptotal <= 0:
return None
ret.append(gptotal)
for i in xrange(top,bottom,step):
p = PeakIO()
p.init_from_dict(sample_peaks[i])
percent = 100.0*gp.overlap_with_other_peaks(p)/gptotal
ret.append(percent)
del p
return ret
def __overlap_fe(self, gold_peaks, sample_peaks, fe_low, fe_up, top,
bottom, step):
ret = ["%d-%d" % (fe_low, fe_up)]
gp = PeakIO()
gp.init_from_dict(gold_peaks)
gp.filter_fc(fe_low, fe_up)
gptotal = gp.total()
if gptotal <= 0:
return None
ret.append(gptotal)
for i in xrange(top, bottom, step):
p = PeakIO()
p.init_from_dict(sample_peaks[i])
percent = 100.0 * gp.overlap_with_other_peaks(p) / gptotal
ret.append(percent)
del p
return ret
