def stackarray(binfile, chr, window, shift): mn = binfile.mapping[chr] m, n = zip(*mn) m = np.array(m, dtype="float") n = np.array(n, dtype="float") w = window / shift m = moving_sum(m, window=w) n = moving_sum(n, window=w) m /= n return m
def linearray(binfile, chr, window, shift): mn = binfile.mapping[chr] m, n = zip(*mn) m = np.array(m, dtype="float") w = window / shift m = moving_sum(m, window=w) return m