Commit d0246c45 authored by Mathieu Giraud's avatar Mathieu Giraud

core/segment.cpp: Karlin-Altschul statistics, e-value computation with...

core/segment.cpp: Karlin-Altschul statistics, e-value computation with sequence sizes in the FineSegmenter
parent b6b66c8d
......@@ -760,8 +760,8 @@ FineSegmenter::FineSegmenter(Sequence seq, Germline *germline, Cost segment_c,
}
/* E-values */
evalue_left = multiplier * segment_cost.toPValue(score_V[0].first);
evalue_right = multiplier * segment_cost.toPValue(score_J[0].first);
evalue_left = multiplier * seq.sequence.size() * germline->rep_5.totalSize() * segment_cost.toPValue(score_V[0].first);
evalue_right = multiplier * seq.sequence.size() * germline->rep_3.totalSize() * segment_cost.toPValue(score_J[0].first);
evalue = evalue_left + evalue_right ;
/* Unsegmentation causes */
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment