Shell Programming and Scripting

Code:


Code:

rs0001      100e-34
rs0003      2.8e-01
rs008        1.9e-90

Code:

rs0001      1.98e-22
rs0004      3.77e-10
rs0003      2.8e-01
rs008        1.9e-90

rs0001      81.98e-22
rs0003      42.8e-01

                 file1              file2             file3
rs0001      100e-34         1.98e-22       81.98e-22
rs0004       NA                3.77e-10       NA
rs0003      2.8e-01          2.8e-01        42.8e-01
rs008        1.9e-90          1.9e-90         NA

awk '{ a[$1]=$2; s[$1]=$1; next } {
  s[$1] = s[$1] " | " $2 - a[$1]; a[$1]=$2
}
END{for(i in s) {print s[i]}}' file*.txt

$ awk '{A[$1] = A[$1] ? A[$1] OFS $2 : $2}END{for(i in A)print i, A[i]}' file*  | sort

$ awk 'BEGIN{printf "\t"}{A[$1] = A[$1] ?  A[$1] OFS $2 : $2}FNR==1{printf FILENAME OFS}END{printf RS; for(i in A)print i, A[i]}' OFS='\t' file* | sort

          file1      file2       file3    
rs0001    100e-34    1.98e-22    81.98e-22
rs0003    2.8e-01    2.8e-01     42.8e-01
rs0004    3.77e-10
rs008     1.9e-90    1.9e-90

awk '{A[$1]=A[$1] ? A[$1] OFS $2 : $2}END{for(i in A)print i, A[i] | "sort"}' file*
rs0001 100e-34 1.98e-22 81.98e-22
rs0003 2.8e-01 2.8e-01 42.8e-01
rs0004 3.77e-10
rs008 1.9e-90 1.9e-90

awk 'FNR==1{j=j OFS FILENAME} NR==FNR{a[$1];} {a[$1]=a[$1]?a[$1] OFS $2:$2} END{gsub(/^ /,X,j); print j;{for(i in a){print i OFS a[i]}}}' check_file_check_file121213*

check_file_check_file1212134 check_file_check_file1212135 check_file_check_file1212136
rs0001 100e-34 1.98e-22 81.98e-22
rs0003 2.8e-01 2.8e-01 42.8e-01
rs0004 3.77e-10
rs008 1.9e-90 1.9e-90

import os
import collections
d=collections.defaultdict(list)
cnt=0
for f in sorted(os.listdir('dir')): 
 with open('dir/'+f) as file:
  for line in file:
   t=line.split(" ")
   if t[0] not in d and cnt==1:
    d[t[0]]=['NA']*cnt    
   d[t[0]].append(t[1]) 
 cnt+=1

for i in sorted(d):
 print(i," ".join(map(lambda x: x.replace("\n",""), d[i])),end=" ")
 if len(d[i])<cnt:
  for i in range(cnt-len(d[i])):
   print('NA',end="")
 print("")