diff options
| author | Peter Carbonetto | 2017-07-07 11:20:56 -0500 |
|---|---|---|
| committer | GitHub | 2017-07-07 11:20:56 -0500 |
| commit | 86e96ede4ff0955bb2d03ac6c1bd7562a3984955 (patch) | |
| tree | 33120540091e7d16b58f389a13949df397535912 /src/lm.cpp | |
| parent | b3747413e6c5c8cd447e979157880676da66a342 (diff) | |
| parent | b9758364059d52e153a9f1b4fcae3bc3f3e68422 (diff) | |
| download | pangemma-86e96ede4ff0955bb2d03ac6c1bd7562a3984955.tar.gz | |
Merge pull request #51 from genenetwork/spacing
Spacing fixes.
Diffstat (limited to 'src/lm.cpp')
| -rw-r--r-- | src/lm.cpp | 6 |
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 + |