Perl - Eingabe abrufen, Bedingungen prüfen

Question

Ich bin davon ausgegangen, dass die extrahierten Regionslinien bei jeder Ausrichtung den CDS-Linien folgen.

Kopieren Sie diesen Code in script.pl:

use strict;
use warnings;

my $input = 1;
my @field = ('CDS','extracted region');
my (%data);
my (%counter);
&zero;

while ( <> ) { 

    ## Omit header.
    next if $. == 1;
    next if $. == 2;

    ## Remove last '\n'.
    chomp;

    ## Split line in tabs.
    my @f = split /\t/;

    ## Is loop over?
    if ( $f[2] =~ /$field[0]/ && $counter{$field[1]} > 1 )
    {
    &comparing;
    &zero;
    }

    ## Count number of $field[0] and $field[1] line
    $counter{$f[2]}++;

    ## Storing data
    @{$data{$f[2]}[$counter{$f[2]}]} = @f;
}

&comparing;

sub zero {
    $data{$field[0]} = [];
    $data{$field[1]} = [];
    $counter{$field[0]} = 0;
    $counter{$field[1]} = 0;
}

sub comparing {
    ## Is same line ($field[0] and $field[1])? if ( $input == 1 )
    if ( $counter{$field[0]} ==  $counter{$field[1]} || $input == 2 )
    {
    &recover;
    &stamp;
    }
}

sub recover {
    my $pos = &input2(0,0)  if ( $input == 2 );
    for my $i ( 1 .. $#{ $data{$field[0]} } ) {
    &input1($i) if ( $input == 1 );
    &input2($i,$pos) if ( $input == 2 );
    }
}

sub input1 {
    #;Extracted interval="376914 -> 377067"
    $data{$field[1]}[$_[0]][8] =~ m/;Extracted interval="(\d+) /;
    $data{$field[0]}[$_[0]][3] = $1;
    $data{$field[1]}[$_[0]][8] =~ m/;Extracted interval="\d+ -> (\d+)"/;
    $data{$field[0]}[$_[0]][4] = $1;
}

sub input2 {
    if ( $_[0] == 0 )
    {
    #;Extracted interval="2010140 <- 2024072"
    $data{$field[1]}[1][8] =~ m/;Extracted interval="(\d+) /;
    $1;
    }
    else
    {
    $data{$field[0]}[$_[0]][3] = $_[1] + $data{$field[0]}[$_[0]][3];
    $data{$field[0]}[$_[0]][4] = $_[1] + $data{$field[0]}[$_[0]][4];
    }
}

sub stamp {
    for my $i ( 1 .. $#{ $data{$field[0]} } ) {
    for my $j ( 0 .. $#{ $data{$field[0]}[$i] } ) {
        print "$data{$field[0]}[$i][$j]\t";
    }
    print "\n";
    }
}

Sie könnten ein Perl-Skript mit input1.txt ausführen:

perl script.pl input1.txt > output1.txt

wenn Sie die Zeile ändern:

my $input = 1;

mit

my $input = 2;

Sie könnten ein Perl-Skript mit input2.txt ausführen:

perl script.pl input2.txt > output2.txt

Ein Perl-Skript kann auch zwei Argumente erhalten: Eingabedatei und Typ [12].

BEARBEITEN

Inhttps://stackoverflow.com/questions/1730333/wie verwende ich Getoptions, um das Standardargument zu erhalten?Es gibt einige Methoden zum Abrufen von Argumenten.

Wenn Sie die Zeile ändern:

my $input = 1;

mit

my $input = 1;
$input = $ARGV[1] if defined $ARGV[1];

Sie könnten ein Perl-Skript mit input1.txt ausführen:

perl script.pl input1.txt > output1.txt

oder

perl script.pl input1.txt 1 > output1.txt

und Sie könnten ein Perl-Skript mit input2.txt ausführen:

perl script.pl input2.txt 2 > output2.txt

Answer 1

