提取A文件中内容在B文件中对应最大数字的行

chenjiali1

一个文件中内容是
A
chr1    990400  990475

另一个文件内容
B
chr1    990400  990425  1.2
chr1    990425  990450  2.5
chr1    990450  990475  2.5

我想在A文件中第二三列间的数目在B文件查询，如chr1    990400  990475 在B中对应的位置内的有
chr1    990400  990425  1.2
chr1    990425  990450  2.5
chr1    990450  990475  2.5
选取其中第三列最大的行，如
chr1    990425  990450  2.5
chr1    990450  990475  2.5
加到原来行的右边 （ 对于这种2个数值完全一样的2个都加过来）
chr1    990400  990475 chr1    990425  990450  2.5
chr1    990400  990475 chr1    990450  990475  2.5
请问各位高手如何实现！

rubyish

我想问下，B文件内容有没有
chr2 ...
chr3 ...

picbhan

学生物的玩编程都伤不起啊。楼主解决了没？

rubyish

哦！我想再问下，A文件内容有没有
chr2 ...
chr3 ...

chenjiali1

回复 4# rubyish


    有呢，谢谢啦！

rubyish

然后:

1. #!/usr/bin/perl
   
2. open my $A, 'A';
   
3. open my $B, 'B';
   
4. my %h;
   
5. while (<$A>) {
   
6.     chomp;
   
7.     my ( $k, $i, $j ) = split;
   
8.     $h{$k}{$_} = [ $_, $i, $j, 0, [''] ];
   
9. }
   
10. 
    
11. while (<$B>) {
    
12.     my ( $k, $n, $x, $t ) = split;
    
13.     for my $K ( keys $h{$k} ) {
    
14.         my $G = $h{$k}{$K};
    
15.         if ( $n >= $G->[1] && $x <= $G->[2] ) {
    
16.             push @{ $G->[4] }, $_ if $t == $G->[3];
    
17.             @{$G}[ 3, 4 ] = ( $t, [$_] ) if $t > $G->[3];
    
18.         }
    
19.     }
    
20. }
    
21. 
    
22. for my $v ( values %h ) {
    
23.     for my $V ( values %$v ) {
    
24.         print map $V->[0] . ( $_ ? " $_" : $/ ), @{ $V->[4] };
    
25.     }
    
26. }

复制代码