#!/usr/bin/perl -w
use strict;

# gravem les dades en un fitxer per crear el file handle
# creem el file handle 5000GENS

if (scalar(@ARGV) < 4) {
print "seleccio.pl nomfitxer ratio pvalue mintumors\n";
exit(1);
}

my $nomfitxer = $ARGV[0];
my $ratio = $ARGV[1];
my $pvalue = $ARGV[2];
my $mintumors = $ARGV[3];

open (GENS,"< $nomfitxer");

my $linia;
$linia = <GENS>;

while (<GENS>) {
chomp;
my @v = split(/\t/);
my $n = 0;
my $i=2;
while ($i < scalar(@v)) {

if (($v[$i] >= $ratio || $v[$i] <= -$ratio) && $v[$i+1] <= $pvalue) {
$n = $n + 1;
}

$i = $i + 2;
}

if ($n >= $mintumors) {

print " $v[0]\t$v[1]";

$i = 2;

while ($i <= scalar(@v)-2) {
print "\t$v[$i]";
$i = $i + 2;
}
print "\n";

}
}

close (GENS);