From c708f435f14e5394aa4e6e8f556aec976d61341f Mon Sep 17 00:00:00 2001
From: Alexander Kabui
Date: Tue, 15 Jun 2021 08:08:38 +0300
Subject: fix script path issues

---
 gn3/computations/biweight.R  | 19 ++++++++++++++++---
 gn3/computations/biweight.py |  8 ++++++--
 2 files changed, 22 insertions(+), 5 deletions(-)

(limited to 'gn3')

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]
-- 
cgit v1.2.3