Commit 40c17d5e authored by Mikaël Salson's avatar Mikaël Salson

KmerMultiSegmenter: TOO_FEW_ZERO (instead of UNSEG_NOISY) iff the e_value is...

KmerMultiSegmenter: TOO_FEW_ZERO (instead of UNSEG_NOISY) iff the e_value is bad globally not either on left or on right
parent e17174a8
...@@ -364,14 +364,17 @@ KmerMultiSegmenter::KmerMultiSegmenter(Sequence seq, MultiGermline *multigermlin ...@@ -364,14 +364,17 @@ KmerMultiSegmenter::KmerMultiSegmenter(Sequence seq, MultiGermline *multigermlin
if (the_kseg->isSegmented()) { if (the_kseg->isSegmented()) {
// the_kseg->evalue also depends on the number of germlines from the *Multi*KmerSegmenter // the_kseg->evalue also depends on the number of germlines from the *Multi*KmerSegmenter
the_kseg->evalue = getNbExpected(); the_kseg->evalue = getNbExpected();
if (the_kseg->evalue > threshold_nb_expected) {
the_kseg->setSegmentationStatus(UNSEG_NOISY);
}
pair <double, double> p = getNbExpectedLeftRight(); pair <double, double> p = getNbExpectedLeftRight();
the_kseg->evalue_left = p.first; the_kseg->evalue_left = p.first;
the_kseg->evalue_right = p.second; the_kseg->evalue_right = p.second;
if (the_kseg->evalue > threshold_nb_expected
&& the_kseg->evalue_left <= threshold_nb_expected
&& the_kseg->evalue_right <= threshold_nb_expected) {
the_kseg->setSegmentationStatus(UNSEG_TOO_FEW_ZERO);
}
if (the_kseg->evalue_left > threshold_nb_expected) { if (the_kseg->evalue_left > threshold_nb_expected) {
the_kseg->setSegmentationStatus(UNSEG_NOISY); // TOO_FEW_V ? the_kseg->setSegmentationStatus(UNSEG_NOISY); // TOO_FEW_V ?
} }
......
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