Ejemplo n.º 1
0
def test_run_idr():
    individual_peaks = glob.glob(os.path.join(peak_dir, 'ctcf_[0-9]_VS_*.gz'))
    individual_peaks = map(_copy_file_to_testdir, individual_peaks)
    pseudo_peaks = glob.glob(os.path.join(peak_dir, 'ctcf_[0-9]_[0-9]_VS_*.gz'))
    pseudo_peaks = map(_copy_file_to_testdir, pseudo_peaks)
    pseudo_peaks.sort()
    pseudo_peaks = idr.bunch(pseudo_peaks)
    pooled_pseudo = glob.glob(os.path.join(peak_dir, 'ctcf_0_pooled_[0-9]_VS_*.gz'))
    pooled_pseudo = map(_copy_file_to_testdir, pooled_pseudo)
    idr_runner = idr.idr_runner(idr_runner_path)
    indiv, pseudo, pool_pseudo = idr.run_idr(individual_peaks, pseudo_peaks,
                                             pooled_pseudo, idr_runner)

    assert all(map(idr.file_exists, indiv))
    assert all(map(idr.file_exists, list(idr.flatten(pooled_pseudo))))
    assert all(map(idr.file_exists, pool_pseudo))
Ejemplo n.º 2
0
def test_post_idr_plot():
    i_files, p_files, pp_files = _get_idr_set()
    idr_plotter = idr.idr_plotter(idr_plotter_path)
    out_files = idr.plot_idr_output(i_files, p_files, pp_files, idr_plotter)
    assert all(map(idr.file_exists, list(idr.flatten(out_files))))