Mercurial > repos > davidvanzessen > shm_csr
comparison merge.r @ 81:b6f9a640e098 draft
Uploaded
author | davidvanzessen |
---|---|
date | Fri, 19 Feb 2021 15:10:54 +0000 |
parents | c33d93683a09 |
children |
comparison
equal
deleted
inserted
replaced
80:a4617f1d1d89 | 81:b6f9a640e098 |
---|---|
1 args <- commandArgs(trailingOnly = TRUE) | |
2 | |
3 input.1 = args[1] | |
4 input.2 = args[2] | |
5 | |
6 fields.1 = args[3] | |
7 fields.2 = args[4] | |
8 | |
9 field.1 = args[5] | |
10 field.2 = args[6] | |
11 | |
12 output = args[7] | |
13 | |
14 dat1 = read.table(input.1, header=T, sep="\t", quote="", stringsAsFactors=F, fill=T, row.names=NULL) | |
15 if(fields.1 != "all"){ | |
16 fields.1 = unlist(strsplit(fields.1, ",")) | |
17 dat1 = dat1[,fields.1] | |
18 } | |
19 dat2 = read.table(input.2, header=T, sep="\t", quote="", stringsAsFactors=F, fill=T, row.names=NULL) | |
20 if(fields.2 != "all"){ | |
21 fields.2 = unlist(strsplit(fields.2, ",")) | |
22 dat2 = dat2[,fields.2] | |
23 } | |
24 | |
25 dat3 = merge(dat1, dat2, by.x=field.1, by.y=field.2) | |
26 | |
27 write.table(dat3, output, sep="\t",quote=F,row.names=F,col.names=T) |