save_dir=img_dirpath) if False: yzer = MotifAnalyzer() motifs_dirpath = yzer.get_filename(dirpath, 'motifs') # Motif finding for enhancers with SRF bg = yzer.get_filename( motifs_dirpath, 'enhancer_at_least_10_tags_15/enhancer_at_least_10_tags_15_regions_for_homer.txt' ) enh_with_srf = datasets[4] yzer.run_homer(enh_with_srf, 'enhancer_with_srf', motifs_dirpath, cpus=6, center=False, reverse=False, preceding=False, size=200, length=[8, 10, 12, 15], bg=bg) # Motif finding for enhancers with SRF and down yzer.run_homer(enh_with_srf[enh_with_srf['balb_nod_notx_1h_fc'] <= -1], 'enhancer_with_srf_down', motifs_dirpath, cpus=6, center=False, reverse=False, preceding=False, size=200, length=[8, 10, 12, 15],
save_dir=img_dirpath, show_plot=True) if True: # Motif finding. yzer = MotifAnalyzer() motifs_dirpath = yzer.get_and_create_path(dirpath, 'motifs') data['id'] = data.index if False: yzer.run_homer(data, 'th_{0}_enhancers'.format(peak), motifs_dirpath, cpus=6, center=True, reverse=False, preceding=False, size=200, length=[8, 10, 12, 15]) bg = yzer.get_filename( motifs_dirpath, 'th_{0}_enhancers/th_{0}_enhancers_regions_for_homer.txt'.format( peak)) multiple = 2 ''' yzer.run_homer(only_th1, 'th1_only_{0}_enhancers'.format(peak), motifs_dirpath, cpus=6, center=True, reverse=False, preceding=False, size=200, length=[8, 10, 12, 15], bg=bg) yzer.run_homer(only_th2, 'th2_only_{0}_enhancers'.format(peak), motifs_dirpath, cpus=6, center=True, reverse=False, preceding=False, size=200, length=[8, 10, 12, 15], bg=bg)