view tools/filters/joinWrapper.pl @ 1:cdcb0ce84a1b

Uploaded
author xuebing
date Fri, 09 Mar 2012 19:45:15 -0500
parents 9071e359b9a3
children
line wrap: on
line source

#! /usr/bin/perl -w

use strict;
use warnings;
use File::Temp "tempfile";

my ($input1, $input2, $field1, $field2, $mode, $OOption, $out_file1) = @ARGV;

die "No arguments\n" unless @ARGV == 7;

my ($fh1, $file1) = tempfile();
my ($fh2, $file2) = tempfile(); 

`sort -k $field1 $input1 > $file1`;
`sort -k $field2 $input2 > $file2`;

my $option = "";
my @fields = ();
my $line = "";

if ($OOption eq "Y") {
  if (defined($fh1)) {
    $line = <$fh1>;
  } else {
    die "Failed to create file $file1\n";
  }
  @fields = split /\t/, $line;
  die "The field you selected does not exist in the input file" if (@fields < $field1);
  my @optionO = ();
  my $i = 0;
  foreach (@fields) {
    ++$i;
    push(@optionO, "1.$i");
  }
  $option = "-o " . join(",", @optionO);
} else {
  $option = "";
}

$ENV{'LC_ALL'} = 'POSIX';

if ($mode eq "V") {
  `join -v 1 $option -1 $field1 -2 $field2 $file1 $file2 | tr " " "\t" > $out_file1`;
} else {
  `join $option -1 $field1 -2 $field2 $file1 $file2 | tr " " "\t" > $out_file1`;
}

`rm $file1 ; rm $file2`;