view SMART/DiffExpAnal/countNumber.pl @ 31:0ab839023fe4

Uploaded

#!/usr/bin/perl -w

use strict;

my $in_file = $ARGV[0];
my $out_file = $ARGV[1];
my $sort_type = $ARGV[2]; # n(umeric) or a(lphanumeric)
my ($line,$ID,$nbOverlaps,%hash);

open(IN, $in_file);
while ($line = <IN>){
	chomp($line);
	$line=~s/\t/|/g;
	my @part=split(/\|/,$line);
	my @split=split(";",$part[$#part]);
	$split[0] =~ m/^(\w+).+$/;
	
	foreach my $i (@split){
		if ($i=~m/nbOverlaps=(.+)/){
			$nbOverlaps=$1;
		}
		if ($i=~m/gene_id=(.+)/){
			$ID=$1;
			$hash{$ID}=$nbOverlaps;
		}
	}
}
close(IN);

open(OUT, ">$out_file");
foreach my $key ( sort keys %hash) {
	print OUT "$key\t$hash{$key}\n";
}
close(OUT);