aboutsummaryrefslogtreecommitdiff
path: root/gn3/computations/biweight.R
diff options
context:
space:
mode:
authorAlexander Kabui2021-06-15 08:08:38 +0300
committerAlexander Kabui2021-06-15 08:08:38 +0300
commitc708f435f14e5394aa4e6e8f556aec976d61341f (patch)
tree5b1f0cf8c2dbc0b2d76d2df34a7bab271b6bc5c3 /gn3/computations/biweight.R
parentf9dec15fb4a1fc24dfb963ed69f2a5cbbd4f31d3 (diff)
downloadgenenetwork3-c708f435f14e5394aa4e6e8f556aec976d61341f.tar.gz
fix script path issues
Diffstat (limited to 'gn3/computations/biweight.R')
-rw-r--r--gn3/computations/biweight.R19
1 files changed, 16 insertions, 3 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))
+ }
+}
+