def test_load_peak(self):
     #peaks_path = 'H3K4me3_H1hesc_VS_K562_all_peak_MAvalues_3peaks.xls'
     peaks_path = 'Ese14_H3K27ac_LICR_Rep2_peaks.xls'
     genome_path = '/home/jiawei/.MotifScan/genome/mm9'
     gene_path = '/home/jiawei/.MotifScan/gene/mm9/refSeq.txt'
     print "loading genes"
     gene_table = peak.load_ref_gene(gene_path)
     print "loading peaks"
     peak_table =  peak.load_peak(peaks_path,genome_path,1000,'macs')
     print "generate randoms"
     peak.generate_random_with_ref2(gene_table,peak_table,genome_path,1)
 def test_generate_random_with_ref(self):
     peak.generate_random_with_ref2(self.gene_table,self.peak_table,self.genome_path,5)