Merge pull request #38 from mfiers/master

Fix bonferroni corr. to prevent p>1
tanghaibao · Jun 4, 2015 · b1f1a4c · b1f1a4c
2 parents 1b5bbc8 + 708f56e
commit b1f1a4c
Showing 1 changed file with 3 additions and 0 deletions.
diff --git a/goatools/multiple_testing.py b/goatools/multiple_testing.py
@@ -36,6 +36,9 @@ class Bonferroni(AbstractCorrection):
     def set_correction(self):
         self.corrected_pvals *= self.n
 
+        # reset all pvals > 1 to 1
+        self.corrected_pvals[self.corrected_pvals > 1] = 1
+
 
 class Sidak(AbstractCorrection):
     """http://en.wikipedia.org/wiki/Bonferroni_correction