diff options
author | Pjotr Prins | 2020-09-28 10:57:58 +0100 |
---|---|---|
committer | Pjotr Prins | 2020-09-28 10:57:58 +0100 |
commit | ca7b513c2c93e70ae3dd888cfc67aaa52b5cc28b (patch) | |
tree | 35648cb5b14fd6fc2e5e9100e74867c839b82f36 | |
parent | 55d19a2dccc0227dcf02eb55733eada1a34f382a (diff) | |
download | pangemma-ca7b513c2c93e70ae3dd888cfc67aaa52b5cc28b.tar.gz |
Show computed variance
-rw-r--r-- | src/debug.cpp | 2 | ||||
-rw-r--r-- | src/gemma_io.cpp | 5 |
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); |