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

algo/core: Debugging KmerStore e-value

parent 68b6ebcb
Pipeline #34408 passed with stages
in 40 seconds
......@@ -381,6 +381,10 @@ double IKmerStore<T>::getProbabilityAtLeastOrAbove(const T kmer, int at_least, i
probability_not_having_system /= (1 - index_load);
}
#ifdef DEBUG_KMS_EVALUE
cerr << "e-value:\tindex_load=" << index_load << ",\tat_least=" << at_least << ",\tlength=" << length <<",\tp-value=" << proba << endl;
#endif
return proba;
}
......
......@@ -614,6 +614,9 @@ KmerMultiSegmenter::KmerMultiSegmenter(Sequence seq, MultiGermline *multigermlin
// E-value multiplier
double multiplier = multi_germline->germlines.size() * nb_reads_for_evalue;
#ifdef DEBUG_KMS_EVALUE
cerr << "multiplier: germline_size=" << multi_germline->germlines.size() <<", nb_reads_for_evalue=" << nb_reads_for_evalue << endl;
#endif
// Iterate over the germlines
for (list<Germline*>::const_iterator it = multigermline->germlines.begin(); it != multigermline->germlines.end(); ++it)
......@@ -690,6 +693,10 @@ void KmerSegmenter::computeSegmentation(int strand, KmerAffect before, KmerAffec
evalue_left = pvalues.first * multiplier ;
evalue_right = pvalues.second * multiplier ;
evalue = evalue_left + evalue_right ;
#ifdef DEBUG_KMS_EVALUE
cerr << "\tmultiplier=" << multiplier << ",\tevalues=[" << evalue_left << ", " << evalue_right << "],\t"
<< "evalue=" << evalue << endl;
#endif
// This can lead to UNSEG_TOO_FEW_ZERO or UNSEG_ONLY_V/J
checkLeftRightEvaluesThreshold(threshold, strand);
......
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