將兩個文件連接在一起（AWK）

Question

我有一個問題。 我想將兩個文件合並在一起。 哪里：

文件1：

 chr**1  10000**   rs200132 A  C  100.000
 chr**2  20000**   rs5000   C  G   80.000

檔案2：

 rs200132  **1:10000**  A   800   200  Nmf 
 rs210111  **1:10000**  G   200   800  VFC 
 rs310000  **1:10000**  C   100   500  tff
 rs50001   **2:20000**  T   500   100  jpp
 rs60000   **2:20000**  A   1000   10  jkl

輸出：

 chr**1  10000**  rs200132  A  A  C   800   200  Nmf
 chr**1  10000**  rs210111  G  A  C   200   800  VFC
 chr**1  10000**  rs310000  C  A  C   100   500  tff
 chr**2  20000**  rs50001   T  C  G   500   100  jpp
 chr**2  20000**  rs60000   A  C  G   1000   10  jkl

然后，從第一個文件開始，在“ chr”之后是標記數字，在第二列中是數字。 相同的標記在第二個文件中，但在第二列中如1：10000。 我想加入這兩個文件，但是對於第一個文件，將有更多的連接行（例如：對於第一個文件的第一行將是第二個文件的三行。）謝謝

Answer 1

perl -lane'
  BEGIN{ $x=pop; %h = map{ $_->[1] => $_ } map [split], <>; @ARGV=$x }
  $F[1] =~ s/.+?://;
  $t = $h{$F[1]};
  print join " ", @$t[0,1], @F[0,2], @$t[3,4], @F[3..5];
' file1 file2

輸出

chr**1 10000** rs200132 A A C 800 200 Nmf
chr**1 10000** rs210111 G A C 200 800 VFC
chr**1 10000** rs310000 C A C 100 500 tff
chr**2 20000** rs50001 T C G 500 100 jpp
chr**2 20000** rs60000 A C G 1000 10 jkl

Answer 2

您可以使用此awk ，

awk 'NR==FNR{a[$2]=$1;b[$2]=$4" "$5;next} {sub(/.*:/,"",$2); $3=$1" "$3" "b[$2]; $1=a[$2];}1' file1 file2

測試：

sat:~# awk 'NR==FNR{a[$2]=$1;b[$2]=$4" "$5;next} {sub(/.*:/,"",$2); $3=$1" "$3" "b[$2]; $1=a[$2];}1' file1 file2
chr**1 10000** rs200132 A A C 800 200 Nmf
chr**1 10000** rs210111 G A C 200 800 VFC
chr**1 10000** rs310000 C A C 100 500 tff
chr**2 20000** rs50001 T C G 500 100 jpp
chr**2 20000** rs60000 A C G 1000 10 jkl

Answer 3

這是使用Perl的另一種方式：

perl -lane '
    if (@ARGV) {
        ($x = $F[0]) =~ s/[^\d]*//;
        $h{$x}{$F[1]} = [ @F[0,1,3,4] ]
    }
    else {
        @t = split(":", $F[1]);
        $r = $h{$t[0]}{$t[1]};
        print join(" ", @$r[0,1], @F[0,2], @$r[2,3], @F[3..5])
    }
' file1 file2 | column -t

結果：

chr1  10000  rs200132  A  A  C  800   200  Nmf
chr1  10000  rs210111  G  A  C  200   800  VFC
chr1  10000  rs310000  C  A  C  100   500  tff
chr2  20000  rs50001   T  C  G  500   100  jpp
chr2  20000  rs60000   A  C  G  1000  10   jkl