Mercurial > repos > yboursin > ebseq_eba
diff GetIg.R @ 0:ddd3dad04441 draft default tip
Uploaded
author | yboursin |
---|---|
date | Thu, 20 Oct 2016 08:59:24 -0400 |
parents | |
children |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/GetIg.R Thu Oct 20 08:59:24 2016 -0400 @@ -0,0 +1,23 @@ +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") +