lecture code viewer

downloads

Code

#!/usr/local/bin/perl

use Data::Dumper;

open(FH,"/home/martink/work/ucsc/hg17/bes/bacend.parsed.txt");
%clonesize = ();
while($line = <FH>) {
    chomp $line;
    ($clone,undef,$chr,$start,$end) = split(" ",$line);
    $size = $end - $start + 1;
    $clonesize{$clone} = $size;
}
close(FH);
open(RX,"rx.txt");
@rx = ();
while($line = <RX>) {
    chomp $line;
    next if $line =~ /^#/;
    push @rx, $line;
}
close(RX);

print Dumper(\@rx);

%sum = ();
%count = ();
for $clone (keys %clonesize) {
    for $rx (@rx) {
	if($clone =~ /$rx/) {
	    $sum{$rx}   += $clonesize{$clone};
	    $count{$rx} ++;
	}
    }
}

for $rx (sort keys %sum) {
    printf("group %10s num %6d avgsize %10.1f\n",$rx,$count{$rx},$sum{$rx}/$count{$rx});
}