2

我有两个文件，文件 1 看起来像：

rs12126768  T   0.0252778071842686
rs908742    A   -0.0175531592475892
rs3107147   A   0.0253753063312283
rs3128309   A   0.0284902703996274
rs16824948  T   0.0696193418800571
rs2643888   A   0.0452601314048646
rs3795272   T   -0.0244976316004002
rs2840531   T   -0.0438473997802763
rs10910060  A   0.0364283566135514
rs1129332   T   -0.0274739770007551

文件 2 看起来像：

1   rs908742    0   2033256
1   rs3795272   0   2284195
1   rs2840531   0   2286408
1   rs1129332   0   2336210
1   rs7512269   0   2418625
1   rs2247308   0   2445406
1   rs1980789   0   2455863
1   rs3122922   0   2456203
1   rs4648453   0   3025712
1   rs6424069   0   3147030

如果文件 1 中的 $1 出现在文件 2 的 $2 中，我希望它打印该行。所以期望的输出是：

rs908742    A   -0.0175531592475892
rs3795272   T   -0.0244976316004002
rs2840531   T   -0.0438473997802763
rs1129332   T   -0.0274739770007551

预先感谢您的任何帮助。

4

2 回答 2

7

试试这个单行：

awk 'NR==FNR{a[$2];next}$1 in a' file2 file1
于 2013-07-03T13:25:50.213 回答
0

我建议这样做的算法是简单的散列。通读 file2，并记住 file2 的 $2 中的所有内容。然后，通读 file1，如果 file1 中的 $1 在内存中，则打印该行。

这是一个可以解决问题的 perl 脚本：

use strict;
open FILE2, "file2.txt";
my %remember;
while(<FILE2>) {
    chomp $_;
    my @fields = split;
    $remember{$fields[1]} = 1;
}
open FILE1, "file1.txt";
while(<FILE1>) {
    chomp $_;
    my @fields = split;
    if(defined $remember{$fields[0]}) {
            print "$_\n";
    }
}
于 2013-07-03T13:27:41.220 回答