急求解！！程序不会写

albertkingking

dme-miR-1-3p        1        TGGAATGTAAAGAAGTATGGAG        t0000002_64648        TGGAATGTAAAGAAGTATGGAG        0        +
dme-miR-1-3p        1        TGGAATGTAAAGAAGTATGGA        t0000036_7701        TGGAATGTAAAGAAGTATGGA        0        +
dme-miR-1-3p        1        TGGAATGTAAAGAAGTATGG        t0000154_3349        TGGAATGTAAAGAAGTATGG        0        +
dme-miR-1-3p        1        TGGAATGTAAAGAAGTATG                     t0009606_145        TGGAATGTAAAGAAGTATG        0        +
dme-miR-2a-3p        1        TATCACAGCCAGCTTTGATGA        t0040809_32        TATCACAGCCAGCTTTGATGA        0        +
dme-miR-2a-3p        1        TATCACAGCCAGCTTTGATGAG        t0077517_15        TATCACAGCCAGCTTTGATGAG        0        +
dme-miR-2a-3p        1        TATCACAGCCAGCTTTGATGAGC        t0115229_9        TATCACAGCCAGCTTTGATGAGC        0        +
dme-miR-2a-3p        3        TCACAGCCAGCTTTGATGAGC        t0040804_32        TCACAGCCAGCTTTGATGAGC        0        +
这是文件内容，想要得到的东西：如果第一列内容相同的话，将第三列中如t0000002_64648下划线后的数字取下进行排序，取数字做大的那一行，比如对dme-miR-1-3p来说保留的那一行是dme-miR-1-3p        1        TGGAATGTAAAGAAGTATGGAG        t0000002_64648        TGGAATGTAAAGAAGTATGGAG        0        +

注意：t0040809_32 t0040804_32，这两个下划线后面的数字是最大的，而且相等，如果这种情况的话两行都保留

albertkingking

求大神帮助啊 ~~~~~小弟万分感谢  

午夜凋零

albertkingking 发表于 2012-07-24 21:44
dme-miR-1-3p        1        TGGAATGTAAAGAAGTATGGAG        t0000002_64648        TGGAATGTAAAGAAGTATGGAG        0        +
dme-miR-1-3p        1        TGGA ...

试试这个。

1. #!/usr/bin/perl
   
2. use warnings;
   
3. use strict;
   
4. 
   
5. my %hash;
   
6. while ( <DATA> ) {
   
7.     chomp;
   
8.     my $head = (split /\s+/, $_)[0];
   
9.     push @{$hash{$head}}, $_;
   
10. }
    
11. 
    
12. foreach my $key ( keys %hash ) {
    
13.     my @array = sort { $b<=>$a } map { m/_(\d+)/; $1 } @{$hash{$key}};
    
14.     my $a = $array[0];
    
15.     foreach ( @{$hash{$key}} ) {
    
16.         print $_.$/ if /_$a/;
    
17.     }
    
18. }
    
19. 
    
20. __DATA__
    
21. dme-miR-1-3p        1        TGGAATGTAAAGAAGTATGGAG        t0000002_64648        TGGAATGTAAAGAAGTATGGAG        0        +
    
22. dme-miR-1-3p        1        TGGAATGTAAAGAAGTATGGA        t0000036_7701        TGGAATGTAAAGAAGTATGGA        0        +
    
23. dme-miR-1-3p        1        TGGAATGTAAAGAAGTATGG        t0000154_3349        TGGAATGTAAAGAAGTATGG        0        +
    
24. dme-miR-1-3p        1        TGGAATGTAAAGAAGTATG                     t0009606_145        TGGAATGTAAAGAAGTATG        0        +
    
25. dme-miR-2a-3p        1        TATCACAGCCAGCTTTGATGA        t0040809_32        TATCACAGCCAGCTTTGATGA        0        +
    
26. dme-miR-2a-3p        1        TATCACAGCCAGCTTTGATGAG        t0077517_15        TATCACAGCCAGCTTTGATGAG        0        +
    
27. dme-miR-2a-3p        1        TATCACAGCCAGCTTTGATGAGC        t0115229_9        TATCACAGCCAGCTTTGATGAGC        0        +
    
28. dme-miR-2a-3p        3        TCACAGCCAGCTTTGATGAGC        t0040804_32

复制代码

cdtits

1. do { @x = split; @y = split /_/, $x[3]; push @{$hash{$x[0]}}, [$y[1], @x]; } for <DATA>;
   
