def test_pooled_caller(self):
peak_caller = idr.spp_peak_caller(spp_path)
control_replicates = idr.tagalign_split(self.control, nfiles=2)
experimental_replicates = idr.tagalign_split(self.experimental,
nfiles=2)
peaks = idr.call_peaks_on_pooled_replicates(control_replicates,
experimental_replicates,
peak_caller)
self.assertTrue(idr.file_exists(peaks))
self.assertTrue(idr.is_peak_file(peaks))
def test_pooled_pseudo_replicate_caller(self):
peak_caller = idr.spp_peak_caller(spp_path)
control_replicates = idr.tagalign_split(self.control, nfiles=2)
experimental_replicates = idr.tagalign_split(self.experimental,
nfiles=2)
with cluster_view(scheduler, queue, jobs) as view:
peaks = idr.call_peaks_on_pooled_pseudoreplicates(control_replicates,
experimental_replicates,
peak_caller, view.map)
self.assertTrue(all(map(idr.file_exists, peaks)))
self.assertTrue(all(map(idr.is_peak_file, peaks)))
def test_spp_caller(self):
peak_caller = idr.spp_peak_caller(spp_path, cores=1)
peaks = peak_caller(self.control, self.experimental)
self.assertTrue(idr.file_exists(peaks))
self.assertTrue(idr.is_peak_file(peaks))
