# HG changeset patch # User ning # Date 1402949864 14400 # Node ID f2327f520c4dd319ec1f1505870bfc9d137a7618 # Parent 97d36e663171141a7f96a993a20a3de51068c62b Deleted selected files diff -r 97d36e663171 -r f2327f520c4d Use/ChoosePatterns.R --- a/Use/ChoosePatterns.R Mon Jun 16 16:17:36 2014 -0400 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,30 +0,0 @@ -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 -