Mercurial > repos > mvdbeek > dapars

diff dapars.py @ 20:71db297a0edd 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:54:18 -0400
parents ed151db39c7e
children eace1635959a
line wrap: on
line diff
--- a/dapars.py	Fri Oct 30 11:36:36 2015 -0400
+++ b/dapars.py	Fri Oct 30 11:54:18 2015 -0400
@@ -354,6 +354,8 @@
     mse_treatment = np.array([mse[1] for mse in mse_list])
     control_breakpoints = list(get_minima(mse_control, local_minimum))
     treatment_breakpoints = list(get_minima(mse_treatment, local_minimum))
+    print control_breakpoints
+    print treatment_breakpoints
     control_abundances = [estimate_abundance(normalized_utr_coverage, bp, num_samples) for bp in control_breakpoints]
     treatment_abundances = [estimate_abundance(normalized_utr_coverage, bp, num_samples) for bp in treatment_breakpoints]
     return control_breakpoints, control_abundances, treatment_breakpoints, treatment_abundances