2. while (($k, $v) = each %hash) {
   
3.     @tmp = reverse sort { $a->[0] <=> $b->[0] } @{$v};
   
4.     $max = $tmp[0][0];
   
5.     do { print join("\t", @{$_}[1..$#$_]), "\n" if ($max == $_->[0]) } for (@tmp);
   
6. }
   

复制代码

yinyuemi

1. my %data;
   
2. my %hash;
   
3. 
   
4. map{split;push @{$data{$_[0]." ".substr($_[3],index($_[3],"_")+1)}},$_}<DATA>;
   
5. map{if(!$hash{$_->[0]}++){print @{$data{$_->[0]." ".$_->[1]}}}}sort{$a->[0] cmp $b->[0] || $b->[1] <=> $a->[1]}map{[split]}keys %data;
   
6. 
   
7. 
   
8. __DATA__
   
9. dme-miR-1-3p        1        TGGAATGTAAAGAAGTATGGAG        t0000002_64648        TGGAATGTAAAGAAGTATGGAG        0        +
   
10. dme-miR-1-3p        1        TGGAATGTAAAGAAGTATGGA        t0000036_7701        TGGAATGTAAAGAAGTATGGA        0        +
    
11. dme-miR-1-3p        1        TGGAATGTAAAGAAGTATGG        t0000154_3349        TGGAATGTAAAGAAGTATGG        0        +
    
12. dme-miR-1-3p        1        TGGAATGTAAAGAAGTATG                     t0009606_145        TGGAATGTAAAGAAGTATG        0        +
    
13. dme-miR-2a-3p        1        TATCACAGCCAGCTTTGATGA        t0040809_32        TATCACAGCCAGCTTTGATGA        0        +
    
14. dme-miR-2a-3p        1        TATCACAGCCAGCTTTGATGAG        t0077517_15        TATCACAGCCAGCTTTGATGAG        0      +
    
15. dme-miR-2a-3p        1        TATCACAGCCAGCTTTGATGAGC        t0115229_9        TATCACAGCCAGCTTTGATGAGC        0     +
    
16. dme-miR-2a-3p        3        TCACAGCCAGCTTTGATGAGC        t0040804_32        TCACAGCCAGCTTTGATGAGC        0        +
    

复制代码

albertkingking

回复 3# 午夜凋零
跪谢 大神~~~~~~


   

albertkingking

回复 4# cdtits
跪谢啦~~~  太牛了{:3_200:}


   

albertkingking

回复 5# yinyuemi
都太牛了 ~~~跪谢{:3_201:}


   

albertkingking

回复 5# yinyuemi
太厉害了 {:3_201:} 跪谢~~~~


   

kk861123

1. #!/usr/bin/perl
   
2. use strict;
   
3. use warnings;
   
4. use Data::Dumper;
   
5. 
   
6. my %tmax_of;
   
7. my %data_of;
   
8. while (<DATA>) {
   
9.    my @fields = (split);
   
10.    my $tmax_n = (split /_/,$fields[3])[1];
    
11.    if (exists $tmax_of{ $fields[0] }) {
    
12.            if ( $tmax_n > $tmax_of{ $fields[0] }) {
    
13.                    $data_of{ $fields[0] } = [$_];
    
14.                    $tmax_of{ $fields[0] } = $tmax_n;
    
15.            }
    
16.            elsif ($tmax_n == $tmax_of{ $fields[0] }) {
    
17.                    push @{ $data_of{ $fields[0] } },$_;
    
18.            }
    
19.    } else {
    
20.        $data_of{ $fields[0] } = [$_];
    
21.            $tmax_of{ $fields[0] } = $tmax_n;
    
22.    }
    
23. }
    
24. 
    
25. print Dumper \%tmax_of;
    
26. print Dumper \%data_of;
    
27. 
    
28. 
    
29. __DATA__
    
30. dme-miR-1-3p        1        TGGAATGTAAAGAAGTATGGAG        t0000002_64648        TGGAATGTAAAGAAGTATGGAG        0        +
    
31. dme-miR-1-3p        1        TGGAATGTAAAGAAGTATGGA        t0000036_7701        TGGAATGTAAAGAAGTATGGA        0        +
    
32. dme-miR-1-3p        1        TGGAATGTAAAGAAGTATGG        t0000154_3349        TGGAATGTAAAGAAGTATGG        0        +
    
33. dme-miR-1-3p        1        TGGAATGTAAAGAAGTATG                     t0009606_145        TGGAATGTAAAGAAGTATG        0        +
    
34. dme-miR-2a-3p        1        TATCACAGCCAGCTTTGATGA        t0040809_32        TATCACAGCCAGCTTTGATGA        0        +
    
35. dme-miR-2a-3p        1        TATCACAGCCAGCTTTGATGAG        t0077517_15        TATCACAGCCAGCTTTGATGAG        0        +
    
36. dme-miR-2a-3p        1        TATCACAGCCAGCTTTGATGAGC        t0115229_9        TATCACAGCCAGCTTTGATGAGC        0        +
    
37. dme-miR-2a-3p        3        TCACAGCCAGCTTTGATGAGC        t0040804_32        TCACAGCCAGCTTTGATGAGC        0        +
    

复制代码