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)