Commit faf2ad94 authored by Mathieu Giraud's avatar Mathieu Giraud

core/dynprog.cpp, core/segment.h: require at least MIN_MATCHES matches

Waiting for a true p-value computation, we just require here that
the score is at least the one of 10 matches.
parent c8d13683
......@@ -102,7 +102,7 @@ int Cost::homo2(char xa, char xb, char y)
double Cost::toPValue(int score)
{
// TODO: compute an actual p-value
return (score <= 0) ? BAD_EVALUE : 1 / (double) score ;
return (score <= MIN_MATCHES * match) ? BAD_EVALUE : 1 / (double) score ;
}
......
......@@ -16,6 +16,7 @@
#define EXTEND_D_ZONE 5
#define MIN_D_LENGTH 5 /* If a D-REGION is smaller than this threshold, it is not output */
#define MIN_MATCHES 10 /* If a V/J-REGION does not give an alignment score with at least this number of matches, the FineSegmenter does not segment the sequence */
#define RATIO_STRAND 2 /* The ratio between the affectations in one
strand and the other, to safely attribute a
......
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