From 57a6af689e85e5cbfe0a825d5b84ed9c451e6ad7 Mon Sep 17 00:00:00 2001
From: Alexander Kabui
Date: Mon, 14 Jun 2021 13:36:40 +0300
Subject: add biweight r script and tests

---
 gn3/computations/biweight.R | 15 +++++++++++++++
 1 file changed, 15 insertions(+)
 create mode 100644 gn3/computations/biweight.R

(limited to 'gn3/computations/biweight.R')

diff --git a/gn3/computations/biweight.R b/gn3/computations/biweight.R
new file mode 100644
index 0000000..d0d8de4
--- /dev/null
+++ b/gn3/computations/biweight.R
@@ -0,0 +1,15 @@
+
+library(WGCNA)
+
+
+myArgs <- commandArgs(trailingOnly = TRUE)
+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)
+    return (list(c(results$bicor)[1],c(results$p)[1]))
+}
+cat(BiweightMidCorrelation(trait_vals,target_vals))
+
+
-- 
cgit v1.2.3