Commit b0659963 authored by Mathieu Giraud's avatar Mathieu Giraud

core/segment.cpp: compute e-value

Fix 6ceffe02. We forgot to re-compute the total e-value.
parent 44d74052
......@@ -427,6 +427,7 @@ void KmerSegmenter::computeSegmentation(int strand, KmerAffect before, KmerAffec
pair <double, double> pvalues = kaa->getLeftRightProbabilityAtLeastOrAbove();
evalue_left = pvalues.first * multiplier ;
evalue_right = pvalues.second * multiplier ;
evalue = evalue_left + evalue_right ;
// We labeled it detected if there were both enough affect_5 and enough affect_3
detected = (max.nb_before_left + max.nb_before_right >= DETECT_THRESHOLD)
......
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