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)
