Plot supplementary figure showing Hah et al error rates against MAX EDGE values when Vespucci is built without knowledge of RefSeq boundaries. ''' from __future__ import division from glasslab.dataanalysis.graphing.seq_grapher import SeqGrapher if __name__ == '__main__': yzer = SeqGrapher() dirpath = 'karmel/Desktop/Projects/GlassLab/Notes_and_Reports/Glass Atlas/NAR_review_data/no_refseq' dirpath = yzer.get_path(dirpath) img_dirpath = yzer.get_and_create_path(dirpath, 'plots') ax = yzer.set_up_plot() title = 'Benchmarking without Foreknowledge of RefSeq' yzer.add_title(title, ax) yzer.add_axis_labels('MAX_EDGE value', 'Error rate defined by Hah et al. (%)') max_edges = [100, 500, 1000, 4000, 5000, 5500, 6000, 10000] error_rates = [ 0.388551822833, 0.372390444765, 0.263807982126, 0.124663089396, 0.121784970634, 0.121807917409, 0.123263849815, 0.142530838464 ] error_pcts = [e * 100 for e in error_rates] yzer.plot(max_edges, error_pcts, '-o') yzer.save_plot_with_dir(save_dir=img_dirpath, title=title) yzer.show_plot()
ctl_y, 'o-', color='black', label='Control (n=8)', linewidth=8) pyplot.plot(low_x, low_y, 'o-', color='blue', label='Low Dose TDB (107 ug, n=4)', linewidth=6) pyplot.plot(med_x, med_y, 'o--', color='green', label='Medium Dose TDB (214 ug, n=4)', linewidth=6) pyplot.plot(high_x, high_y, 'o--', color='red', label='High Dose TDB (321 ug, n=4)', linewidth=4) pyplot.legend(loc='lower left') pyplot.xticks(x_vals, x_labels) yzer.ylim(ax, 0, 100) yzer.xlim(ax, 0, max(x_vals)) yzer.save_plot_with_dir(img_dirpath, None, title) yzer.show_plot()