def test_load_motif(self): #motif_db_path = '/home/jiawei/MAmotif/motifscan_pkg/MotifScan/test/test_motif_db' motif_path = '/mnt/MAmotif/motifscan_pkg/MotifScan/motif/hg19/hg19_jaspar_2014_all_1e-4.txt' motif_list_path = '' motif_table = motif.load_motif(motif_list_path, motif_path)
def test_motif_scan(self): peak_path = 'Ese14_H3K27ac_LICR_Rep2_peaks.xls' p = peak.load_peak(peak_path,'/home/jiawei/.MotifScan/genome/hg19',1000,'macs') m = motif.load_motif('/mnt/MAmotif/motifscan_pkg/MotifScan/motif_list/jaspar_2014_vertebrates_motif_list.txt','/mnt/MAmotif/motifscan_pkg/MotifScan/motif/mm9/mm9_jaspar_2014_all_1e-4.txt') core.motif_scan(p,m,np.array([0.29,0.21,0.21,0.29]),'./tmp')