Mercurial > repos > yboursin > ebseq_eba

view GetIg.R @ 0:ddd3dad04441 draft default tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author yboursin
date Thu, 20 Oct 2016 08:59:24 -0400
parents
children
line wrap: on
line source

sink(file="/tmp/none")
sink("/dev/null")
options(warn=-1)
options(echo=F) 

invisible("EBSeq")
suppressMessages(library("EBSeq"))

args <- commandArgs(trailingOnly = T)
inputfile <- args[1]
outputfile <- args[2]

print(args)

a1=read.csv(inputfile,stringsAsFactors=F,header=F, sep="\t")
Ng=GetNg(a1[[1]],a1[[2]])
Ig=Ng$IsoformNgTrun




write.table(Ig,file=outputfile,quote=F,col.names=F,row.names=F,sep = "\t")