aboutsummaryrefslogtreecommitdiff
path: root/gn3/computations/biweight.py
blob: c17de8e79e09236e739affb565094d53d245c6b0 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
"""module contains script to call biweight mid\
correlation in R"""

import subprocess
from typing import List


def call_biweight_script(trait_vals: List,
                         target_vals: List,
                         path_to_script: str = "./biweight_R",
                         command: str = "Rscript"
                         ):
    '''biweight function'''
    args_1 = ' '.join(str(trait_val) for trait_val in trait_vals)
    args_2 = ' '.join(str(target_val) for target_val in target_vals)
    cmd = [command, path_to_script] + [args_1] + [args_2]

    results = subprocess.check_output(cmd, universal_newlines=True)

    return tuple([float(y) for y in results.split()])