aboutsummaryrefslogtreecommitdiff
path: root/src/lm.cpp
diff options
context:
space:
mode:
authorPjotr Prins2017-07-07 06:54:26 +0000
committerPjotr Prins2017-07-07 06:54:26 +0000
commitb9758364059d52e153a9f1b4fcae3bc3f3e68422 (patch)
treecc3b526c1621ca452ded085114d7c40559c09887 /src/lm.cpp
parentdd72b87354d1d3f6d3aa42ed0123a23880e9cb15 (diff)
downloadpangemma-b9758364059d52e153a9f1b4fcae3bc3f3e68422.tar.gz
Fix spacing
Diffstat (limited to 'src/lm.cpp')
-rw-r--r--src/lm.cpp6
1 files changed, 3 insertions, 3 deletions
diff --git a/src/lm.cpp b/src/lm.cpp
index d3ad7f3..94729db 100644
--- a/src/lm.cpp
+++ b/src/lm.cpp
@@ -412,13 +412,13 @@ void LM::Analyzebgen (const gsl_matrix *W, const gsl_vector *y) {
string chr;
std::cout << "Warning: WJA hard coded SNP missingness " <<
"threshold of 10%" << std::endl;
-
+
// Start reading genotypes and analyze.
for (size_t t=0; t<indicator_snp.size(); ++t) {
if (t%d_pace==0 || t==(ns_total-1)) {
ProgressBar ("Reading SNPs ", t, ns_total-1);
}
-
+
// Read SNP header.
id.clear();
rs.clear();
@@ -500,7 +500,7 @@ void LM::Analyzebgen (const gsl_matrix *W, const gsl_vector *y) {
static_cast<double>(unzipped_data[i*3+1])/32768.0;
bgen_geno_prob_BB=
static_cast<double>(unzipped_data[i*3+2])/32768.0;
-
+
// WJA
bgen_geno_prob_non_miss=
bgen_geno_prob_AA +