comparison aurora_wgcna_render.R @ 14:4ec9d51ed4e0 draft

Uploaded
author spficklin
date Fri, 22 Nov 2019 21:42:39 +0000
parents 083886b3f3db
children
comparison
equal deleted inserted replaced
13:2c8beaaa2374 14:4ec9d51ed4e0
9 }) 9 })
10 10
11 # We need to not crash galaxy with an UTF8 error on German LC settings. 11 # We need to not crash galaxy with an UTF8 error on German LC settings.
12 loc = Sys.setlocale("LC_MESSAGES", "en_US.UTF-8") 12 loc = Sys.setlocale("LC_MESSAGES", "en_US.UTF-8")
13 13
14 # suppress warning
15 options(warn = -1) 14 options(warn = -1)
16
17 options(stringsAsFactors=FALSE, useFancyQuotes=FALSE)
18
19 args = commandArgs(trailingOnly=TRUE) 15 args = commandArgs(trailingOnly=TRUE)
20 16
21 suppressPackageStartupMessages({ 17 suppressPackageStartupMessages({
18 library(rmarkdown)
22 library(getopt) 19 library(getopt)
23 library(tools)
24 library(WGCNA) 20 library(WGCNA)
21 library(caret)
22 library(ggplot2)
23 library(reshape2)
24 # Uncomment for HTML-based reports.
25 #library(DT)
26 #library(htmltools)
25 }) 27 })
26 28
27 # Load libraries
28 library(rmarkdown)
29 library(WGCNA)
30 library(DT)
31 library(htmltools)
32 library(ggplot2)
33 library(reshape2)
34 library(caret)
35 29
36 # ------------------------------------------------------------------ 30 # ------------------------------------------------------------------
37 # Handle arguments from command line 31 # Handle arguments from command line
38 # ------------------------------------------------------------------ 32 # ------------------------------------------------------------------
39 33
42 # Input Files 36 # Input Files
43 'trait_data', 't', 2, 'character', 37 'trait_data', 't', 2, 'character',
44 'expression_data', 'e', 1, 'character', 38 'expression_data', 'e', 1, 'character',
45 39
46 # Input Arguments 40 # Input Arguments
47 'missing_value', 'i', 1, 'character', 41 'missing_value1', 'i', 1, 'character',
42 'missing_value2', 'o', 2, 'character',
48 'sname_col', 'c', 2, 'integer', 43 'sname_col', 'c', 2, 'integer',
49 'min_cluster_size', 's', 1, 'integer', 44 'min_cluster_size', 's', 1, 'integer',
50 'height_cut', 'h', 2, 'double', 45 'height_cut', 'h', 2, 'double',
51 'power', 'p', 2, 'double', 46 'power', 'p', 2, 'double',
52 'block_size', 'b', 1, 'integer', 47 'block_size', 'b', 1, 'integer',
59 'network_edges_file', 'w', 1, 'character', 54 'network_edges_file', 'w', 1, 'character',
60 'gene_association_file', 'g', 2, 'character', 55 'gene_association_file', 'g', 2, 'character',
61 'module_association_file', 'm', 2, 'character', 56 'module_association_file', 'm', 2, 'character',
62 'module_association_report', 'q', 2, 'character', 57 'module_association_report', 'q', 2, 'character',
63 'network_construction_report', 'r', 1, 'character', 58 'network_construction_report', 'r', 1, 'character',
59 'updated_trait_matrix', 'z', 2, 'character',
64 'r_data', 'a', 1, 'character', 60 'r_data', 'a', 1, 'character',
65 'render_log_file', 'l', 1, 'character' 61 'render_log_file', 'l', 1, 'character',
62 'filtered_GEM', 'd', 1, 'character'
66 ), 63 ),
67 byrow=TRUE, ncol=4) 64 byrow=TRUE, ncol=4)
68 65
69 opt = getopt(spec) 66 opt = getopt(spec)
70 67