Mercurial > repos > yboursin > ebseq_eba
annotate GetIg.R @ 0:ddd3dad04441 draft default tip
Uploaded
| author | yboursin |
|---|---|
| date | Thu, 20 Oct 2016 08:59:24 -0400 |
| parents | |
| children |
| rev | line source |
|---|---|
| 0 | 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 |