Mercurial > repos > mvdbeek > dapars

comparison dapars.py @ 19:ed151db39c7e draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository https://github.com/mvdbeek/dapars commit b1b007c561ea6c9db145c88b6b128d66ecd05e24-dirty
author mvdbeek
date Fri, 30 Oct 2015 11:36:36 -0400
parents 917a2f7ab841
children 71db297a0edd
comparison
equal deleted inserted replaced
18:e041f9f3271f 19:ed151db39c7e
325 ax2.set_title("coverage plot") 325 ax2.set_title("coverage plot")
326 ax2.set_xlabel("nt after UTR start") 326 ax2.set_xlabel("nt after UTR start")
327 ax2.set_ylabel("normalized nucleotide coverage") 327 ax2.set_ylabel("normalized nucleotide coverage")
328 mse_control = [ condition[0] for condition in mse_list] 328 mse_control = [ condition[0] for condition in mse_list]
329 mse_treatment = [ condition[1] for condition in mse_list] 329 mse_treatment = [ condition[1] for condition in mse_list]
330 minima_control = get_minima(np.array(mse_control)) 330 #minima_control = get_minima(np.array(mse_control))
331 minima_treatment = get_minima(np.array(mse_treatment)) 331 #minima_treatment = get_minima(np.array(mse_treatment))
332 control = normalized_utr_coverage[:num_control] 332 control = normalized_utr_coverage[:num_control]
333 treatment = normalized_utr_coverage[num_control:] 333 treatment = normalized_utr_coverage[num_control:]
334 ax1.plot(mse_control, "b-") 334 ax1.plot(mse_control, "b-")
335 ax1.plot(mse_treatment, "r-") 335 ax1.plot(mse_treatment, "r-")
336 [ax2.plot(cov, "b-") for cov in control] 336 [ax2.plot(cov, "b-") for cov in control]
337 [ax2.plot(cov, "r-") for cov in treatment] 337 [ax2.plot(cov, "r-") for cov in treatment]
338 [ax2.axvline(val, color="b", alpha=0.25) for val in minima_control] 338 #[ax2.axvline(val, color="b", alpha=0.25) for val in minima_control]
339 ax2.axvline(mse_control.index(min(mse_control)), color="b", alpha=1) 339 ax2.axvline(mse_control.index(min(mse_control)), color="b", alpha=1)
340 [ax2.axvline(val, color="r", alpha=0.25) for val in minima_treatment] 340 #[ax2.axvline(val, color="r", alpha=0.25) for val in minima_treatment]
341 ax2.axvline(mse_treatment.index(min(mse_treatment)), color="r", alpha=1) 341 ax2.axvline(mse_treatment.index(min(mse_treatment)), color="r", alpha=1)
342 fig.add_subplot(ax1) 342 fig.add_subplot(ax1)
343 fig.add_subplot(ax2) 343 fig.add_subplot(ax2)
344 gs.tight_layout(fig) 344 gs.tight_layout(fig)
345 fig.savefig(os.path.join(plot_path, "{utr}.svg".format(utr=utr))) 345 fig.savefig(os.path.join(plot_path, "{utr}.svg".format(utr=utr)))