diff options
| -rw-r--r-- | gn3/computations/biweight.R | 31 |
1 files changed, 12 insertions, 19 deletions
diff --git a/gn3/computations/biweight.R b/gn3/computations/biweight.R index 650c5fd..8730536 100644 --- a/gn3/computations/biweight.R +++ b/gn3/computations/biweight.R @@ -1,28 +1,21 @@ -# library(WGCNA) +library(WGCNA) +FetchArgs <- function(){ + myArgs <- commandArgs(trailingOnly = TRUE) + trait_vals <- as.numeric(unlist(strsplit(myArgs[1], split=" "))) + target_vals <- as.numeric(unlist(strsplit(myArgs[2], split=" "))) -myArgs <- commandArgs(trailingOnly = TRUE) -trait_vals <- as.numeric(unlist(strsplit(myArgs[1], split=" "))) -target_vals <- as.numeric(unlist(strsplit(myArgs[2], split=" "))) + return(list(trait_vals= c(trait_vals),target_vals = c(target_vals))) -BiweightMidCorrelation <- function(trait_val,target_val){ - results <- bicorAndPvalue(trait_val,target_val) - return (list(c(results$bicor)[1],c(results$p)[1])) } +BiweightMidCorrelation <- function(trait_val,target_val){ + results <- bicorAndPvalue(c(trait_val),c(target_val)) + return ((c(c(results$bicor)[1],c(results$p)[1]))) +} +results <- (BiweightMidCorrelation(FetchArgs()[1],FetchArgs()[2])) - - -# 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)) - } -} - +cat(results)
\ No newline at end of file |