aboutsummaryrefslogtreecommitdiff
path: root/gn3/api/rqtl.py
blob: 82cf34f3d59d072edc7509e45f7bf8d598f68710 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
import os

from flask import Blueprint
from flask import current_app
from flask import jsonify
from flask import request

from gn3.computations.rqtl import generate_rqtl_cmd
from gn3.computations.gemma import do_paths_exist

rqtl = Blueprint("rqtl", __name__)

@rqtl.route("/compute", methods=["POST"])
def compute():
    working_dir = os.path.join(current_app.config.get("TMPDIR"))

    genofile = request.form['geno_file']
    phenofile = request.form['pheno_file']

    if not do_paths_exist([genofile, phenofile]):
        raise FileNotFoundError

    kwarg_list = ["addcovar", "model", "method", "interval", "nperm", "scale", "control_marker"]

    rqtl_kwargs = {"geno": genofile, "pheno": phenofile}
    for kwarg in kwarg_list:
        if kwarg in request.form:
            rqtl_kwargs[kwarg] = request.form[kwarg]

    results = generate_rqtl_cmd(
        rqtl_wrapper_cmd = current_app.config.get("RQTL_WRAPPER_CMD"),
        output_dir = current_app.config.get('TMPDIR'),
        rqtl_wrapper_kwargs = rqtl_kwargs
    )

    return jsonify(results)