文件匹配相同ID根据比较结果输出ID和匹配项——求大神赐教...

little_joe

感谢感谢！ 回复 6# jason680


   

little_joe

恩恩，数据的问题， 感谢感谢！回复 7# RE_HASH


   

sunzhiguolu

1. #!/usr/bin/perl
   
2. use strict;
   
3. use warnings;
   
4. 
   
5. my (%hStat, $sID);
   
6. while (<>){
   
7.     my @aCols = split;
   
8.     local $_ = $aCols[0];
   
9.     if (m/\A>/ and not defined $aCols[1]){
   
10.         $sID = $_;
    
11.         next;
    
12.     }
    
13.     if (!defined $aCols[1]){
    
14.         push (@{$hStat{$sID}}, $_);
    
15.         next;
    
16.     }
    
17.     next if (!/\A>/);
    
18.     if (exists $hStat{$_}){
    
19.         my $sPrint = 1;
    
20.         foreach my $sV (grep {$_ >= $aCols[1] and $_ <= $aCols[-1]} @{$hStat{$_}}){
    
21.             print "$sID\n" if ($sPrint++ == 1);
    
22.             print "$sV\n";
    
23.         }
    
24.     }
    
25. }
    

复制代码

little_joe

谢谢您的帮助，我试了一下发现输出结果不太匹配，例如文件输出为 但是事实上文件搜索并没有 匹配到的数字，请问这个是怎么回事呢，还望大神能够给解答一下，感谢感谢！文件如下，期待大神回信》》》 文件.zip (17.31 KB, 下载次数: 13)

2016-05-24 11:26 上传

点击文件名下载附件

回复 13# sunzhiguolu


   

104359176

现在怎么了，大神，跪拜这些词泛滥，难道为了获取帮助，连节操都失去了吗？大家都是从获取别人无私的帮助中过来的，不必如此。

sunzhiguolu

回复 14# little_joe
再试下,
perl ba.pl b.txt a.txt

1. #!/usr/bin/perl
   
2. use strict;
   
3. use warnings;
   
4. 
   
5. my ($sID, %hStat, $sPrint);
   
6. while (<>){
   
7.     my @aCols = split;
   
8.     if (@ARGV){
   
9.         next if (/\A\s+/);
   
10.         next if (exists $hStat{$aCols[0]});
    
11.         $hStat{$aCols[0]} = {'Min' => $aCols[1], 'Max' => $aCols[-1]};
    
12.         next;
    
13.     }
    
14. 
    
15.     if (/\A(>\S+)/){
    
16.         if (!exists $hStat{$1}){
    
17.             $sPrint = 0;
    
18.             next;
    
19.         }
    
20.         ($sPrint, $sID) = (1, $1);
    
21.         next;
    
22.     }
    
23.     if ($sPrint and $aCols[0] >= $hStat{$sID}{'Min'} && $aCols[0] <= $hStat{$sID}{'Max'}){
    
24.         print $sID, "\n" if ($sPrint++ == 1);
    
25.         print;
    
26.     }
    
27. }
    

复制代码

>NP_414681.1
683
690
691
692
694
>NP_414854.1
395
407

   

sunzhiguolu

回复 14# little_joe
13 楼的帖子, 代码的确有问题, 感谢提醒, 谢谢... (把有问题的部分已经修改了)
perl ab.pl a.txt b.txt

1. #!/usr/bin/perl
   
2. use strict;
   
3. use warnings;
   
4. 
   
5. my ($sID, %hStat);
   
6. while (<>){
   
7.     my @aCols = split;
   
8.     if (/\A>/ and !defined ($aCols[1])){
   
9.         $sID = $aCols[0];
   
10.         next;
    
11.     }
    
12.     if (!defined $aCols[1]){
    
13.         push (@{$hStat{$sID}}, $_);
    
14.         next;
    
15.     }
    
16. 
    
17.     next if (!/\A>/);
    
18.     my $sPrint = 0;
    
19.     if (exists $hStat{$aCols[0]}){
    
20.         foreach (grep {$_ >= $aCols[1] and $_ <= $aCols[-1]} @{$hStat{$aCols[0]}}){
    
21.             print $aCols[0], "\n" if (!$sPrint++);
    
22.             print;
    
23.         }
    
24.     }
    
25. }
    

复制代码

>NP_414681.1
683
690
691
692
694
>NP_414854.1
395
407

little_joe

这次结果对的，谢谢O(∩_∩)O谢谢回复 17# sunzhiguolu


   

little_joe

谢谢提醒，再次感谢各位的帮助和提醒回复 15# 104359176


   

sunzhiguolu

回复 6# jason680
通过这个问题, 终于找到问题所在了, 非常感谢大神指点...