aboutsummaryrefslogtreecommitdiff
path: root/src
diff options
context:
space:
mode:
authorPjotr Prins2020-09-28 10:57:58 +0100
committerPjotr Prins2020-09-28 10:57:58 +0100
commitca7b513c2c93e70ae3dd888cfc67aaa52b5cc28b (patch)
tree35648cb5b14fd6fc2e5e9100e74867c839b82f36 /src
parent55d19a2dccc0227dcf02eb55733eada1a34f382a (diff)
downloadpangemma-ca7b513c2c93e70ae3dd888cfc67aaa52b5cc28b.tar.gz
Show computed variance
Diffstat (limited to 'src')
-rw-r--r--src/debug.cpp2
-rw-r--r--src/gemma_io.cpp5
2 files changed, 5 insertions, 2 deletions
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);