diff options
author | Pjotr Prins | 2017-10-06 07:54:01 +0000 |
---|---|---|
committer | Pjotr Prins | 2017-10-06 07:54:01 +0000 |
commit | df161be507ac0ad1d67a6528ebc664acec89fc9c (patch) | |
tree | 1ac1b8d216f024d8a0eac40d4a67f9f41e9ea86a /src/gemma.cpp | |
parent | 45aed16d2a5a39bb02beec8fde16d86185cd413f (diff) | |
download | pangemma-df161be507ac0ad1d67a6528ebc664acec89fc9c.tar.gz |
Oops: fixed logic problem after removing Oxford data
Diffstat (limited to 'src/gemma.cpp')
-rw-r--r-- | src/gemma.cpp | 19 |
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(); |