From ca7b513c2c93e70ae3dd888cfc67aaa52b5cc28b Mon Sep 17 00:00:00 2001
From: Pjotr Prins
Date: Mon, 28 Sep 2020 10:57:58 +0100
Subject: Show computed variance

---
 src/debug.cpp    | 2 +-
 src/gemma_io.cpp | 5 ++++-
 2 files changed, 5 insertions(+), 2 deletions(-)

(limited to 'src')

diff --git a/src/debug.cpp b/src/debug.cpp
index 3b92131..8ef8797 100644
--- a/src/debug.cpp
+++ b/src/debug.cpp
@@ -163,7 +163,7 @@ void write(const char *s, const char *msg) {
 
 void write(const double d, const char *msg) {
   if (!is_debug_data_mode()) return;
-  cout << std::setprecision(3) << d << ": " << msg << endl;
+  cout << std::setprecision(6) << d << ": " << msg << endl;
 }
 
 void write(const gsl_vector *v, const char *msg) {
diff --git a/src/gemma_io.cpp b/src/gemma_io.cpp
index 9c10e62..55f0aad 100644
--- a/src/gemma_io.cpp
+++ b/src/gemma_io.cpp
@@ -1510,7 +1510,10 @@ bool BimbamKin(const string file_geno, const set<string> ksnps,
       gsl_vector_scale(geno, 1.0 / sqrt(geno_var));
     }
 
-    if (ns_test<1) write(geno,"geno z-scored");
+    if (ns_test<1) {
+      write(geno_var,"geno var");
+      write(geno,"geno z-scored");
+    }
 
     // set the SNP column ns_test
     gsl_vector_view Xlarge_col = gsl_matrix_column(Xlarge, ns_test % msize);
-- 
cgit v1.2.3