comparison GetIg.R @ 3:f2a1b4f9ccd5 draft

Uploaded
author ning
date Fri, 15 Mar 2013 18:20:25 -0400
parents
children
comparison
equal deleted inserted replaced
2:186aff8bcca3 3:f2a1b4f9ccd5
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 outputfile <- args[2]
12
13 print(args)
14
15 a1=read.csv(inputfile,stringsAsFactors=F,header=F, sep="\t")
16 Ng=GetNg(a1[[1]],a1[[2]])
17 Ig=Ng$IsoformNgTrun
18
19
20
21
22 write.table(Ig,file=outputfile,quote=F,col.names=F,row.names=F,sep = "\t")
23