Sha256: 26b83d3336bc9ae6dfec57ecacbb84b126863a04291f30b860407b43fc024325
Contents?: true
Size: 767 Bytes
Versions: 7
Compression:
Stored size: 767 Bytes
Contents
require 'spec_helper'
require 'ms/error_rate/qvalue'
Hit = Struct.new(:score, :charge)
describe 'calculating q-values' do
before do
scores = [14,15,13,12,11]
qvals_expected = [0.5 ,0.0, 2.0/3.0, 3.0/4, 4.0/5]
@target_hits = scores.zip(Array.new(scores.size, 2)).map {|pair| Hit.new(*pair) }
@decoy_hits = scores.zip(Array.new(scores.size, 2)).map {|pair| Hit.new(pair.first-0.5, pair.last) }
@qval_by_hit = {}
@target_hits.zip(qvals_expected) {|hit, qval| @qval_by_hit[hit] = qval }
end
it 'can calculate qvalues on target deccoy sets' do
pairs = Ms::ErrorRate::Qvalue.target_decoy_qvalues(@target_hits, @decoy_hits)
pairs.each do |hit, qval|
@qval_by_hit[hit].should.be.close(qval, 0.00000001)
end
end
end
Version data entries
7 entries across 7 versions & 1 rubygems