def test_no_collapse_real_data(self):
positions = pd.read_csv("testFiles/testAssignRegions.csv")
cooler_file = cooler.Cooler(
"testFiles/test3_realdata.mcool::resolutions/50000")
result = HT.do_pileup_iccf(cooler_file,
positions,
proc=1,
collapse=False)
expected = np.load("testFiles/real_data_iccf_pileup_not_collapsed.npy")
self.assertTrue(np.allclose(result, expected, equal_nan=True))
def test_collapse_real_data(self):
"""tests whether pileup works on real data."""
positions = pd.read_csv("testFiles/testAssignRegions.csv")
cooler_file = cooler.Cooler(
"testFiles/test3_realdata.mcool::resolutions/50000")
result = HT.do_pileup_iccf(cooler_file,
positions,
proc=1,
collapse=True)
expected = np.load("testFiles/real_data_iccf_pileup_collapsed.npy")
self.assertTrue(np.allclose(result, expected))
def test_collapse(self):
"""Tests pileup of synthetic Hi-C data, with collapsing results."""
position_frame = pd.read_csv("testFiles/posPileups.csv")
arms = pd.DataFrame({
"chrom": "chrSyn",
"start": 0,
"end": 4990000
},
index=[0])
assigned = HT.assign_regions(50000, 10000, position_frame["chrom"],
position_frame["pos"], arms)
cooler_file = cooler.Cooler(
"testFiles/test2.mcool::/resolutions/10000")
result = HT.do_pileup_iccf(cooler_file,
assigned,
proc=1,
collapse=True)
expected = np.load("testFiles/test_pileups_iccf_collapse.npy")
self.assertTrue(np.allclose(result, expected))
