|
40
|
1 sink(file="/tmp/none")
|
|
|
2 sink("/dev/null")
|
|
|
3 options(warn=-1)
|
|
|
4 options(echo=F)
|
|
|
5
|
|
|
6 invisible("EBSeq")
|
|
|
7 suppressMessages(library("EBSeq"))
|
|
|
8
|
|
|
9 args <- commandArgs()
|
|
|
10 inputfile <- args[6]
|
|
|
11 outputfile <- args[7]
|
|
|
12 #PairwisePlots <-args[6]
|
|
|
13
|
|
|
14 print(args)
|
|
|
15
|
|
|
16 Conds=strsplit(inputfile,split=",")[[1]]
|
|
|
17
|
|
|
18
|
|
|
19 Out=GetPatterns(Conds)
|
|
|
20
|
|
|
21
|
|
|
22 #write.table(Out,file=outputfile,quote=F,col.names=T,row.names=T,sep = "\t")
|
|
|
23 write.table(Out,file=outputfile,quote=F,col.names=NA,row.names=T,sep = "\t") ##cms - changed on 11/11/13
|
|
|
24
|
