Commit a2471bbd authored by Mikaël Salson's avatar Mikaël Salson

KmerSegmenter: Add getNbExpectedLeftRight()

parent afaae79d
......@@ -384,6 +384,10 @@ KmerMultiSegmenter::KmerMultiSegmenter(Sequence seq, MultiGermline *multigermlin
double KmerMultiSegmenter::getNbExpected() const {
double proba = the_kseg->getKmerAffectAnalyser()->getProbabilityAtLeastOrAbove(the_kseg->score);
return multi_germline->germlines.size() * proba;
pair<double,double> KmerMultiSegmenter::getNbExpectedLeftRight() const {
pair <double, double> p = the_kseg->getKmerAffectAnalyser()->getLeftRightProbabilityAtLeastOrAbove();
return pair<double, double>(p.first * multi_germline->germlines.size(), p.second * multi_germline->germlines.size());
}
KmerMultiSegmenter::~KmerMultiSegmenter() {
......
......@@ -214,6 +214,12 @@ class KmerMultiSegmenter
*/
double getNbExpected() const;
/**
* @return expected number of Segmenter that would have yield the maximum score by chance
* on the left part of the read and on the right part of the read respectively.
*/
pair<double,double> getNbExpectedLeftRight() const;
~KmerMultiSegmenter();
KmerSegmenter *the_kseg;
......
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