aboutsummaryrefslogtreecommitdiff
path: root/scripts/frequency
diff options
context:
space:
mode:
Diffstat (limited to 'scripts/frequency')
-rwxr-xr-xscripts/frequency51
1 files changed, 51 insertions, 0 deletions
diff --git a/scripts/frequency b/scripts/frequency
new file mode 100755
index 00000000..1f9c66e6
--- /dev/null
+++ b/scripts/frequency
@@ -0,0 +1,51 @@
+#!/usr/bin/perl -w
+
+use strict;
+
+my $bins = 200;
+open(FH, $ARGV[0]);
+
+my $min = 99999999999999999999.0;
+my $max = 0.0;
+my $line;
+
+while ( $line = <FH> ) {
+
+ chomp($line);
+ if ( $line < $min ) {
+ $min = $line;
+ }
+ if ( $line > $max ) {
+ $max = $line;
+ }
+
+}
+
+printf("%f -> %f\n", $min, $max);
+
+seek(FH, 0, 0);
+
+my $binsize = ($max+1-$min)/$bins;
+my @hist;
+my $i;
+
+for ( $i=0; $i<$bins; $i++ ) {
+ $hist[$i] = 0;
+}
+
+while ( $line = <FH> ) {
+
+ my $bin;
+
+ chomp($line);
+
+ $bin = int(($line-$min)/$binsize);
+ $hist[$bin]++;
+
+}
+
+for ( $i=0; $i<$bins; $i++ ) {
+ printf("%f\t%i\n", $min+(($i+0.5)*$binsize), $hist[$i]);
+}
+
+close(FH);