问在perl脚本中调用awk

EN

#!/usr/bin/env perl

use strict;
use warnings;
use autodie;  # avoid a bunch of `or die` clauses

# First, load the criteria for splitting into the output files
my %files = ();

# extra block level wrapping the $ranges file access; file
# is automatically closed at the end of the block
{ 
  open my $ranges, '<', 'sim.snps';
  while (<$ranges>)  {
    (undef, my $key, undef, my ($min, $max)) = split;
    $files{$key} = { min => $min, max => $max };

    # go ahead and open the output file while we're here
    open $files{$file}{fh}, '>', "$key.legend";
  }
} # $ranges filehandle closed here

# another file-access block
{
  # open the data file
  open my $data, '<', '../phasing_and_imputation/1000GP_Phase3_chr1.legend';

  while (<$data>) {
    # split the data into fields
    my @f = split;

    # loop over the output files and write the relevant parts of this line
    # to the ones that want it
    while (my ($file, $data) = each %files) {
      if ($f[1] >= $data->{min}  && $f[1] <= $data->{max}) {
        print { $data->{fh} } join(' ', @f[0..3]), "\n";
      }
    }
  }
} # data file closed here

# close the output files
foreach my $data (values %files) {
  close $data->{fh};
}