diff options
author | Alexander Kabui | 2021-06-15 08:08:38 +0300 |
---|---|---|
committer | Alexander Kabui | 2021-06-15 08:08:38 +0300 |
commit | c708f435f14e5394aa4e6e8f556aec976d61341f (patch) | |
tree | 5b1f0cf8c2dbc0b2d76d2df34a7bab271b6bc5c3 /gn3/computations/biweight.R | |
parent | f9dec15fb4a1fc24dfb963ed69f2a5cbbd4f31d3 (diff) | |
download | genenetwork3-c708f435f14e5394aa4e6e8f556aec976d61341f.tar.gz |
fix script path issues
Diffstat (limited to 'gn3/computations/biweight.R')
-rw-r--r-- | gn3/computations/biweight.R | 19 |
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)) + } +} + |