diff options
| author | zsloan | 2020-03-19 11:46:49 -0500 |
|---|---|---|
| committer | zsloan | 2020-03-19 11:46:49 -0500 |
| commit | f97da31add711dcc0e75797b206821ca2d05c105 (patch) | |
| tree | 28514200563f2fcd10c67cf16ab057b08a630977 | |
| parent | f0f97202cf77a4459dd7ca525e94e34a56ddb147 (diff) | |
| download | genenetwork2-f97da31add711dcc0e75797b206821ca2d05c105.tar.gz | |
Made a minor change just to help qtlreaper deal with genotype files with positions in bases instead of megabases
| -rw-r--r-- | wqflask/wqflask/marker_regression/qtlreaper_mapping.py | 5 |
1 files changed, 4 insertions, 1 deletions
diff --git a/wqflask/wqflask/marker_regression/qtlreaper_mapping.py b/wqflask/wqflask/marker_regression/qtlreaper_mapping.py index 41764c1c..0c560582 100644 --- a/wqflask/wqflask/marker_regression/qtlreaper_mapping.py +++ b/wqflask/wqflask/marker_regression/qtlreaper_mapping.py @@ -97,7 +97,10 @@ def parse_reaper_output(gwa_filename, permu_filename, bootstrap_filename): except: marker['chr'] = line.split("\t")[2] marker['cM'] = float(line.split("\t")[3]) - marker['Mb'] = float(line.split("\t")[4]) + if float(line.split("\t")[4]) > 1000: + marker['Mb'] = float(line.split("\t")[4])/1000000 + else: + marker['Mb'] = float(line.split("\t")[4]) if float(line.split("\t")[7]) != 1: marker['p_value'] = float(line.split("\t")[7]) marker['lrs_value'] = float(line.split("\t")[5]) |