Commit 8f7219b8 authored by Mathieu Giraud's avatar Mathieu Giraud
Browse files

vidjil.cpp: FineSegmenter in -c clones also considers e-value

Follows what is done for D in e6ffb915, and what was already done for V/J in -c segment.
Even if it may be counter-intuitive, the results of segmentation thus actually depends on the number of reads:
'-e 1' should mean that we expect one false positive segmentation in the output.
Even if some e-values slightly change, this commit does not affect the results of the tests.

We should investigate what happens when a FineSegmenter fails here.
Moreover, is it sane to have the same '-e' paramater for fixing both e-values of Kmer and FineSegmenter?
Or should we rather have another paramater?
parent e51fc02e
...@@ -1410,7 +1410,7 @@ int main (int argc, char **argv) ...@@ -1410,7 +1410,7 @@ int main (int argc, char **argv)
// FineSegmenter // FineSegmenter
FineSegmenter seg(representative, segmented_germline, segment_cost); FineSegmenter seg(representative, segmented_germline, segment_cost, expected_value, sort_clones.size());
if (segmented_germline->seg_method == SEG_METHOD_543) if (segmented_germline->seg_method == SEG_METHOD_543)
seg.FineSegmentD(segmented_germline, several_D, expected_value_D, sort_clones.size()); seg.FineSegmentD(segmented_germline, several_D, expected_value_D, sort_clones.size());
......
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