show_2x_range=(not show_points),
show_legend=(not show_points),
plot_regression=(not show_points),
show_count=(not show_points),
show_correlation=False,
show_plot=False,
ax=ax)
#yzer.xlim(ax, min(data[xcolname]), max(data[xcolname]))
#yzer.ylim(ax, min(data[ycolname]), max(data[ycolname]))
yzer.save_plot(
yzer.get_filename(
img_dirpath, '{0}_in_{1}_vs_KLA-Dex_by_{2}.png'.format(
main, basal_cond, compare)))
yzer.show_plot()
if False:
for main, compare, basal_cond in (('p65', 'GR', 'KLA'), ('GR', 'p65',
'Dex')):
data = yzer.import_file(
yzer.get_filename(dirpath, 'motifs', 'from_peaks',
'{0}_kla_dex_vectors.txt'.format(main)))
data = data.fillna(1)
data = data.groupby(['id', 'chr_name'], as_index=False).mean()
data['ratio_main'] = data['tag_count'] / data['tag_count_2']
data['ratio_compare'] = data['tag_count_3'] / data['tag_count_4']
ax = yzer.scatterplot(
data,
key = 'p65_kla_tag_count'
datasets = [down[key],nc[key],up[key]]
datasets = [d['p65_kla_dex_tag_count'] - d[key] for d in [down, nc, up]]
#title = 'Tags in p65 peaks in KLA 1h + DMSO 2h: Distal'
title = 'Diff in tags in p65 peaks in KLA 1h + Dex 2h vs KLA 1h + DMSO 2h: RefSeq'
ax = grapher.boxplot(datasets,
['Down in Dex 2h','No change in Dex 2h','Up in Dex 2h',],
title=title,
xlabel='Condition',
ylabel='Total tags in all peaks overlapping transcript',
show_outliers=False, show_plot=False)
grapher.save_plot(grapher.get_filename(base_dirpath, 'boxplots', 'dex',
title.replace(' ','_')))
grapher.show_plot()
for sub in datasets: print sub.mean()
# Boxplots for gr_kla_dex peaks by lfc in Dex
if False:
#data = data[data['distal'] == 't']
#data = data[data['has_refseq'] == 1]
trans = data[(data['kla_1_lfc'] >= 1) & (data['dex_over_kla_1_lfc'] <= -.58)]
rest = data[(data['kla_1_lfc'] < 1) | (data['dex_over_kla_1_lfc'] > -.58)]
key = 'gr_dex_tag_count'
datasets = [rest[key],trans[key]]
datasets = [d['gr_kla_dex_tag_count'] - d[key] for d in [rest, trans]]
title = 'Tags in p65 peaks in KLA 1h + Dex 2h: Distal'
