UNIX for Dummies Questions & Answers

Hi,

I would like to know how can I merge files based on their coordinates, but mantaining the score of each file in the output file like:

Note: 1st column is for chromosome, 2nd for start, 3rd for end of segment, 4th for score

file1:

Code:

1 200 300 20
1 400 500 30

file2:

Code:

1 200 350 30
1 400 600 10

file3:

Code:

1 200 300 40

Desired output file:
(Note: 1st column is for chromosome, 2nd for start, 3rd for end of segment, 4th-to-last column is for score of each file)

Code:

1 200 300 20 30 40
1 400 500 30 10 0

Many thanks!

try:

Code:

awk '
NR==FNR {a[$1" "$2]=$1" "$2" "$3;}
{if (!ff[FILENAME]) {f[fc++]=FILENAME;ff[FILENAME]=FILENAME}}
{b[$1" "$2]=$1" "$2;c[$1" "$2" "FILENAME]=$NF;}
END {
  for (i=0; i<fc; i++) for (j in b) a[j]=a[j]" "((c[j" "f[i]])?(c[j" "f[i]]):0);
  for (i in a) print a[i];
}' file[1-3]