def test_genome_wide_norm_without_diag(self):
"""Tests functionality to extract ICCF values
genome-wide with median normalization from synthetic Hi-C data
with blanking the main diagonal."""
arms = pd.DataFrame({
"chrom": "chrSyn",
"start": 0,
"end": 4990000
},
index=[0])
cooler_file = cooler.Cooler(
"testFiles/test2.mcool::/resolutions/10000")
pairing_score = HT.get_pairing_score(cooler_file,
50000,
arms=arms,
norm=True,
blank_diag=True)
expected = pd.read_csv(
"testFiles/test_pairingScore_ICCF_genomeWide_withoutDiag_norm.csv",
dtype={
name: pairing_score.dtypes[name]
for name in pairing_score.dtypes.index
},
index_col=0,
)
assert_frame_equal(pairing_score, expected)
def test_specific_regions_without_diag(self): # Simulated cooler
"""Tests functionality to extract ICCF values
at specific regions from synthetic Hi-C data
with blanking the main diagonal."""
position_frame = pd.read_csv("testFiles/posPileups.csv")
position_frame.loc[:, "mid"] = position_frame["pos"]
arms = pd.DataFrame({
"chrom": "chrSyn",
"start": 0,
"end": 4990000
},
index=[0])
cooler_file = cooler.Cooler(
"testFiles/test2.mcool::/resolutions/10000")
pairing_score = HT.get_pairing_score(
cooler_file,
50000,
regions=position_frame,
arms=arms,
norm=False,
blank_diag=True,
)
expected = pd.read_csv(
"testFiles/test_pairingScore_ICCF_specificRegions_withoutDiag.csv")
assert_frame_equal(pairing_score, expected)