我正在处理包含制表符分隔序列的输入文件。序列组通过换行符分隔。该文件看起来像:如何计算散列中存在的键的数量?
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";
#....
感谢您的帮助。真的很感激它。
这太好了。谢谢。 – Rob