diff options
| author | Peter Carbonetto | 2017-10-09 16:23:05 -0500 |
|---|---|---|
| committer | GitHub | 2017-10-09 16:23:05 -0500 |
| commit | ab2a1830659ac970d97c107f145105fb410eb219 (patch) | |
| tree | 14d3cab14fca4b68eb0916992344b8248a380a61 /src/mvlmm.cpp | |
| parent | 1e99b37e42edb7981b44f1409cc3e55a47057878 (diff) | |
| parent | 86323ccaf26ad0a3b706a67a0014dd04b9965823 (diff) | |
| download | pangemma-ab2a1830659ac970d97c107f145105fb410eb219.tar.gz | |
Merge pull request #92 from genenetwork/gemma-master
Fixed for a number of outstanding issues.
Diffstat (limited to 'src/mvlmm.cpp')
| -rw-r--r-- | src/mvlmm.cpp | 5 |
1 files changed, 5 insertions, 0 deletions
diff --git a/src/mvlmm.cpp b/src/mvlmm.cpp index be9fd78..c5efb6e 100644 --- a/src/mvlmm.cpp +++ b/src/mvlmm.cpp @@ -2953,6 +2953,7 @@ double PCRT(const size_t mode, const size_t d_size, const double p_value, // WJA added. void MVLMM::Analyzebgen(const gsl_matrix *U, const gsl_vector *eval, const gsl_matrix *UtW, const gsl_matrix *UtY) { + debug_msg("entering"); string file_bgen = file_oxford + ".bgen"; ifstream infile(file_bgen.c_str(), ios::binary); if (!infile) { @@ -3501,6 +3502,7 @@ void MVLMM::Analyzebgen(const gsl_matrix *U, const gsl_vector *eval, void MVLMM::AnalyzeBimbam(const gsl_matrix *U, const gsl_vector *eval, const gsl_matrix *UtW, const gsl_matrix *UtY) { + debug_msg("entering"); igzstream infile(file_geno.c_str(), igzstream::in); if (!infile) { cout << "error reading genotype file:" << file_geno << endl; @@ -3938,6 +3940,7 @@ void MVLMM::AnalyzeBimbam(const gsl_matrix *U, const gsl_vector *eval, void MVLMM::AnalyzePlink(const gsl_matrix *U, const gsl_vector *eval, const gsl_matrix *UtW, const gsl_matrix *UtY) { + debug_msg("entering"); string file_bed = file_bfile + ".bed"; ifstream infile(file_bed.c_str(), ios::binary); if (!infile) { @@ -4475,6 +4478,7 @@ void CalcMvLmmVgVeBeta(const gsl_vector *eval, const gsl_matrix *UtW, void MVLMM::AnalyzeBimbamGXE(const gsl_matrix *U, const gsl_vector *eval, const gsl_matrix *UtW, const gsl_matrix *UtY, const gsl_vector *env) { + debug_msg("entering"); igzstream infile(file_geno.c_str(), igzstream::in); if (!infile) { cout << "error reading genotype file:" << file_geno << endl; @@ -4920,6 +4924,7 @@ void MVLMM::AnalyzeBimbamGXE(const gsl_matrix *U, const gsl_vector *eval, void MVLMM::AnalyzePlinkGXE(const gsl_matrix *U, const gsl_vector *eval, const gsl_matrix *UtW, const gsl_matrix *UtY, const gsl_vector *env) { + debug_msg("entering"); string file_bed = file_bfile + ".bed"; ifstream infile(file_bed.c_str(), ios::binary); if (!infile) { |