Ich bin davon ausgegangen, dass die extrahierten Regionslinien bei jeder Ausrichtung den CDS-Linien folgen.

Kopieren Sie diesen Code in script.pl:

use strict;
use warnings;

my $input = 1;
my @field = ('CDS','extracted region');
my (%data);
my (%counter);
&zero;

while ( <> ) { 

    ## Omit header.
    next if $. == 1;
    next if $. == 2;

    ## Remove last '\n'.
    chomp;

    ## Split line in tabs.
    my @f = split /\t/;

    ## Is loop over?
    if ( $f[2] =~ /$field[0]/ && $counter{$field[1]} > 1 )
    {
    &comparing;
    &zero;
    }

    ## Count number of $field[0] and $field[1] line
    $counter{$f[2]}++;

    ## Storing data
    @{$data{$f[2]}[$counter{$f[2]}]} = @f;
}

&comparing;

sub zero {
    $data{$field[0]} = [];
    $data{$field[1]} = [];
    $counter{$field[0]} = 0;
    $counter{$field[1]} = 0;
}

sub comparing {
    ## Is same line ($field[0] and $field[1])? if ( $input == 1 )
    if ( $counter{$field[0]} ==  $counter{$field[1]} || $input == 2 )
    {
    &recover;
    &stamp;
    }
}

sub recover {
    my $pos = &input2(0,0)  if ( $input == 2 );
    for my $i ( 1 .. $#{ $data{$field[0]} } ) {
    &input1($i) if ( $input == 1 );
    &input2($i,$pos) if ( $input == 2 );
    }
}

sub input1 {
    #;Extracted interval="376914 -> 377067"
    $data{$field[1]}[$_[0]][8] =~ m/;Extracted interval="(\d+) /;
    $data{$field[0]}[$_[0]][3] = $1;
    $data{$field[1]}[$_[0]][8] =~ m/;Extracted interval="\d+ -> (\d+)"/;
    $data{$field[0]}[$_[0]][4] = $1;
}

sub input2 {
    if ( $_[0] == 0 )
    {
    #;Extracted interval="2010140 <- 2024072"
    $data{$field[1]}[1][8] =~ m/;Extracted interval="(\d+) /;
    $1;
    }
    else
    {
    $data{$field[0]}[$_[0]][3] = $_[1] + $data{$field[0]}[$_[0]][3];
    $data{$field[0]}[$_[0]][4] = $_[1] + $data{$field[0]}[$_[0]][4];
    }
}

sub stamp {
    for my $i ( 1 .. $#{ $data{$field[0]} } ) {
    for my $j ( 0 .. $#{ $data{$field[0]}[$i] } ) {
        print "$data{$field[0]}[$i][$j]\t";
    }
    print "\n";
    }
}

Sie könnten ein Perl-Skript mit input1.txt ausführen:

perl script.pl input1.txt > output1.txt

wenn Sie die Zeile ändern:

my $input = 1;

mit

my $input = 2;

Sie könnten ein Perl-Skript mit input2.txt ausführen:

perl script.pl input2.txt > output2.txt

Ein Perl-Skript kann auch zwei Argumente erhalten: Eingabedatei und Typ [12].

BEARBEITEN

Inhttps://stackoverflow.com/questions/1730333/wie verwende ich Getoptions, um das Standardargument zu erhalten?Es gibt einige Methoden zum Abrufen von Argumenten.

Wenn Sie die Zeile ändern:

my $input = 1;

mit

my $input = 1;
$input = $ARGV[1] if defined $ARGV[1];

Sie könnten ein Perl-Skript mit input1.txt ausführen:

perl script.pl input1.txt > output1.txt

oder

perl script.pl input1.txt 1 > output1.txt

und Sie könnten ein Perl-Skript mit input2.txt ausführen:

perl script.pl input2.txt 2 > output2.txt

Perl - Eingabe abrufen, Bedingungen prüfen

Antwort1

verwandte Informationen