Class: Ms::Quant::Qspec
- Inherits:
-
Object
- Object
- Ms::Quant::Qspec
- Defined in:
- lib/ms/quant/qspec.rb
Defined Under Namespace
Classes: Results
Constant Summary collapse
- NBURNIN =
personal communication with Hyungwon Choi: “We typically use nburn=2000, niter=10000, which is quite sufficient to guarantee the reproducibility of results using the same data.”
2000
- NITER =
10000
- INIT_HEADER =
%w(protid protLen)
- DELIMITER =
"\t"
- SUBMITTED_TO_QSPEC =
'submitted_to_qspec.txt'
Class Method Summary collapse
-
.conditions_to_ints(conditions) ⇒ Object
takes an ordered list of conditions [‘cond1’, ‘cond1’, ‘cond2’, ‘cond2’] and returns an array of ints [0,0,0,1,1,1…].
-
.executable(conditions) ⇒ Object
returns the right executable based on the array of conditions.
-
.results_array(resultsfile) ⇒ Object
returns an array of Results structs which is each row of the returned file works with V2 of QSpec.
Instance Method Summary collapse
- #conditions ⇒ Object
-
#initialize(protname_length_pairs, condition_to_count_array) ⇒ Qspec
constructor
protname_length_pairs is an array of doublets: [protname, length] condition_to_count_array is an array doublets: [condition, array_of_counts].
-
#run(normalize = true, opts = {}) ⇒ Object
returns an array of Qspec::Results objects (each object can be considered a row of data).
-
#write(filename) ⇒ Object
writes a qspec formatted file to filename.
Constructor Details
#initialize(protname_length_pairs, condition_to_count_array) ⇒ Qspec
protname_length_pairs is an array of doublets: [protname, length] condition_to_count_array is an array doublets: [condition, array_of_counts]
54 55 56 57 |
# File 'lib/ms/quant/qspec.rb', line 54 def initialize(protname_length_pairs, condition_to_count_array) @protname_length_pairs = protname_length_pairs @condition_to_count_array = condition_to_count_array end |
Class Method Details
.conditions_to_ints(conditions) ⇒ Object
takes an ordered list of conditions [‘cond1’, ‘cond1’, ‘cond2’, ‘cond2’] and returns an array of ints [0,0,0,1,1,1…]
18 19 20 21 22 23 24 25 26 27 28 29 |
# File 'lib/ms/quant/qspec.rb', line 18 def self.conditions_to_ints(conditions) i = 0 current_condition = conditions.first conditions.map do |cond| if current_condition == cond ; i else i += 1 current_condition = cond i end end end |
.executable(conditions) ⇒ Object
returns the right executable based on the array of conditions
47 48 49 50 |
# File 'lib/ms/quant/qspec.rb', line 47 def self.executable(conditions) biggest_size = conditions.group_by {|v| v }.values.map(&:size).max (biggest_size >= 3) ? 'qspecgp' : 'qspec' end |
.results_array(resultsfile) ⇒ Object
returns an array of Results structs which is each row of the returned file works with V2 of QSpec
33 34 35 36 37 38 39 40 41 42 43 44 |
# File 'lib/ms/quant/qspec.rb', line 33 def self.results_array(resultsfile) rows = IO.readlines(resultsfile).map {|line| line.chomp.split("\t") } headers = rows.shift start_bayes = headers.index {|v| v =~ /BayesFactor/i } rows.map do |row| data = [row[0]] data.push( row[1...start_bayes].map(&:to_f) ) data.push( *row[start_bayes,4].map(&:to_f) ) data.push( row[start_bayes+4] ) Results.new(*data) end end |
Instance Method Details
#conditions ⇒ Object
59 60 61 |
# File 'lib/ms/quant/qspec.rb', line 59 def conditions @condition_to_count_array.map(&:first) end |
#run(normalize = true, opts = {}) ⇒ Object
returns an array of Qspec::Results objects (each object can be considered a row of data)
79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 |
# File 'lib/ms/quant/qspec.rb', line 79 def run(normalize=true, opts={}) puts "normalize: #{normalize}" if $VERBOSE tfile = Tempfile.new("qspec") write(tfile.path) if opts[:keep] local_file = File.join(Dir.pwd,File.basename(tfile.path)) FileUtils.cp(tfile.path, local_file, :verbose => $VERBOSE) puts "(copy of) file submitted to qspec: #{local_file}" if $VERBOSE end qspec_exe = self.class.executable(conditions) cmd = [qspec_exe, tfile.path, NBURNIN, NITER, (normalize ? 1 : 0)].join(' ') if $VERBOSE puts "running #{cmd}" if $VERBOSE else cmd << " 2>&1" end reply = `#{cmd}` puts reply if $VERBOSE outfile = tfile.path + '_' + qspec_exe results = self.class.results_array(outfile) if opts[:keep] local_outfile = File.join(Dir.pwd, File.basename(outfile)) FileUtils.cp(outfile, local_outfile, :verbose => $VERBOSE) puts "(copy of) file returned from qspec: #{outfile}" end tfile.unlink results end |
#write(filename) ⇒ Object
writes a qspec formatted file to filename
64 65 66 67 68 69 70 71 72 73 74 75 |
# File 'lib/ms/quant/qspec.rb', line 64 def write(filename) ints = Ms::Quant::Qspec.conditions_to_ints(conditions) header_cats = INIT_HEADER + ints rows = @protname_length_pairs.map {|pair| pair.map.to_a } @condition_to_count_array.each do |cond,counts| rows.zip(counts) {|row,cnt| row << cnt } end File.open(filename,'w') do |out| out.puts header_cats.join(DELIMITER) rows.each {|row| out.puts row.join(DELIMITER) } end end |