lib/lederhosen/tasks/cluster.rb in lederhosen-0.1.8 vs lib/lederhosen/tasks/cluster.rb in lederhosen-0.1.9

- old
+ new

@@ -6,23 +6,35 @@ class CLI desc "cluster", "cluster a fasta file using UCLUST" - method_option :input, :type => :string, :required => true - method_option :output, :type => :string, :required => true - method_option :identity, :type => :numeric, :required => true + method_option :input, :type => :string, :required => true + method_option :output, :type => :string, :required => true + method_option :identity, :type => :numeric, :required => true + method_option :stepwords, :type => :numeric, :default => 8 + method_option :wordlen, :type => :numeric, :default => 8 + method_option :maxaccepts, :type => :numeric, :default => 1 + method_option :maxrejects, :type => :numeric, :default => 8 def cluster - identity = options[:identity] - output = options[:output] - input = options[:input] + identity = options[:identity] + output = options[:output] + input = options[:input] + stepwords = options[:stepwords] + maxaccepts = options[:maxaccepts] + maxrejects = options[:maxrejects] + wordlen = options[:wordlen] cmd = [ 'uclust', "--input #{input}", "--uc #{output}", "--id #{identity}", + "--stepwords #{stepwords}", + "--maxaccepts #{maxaccepts}", + "--maxrejects #{maxrejects}", + "--w #{wordlen}" ].join(' ') @shell.mute { run cmd } end