2016-07-22 41 views
0

我正在处理包含制表符分隔序列的输入文件。序列组通过换行符分隔。该文件看起来像:如何计算散列中存在的键的数量?

TAGC  TAGC  TAGC  HELP 
    TAGC  TAGC  TAGC 
    TAGC  HELP 
    TAGC 

这里是我的代码有:

use strict; 
    use warnings; 

    open(INFILE, "<", "/path/to/infile.txt") or die $!; 

    my %hash = ( 
      TAGC => 'THIS_EXISTS', 
      GCTA => 'THIS_DOESNT_EXIST', 
    ); 

    while (my $line = <INFILE>){ 
      chomp $line; 
      my $hash; 
      my @elements = split "\t", $line; 
      open my $out, '>', "/path/to/outfile.txt" or die $!; 
      foreach my $sequence(@elements){ 
        if (exists $hash{$sequence}){ 
         print $out ">$sequence\n$hash{$sequence}\n"; 
        } 
        else 
        } 
         $count++; 
         print "Doesn't exist ", $count, "\n"; 
        } 
      } 
    } 

我怎么知道有多少序列存在之前,我打印?我需要将该信息放入输出文件的名称中。

理想情况下,我会在文件名中包含一个变量。不幸的是,我不能仅仅采用@elements的标量,因为有些序列不会被打印出来。当我尝试将存在的键推入数组中,然后打印该数组的标量时,我仍然得不到所需的结果。这是我试过的(所有变量都需要全局性的):

open my $out, '>', "/path/to/file.$number.txt" or die $!;  
    foreach my $sequence(@elements){ 
      if (exists $hash{$sequence}){ 
        push(@Array, $hash{$sequence}, "\n"); 
        my $number = @Array; 
        print $out ">$sequence\n$hash{$sequence}\n"; 
      #.... 

感谢您的帮助。真的很感激它。

回答

2
my $sequences = grep exists $hash{$_}, @elements; 
open my $out, '>', "/path/to/outfile_containing_$sequences.txt" or die $!; 

在列表上下文中,grep按标准过滤列表;在标量上下文中,它返回满足条件的元素数。

+0

这太好了。谢谢。 – Rob

0

最简单的方法是跟踪您在一个变量中打印多少个键,并且一旦循环结束,只需使用计算的数字重命名该文件即可。 Perl带有一个内置函数来执行此操作。该代码将是这样的:

use strict; 
use warnings; 

open(INFILE, "<", "/path/to/infile.txt") or die $!; 

my %hash = ( 
     TAGC => 'THIS_EXISTS', 
     GCTA => 'THIS_DOESNT_EXIST', 
); 
my $ammt; 

while (my $line = <INFILE>){ 
     chomp $line; 
     my $hash; 
     my @elements = split "\t", $line; 
     open my $out, '>', "/path/to/outfile.txt" or die $!; 
     foreach my $sequence(@elements){ 
       if (exists $hash{$sequence}){ 
        print $out ">$sequence\n$hash{$sequence}\n"; 
        $ammt++; 
       } 
       else 
       } 
       print "Doesn't exist ", $count, "\n"; 
       } 
     } 
} 

rename "/path/to/outfile.txt", "/path/to/outfile${ammt}.txt" or die $!; 

我删除了$count变量,因为它不是在你的代码中声明(严会抱怨这一点)。 Here'srename的官方文档。既然它返回True或False,你可以检查它是否成功。

顺便说一句,请注意:

push(@Array, $hash{$sequence}, "\n"); 

被存储两个项目($hash{$sequence}\n),使计数将在两次因为它应该是。