aboutsummaryrefslogtreecommitdiff
path: root/src/lm.cpp
diff options
context:
space:
mode:
Diffstat (limited to 'src/lm.cpp')
-rw-r--r--src/lm.cpp4
1 files changed, 4 insertions, 0 deletions
diff --git a/src/lm.cpp b/src/lm.cpp
index f8fc43d..83558e4 100644
--- a/src/lm.cpp
+++ b/src/lm.cpp
@@ -290,6 +290,7 @@ void LmCalcP(const size_t test_mode, const double yPwy, const double xPwy,
}
void LM::AnalyzeGene(const gsl_matrix *W, const gsl_vector *x) {
+ debug_msg("entering");
ifstream infile(file_gene.c_str(), ifstream::in);
if (!infile) {
cout << "error reading gene expression file:" << file_gene << endl;
@@ -382,6 +383,7 @@ void LM::AnalyzeGene(const gsl_matrix *W, const gsl_vector *x) {
// WJA added
void LM::Analyzebgen(const gsl_matrix *W, const gsl_vector *y) {
+ debug_msg("entering");
string file_bgen = file_oxford + ".bgen";
ifstream infile(file_bgen.c_str(), ios::binary);
if (!infile) {
@@ -606,6 +608,7 @@ void LM::Analyzebgen(const gsl_matrix *W, const gsl_vector *y) {
}
void LM::AnalyzeBimbam(const gsl_matrix *W, const gsl_vector *y) {
+ debug_msg("entering");
igzstream infile(file_geno.c_str(), igzstream::in);
if (!infile) {
cout << "error reading genotype file:" << file_geno << endl;
@@ -720,6 +723,7 @@ void LM::AnalyzeBimbam(const gsl_matrix *W, const gsl_vector *y) {
}
void LM::AnalyzePlink(const gsl_matrix *W, const gsl_vector *y) {
+ debug_msg("entering");
string file_bed = file_bfile + ".bed";
ifstream infile(file_bed.c_str(), ios::binary);
if (!infile) {