add different thresholds for each method

lib/convergence_detection.ml

@@ -63,6 +63,15 @@ let plot_merge_results ~res_by_tools ~tree ~faa ~tsv : svg workflow =
     string opt
   ) |> seq ~sep:","
   in
+  let meths_t = List.map res_by_tools ~f:(fun res ->
+    let opt = match res with
+      | `Pcoc _ -> "PCOC:0.99"
+      | `Pcoc_gamma _ -> "PCOC_gamma:0.99"
+      | `Diffsel _ -> "Diffsel:0.9"
+    in
+    string opt
+  ) |> seq ~sep:","
+  in
   let package_diffsel_script_utils = tmp // "diffsel_script_utils.py" in
   let package_plot_data = tmp // "plot_data.py" in
   let script_plot_convergent_sites = tmp // "plot_convergent_sites.py" in
@@ -85,6 +94,7 @@ let plot_merge_results ~res_by_tools ~tree ~faa ~tsv : svg workflow =
           opt "-tree" dep tree ;
           opt "-out" ident out ;
           opt "-meth" ident meths ;
+          opt "-t" ident meths_t ;
           flag string "--all_sites" plot_all_sites ;
         ]
       ]