Mercurial > repos > davidvanzessen > extract_duplicates
view extract_duplicates.r @ 0:02cf2dd19564 draft
Uploaded
author | davidvanzessen |
---|---|
date | Fri, 21 Aug 2015 10:49:14 -0400 |
parents | |
children | a3c4e3e62e10 |
line wrap: on
line source
args <- commandArgs(trailingOnly = TRUE) input=args[1] column=as.numeric(args[2]) header=(args[3] == "yes") out_file=args[4] dat = read.table(input, header=header, sep="\t", fill=T, stringsAsFactors=F) duplicates = dat[duplicated(dat[,column]),column] dat = dat[dat[,column] %in% duplicates,] dat = dat[order(dat[,column]),] write.table(dat, out_file, sep="\t", row.names=F, col.names=header, quote=F)