annotate GetIg.R @ 12:1108020e2d84 draft

Uploaded
author ning
date Fri, 15 Mar 2013 18:23:14 -0400
parents f2a1b4f9ccd5
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
3
f2a1b4f9ccd5 Uploaded
ning
parents:
diff changeset
1 sink(file="/tmp/none")
f2a1b4f9ccd5 Uploaded
ning
parents:
diff changeset
2 sink("/dev/null")
f2a1b4f9ccd5 Uploaded
ning
parents:
diff changeset
3 options(warn=-1)
f2a1b4f9ccd5 Uploaded
ning
parents:
diff changeset
4 options(echo=F)
f2a1b4f9ccd5 Uploaded
ning
parents:
diff changeset
5
f2a1b4f9ccd5 Uploaded
ning
parents:
diff changeset
6 invisible("EBSeq")
f2a1b4f9ccd5 Uploaded
ning
parents:
diff changeset
7 suppressMessages(library("EBSeq"))
f2a1b4f9ccd5 Uploaded
ning
parents:
diff changeset
8
f2a1b4f9ccd5 Uploaded
ning
parents:
diff changeset
9 args <- commandArgs(trailingOnly = T)
f2a1b4f9ccd5 Uploaded
ning
parents:
diff changeset
10 inputfile <- args[1]
f2a1b4f9ccd5 Uploaded
ning
parents:
diff changeset
11 outputfile <- args[2]
f2a1b4f9ccd5 Uploaded
ning
parents:
diff changeset
12
f2a1b4f9ccd5 Uploaded
ning
parents:
diff changeset
13 print(args)
f2a1b4f9ccd5 Uploaded
ning
parents:
diff changeset
14
f2a1b4f9ccd5 Uploaded
ning
parents:
diff changeset
15 a1=read.csv(inputfile,stringsAsFactors=F,header=F, sep="\t")
f2a1b4f9ccd5 Uploaded
ning
parents:
diff changeset
16 Ng=GetNg(a1[[1]],a1[[2]])
f2a1b4f9ccd5 Uploaded
ning
parents:
diff changeset
17 Ig=Ng$IsoformNgTrun
f2a1b4f9ccd5 Uploaded
ning
parents:
diff changeset
18
f2a1b4f9ccd5 Uploaded
ning
parents:
diff changeset
19
f2a1b4f9ccd5 Uploaded
ning
parents:
diff changeset
20
f2a1b4f9ccd5 Uploaded
ning
parents:
diff changeset
21
f2a1b4f9ccd5 Uploaded
ning
parents:
diff changeset
22 write.table(Ig,file=outputfile,quote=F,col.names=F,row.names=F,sep = "\t")
f2a1b4f9ccd5 Uploaded
ning
parents:
diff changeset
23