about summary refs log tree commit diff
diff options
context:
space:
mode:
-rw-r--r--src/gemma.cpp4
-rw-r--r--src/gemma_io.cpp2
2 files changed, 3 insertions, 3 deletions
diff --git a/src/gemma.cpp b/src/gemma.cpp
index 7b955cd..1938e56 100644
--- a/src/gemma.cpp
+++ b/src/gemma.cpp
@@ -1129,7 +1129,7 @@ void GEMMA::Assign(int argc, char **argv, PARAM &cPar) {
         break;
       }
       if (argv[i + 1] == NULL || argv[i + 1][0] == '-') {
-        cPar.a_mode = M_KIN;
+        cPar.a_mode = M_KIN; // default
         continue;
       }
       ++i;
@@ -1913,7 +1913,7 @@ void GEMMA::BatchRun(PARAM &cPar) {
 
     if (cPar.a_mode == M_KIN) {
       cPar.WriteMatrix(G, "cXX");
-    } else {
+    } else { // M_KIN2
       cPar.WriteMatrix(G, "sXX");
     }
 
diff --git a/src/gemma_io.cpp b/src/gemma_io.cpp
index 3302bad..b4f8d39 100644
--- a/src/gemma_io.cpp
+++ b/src/gemma_io.cpp
@@ -1671,7 +1671,7 @@ bool PlinkKin(const string &file_bed, vector<int> &indicator_snp,
 
     gsl_vector_add_constant(geno, -1.0 * geno_mean);
 
-    if (k_mode == 2 && geno_var != 0) {
+    if (k_mode == 2 && geno_var != 0) { // k_mode is a_mode-20
       gsl_vector_scale(geno, 1.0 / sqrt(geno_var));
     }
     gsl_vector_view Xlarge_col = gsl_matrix_column(Xlarge, ns_test % msize);