HAP-240 avoid division by zero

Illumina · Oct 5, 2016 · 03a8e0a · 03a8e0a
1 parent 32bda78
commit 03a8e0a
Showing 1 changed file with 4 additions and 0 deletions.
diff --git a/src/python/Tools/ci.py b/src/python/Tools/ci.py
@@ -23,6 +23,10 @@ def jeffreysCI(x, n, alpha=0.05):
     Brown, Cai and DasGupta: Interval Estimation for a Binomial Proportion.
     2001, doi:10.1214/ss/1009213286'''
 
+    # HAP-240 avoid division by zero
+    if n == 0:
+        return 0.0, 0.0, 1.0
+
     p = x / n
     beta = stats.distributions.beta(x+0.5, n-x+0.5)