diff options
Diffstat (limited to 'gn3/computations/biweight.R')
| -rw-r--r-- | gn3/computations/biweight.R | 15 |
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)) + + |