about summary refs log tree commit diff
path: root/src/gemma.cpp
diff options
context:
space:
mode:
authorPjotr Prins2017-10-06 07:54:01 +0000
committerPjotr Prins2017-10-06 07:54:01 +0000
commitdf161be507ac0ad1d67a6528ebc664acec89fc9c (patch)
tree1ac1b8d216f024d8a0eac40d4a67f9f41e9ea86a /src/gemma.cpp
parent45aed16d2a5a39bb02beec8fde16d86185cd413f (diff)
downloadpangemma-df161be507ac0ad1d67a6528ebc664acec89fc9c.tar.gz
Oops: fixed logic problem after removing Oxford data
Diffstat (limited to 'src/gemma.cpp')
-rw-r--r--src/gemma.cpp19
1 files changed, 12 insertions, 7 deletions
diff --git a/src/gemma.cpp b/src/gemma.cpp
index 5fbd86c..c890055 100644
--- a/src/gemma.cpp
+++ b/src/gemma.cpp
@@ -2736,20 +2736,25 @@ void GEMMA::BatchRun(PARAM &cPar) {
           gsl_vector_view UtY_col = gsl_matrix_column(UtY, 0);
 
           if (!cPar.file_bfile.empty()) {
+            // PLINK analysis
             if (cPar.file_gxe.empty()) {
               cLmm.AnalyzePlink(U, eval, UtW, &UtY_col.vector, W,
                                 &Y_col.vector);
-            } else {
+            }
+            else {
               cLmm.AnalyzePlinkGXE(U, eval, UtW, &UtY_col.vector, W,
                                    &Y_col.vector, env);
             }
           }
-          if (cPar.file_gxe.empty()) {
-            cLmm.AnalyzeBimbam(U, eval, UtW, &UtY_col.vector, W,
-                               &Y_col.vector, cPar.setGWASnps);
-          } else {
-            cLmm.AnalyzeBimbamGXE(U, eval, UtW, &UtY_col.vector, W,
-                                  &Y_col.vector, env);
+          else {
+            // BIMBAM analysis
+            if (cPar.file_gxe.empty()) {
+              cLmm.AnalyzeBimbam(U, eval, UtW, &UtY_col.vector, W,
+                                 &Y_col.vector, cPar.setGWASnps);
+            } else {
+              cLmm.AnalyzeBimbamGXE(U, eval, UtW, &UtY_col.vector, W,
+                                    &Y_col.vector, env);
+            }
           }
 
           cLmm.WriteFiles();