Commit e17174a8 authored by Mikaël Salson's avatar Mikaël Salson
Browse files

KmerMultiSegmenter: getNbExpected uses the e-value from left and right regions

parent a2471bbd
...@@ -368,7 +368,7 @@ KmerMultiSegmenter::KmerMultiSegmenter(Sequence seq, MultiGermline *multigermlin ...@@ -368,7 +368,7 @@ KmerMultiSegmenter::KmerMultiSegmenter(Sequence seq, MultiGermline *multigermlin
the_kseg->setSegmentationStatus(UNSEG_NOISY); the_kseg->setSegmentationStatus(UNSEG_NOISY);
} }
pair <double, double> p = the_kseg->getKmerAffectAnalyser()->getLeftRightProbabilityAtLeastOrAbove(); 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;
...@@ -382,8 +382,9 @@ KmerMultiSegmenter::KmerMultiSegmenter(Sequence seq, MultiGermline *multigermlin ...@@ -382,8 +382,9 @@ KmerMultiSegmenter::KmerMultiSegmenter(Sequence seq, MultiGermline *multigermlin
} }
double KmerMultiSegmenter::getNbExpected() const { double KmerMultiSegmenter::getNbExpected() const {
double proba = the_kseg->getKmerAffectAnalyser()->getProbabilityAtLeastOrAbove(the_kseg->score); pair <double, double> p = getNbExpectedLeftRight();
return multi_germline->germlines.size() * proba; return (p.first + p.second);
}
pair<double,double> KmerMultiSegmenter::getNbExpectedLeftRight() const { pair<double,double> KmerMultiSegmenter::getNbExpectedLeftRight() const {
pair <double, double> p = the_kseg->getKmerAffectAnalyser()->getLeftRightProbabilityAtLeastOrAbove(); pair <double, double> p = the_kseg->getKmerAffectAnalyser()->getLeftRightProbabilityAtLeastOrAbove();
......
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