# HG changeset patch # User mvdbeek # Date 1446220458 14400 # Node ID 71db297a0edd426327426c0d78e2040aef93e721 # Parent ed151db39c7e1d5d7fbbc3d2ffa639da9a0f2af9 planemo upload for repository https://github.com/mvdbeek/dapars commit b1b007c561ea6c9db145c88b6b128d66ecd05e24-dirty diff -r ed151db39c7e -r 71db297a0edd dapars.py --- 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