|
15
|
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(trailingOnly = T)
|
|
|
10 inputfile <- args[1]
|
|
|
11 Idx <- args[2]
|
|
|
12 outputfile <- args[3]
|
|
|
13
|
|
|
14 print(args)
|
|
|
15
|
|
|
16 ReadIn=read.table(inputfile,stringsAsFactors=F,header=T, sep="\t")
|
|
|
17
|
|
|
18 IndexIn=strsplit(Idx,split=",")[[1]]
|
|
|
19 Index=as.numeric(IndexIn)
|
|
|
20
|
|
|
21 Mat=data.matrix(ReadIn)
|
|
|
22
|
|
|
23 Out=Mat[Index,]
|
|
|
24
|
|
|
25
|
|
|
26 write.table(Out,file=outputfile,quote=F,col.names=T,row.names=T,sep = "\t")
|
|
|
27
|
