diff options
| author | Alexander_Kabui | 2024-09-10 18:01:20 +0300 |
|---|---|---|
| committer | Alexander_Kabui | 2024-09-10 18:01:20 +0300 |
| commit | 613e8e973febdef73cffd6ca85a96d66893ffff0 (patch) | |
| tree | e88aad343d51570e15e05998274b8f152b06160e /gn2/wqflask | |
| parent | 16207edb77eef75903a007df6b2bd84cb4d5b943 (diff) | |
| parent | d8e1e010a268616dbf66656f0af65848c6869f7f (diff) | |
| download | genenetwork2-613e8e973febdef73cffd6ca85a96d66893ffff0.tar.gz | |
Merge changes
Diffstat (limited to 'gn2/wqflask')
| -rw-r--r-- | gn2/wqflask/marker_regression/gemma_mapping.py | 5 |
1 files changed, 4 insertions, 1 deletions
diff --git a/gn2/wqflask/marker_regression/gemma_mapping.py b/gn2/wqflask/marker_regression/gemma_mapping.py index d8851486..aa059f93 100644 --- a/gn2/wqflask/marker_regression/gemma_mapping.py +++ b/gn2/wqflask/marker_regression/gemma_mapping.py @@ -3,6 +3,7 @@ import math import string import random import json +import logging from gn2.base import webqtlConfig from gn2.base.trait import create_trait @@ -73,6 +74,7 @@ def run_gemma(this_trait, this_dataset, samples, vals, covariates, use_loco, f"{pheno_filepath} -a " f"{snps_filepath} -gk > " f"{k_json_output_filepath}") + logging.debug("generate_k_command: %s", generate_k_command) os.system(generate_k_command) gemma_command = (f"{GEMMA_WRAPPER_COMMAND} --json --loco " @@ -102,7 +104,7 @@ def run_gemma(this_trait, this_dataset, samples, vals, covariates, use_loco, f"{flat_files('genotype/bimbam')}/" f"{genofile_name}_snps.txt -gk > " f"{TEMPDIR}/gn2/{k_output_filename}.json") - + logging.debug("generate_k_command: %s", generate_k_command) os.system(generate_k_command) gemma_command = (f"{GEMMA_WRAPPER_COMMAND} --json --input " @@ -121,6 +123,7 @@ def run_gemma(this_trait, this_dataset, samples, vals, covariates, use_loco, else: gemma_command += f" > {TEMPDIR}/gn2/{gwa_output_filename}.json" + logging.debug("gemma_command: %s", gemma_command) os.system(gemma_command) else: gwa_output_filename = output_files |