Mercurial > repos > davidvanzessen > mutation_analysis

view filter.r @ 3:a0b27058dcac draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author davidvanzessen
date Wed, 17 Sep 2014 07:25:17 -0400
parents 74d2bc479bee
children
line wrap: on
line source

args <- commandArgs(trailingOnly = TRUE)

input = args[1]
filterf = args[2]
output = args[3]


dat = read.table(input, header=T, sep="\t", fill=T, stringsAsFactors=F)
filterdat = read.table(filterf, header=T, sep="\t", fill=T, stringsAsFactors=F)

dat = dat[dat$Sequence.ID %in% filterdat$ID,]

write.table(x=dat, file=output, sep="\t",quote=F,row.names=F,col.names=T)