From 8bd32ac3ed3d0bf2f7cb01a6d1d9b6af4c6a33c0 Mon Sep 17 00:00:00 2001 From: msbentsen Date: Fri, 12 Apr 2019 10:10:48 +0200 Subject: [PATCH] multi-norm for bindetect --- tobias/footprinting/BINDetect.py | 16 ++++++++++++++++ 1 file changed, 16 insertions(+) diff --git a/tobias/footprinting/BINDetect.py b/tobias/footprinting/BINDetect.py index 09c8ca8..4e4fba3 100644 --- a/tobias/footprinting/BINDetect.py +++ b/tobias/footprinting/BINDetect.py @@ -463,6 +463,22 @@ def run_bindetect(args): figures.append((ax, fig)) + + + #Quantile normalization + ref = args.cond_names[0] + for cond in args.cond_names[1:]: + + ref_quantiles = scipy.stats.lognorm(5, *all_log_params[ref]) + cond_quantiles = scipy.stats.lognorm(5, *all_log_params[cond]) + + print(ref_quantiles) + print(cond_quantiles) + + plt.plot(ref_quantiles, cond_quantiles) + #all_log_params[bigwig] + + #Only plot if args.debug is True if args.debug: #Set x-max of all plots equal