def get_duoPeaks(counts, rcounts): """ Given two lists of counts, returns another list identifying the points where both counts are at a peak or plateau """ seg1 = get_peaks(counts) segRev = get_peaks(rcounts, True) segList = list(set(seg1) & set(segRev)) segList.sort() return segList
def get_sumPeaks(counts,rcounts): """ Given two lists of counts, returns another list identifying the points where the sum of the counts is at a peak or plateau """ return get_peaks([counts[i]+rcounts[i+1] for i in range(len(counts)-1)])