aboutsummaryrefslogtreecommitdiff
path: root/gn3/computations
diff options
context:
space:
mode:
Diffstat (limited to 'gn3/computations')
-rw-r--r--gn3/computations/biweight.R19
-rw-r--r--gn3/computations/biweight.py8
2 files changed, 22 insertions, 5 deletions
diff --git a/gn3/computations/biweight.R b/gn3/computations/biweight.R
index 05f7145..650c5fd 100644
--- a/gn3/computations/biweight.R
+++ b/gn3/computations/biweight.R
@@ -1,5 +1,5 @@
-library(WGCNA)
+# library(WGCNA)
myArgs <- commandArgs(trailingOnly = TRUE)
@@ -7,9 +7,22 @@ trait_vals <- as.numeric(unlist(strsplit(myArgs[1], split=" ")))
target_vals <- as.numeric(unlist(strsplit(myArgs[2], split=" ")))
BiweightMidCorrelation <- function(trait_val,target_val){
- results <- bicorAndPvalue(x,y)
+ results <- bicorAndPvalue(trait_val,target_val)
return (list(c(results$bicor)[1],c(results$p)[1]))
}
-cat(BiweightMidCorrelation(trait_vals,target_vals))
+
+
+
+
+# the idea is that you get the entire dataset in any format
+# and then do ther correlation
+
+ComputeAll <-function(trait_val,target_dataset) {
+ for target_val in target_dataset {
+ results = BiweightMidCorrelation(trait_val,target_val)
+ cat(BiweightMidCorrelation(trait_vals,target_vals))
+ }
+}
+
diff --git a/gn3/computations/biweight.py b/gn3/computations/biweight.py
index c17de8e..e598a5b 100644
--- a/gn3/computations/biweight.py
+++ b/gn3/computations/biweight.py
@@ -2,17 +2,21 @@
"""module contains script to call biweight mid\
correlation in R"""
-
import subprocess
+import os
+from pathlib import Path
from typing import List
+FILE_PATH = os.path.join(Path(__file__).parent.absolute(), "biweight.R")
+
def call_biweight_script(trait_vals: List,
target_vals: List,
- path_to_script: str = "./biweight_R",
+ path_to_script: str = FILE_PATH,
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]