diff options
| author | Alexander Kabui | 2021-06-16 13:41:28 +0300 |
|---|---|---|
| committer | Alexander Kabui | 2021-06-16 13:41:28 +0300 |
| commit | c0931bda121c1b39e1ae2adc7bfa49c2ca207a4d (patch) | |
| tree | 7dd793ee03bdd7d0f85130623d40f4606db99214 /gn3/computations/biweight.R | |
| parent | ceaea7930b0205993b9357bd1317fefe0c51b3b6 (diff) | |
| download | genenetwork3-c0931bda121c1b39e1ae2adc7bfa49c2ca207a4d.tar.gz | |
refactor script and add tests
Diffstat (limited to 'gn3/computations/biweight.R')
| -rw-r--r-- | gn3/computations/biweight.R | 14 |
1 files changed, 8 insertions, 6 deletions
diff --git a/gn3/computations/biweight.R b/gn3/computations/biweight.R index 8730536..e424360 100644 --- a/gn3/computations/biweight.R +++ b/gn3/computations/biweight.R @@ -1,10 +1,11 @@ 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=" "))) +arg_values <- commandArgs(trailingOnly = TRUE) +ParseArgs <- function(args){ + + trait_vals <- as.numeric(unlist(strsplit(args[1], split=" "))) + target_vals <- as.numeric(unlist(strsplit(args[2], split=" "))) return(list(trait_vals= c(trait_vals),target_vals = c(target_vals))) @@ -13,9 +14,10 @@ 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])) +parsed_values <- ParseArgs(arg_values) -cat(results)
\ No newline at end of file +cat((BiweightMidCorrelation(parsed_values[1],parsed_values[2])))
\ No newline at end of file |