Fix for a rare overflow. [Found by @ipiras]

andreyshabalin · May 9, 2022 · 7e5c1cb · 7e5c1cb
1 parent 681b367
commit 7e5c1cb
Showing 1 changed file with 1 addition and 0 deletions.
diff --git a/R/Matrix_eQTL_engine.R b/R/Matrix_eQTL_engine.R
@@ -1586,6 +1586,7 @@ Matrix_eQTL_main = function(
             threshfun = function(pv){
                 thr = qt(pv/2, dfFull, lower.tail = FALSE);
                 thr = thr^2;
+                thr[is.infinite(thr)] = .Machine$double.xmax;
                 thr = sqrt(  thr / (dfFull + thr) );
                 thr[pv >= 1] = 0;
                 thr[pv <= 0] = 1;