Merge pull request #1 from gesinecauer/torm_bug

torm is a filter that removes beads with 0 reads. The lower half of c…
gesinecauer · Mar 13, 2018 · 509a8f5 · 509a8f5
2 parents cd6c547 + 0c5bb6c
commit 509a8f5
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/pastis/algorithms.py b/pastis/algorithms.py
@@ -288,7 +288,7 @@ def run_pm2(directory):
               verbose=options["verbose"])
     X = pm2.fit(counts)
 
-    torm = np.array((counts.sum(axis=0) == 0)).flatten()
+    torm = np.array(((counts + counts.transpose()).sum(axis=0) == 0)).flatten()
 
     X[torm] = np.nan