From 254113f4e02d5b423af302e9facb0ed580749ef5 Mon Sep 17 00:00:00 2001 From: Pjotr Prins Date: Sun, 30 Nov 2025 08:32:33 +0100 Subject: K: reads geno mdb --- src/param.cpp | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) (limited to 'src/param.cpp') diff --git a/src/param.cpp b/src/param.cpp index 05e4a7b..3b74e53 100644 --- a/src/param.cpp +++ b/src/param.cpp @@ -1294,7 +1294,7 @@ void PARAM::ReadBIMBAMGenotypes(gsl_matrix *UtX, gsl_matrix *K, const bool calc_ } void PARAM::MdbCalcKin() { - error = !mdb_calc_kin(file_geno, is_loco(), indicator_snp, a_mode - 20); + error = !mdb_calc_kin(file_geno, is_loco(), a_mode - 20, indicator_snp); return; } @@ -1310,6 +1310,7 @@ void PARAM::CalcKin(gsl_matrix *matrix_kin) { error = true; } } else { + // indicator_snp is an array of bool error = !BimbamKin(file_geno, setKSnps, indicator_snp, a_mode - 20, d_pace, matrix_kin, ni_max == 0); } -- cgit 1.4.1