aboutsummaryrefslogtreecommitdiff
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);