sub dataextractor($) { my $input = shift; my $output = shift; my $order = shift; my (@data); open(OUT,">".$output) or die "Can't open the file $output\n"; open(IN,"<".$input) or die "Can't open the file $input\n"; while () { if ($. eq 1) { print OUT $_; next; } s/\r?\n$//; @data = split (',',$_); if ($data[1] eq $site_peakhr{$data[$order]}) { print OUT $_."\n"; } }#while close(IN); close(OUT); return($output); }