Adding loco support for mdb

1 files changed, 3 insertions, 1 deletions

diff --git a/src/lmm.cpp b/src/lmm.cpp
index 375013c..7a80b8a 100644
--- a/src/lmm.cpp
+++ b/src/lmm.cpp
@@ -2155,8 +2155,10 @@ void LMM::mdb_analyze(std::function< SnpNameValues2(size_t) >& fetch_snp,
 
 void LMM::mdb_calc_gwa(const gsl_matrix *U, const gsl_vector *eval,
                           const gsl_matrix *UtW, const gsl_vector *Uty,
-                          const gsl_matrix *W, const gsl_vector *y) {
+                       const gsl_matrix *W, const gsl_vector *y, const string loco) {
   checkpoint("mdb-calc-gwa",file_geno);
+  bool is_loco = !loco.empty();
+  cout << loco << " CHR!!!!" << endl;
 
   // const auto num_snps = indicator_snp.size();
   // enforce_msg(num_snps > 0,"Zero SNPs to process - data corrupt?");