Commit d06c4dbd authored by Philippe Veber's avatar Philippe Veber
Browse files

update wrt biocaml

parent 9fba8752
......@@ -220,7 +220,7 @@ let oracle d =
Convergence_detection.oracle ~n_h0 ~n_ha
let multinomial_benchmark d =
Utils.rp_curve
Utils.recall_precision_curve
~oracle:(oracle d)
~labels:[|"LRT";"LRTsim";"sparse";"sparse_sim"|]
~results:[
......
......@@ -31,7 +31,7 @@ let%pworkflow fasta_cappend fa1 fa2 =
|> List.map ~f:(fun (d, s) -> sprintf ">%s\n%s" d s)
|> Out_channel.write_lines [%dest]
let%pworkflow rp_curve ~labels ~oracle ~results : Bistro.pdf Bistro.pworkflow =
let%pworkflow recall_precision_curve ~labels ~oracle ~results : Bistro.pdf Bistro.pworkflow =
let n = Array.length labels in
let colors = Array.sub [|"blue";"green";"red";"tan1";"royalblue4";"darkorchid"|] ~pos:0 ~len:n in
let labels = [%param labels] in
......@@ -49,7 +49,8 @@ let%pworkflow rp_curve ~labels ~oracle ~results : Bistro.pdf Bistro.pworkflow =
let results = List.map ~f:load_results [%eval Bistro.Workflow.eval_paths results] in
let result_curves =
List.map results ~f:(fun scores ->
Biocaml_unix.Bin_pred.rp_curve ~scores ~labels:oracle
Biocaml_unix.Bin_pred.recall_precision_curve ~scores ~labels:oracle
|> fst
|> Array.unzip
)
in
......
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