annotate extract_duplicates.r @ 1:a3c4e3e62e10 draft

Uploaded
author davidvanzessen
date Tue, 01 Sep 2015 07:49:47 -0400
parents 02cf2dd19564
children 1f1640608245
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
02cf2dd19564 Uploaded
davidvanzessen
parents:
diff changeset
1 args <- commandArgs(trailingOnly = TRUE)
02cf2dd19564 Uploaded
davidvanzessen
parents:
diff changeset
2
02cf2dd19564 Uploaded
davidvanzessen
parents:
diff changeset
3 input=args[1]
02cf2dd19564 Uploaded
davidvanzessen
parents:
diff changeset
4 column=as.numeric(args[2])
02cf2dd19564 Uploaded
davidvanzessen
parents:
diff changeset
5 header=(args[3] == "yes")
02cf2dd19564 Uploaded
davidvanzessen
parents:
diff changeset
6 out_file=args[4]
02cf2dd19564 Uploaded
davidvanzessen
parents:
diff changeset
7
1
a3c4e3e62e10 Uploaded
davidvanzessen
parents: 0
diff changeset
8 dat = read.table(input, header=header, sep="\t", fill=T, stringsAsFactors=F, quote="")
0
02cf2dd19564 Uploaded
davidvanzessen
parents:
diff changeset
9
02cf2dd19564 Uploaded
davidvanzessen
parents:
diff changeset
10 duplicates = dat[duplicated(dat[,column]),column]
02cf2dd19564 Uploaded
davidvanzessen
parents:
diff changeset
11
02cf2dd19564 Uploaded
davidvanzessen
parents:
diff changeset
12 dat = dat[dat[,column] %in% duplicates,]
02cf2dd19564 Uploaded
davidvanzessen
parents:
diff changeset
13
02cf2dd19564 Uploaded
davidvanzessen
parents:
diff changeset
14 dat = dat[order(dat[,column]),]
02cf2dd19564 Uploaded
davidvanzessen
parents:
diff changeset
15
02cf2dd19564 Uploaded
davidvanzessen
parents:
diff changeset
16 write.table(dat, out_file, sep="\t", row.names=F, col.names=header, quote=F)