diff options
author | Pjotr Prins | 2018-02-10 08:22:41 +0000 |
---|---|---|
committer | Pjotr Prins | 2018-02-10 08:22:41 +0000 |
commit | b3d9c748008d401ca5b5fa5d60a6a1eff505ee3e (patch) | |
tree | 2e4ffdecc6c85c550c4170c9c440a0328aecb28d /wqflask/wqflask | |
parent | 97ccc478543a7bb08fe0cf178e320dd5811d8bfa (diff) | |
download | genenetwork2-b3d9c748008d401ca5b5fa5d60a6a1eff505ee3e.tar.gz |
Fixed finding PATHs - one return statement!
Diffstat (limited to 'wqflask/wqflask')
-rw-r--r-- | wqflask/wqflask/marker_regression/gemma_mapping.py | 7 |
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 |