vidjil.cpp: safer D e-value computation, multiply by number of clones or reads

e6ffb915 · Mathieu Giraud · 47a7c984 · e6ffb915
Commit e6ffb915 authored Aug 03, 2016 by Mathieu Giraud
Hide whitespace changes
Inline Side-by-side

Showing with 2 additions and 2 deletions

algo/vidjil.cpp algo/vidjil.cpp +2 -2

No files found.
--- a/algo/vidjil.cpp
+++ b/algo/vidjil.cpp
@@ -1386,7 +1386,7 @@ int main (int argc, char **argv)
        FineSegmenter seg(representative, segmented_germline, segment_cost);
 	
        if (segmented_germline->seg_method == SEG_METHOD_543)
-	  seg.FineSegmentD(segmented_germline, several_D, expected_value_D);
+	  seg.FineSegmentD(segmented_germline, several_D, expected_value_D, sort_clones.size());

        if (detect_CDR3)
          seg.findCDR3();
@@ -1603,7 +1603,7 @@ int main (int argc, char **argv)
            if (s.isSegmented()) 
              {
                if (germline->seg_method == SEG_METHOD_543)
-                  s.FineSegmentD(germline, several_D, expected_value_D);
+                  s.FineSegmentD(germline, several_D, expected_value_D, nb_reads_for_evalue);

                if (detect_CDR3)
                  s.findCDR3();