# HG changeset patch
# User ning
# Date 1363386025 14400
# Node ID f2a1b4f9ccd5b1c9444d3bce3c7a6a79db4c1a78
# Parent  186aff8bcca3ab207ea297a224bc1d82e27a27c0
Uploaded

diff -r 186aff8bcca3 -r f2a1b4f9ccd5 GetIg.R
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/GetIg.R	Fri Mar 15 18:20:25 2013 -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")
+