fixed encoding bug in read_454

[biopieces.git] / bp_bin / analyze_vals

diff --git a/bp_bin/analyze_vals b/bp_bin/analyze_vals
index aa9d12ae536b280e196598fd54a03cbc03caad2b..951f321f808e7b4e907d11d791e29c37cb1d64c2 100755 (executable)
--- a/bp_bin/analyze_vals
+++ b/bp_bin/analyze_vals
@@ -31,8 +31,10 @@
 require 'maasha/biopieces'
 
 casts = []
-casts << {long: 'keys',    short: 'k', type: 'list', mandatory: false, default: nil, allowed: nil, disallowed: nil}
-casts << {long: 'no_keys', short: 'K', type: 'list', mandatory: false, default: nil, allowed: nil, disallowed: nil}
+casts << {long: 'keys',      short: 'k', type: 'list', mandatory: false, default: nil, allowed: nil, disallowed: nil}
+casts << {long: 'no_keys',   short: 'K', type: 'list', mandatory: false, default: nil, allowed: nil, disallowed: nil}
+casts << {long: 'no_stream', short: 'x', type: 'flag', mandatory: false, default: nil, allowed: nil, disallowed: nil}
+casts << {long: 'data_out',  short: 'o', type: 'file', mandatory: false, default: nil, allowed: nil, disallowed: nil}
 
 options = Biopieces.options_parse(ARGV, casts)
 
@@ -91,6 +93,12 @@ Biopieces.open(options[:stream_in], options[:stream_out]) do |input, output|
 
       stats[key][:count] += 1
     end
+
+    output.puts record unless options[:no_stream]
+  end
+
+  if options[:data_out]
+    data_out = File.open(options[:data_out], 'w')
   end
 
   stats.each do |key, value|
@@ -98,13 +106,17 @@ Biopieces.open(options[:stream_in], options[:stream_out]) do |input, output|
       :KEY   => key,
       :TYPE  => value[:type].to_s.capitalize,
       :COUNT => value[:count],
-      :MIN   => value[:min],
-      :MAX   => value[:max],
-      :SUM   => value[:sum],
-      :MEAN  => (value[:sum] / value[:count].to_f).round(1)
+      :MIN   => value[:min].is_a?(Float) ? "%0.2f" % value[:min] : value[:min],
+      :MAX   => value[:max].is_a?(Float) ? "%0.2f" % value[:max] : value[:max],
+      :SUM   => value[:sum].is_a?(Float) ? "%0.2f" % value[:sum] : value[:sum],
+      :MEAN  => "%0.2f" % (value[:sum] / value[:count].to_f)
     }
 
-    output.puts stat_record
+    if options[:data_out]
+      data_out.puts stat_record
+    else
+      output.puts stat_record
+    end
   end
 end