about summary refs log tree commit diff
diff options
context:
space:
mode:
authorzsloan2022-06-01 20:43:00 +0000
committerzsloan2022-06-01 20:43:00 +0000
commitcba418b1953c36558f0bb05da98165d52cbfb52b (patch)
tree37b110e42a76005043f744aa1425c9485cb6caaa
parent508a24488190ee5ca22294fd4331d857efcf170b (diff)
downloadgenenetwork2-cba418b1953c36558f0bb05da98165d52cbfb52b.tar.gz
Fix the way R/qtl is called in the GN2 REST API
-rw-r--r--wqflask/wqflask/api/mapping.py18
1 files changed, 8 insertions, 10 deletions
diff --git a/wqflask/wqflask/api/mapping.py b/wqflask/wqflask/api/mapping.py
index f8b0d8bd..d1cc368e 100644
--- a/wqflask/wqflask/api/mapping.py
+++ b/wqflask/wqflask/api/mapping.py
@@ -52,17 +52,15 @@ def do_mapping_for_api(start_vars):
     elif mapping_params['mapping_method'] == "rqtl":
         header_row = ["name", "chr", "cM", "lod_score"]
         if mapping_params['num_perm'] > 0:
-            _sperm_output, _suggestive, _significant, result_markers = rqtl_mapping.run_rqtl_geno(vals, dataset, mapping_params['rqtl_method'], mapping_params['rqtl_model'],
-                                                                                                  mapping_params['perm_check'], mapping_params[
-                                                                                                      'num_perm'],
-                                                                                                  mapping_params['do_control'], mapping_params[
-                'control_marker'],
-                mapping_params['manhattan_plot'], mapping_params['pair_scan'])
+            _sperm_output, _suggestive, _significant, result_markers = rqtl_mapping.run_rqtl(this_trait.name, vals, samples, dataset, None, "Mb", mapping_params['rqtl_model'],
+                                                                                             mapping_params['rqtl_method'], mapping_params['num_perm'], None,
+                                                                                             mapping_params['do_control'], mapping_params['control_marker'],
+                                                                                             mapping_params['manhattan_plot'], None)
         else:
-            result_markers = rqtl_mapping.run_rqtl_geno(vals, dataset, mapping_params['rqtl_method'], mapping_params['rqtl_model'],
-                                                        mapping_params['perm_check'], mapping_params['num_perm'],
-                                                        mapping_params['do_control'], mapping_params['control_marker'],
-                                                        mapping_params['manhattan_plot'], mapping_params['pair_scan'])
+            result_markers = rqtl_mapping.run_rqtl(this_trait.name, vals, samples, dataset, None, "Mb", mapping_params['rqtl_model'],
+                                                   mapping_params['rqtl_method'], mapping_params['num_perm'], None,
+                                                   mapping_params['do_control'], mapping_params['control_marker'],
+                                                   mapping_params['manhattan_plot'], None)
 
     if mapping_params['limit_to']:
         result_markers = result_markers[:mapping_params['limit_to']]