☰

Perl-最小値と最大値

tag-icon • tag-icon perl
Perl-最小値と最大値

入力ファイル

Xm_ABL1 Geneious    extracted region    1   168 .   +   .   Name=Extracted region from gi|371443098|gb|JH556762.1|;Extracted interval="3512970000000 -> 3514640000000"
Xm_ABL1 Geneious    extracted region    169 334 .   +   .   Name=Extracted region from gi|371443098|gb|JH556762.1|;Extracted interval="3717850000000 -> 3719500000000"

Perlコードの一部

 if ($array[1] =~ /extracted region/){
            die "No CDS record for $key!\n" unless $metadata->{$key};
    (my $label = $array[7]) =~ s/.*region from (.*)\|;.*/$1/;
    $label =~ s/\|/_/g;
    $group->{$label} ||= { 
            pos1 => 1e10,
            pos2 => 0,
            metadata => $metadata->{$key},
            sequences => [],
    };
    (my $pos1, my $arr, my $pos2) = ($array[7]=~/.*interval=\"(\d+) (<?->?) (\d+)\"$/gm);
    # capture hi/lo values for group
    $group->{$label}->{pos1} = $pos1 if $pos1 < $group->{$label}->{pos1};
    $group->{$label}->{pos2} = $pos2 if $pos2 > $group->{$label}->{pos2};
    # push this sequence onto the group's array
    push(@{ $group->{$label}->{sequences} }, [ $pos1, $pos2, $arrow->{$arr} ]);
}

コード $pos1=3512970000000,3717850000000 & $pos2=3514640000000,3719500000000。私のコードは、最小値と最大値を検索する新しい行を出力します ($pos1 が 10,000 未満の場合)。ただし、値が 10,0000 より大きい場合は、pos1 の最小値の出力でエラーが発生します。$pos1 の最小値と $pos2 の最大値を検索するためのデバッグについて、ご協力いただければ幸いです。

答え1

$group->{$label}->{pos1}常に定義されているので、6行目で初期化しました。

pos1 => 1e10,

つまり、 まで正しく動作します。最小値(サンプル入力のように)$pos1 <= 1e10を印刷したい場合は、初期化する必要があります(行 6):$group->{$label}->{pos1} = 3512970000000$group->{$label}->{pos1} = -1

pos1 => -1,

13行目を次のように変更します。

$group->{$label}->{pos1} = $pos1 if (($group->{$label}->{pos1} < 0) || ($pos1 < $group->{$label}->{pos1}));

答え2

13行目のコードを変更して、defiendを使用する

 $group->{$label}->{pos1} = $pos1 if ((!defined $group->{$label}->{pos1})
 ||  ($pos1 < $group->{$label}->{pos1}));

関連情報