add biweight r script and tests

1 files changed, 15 insertions, 0 deletions

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))
+
+