e8a8edec508fcf886d5c6329ff0729b1cff845e8,hicexplorer/utilities.py,,expected_interactions_norm,#Any#Any#Any#,239

Before Change

    for i, distance_ in enumerate(distance):
        expected_interactions[distance_] += pSubmatrix.data[i]

    count_times_i = np.arange(float(len(expected_interactions)))
    pChromosome_count = np.int(pChromosome_count)
    pLength_chromosome = np.int(pLength_chromosome)
    count_times_i *= pChromosome_count
    count_times_i -= pLength_chromosome
    count_times_i *= np.int(-1)
    count_times_i *= np.sqrt(count_times_i)
    expected_interactions /= count_times_i
    log.debug("exp_obs_matrix_lieberman {}".format(expected_interactions))

After Change

    row, col = pSubmatrix.nonzero()
    distance = np.absolute(row - col)

    occurences = np.zeros(pSubmatrix.shape[0])
    for i, distance_ in enumerate(distance):
        expected_interactions[distance_] += pSubmatrix.data[i]
        occurences[distance_] += 1
    expected_interactions /= occurences
    return expected_interactions

def exp_obs_matrix_lieberman(pSubmatrix, pLength_chromosome, pChromosome_count):

Instances