about summary refs log tree commit diff
path: root/src/param.cpp
diff options
context:
space:
mode:
Diffstat (limited to 'src/param.cpp')
-rw-r--r--src/param.cpp17
1 files changed, 7 insertions, 10 deletions
diff --git a/src/param.cpp b/src/param.cpp
index 017d588..a720801 100644
--- a/src/param.cpp
+++ b/src/param.cpp
@@ -1293,28 +1293,25 @@ void PARAM::ReadBIMBAMGenotypes(gsl_matrix *UtX, gsl_matrix *K, const bool calc_
   return;
 }
 
+void PARAM::MdbCalcKin() {
+  error = !mdb_calc_kin(file_geno, setKSnps, indicator_snp, a_mode - 20);
+  return;
+}
+
 void PARAM::CalcKin(gsl_matrix *matrix_kin) {
   checkpoint_nofile("PARAM::CalcKin");
   string file_str;
-
   gsl_matrix_set_zero(matrix_kin);
 
   if (!file_bfile.empty()) {
     file_str = file_bfile + ".bed";
-    // enforce_msg(loco.empty(), "FIXME: LOCO nyi");
     if (PlinkKin(file_str, indicator_snp, a_mode - 20, d_pace, matrix_kin) ==
         false) {
       error = true;
     }
   } else {
-    file_str = file_geno;
-    if (is_mdb)
-      error = !BimbamKin(file_str, setKSnps, indicator_snp, a_mode - 20, d_pace,
-                      matrix_kin, ni_max == 0);
-    else
-      error = !BimbamKin(file_str, setKSnps, indicator_snp, a_mode - 20, d_pace,
-                         matrix_kin, ni_max == 0);
-
+    error = !BimbamKin(file_geno, setKSnps, indicator_snp, a_mode - 20, d_pace,
+                       matrix_kin, ni_max == 0);
   }
 
   return;