数据的多条件过滤

gooderpan

有一个文件file（空格分隔）：
>NM_198493|m.1  NM_198493|g.1  type:5prime_partial  len:287  gc:universal  NM_198493:2-862(+)
>NM_198493|m.2  NM_198493|g.2  type:complete  len:153  gc:universal  NM_198493:1865-2323(+)
>NM_198493|m.3  NM_198493|g.3  type:complete  len:118  gc:universal  NM_198493:1954-2307(+)
>NM_173852|m.4  NM_173852|g.4  type:complete  len:163  gc:universal  NM_173852:27-515(+)
>NM_001135254|m.9  NM_001135254|g.9  type:complete  len:506 gc:universal  NM_001135254:599-2116(+)
>NM_001135254|m.10  NM_001135254|g.10  type:5prime_partial  len:142 gc:universal  NM_001135254:2-427(+)

如果NM编号相同的情况下选出len最大的，想要的结果如下：
>NM_198493|m.1  NM_198493|g.1  type:5prime_partial  len:287  gc:universal  NM_198493:2-862(+)
>NM_173852|m.4  NM_173852|g.4  type:complete  len:163  gc:universal  NM_173852:27-515(+)
>NM_001135254|m.9  NM_001135254|g.9  type:complete  len:506 gc:universal  NM_001135254:599-2116(+)

请问大神，该怎么弄呀？？？

sunzhiguolu

1. #!/usr/bin/perl
   
2. use strict;
   
3. use warnings;
   
4. 
   
5. my @aData = "";
   
6. while (<DATA>){
   
7.     my ($id, $leng) = /(?|>(\w+)|len:(\d+))/g;
   
8.     if ($aData[0] ne $id){
   
9.         print $aData[-1] if ($aData[0] ne "");
   
10.         @aData = ($id, $leng, $_);
    
11.         next;
    
12.     }
    
13.     @aData = ($id, $leng, $_) if ($aData[1] < $leng);
    
14. }
    
15. 
    
16. print $aData[-1];
    
17. 
    
18. __DATA__
    
19. >NM_198493|m.1  NM_198493|g.1  type:5prime_partial  len:287  gc:universal  NM_198493:2-862(+)
    
20. >NM_198493|m.2  NM_198493|g.2  type:complete  len:153  gc:universal  NM_198493:1865-2323(+)
    
21. >NM_198493|m.3  NM_198493|g.3  type:complete  len:118  gc:universal  NM_198493:1954-2307(+)
    
22. >NM_173852|m.4  NM_173852|g.4  type:complete  len:163  gc:universal  NM_173852:27-515(+)
    
23. >NM_001135254|m.9  NM_001135254|g.9  type:complete  len:506 gc:universal  NM_001135254:599-2116(+)
    
24. >NM_001135254|m.10  NM_001135254|g.10  type:5prime_partial  len:142 gc:universal  NM_001135254:2-427(+)

复制代码

华小飞_Perl

1. #!/usr/bin/perl
   
2. 
   
3. use warnings;
   
4. use strict;
   
5. use List::Util qw(max);
   
6. 
   
7. my %len_max;
   
8. while (<>) {
   
9.         chomp;
   
10.         if (/>(\w+).*?len:(\d+)/) {
    
11.                 $len_max{$1}{$2} = $_;
    
12.         }
    
13. }
    
14. 
    
15. for my $name ( keys %len_max ) {
    
16.         my $max = max keys %{ $len_max{$name} };
    
17.         print $len_max{$name}{$max}, "\n";
    
18. }

复制代码