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