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'