Class: Bioinform::CLI::ConvertMotif
- Inherits:
-
Object
- Object
- Bioinform::CLI::ConvertMotif
- Defined in:
- lib/bioinform/cli/convert_motif.rb
Class Method Summary collapse
Instance Method Summary collapse
- #arguments ⇒ Object
- #main(argv) ⇒ Object
- #option_parser ⇒ Object
- #options ⇒ Object
- #parse!(argv, default_options = {}) ⇒ Object
Class Method Details
.main(argv) ⇒ Object
136 137 138 |
# File 'lib/bioinform/cli/convert_motif.rb', line 136 def self.main(argv) self.new.main(argv) end |
Instance Method Details
#arguments ⇒ Object
11 12 13 |
# File 'lib/bioinform/cli/convert_motif.rb', line 11 def arguments @arguments ||= [] end |
#main(argv) ⇒ Object
18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 |
# File 'lib/bioinform/cli/convert_motif.rb', line 18 def main(argv) parse!(argv, filename_format: './{name}.{ext}') motif_files = arguments motif_files += $stdin.read.shellsplit unless $stdin.tty? if motif_files.empty? puts option_parser.help() return end output_motifs = [] motifs = motif_files.map do |filename| input = File.read(filename) case [:model_from] when 'pwm' MotifModel::PWM.from_string(input) when 'pcm' MotifModel::PCM.from_string(input) when 'ppm' MotifModel::PPM.from_string(input) else raise "Unknown value of model-from parameter: `#{[:model_from]}`" end end pcm2pwm_converter = ConversionAlgorithms::PCM2PWMConverter.new(pseudocount: :log, background: Background::Uniform) pcm2ppm_converter = ConversionAlgorithms::PCM2PPMConverter.new ppm2pcm_converter = ConversionAlgorithms::PPM2PCMConverter.new(count: 100) motifs.each do |motif| begin case [:model_to] when 'pwm' if MotifModel.acts_as_pcm?(motif) output_motifs << pcm2pwm_converter.convert(motif) elsif MotifModel.acts_as_ppm?(motif) output_motifs << pcm2pwm_converter.convert(ppm2pcm_converter.convert(motif)) elsif MotifModel.acts_as_pwm?(motif) output_motifs << motif else raise "Can't be here" end when 'pcm' if MotifModel.acts_as_pcm?(motif) output_motifs << motif elsif MotifModel.acts_as_ppm?(motif) output_motifs << ppm2pcm_converter.convert(motif) elsif MotifModel.acts_as_pwm?(motif) raise 'Not yet implemented' else raise "Can't be here" end when 'ppm' if MotifModel.acts_as_pcm?(motif) output_motifs << pcm2ppm_converter.convert(motif) elsif MotifModel.acts_as_ppm?(motif) output_motifs << motif elsif MotifModel.acts_as_pwm?(motif) raise 'Not yet implemented' else raise "Can't be here" end else raise "Unknown value of model-to parameter: `#{[:model_to]}`" end rescue $stderr.puts "One can't convert from #{[:model_from]} data-model to #{[:model_to]} data-model" raise end end puts output_motifs.join("\n\n") rescue => e $stderr.puts "Error! Conversion wasn't performed" $stderr.puts e $stderr.puts e.backtrace end |
#option_parser ⇒ Object
92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 |
# File 'lib/bioinform/cli/convert_motif.rb', line 92 def option_parser @option_parser ||= OptionParser.new do |cli| = <<-BANNER Usage: convert_motif [options] <motif-files>... ls | convert_motif [options] convert_motif - tool for converting motifs from different input formats to different output formats. It can change both formatting style and motif models. Resulting model is sent to stdout (this can be overriden with --save option). BANNER cli.version = ::Bioinform::VERSION cli.summary_indent = '' cli. = strip_doc() cli.separator "" cli.separator "Options:" cli.on('--parser PARSER', 'Parser for input motif.'){|parser| [:parser] = parser} cli.on('--formatter FORMATTER', 'Formatter for output motif.'){|formatter| [:formatter] = formatter} cli.on('--from MODEL_OF_INPUT', 'Specify motif model of input.', 'It can be overriden by --parser option', '(when parser implies certain input model)'){|model_from| [:model_from] = model_from} cli.on('--to MODEL_OF_OUTPUT', 'Specify motif model to convert to.', 'It can be overriden by --formatter option', '(when formatter implies certain output model)'){|model_to| [:model_to] = model_to} cli.on('--algorithm ALGORITHM', 'Conversion algorithm to transform model.'){|conversion_algorithm| [:conversion_algorithm] = conversion_algorithm} cli.on('--save [FILENAME_FORMAT]', 'Save resulting motifs to according files.', 'filename format by default is ' + [:filename_format], 'one can specify output folder here', '{name} is a placeholder for motif name', '{ext} is a placeholder for motif model') {|filename_format| [:filename_format] = filename_format if filename_format} cli.on('--[no-]force', 'Overwrite existing files.'){|force| [:force] = force} cli.on('--[no-]silent', 'Suppress error messages and notifications.'){|silent| [:silent] = silent} end end |
#options ⇒ Object
14 15 16 |
# File 'lib/bioinform/cli/convert_motif.rb', line 14 def @options ||= {} end |
#parse!(argv, default_options = {}) ⇒ Object
129 130 131 132 133 |
# File 'lib/bioinform/cli/convert_motif.rb', line 129 def parse!(argv, = {}) @options = .dup option_parser.parse!(argv) @arguments = argv end |