about summary refs log tree commit diff
path: root/wqflask/wqflask
diff options
context:
space:
mode:
authorPjotr Prins2018-02-10 08:22:41 +0000
committerPjotr Prins2018-02-10 08:22:41 +0000
commitb3d9c748008d401ca5b5fa5d60a6a1eff505ee3e (patch)
tree2e4ffdecc6c85c550c4170c9c440a0328aecb28d /wqflask/wqflask
parent97ccc478543a7bb08fe0cf178e320dd5811d8bfa (diff)
downloadgenenetwork2-b3d9c748008d401ca5b5fa5d60a6a1eff505ee3e.tar.gz
Fixed finding PATHs - one return statement!
Diffstat (limited to 'wqflask/wqflask')
-rw-r--r--wqflask/wqflask/marker_regression/gemma_mapping.py7
1 files changed, 4 insertions, 3 deletions
diff --git a/wqflask/wqflask/marker_regression/gemma_mapping.py b/wqflask/wqflask/marker_regression/gemma_mapping.py
index a24e43d4..c004983e 100644
--- a/wqflask/wqflask/marker_regression/gemma_mapping.py
+++ b/wqflask/wqflask/marker_regression/gemma_mapping.py
@@ -3,7 +3,7 @@ import os, math, string, random, json
 from base import webqtlConfig
 from base.trait import GeneralTrait
 from base.data_set import create_dataset
-from utility.tools import flat_files, GEMMA_COMMAND, GEMMA_WRAPPER_COMMAND, TEMPDIR
+from utility.tools import flat_files, GEMMA_COMMAND, GEMMA_WRAPPER_COMMAND, TEMPDIR, assert_bin
 
 import utility.logger
 logger = utility.logger.getLogger(__name__ )
@@ -11,6 +11,7 @@ logger = utility.logger.getLogger(__name__ )
 def run_gemma(this_dataset, samples, vals, covariates, method, use_loco):
     """Generates p-values for each marker using GEMMA"""
 
+    assert_bin(GEMMA_COMMAND);
     if this_dataset.group.genofile != None:
         genofile_name = this_dataset.group.genofile[:-5]
     else:
@@ -27,7 +28,7 @@ def run_gemma(this_dataset, samples, vals, covariates, method, use_loco):
         if i < (len(this_chromosomes) - 1):
             chr_list_string += this_chromosomes[i+1].name + ","
         else:
-            chr_list_string += this_chromosomes[i+1].name  
+            chr_list_string += this_chromosomes[i+1].name
 
     if covariates != "":
         gen_covariates_file(this_dataset, covariates)
@@ -247,4 +248,4 @@ def parse_loco_output(this_dataset, gwa_output_filename):
                     included_markers.append(line.split("\t")[1])
                     p_values.append(float(line.split("\t")[10]))
 
-    return marker_obs
\ No newline at end of file
+    return marker_obs