aboutsummaryrefslogtreecommitdiff
path: root/gn3/computations/ctl.py
blob: 1a71056b72107c5e8e139a470f2fd65f8d6c93ae (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
"""module contains code to process ctl analysis data"""
import json
from gn3.commands import run_cmd

from gn3.computations.wgcna import dump_wgcna_data
from gn3.computations.wgcna import compose_wgcna_cmd
from gn3.computations.wgcna import process_image

from gn3.settings import TMPDIR


def call_ctl_script(data):
    """function to call ctl script"""
    data["imgDir"] = TMPDIR
    temp_file_name = dump_wgcna_data(data)
    cmd = compose_wgcna_cmd("ctl_analysis.R", temp_file_name)

    try:
        cmd_results = run_cmd(cmd)
        with open(temp_file_name, "r") as outputfile:
            if cmd_results["code"] != 0:
                return cmd_results

            output_file_data = json.load(outputfile)

            output_file_data["image_data"] = process_image(
                output_file_data["image_loc"]).decode("ascii")

            output_file_data["ctl_plots"] = [process_image(ctl_plot).decode("ascii") for
                                             ctl_plot in output_file_data["ctl_plots"]]

            return output_file_data

    except Exception as e:
        return str(e)