def test_peak_result_test(self):
p1 = pd.read_pickle('MAmotif_test_input/h1hesc_peak_result.pkl')
r1 = pd.read_pickle('MAmotif_test_input/h1hesc_rnd_result.pkl')
e1 = pd.read_pickle('MAmotif_test_input/h1hesc_enrich_result.pkl')
p2 = pd.read_pickle('MAmotif_test_input/k562_peak_result.pkl')
r2 = pd.read_pickle('MAmotif_test_input/k562_rnd_result.pkl')
e2 = pd.read_pickle('MAmotif_test_input/k562_enrich_result.pkl')
p3 = pd.read_pickle('MAmotif_test_input/common_peak_result.pkl')
r3 = pd.read_pickle('MAmotif_test_input/common_rnd_result.pkl')
e3 = pd.read_pickle('MAmotif_test_input/common_enrich_result.pkl')
[p,r,e] = motifscan_general.merge_two_results(p1,p3,r1,r3,e1)
core.peak_result_test(p2)
def test_peak_test(self):
common_motif_dir = '../MAmotif_H3K4me3_H1hesc_VS_K562/motifscan_on_common_peaks'
a_motif_dir = '../MAmotif_H3K4me3_H1hesc_VS_K562/motifscan_on_H1hesc_H3k4me3_Broad_Rep1_P100_unique'
print 'Loading motifscan results...'
common_peak_result = pd.read_pickle("%s/peak_result.pkl"%common_motif_dir)
common_rnd_result = pd.read_pickle("%s/rnd_result.pkl"%common_motif_dir)
common_enrich_result = pd.read_pickle("%s/enrich_result.pkl"%common_motif_dir)
a_peak_result = pd.read_pickle("%s/peak_result.pkl"%a_motif_dir)
a_rnd_result = pd.read_pickle("%s/rnd_result.pkl"%a_motif_dir)
print 'Merging peaks...'
[positive_peak_result,positive_rnd_result,positive_enrich_result] = \
motifscan_general.merge_two_results(common_peak_result,a_peak_result,common_rnd_result,a_rnd_result,common_enrich_result)
core.peak_result_test(positive_peak_result, positive_rnd_result,plot_output_dir='enrichment_plot')
