e8a8edec508fcf886d5c6329ff0729b1cff845e8,hicexplorer/utilities.py,,exp_obs_matrix_norm,#Any#Any#Any#,287

Before Change

    distance = np.ceil(np.absolute(row - col) / 2).astype(np.int32)

    if len(pSubmatrix.data) > 0:
        data_type = type(pSubmatrix.data[0])

        expected = expected_interactions_in_distance_[distance]
        pSubmatrix.data = pSubmatrix.data.astype(np.float32)
        pSubmatrix.data /= expected

After Change

    expected_interactions_in_distance = expected_interactions_norm(pLength_chromosome, pChromosome_count, pSubmatrix)
    
    row_sums = np.array(pSubmatrix.sum(axis=1).T).flatten()
    total_interactions = pSubmatrix.sum()

    row, col = pSubmatrix.nonzero()
    // data = pSubmatrix.data.tolist()
    for i in range(len(row)):
        expected = expected_interactions_in_distance[np.absolute(row[i]-col[i])] 
        expected /= row_sums[row[i]] * row_sums[col[i]] / total_interactions
        pSubmatrix.data[i] /= expected
    return pSubmatrix

def toString(s):

Instances