./sircel/Sircel_master.py

pachterlab · Sep 26, 2017 · 4e7b999 · 4e7b999
1 parent e565d61
commit 4e7b999
Showing 1 changed file with 3 additions and 2 deletions.
diff --git a/sircel/Split_reads.py b/sircel/Split_reads.py
@@ -433,8 +433,9 @@ def get_threshold(grad, second_grad, lmax, num_cells, unique_paths_sorted):
 		threshold = lmax[-1]
 	except IndexError:
 		return len(unique_paths_sorted)
-	for i in lmax:
-		print(i, grad[i], second_grad[i])
+	for i in lmax:		
+		print(i, grad[i - LOCAL_WINDOW_LEN])
+
 
 		if(grad[i] > grad[threshold]):
 			threshold = i