comparison filter.r @ 0:74d2bc479bee draft

Uploaded
author davidvanzessen
date Mon, 18 Aug 2014 04:04:37 -0400
parents
children
comparison
equal deleted inserted replaced
-1:000000000000 0:74d2bc479bee
1 args <- commandArgs(trailingOnly = TRUE)
2
3 input = args[1]
4 filterf = args[2]
5 output = args[3]
6
7
8 dat = read.table(input, header=T, sep="\t", fill=T, stringsAsFactors=F)
9 filterdat = read.table(filterf, header=T, sep="\t", fill=T, stringsAsFactors=F)
10
11 dat = dat[dat$Sequence.ID %in% filterdat$ID,]
12
13 write.table(x=dat, file=output, sep="\t",quote=F,row.names=F,col.names=T)