def filter_peaks_based_on_peak_shape_worker2(max_peaks, i_list, rt_list,
max_peaks_all):
max_fold_change_score = -1
best_rt = -1
best_i = -1
for rt, i in max_peaks:
rt_left, rt_right = chrom.get_peak_boundary(rt_list, i_list, rt)
i_apex = float(i)
i_left = peak_groups.get_intensity_for_closest_rt(
rt_left, rt_list, i_list)
i_right = peak_groups.get_intensity_for_closest_rt(
rt_right, rt_list, i_list)
fold_change_left = i_apex / (i_left + 1.0)
fold_change_right = i_apex / (i_right + 1.0)
fold_change_score = fold_change_left * fold_change_right
if fold_change_score > max_fold_change_score:
max_fold_change_score = fold_change_score
best_rt = rt
best_i = i
max_peaks_all.append((best_rt, best_i))
return max_peaks_all
def filter_peaks_based_on_peak_shape_worker2(max_peaks, i_list, rt_list, max_peaks_all):
max_fold_change_score = -1
best_rt = -1
best_i = -1
for rt, i in max_peaks:
rt_left, rt_right = chrom.get_peak_boundary(rt_list, i_list, rt)
i_apex = float(i)
i_left = peak_groups.get_intensity_for_closest_rt(rt_left, rt_list, i_list)
i_right = peak_groups.get_intensity_for_closest_rt(rt_right, rt_list, i_list)
fold_change_left = i_apex / (i_left + 1.0)
fold_change_right = i_apex / (i_right + 1.0)
fold_change_score = fold_change_left * fold_change_right
if fold_change_score > max_fold_change_score:
max_fold_change_score = fold_change_score
best_rt = rt
best_i = i
max_peaks_all.append((best_rt, best_i))
return max_peaks_all
def filter_peaks_based_on_peak_shape_worker(max_peaks, i_list, rt_list, max_peaks_all):
for rt, i in max_peaks:
rt_left, rt_right = chrom.get_peak_boundary(rt_list, i_list, rt)
i_apex = float(i)
i_left = peak_groups.get_intensity_for_closest_rt(rt_left, rt_list, i_list)
i_right = peak_groups.get_intensity_for_closest_rt(rt_right, rt_list, i_list)
fold_change_left = i_apex / (i_left + 1.0)
fold_change_right = i_apex / (i_right + 1.0)
if fold_change_left >= parameters.PEAK_SHAPE_FOLD_VARIATION_CRUDE and fold_change_right >= parameters.PEAK_SHAPE_FOLD_VARIATION_CRUDE:
max_peaks_all.append((rt, i))
return max_peaks_all
def filter_peaks_based_on_peak_shape_worker(max_peaks, i_list, rt_list,
max_peaks_all):
for rt, i in max_peaks:
rt_left, rt_right = chrom.get_peak_boundary(rt_list, i_list, rt)
i_apex = float(i)
i_left = peak_groups.get_intensity_for_closest_rt(
rt_left, rt_list, i_list)
i_right = peak_groups.get_intensity_for_closest_rt(
rt_right, rt_list, i_list)
fold_change_left = i_apex / (i_left + 1.0)
fold_change_right = i_apex / (i_right + 1.0)
if fold_change_left >= parameters.PEAK_SHAPE_FOLD_VARIATION_CRUDE and fold_change_right >= parameters.PEAK_SHAPE_FOLD_VARIATION_CRUDE:
max_peaks_all.append((rt, i))
return max_peaks_all
