From cf4b9ce2445912db7af04d349ca19a1bda7b9de6 Mon Sep 17 00:00:00 2001
From: Pjotr Prins
Date: Wed, 30 Sep 2020 08:55:06 +0100
Subject: Align R2 filter with docs

---
 src/gemma_io.cpp | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/src/gemma_io.cpp b/src/gemma_io.cpp
index b4f8d39..4bcba62 100644
--- a/src/gemma_io.cpp
+++ b/src/gemma_io.cpp
@@ -824,7 +824,7 @@ bool ReadFile_geno(const string &file_geno, const set<string> &setSnps,
     gsl_blas_ddot(genotype, genotype, &v_x);
     gsl_blas_ddot(Wtx, WtWiWtx, &v_w);
 
-    if (W->size2 != 1 && v_w / v_x >= r2_level) {
+    if (W->size2 != 1 && v_w / v_x > r2_level) {
       indicator_snp.push_back(0);
       continue;
     }
-- 
cgit v1.2.3