-1

我有一个很困扰我的问题...我有一个包含两列的文件（感谢您在上一个问题中的帮助），例如：

14430001        0.040
14430002        0.000
14430003        0.990
14430004        1.000
14430005        0.050
14430006        0.490
....................

第一列是坐标第二个概率。我试图找到概率 >=0.990 并且大小超过 100 的块。作为输出我想是这样的：

14430001        14430250
14431100        14431328
18750003        18750345
.......................

其中第一列是每个块的起点坐标，第二列是它的终点。

我写了这个脚本：

use strict;
#use warnings;
use POSIX;

my $scores_file = $ARGV[0];

#finds the highly conserved subsequences


open my $scores_info, $scores_file or die "Could not open $scores_file: $!";
#open(my $fh, '>', $coords_file) or die;
my $count = 0;
my $cons = "";
my $newcons = "";
while( my $sline = <$scores_info>)  {
      my @data = split('\t', $sline);
      my $coord = $data[0];
      my $prob = $data[1];
     if ($data[1] >= 0.990) {
      #$cons = "$cons + '\n' + $sline + '\n'";
    $cons = join("\n", $cons, $sline);
    # print $cons;
     $count++;
    if($count >= 100) {

    $newcons = join("\n", $newcons, $cons);
    my @array = split /'\n'/, $newcons;
    print @array;
            }
}
 else {
   $cons = "";
   $count = 0;
   }

}

它给了我概率> = 0.990（第一个如果有效）的线条，但坐标是错误的。当我尝试将它打印到它堆叠的文件中时，所以我只有一个样本来检查它。如果我的解释没有帮助，我很抱歉，但我是编程新手。

拜托，我需要你的帮助......非常感谢你提前！！！

4

1 回答 1

2

您似乎使用了太多变量。此外，在拆分数组并将其部分分配给变量后，使用新变量而不是原始数组。

sub output {
    my ($from, $to) = @_;
    print "$from\t$to\n";
}

my $threshold = 0.980;   # Or is it 0.990?
my $count = 0;
my ($start, $last);
while (my $sline = <$scores_info>) {
    my ($coord, $prob) = split /\t/, $sline;
    if ($prob >= $threshold) {
        $count++;
        defined $start or $start = $coord;
        $last = $coord;
    } else {
        output($start, $last) if $count > 100;
        undef $start;
        $count = 0;
    }
}
output($start, $last) if $count > 100;

（未经测试）

于 2013-08-25T23:16:46.220 回答