# HG changeset patch
# User ning
# Date 1402949803 14400
# Node ID fdbabf28d095e9942e799389413cd5b83c1e8a08
# Parent  67e9a98ad324d3f035c6e60de43ffebc03126164
Uploaded

diff -r 67e9a98ad324 -r fdbabf28d095 Use/ChoosePatterns.R
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/Use/ChoosePatterns.R	Mon Jun 16 16:16:43 2014 -0400
@@ -0,0 +1,30 @@
+sink(file="/tmp/none")
+sink("/dev/null")
+options(warn=-1)
+options(echo=F) 
+
+#invisible("EBSeq")
+#suppressMessages(library("EBSeq"))  ##cms - changed on 11/11/13. Don't need EBSeq package.
+
+args <- commandArgs(trailingOnly = T)
+inputfile <- args[1]
+Idx <- args[2]
+outputfile <- args[3]
+
+print(args)
+
+ReadIn=read.table(inputfile,stringsAsFactors=F,header=T, sep="\t")
+
+IndexIn=strsplit(Idx,split=",")[[1]]
+Index=as.numeric(IndexIn)
+
+#Mat=data.matrix(ReadIn)
+Mat=data.matrix(ReadIn[,-1])  ##cms - changed on 11/11/13. Otherwise pattern IDs are replaced with NA.
+rownames(Mat)=ReadIn[,1]  ##cms - also added on 11/11/13 to add pattern IDs back in.
+
+Out=Mat[Index,]
+
+
+#write.table(Out,file=outputfile,quote=F,col.names=T,row.names=T,sep = "\t")
+write.table(Out,file=outputfile,quote=F,col.names=NA,row.names=T,sep = "\t")  ##cms - changed on 11/11/13
+