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

segment: the probability of not having a system

depends on 1 - index_load, not the opposite
parent 0b8a6e8e
......@@ -372,7 +372,7 @@ double KmerSegmenter::getProbabilityAtLeastOrAbove(int at_least) const {
for (int i=at_least; i<=n; i++) {
proba += nChoosek(n, i) * probability_having_system * probability_not_having_system;
probability_having_system *= index_load;
probability_not_having_system /= index_load;
probability_not_having_system /= (1 - index_load);
}
return proba;
......
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