annotate filter.r @ 3:a0b27058dcac draft

Uploaded
author davidvanzessen
date Wed, 17 Sep 2014 07:25:17 -0400
parents 74d2bc479bee
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
74d2bc479bee Uploaded
davidvanzessen
parents:
diff changeset
1 args <- commandArgs(trailingOnly = TRUE)
74d2bc479bee Uploaded
davidvanzessen
parents:
diff changeset
2
74d2bc479bee Uploaded
davidvanzessen
parents:
diff changeset
3 input = args[1]
74d2bc479bee Uploaded
davidvanzessen
parents:
diff changeset
4 filterf = args[2]
74d2bc479bee Uploaded
davidvanzessen
parents:
diff changeset
5 output = args[3]
74d2bc479bee Uploaded
davidvanzessen
parents:
diff changeset
6
74d2bc479bee Uploaded
davidvanzessen
parents:
diff changeset
7
74d2bc479bee Uploaded
davidvanzessen
parents:
diff changeset
8 dat = read.table(input, header=T, sep="\t", fill=T, stringsAsFactors=F)
74d2bc479bee Uploaded
davidvanzessen
parents:
diff changeset
9 filterdat = read.table(filterf, header=T, sep="\t", fill=T, stringsAsFactors=F)
74d2bc479bee Uploaded
davidvanzessen
parents:
diff changeset
10
74d2bc479bee Uploaded
davidvanzessen
parents:
diff changeset
11 dat = dat[dat$Sequence.ID %in% filterdat$ID,]
74d2bc479bee Uploaded
davidvanzessen
parents:
diff changeset
12
74d2bc479bee Uploaded
davidvanzessen
parents:
diff changeset
13 write.table(x=dat, file=output, sep="\t",quote=F,row.names=F,col.names=T)