diff options
| author | zsloan | 2022-08-20 19:20:20 +0000 |
|---|---|---|
| committer | zsloan | 2022-08-20 19:20:20 +0000 |
| commit | e15d5649b3a98aa39b0d8c095f4f9a44c99e4acd (patch) | |
| tree | 210647afd86555df495a5ef39c0ddd8faf9ba1b2 | |
| parent | c33ac9566645b08370a8a85f8a3df5d0c8c66f28 (diff) | |
| download | genenetwork2-e15d5649b3a98aa39b0d8c095f4f9a44c99e4acd.tar.gz | |
Divide GEMMA additive effect by 2
| -rw-r--r-- | wqflask/wqflask/marker_regression/gemma_mapping.py | 2 |
1 files changed, 1 insertions, 1 deletions
diff --git a/wqflask/wqflask/marker_regression/gemma_mapping.py b/wqflask/wqflask/marker_regression/gemma_mapping.py index ec12f9b7..6b525733 100644 --- a/wqflask/wqflask/marker_regression/gemma_mapping.py +++ b/wqflask/wqflask/marker_regression/gemma_mapping.py @@ -220,7 +220,7 @@ def parse_loco_output(this_dataset, gwa_output_filename, loco="True"): marker['chr'] = line.split("\t")[0] marker['Mb'] = float(line.split("\t")[2]) / 1000000 marker['p_value'] = float(line.split("\t")[10]) - marker['additive'] = -(float(line.split("\t")[7])) + marker['additive'] = -(float(line.split("\t")[7])/2) if math.isnan(marker['p_value']) or (marker['p_value'] <= 0): marker['lod_score'] = 0 else: |