如何用perl處理？

bioinfor

有一組數(shù)據(jù)，例如下：
ACACA        matches        ACACA
ACAS2L        previous symbol of        ACSS1
AD-020        alias of        TMEM167B
ADCY3      matches        ADCY3
ADCY3        previous symbol of        ADCY8
AIP        alias of        AURKAIP1
AIP        matches        AK2

在第一列如果有重復(fù)，則選擇行里含有matches的項(xiàng)，而刪除不含matches的項(xiàng)，如ADCY3有重復(fù)，則保留ADCY3      matches        ADCY3；AIP有重復(fù)，則保留AIP        matches        AK2
行，上例的最終處理結(jié)果應(yīng)該為：

ACACA        matches        ACACA
ACAS2L        previous symbol of        ACSS1
AD-020        alias of        TMEM167B
ADCY3      matches        ADCY3
AIP        matches        AK2


請教這一程序如何編寫？

附件為data。謝謝！

bernana

回復(fù) 1# bioinfor

1. 
   
2. #!/usr/bin/perl
   
3. use strict;
   
4. use warnings;
   
5. use Data::Dumper;
   
6. my %hash;
   
7. open(IN,"$ARGV[0] " )||die "can't open $ARGV[0]";
   
8. while (<IN>) {
   
9.         chomp;
   
10.         my @list=split (/\s+/,$_,2);
    
11.         if (!exists $hash{$list[0]}) {
    
12.                 $hash{$list[0]}=$_;
    
13.         }elsif($list[1]=~/matches/){
    
14.                 $hash{$list[0]}=$_;
    
15.         }
    
16. }
    
17. close IN;
    
18. 
    
19. print Dumper %hash;
    
20. 
    
21. __DATA__
    
22. ACACA        matches        ACACA
    
23. ACAS2L        previous symbol of        ACSS1
    
24. AD-020        alias of        TMEM167B
    
25. ADCY3      matches        ADCY3
    
26. ADCY3        previous symbol of        ADCY8
    
27. AIP        alias of        AURKAIP1
    
28. AIP        matches        AK2
    

復(fù)制代碼

bioinfor

運(yùn)行有錯(cuò)誤啊

yinyuemi

回復(fù) 1# bioinfor


    shell 行么？

awk 'NR==FNR{count[$1]++;next}(count[$1]==1 || /match/)'  urfile urfile

bernana

運(yùn)行有錯(cuò)誤啊
bioinfor 發(fā)表于 2011-01-30 11:54

    啥錯(cuò)誤？貼出來呀。我這里沒問題

zzy7186

try awk :

1. 
   
2. awk '{a[$1]++;if(a[$1] == 1){b[$1] = $0}else{if($2~/\<matches\>/)b[$1]=$0;}}END{for(i in b)print b[i];}' urlfile
   

復(fù)制代碼

zzy7186

perl try :

1. 
   
2. 
   
3. #!/usr/bin/perl
   
4. use warnings;
   
5. use strict;
   
6. my %result;
   
7. while (<>){
   
8.         my $F1 = (split /\s+/,$_)[0];
   
9.         #print $F1,"\n";
   
10.         if($result{$F1}){
    
11.                 $result{$F1} = $_ if ($_ =~ /\bmatches\b/);
    
12.         }else{
    
13.                 $result{$F1} = $_;
    
14.         }
    
15. }
    
16. foreach (values %result){
    
17.         print;
    
18. }
    

復(fù)制代碼