Mercurial > repos > ppericard > viscorvar
comparison viscorvar_networkvar.R @ 2:c8533e9298e5 draft
"planemo upload for repository https://gitlab.com/bilille/galaxy-viscorvar commit 8cb5630238352459037b3647eebfabb5554566f6-dirty"
author | ppericard |
---|---|
date | Fri, 23 Oct 2020 10:15:56 +0000 |
parents | |
children | 88c1fd2ac110 |
comparison
equal
deleted
inserted
replaced
1:e93350dc99f1 | 2:c8533e9298e5 |
---|---|
1 #!/usr/bin/env Rscript | |
2 | |
3 # Setup R error handling to go to stderr | |
4 options( show.error.messages=F, error = function () { cat( geterrmessage(), file=stderr() ); q( "no", 1, F ) } ) | |
5 | |
6 # we need that to not crash galaxy with an UTF8 error on German LC settings. | |
7 loc <- Sys.setlocale("LC_MESSAGES", "en_US.UTF-8") | |
8 | |
9 ## Get parameters ## | |
10 suppressPackageStartupMessages(require(argparse)) | |
11 | |
12 parser <- ArgumentParser(description='Run the networkVar function') | |
13 | |
14 parser$add_argument('--mat_similarity_rdata', dest='mat_similarity_rdata', required=TRUE, help="matSimilarity RData file") | |
15 parser$add_argument('--var_list_file', dest='var_list_file', required=TRUE, help="Variables list file") | |
16 # parser$add_argument('--interest_var_file', dest='interest_var_file', required=FALSE, help="Variables of interest file") | |
17 parser$add_argument('--response_var', dest='response_var', required=TRUE, help="Response variables") | |
18 parser$add_argument('--output_graph', dest='output_graph', required=TRUE, help="Output graphml") | |
19 | |
20 args <- parser$parse_args() | |
21 | |
22 ## Print parameters | |
23 print("matSimilarity RData file:") | |
24 print(args$mat_similarity_rdata) | |
25 print("Variables list file:") | |
26 print(args$var_list_file) | |
27 # print("Variables of interest:") | |
28 # print(args$interest_var_file) | |
29 print("Response variables:") | |
30 print(args$response_var) | |
31 print("Output graphml:") | |
32 print(args$output_graph) | |
33 | |
34 ## Loading libraries | |
35 # suppressPackageStartupMessages(require(mixOmics)) | |
36 suppressPackageStartupMessages(require(igraph)) | |
37 | |
38 # R script call | |
39 source_local <- function(fname) | |
40 { | |
41 argv <- commandArgs(trailingOnly = FALSE) | |
42 base_dir <- dirname(substring(argv[grep("--file=", argv)], 8)) | |
43 source(paste(base_dir, fname, sep="/")) | |
44 } | |
45 | |
46 ## Loading local functions | |
47 suppressPackageStartupMessages(require(visCorVar)) | |
48 | |
49 # Loading input Rdata file | |
50 # loads res_compute_mat_similarity object | |
51 load(args$mat_similarity_rdata) | |
52 | |
53 # | |
54 names_block_variables = as.character(read.table(args$var_list_file, header=FALSE)[[1]]) | |
55 names_response_variables = strsplit(args$response_var, ",")[[1]] | |
56 | |
57 # interest_var_vec = NULL | |
58 # if (args$interest_var_file != 'None') { | |
59 # interest_var_vec = as.character(read.table(args$interest_var_file, header=FALSE)[[1]]) | |
60 # } | |
61 | |
62 print("names_block_variables:") | |
63 print(names_block_variables) | |
64 # print("interest_var_vec:") | |
65 # print(interest_var_vec) | |
66 print("names_response_variables:") | |
67 print(names_response_variables) | |
68 | |
69 res_networkVar = networkVar(res_compute_mat_similarity = res_compute_mat_similarity, | |
70 names_block_variables = names_block_variables, | |
71 names_response_variables = names_response_variables, | |
72 cutoff = 0) | |
73 | |
74 print("networkVar() completed") | |
75 print("res_networkVar:") | |
76 print(res_networkVar) | |
77 | |
78 write.graph(res_networkVar$gR, file = args$output_graph, format = "graphml") |