Avoiding taking log of 0

nowling-lab · Nov 7, 2021 · eba3ec4 · eba3ec4
1 parent ba6909a
commit eba3ec4
Showing 1 changed file with 3 additions and 0 deletions.
diff --git a/bin/asaph_detect_and_localize b/bin/asaph_detect_and_localize
@@ -282,6 +282,9 @@ if __name__ == "__main__":
         df = read_snp_table(args.pca_associations_tsv,
                             args.component,
                             args.chromosome)
+
+        # avoid domain errors from trying to take the log of 0
+        df["pvalue"] = np.maximum(df["pvalue"], np.power(10., -300.))
 
         df = mark_significant_snps(df)