Mercurial > repos > xuebing > sharplabtool

comparison Start_GenometriCorr.R @ 11:b7f1d9f8f3bc

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author xuebing
date Sat, 10 Mar 2012 07:59:27 -0500
parents
children
comparison
equal deleted inserted replaced
10:1558594a3c2e 11:b7f1d9f8f3bc
1 # Start_GenometriCorr.R
2
3 ###################################################
4 # #
5 # command-line interface to GenometriCorr #
6 # functions, for use with Galaxy. #
7 # #
8 ###################################################
9
10 capture.output <- function (result, pdffile, output_options)
11 {
12 if(output_options != "stats")
13 {
14 pdf(file=pdffile, width=10, height=19, paper="special")
15
16 if (output_options != "vis") #need to do a plot
17 {
18 mymat <- matrix(ncol=3, nrow=4)
19 mymat[1,1] <- 1
20 mymat[1,2] <- 2
21 mymat[1,3] <- 3
22 mymat[2,1] <- 4
23 mymat[2,2] <- 5
24 mymat[2,3] <- 6
25 mymat[3,1] <- 7
26 mymat[3,2] <- 8
27 mymat[3,3] <- 9
28 mymat[4,1] <- 10
29 mymat[4,2] <- 11
30 mymat[4,3] <- 12
31
32 layout(mymat, heights=c(0.2,0.2,0.2,0.2))
33 plot(result, pdffile, make.new=FALSE)
34 }
35 if (output_options != "plot") #need to do the bigger graphic
36 {
37 mymat <- matrix(ncol=2, nrow=8)
38 mymat[1,1] <- 2
39 mymat[1,2] <- 3
40 mymat[2,1] <- 4
41 mymat[2,2] <- 4
42 mymat[3,1] <- 1
43 mymat[3,2] <- 1
44 mymat[4,1] <- 5
45 mymat[4,2] <- 6
46 mymat[5,1] <- 7
47 mymat[5,2] <- 7
48 mymat[6,1] <- 8
49 mymat[6,2] <- 9
50 mymat[7,1] <- 10
51 mymat[7,2] <- 10
52 mymat[8,1] <- 11
53 mymat[8,2] <- 12
54 layoutresults <- 3
55
56 layout(mymat, heights=c(0.05,0.05,0.15,0.15,0.15,0.15,0.15,0.15))
57 visualize(result, pdffile, make.new=FALSE)
58 }
59 dev.off()
60 }
61
62 if (output_options == "stats")
63 {
64 show(result)
65 }
66 }
67
68
69
70 # Reads the command line arguments
71 args <- commandArgs(trailingOnly=T)
72
73 suppressPackageStartupMessages(library('GenometriCorr', warn.conflicts=F, verbose=F))
74 suppressPackageStartupMessages(library('graphics', warn.conflicts=F, verbose=F))
75 suppressPackageStartupMessages(library('gdata', warn.conflicts=F, verbose=F))
76 suppressPackageStartupMessages(library('gplots', warn.conflicts=F, verbose=F))
77 suppressPackageStartupMessages(library('gtools', warn.conflicts=F, verbose=F))
78 suppressPackageStartupMessages(library('caTools', warn.conflicts=F, verbose=F))
79 suppressPackageStartupMessages(library('grid', warn.conflicts=F, verbose=F))
80
81
82
83 # Variables
84 query_file <- ""
85 reference_file <- ""
86 config_file <- ""
87 output_options <- ""
88
89 # Parse the command line arguments
90
91 config_file <- args[1]
92 query_file <- as.character(args[2])
93 reference_file <- as.character(args[3])
94 output_options <- args[4]
95 pdffile <- args[5]
96
97 conf<-new("GenometriCorrConfig",config_file)
98
99 print('OK')
100
101 result<-suppressWarnings(suppressPackageStartupMessages(GenometriCorr:::run.config(conf,query=query_file,reference=reference_file)))
102 print('OK2')
103
104 hideoutput <- capture.output(result, pdffile=args[5], output_options)
105