42dcbde0e602cba7adb911a2072df3887afd61f4,hicexplorer/utilities.py,,exp_obs_matrix_lieberman,#Any#Any#Any#,215

Before Change

    expected_interactions_in_distance_ = expected_interactions_in_distance(pLength_chromosome, pChromosome_count, pSubmatrix)
    row, col = pSubmatrix.nonzero()
    distance = np.ceil(np.absolute(row - col) / 2).astype(np.int32)
    for i in range(len(pSubmatrix.data)):
        try:
            if expected_interactions_in_distance_[distance[i]] == 0:
                pSubmatrix.data[i] = 0.0
            else:
                pSubmatrix.data[i] = pSubmatrix.data[i] / expected_interactions_in_distance_[distance[i]]
        except Exception:
            log.debug("pSubmatrix.data[i]: {}".format(pSubmatrix.data[i]))
            log.debug("distance[i]: {}".format(distance[i]))
            log.debug("expected_interactions_in_distance_[distance[i]]: {} ".format(expected_interactions_in_distance_[distance[i]]))
            exit(1)

    return pSubmatrix


def toString(s):

After Change

        expected = expected_interactions_in_distance_[distance]
        // log.info("len(expected); {}".format(len(expected)))
        // log.info("len(data): {}".format(len(pSubmatrix.data)))
        pSubmatrix.data = pSubmatrix.data.astype(np.float32)
        pSubmatrix.data /= expected
        // data2 = deepcopy(data)
        pSubmatrix.data = convertInfsToZeros_ArrayFloat(pSubmatrix.data).astype(data_type)


    // for i in range(len(pSubmatrix.data)):
    //     try:

Instances