Mercurial > repos > ning > ebseq
annotate GetIg.R @ 8:03aa1c6fe091 draft
Uploaded
| author | ning |
|---|---|
| date | Fri, 15 Mar 2013 18:21:47 -0400 |
| parents | f2a1b4f9ccd5 |
| children |
| rev | line source |
|---|---|
| 3 | 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 |