multi-norm for bindetect

loosolab · Apr 12, 2019 · 8bd32ac · 8bd32ac
1 parent afe25f0
commit 8bd32ac
Showing 1 changed file with 16 additions and 0 deletions.
diff --git 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