Mercurial > repos > mvdbeek > dapars
comparison dapars.py @ 20:71db297a0edd draft
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 |
comparison
equal
deleted
inserted
replaced
| 19:ed151db39c7e | 20:71db297a0edd |
|---|---|
| 352 """ | 352 """ |
| 353 mse_control = np.array([mse[0] for mse in mse_list]) | 353 mse_control = np.array([mse[0] for mse in mse_list]) |
| 354 mse_treatment = np.array([mse[1] for mse in mse_list]) | 354 mse_treatment = np.array([mse[1] for mse in mse_list]) |
| 355 control_breakpoints = list(get_minima(mse_control, local_minimum)) | 355 control_breakpoints = list(get_minima(mse_control, local_minimum)) |
| 356 treatment_breakpoints = list(get_minima(mse_treatment, local_minimum)) | 356 treatment_breakpoints = list(get_minima(mse_treatment, local_minimum)) |
| 357 print control_breakpoints | |
| 358 print treatment_breakpoints | |
| 357 control_abundances = [estimate_abundance(normalized_utr_coverage, bp, num_samples) for bp in control_breakpoints] | 359 control_abundances = [estimate_abundance(normalized_utr_coverage, bp, num_samples) for bp in control_breakpoints] |
| 358 treatment_abundances = [estimate_abundance(normalized_utr_coverage, bp, num_samples) for bp in treatment_breakpoints] | 360 treatment_abundances = [estimate_abundance(normalized_utr_coverage, bp, num_samples) for bp in treatment_breakpoints] |
| 359 return control_breakpoints, control_abundances, treatment_breakpoints, treatment_abundances | 361 return control_breakpoints, control_abundances, treatment_breakpoints, treatment_abundances |
| 360 | 362 |
| 361 def get_minima(a, local_minimum=False): | 363 def get_minima(a, local_minimum=False): |