inprogress
authorcmzmasek@gmail.com <cmzmasek@gmail.com@ca865154-3058-d1c3-3e42-d8f55a55bdbd>
Sat, 5 Jan 2013 00:42:04 +0000 (00:42 +0000)
committercmzmasek@gmail.com <cmzmasek@gmail.com@ca865154-3058-d1c3-3e42-d8f55a55bdbd>
Sat, 5 Jan 2013 00:42:04 +0000 (00:42 +0000)
forester/ruby/evoruby/lib/evo/tool/hmmscan_summary.rb

index a1baafb..3a9d662 100644 (file)
@@ -372,13 +372,14 @@ module Evoruby
         
         puts r.model
          puts r.fs_e_value
+          puts fs_e_value_threshold
         
         if r.model == hmm_for_protein_output
-          if r.fs_e_value > fs_e_value_threshold
+          if i_e_value_threshold >= 0.0 && r.fs_e_value > fs_e_value_threshold
             return
           end
         end
-        if r.i_e_value <= i_e_value_threshold
+        if i_e_value_threshold < 0 || r.i_e_value <= i_e_value_threshold
           hmmscan_results_per_protein_filtered << r
           if r.model == hmm_for_protein_output
             dc += 1