Commit c8363867 authored by Carine Rey's avatar Carine Rey
Browse files

add in the repo the figure for the paper

parent 9c41e5c0
......@@ -16,6 +16,7 @@ type t_choices = {
t_choices_recall09: text_file workflow ;
t_choices_plot: text_file workflow ;
t_choices_condensed_plot: text_file workflow ;
rp_plot: text_file workflow ;
tree_prefix: string;
}
......@@ -272,8 +273,9 @@ let get_t_choices ~tree_prefix ~(dataset_results_l: dataset_res list) : t_choice
let t_choices_complete = t_choices_dir / selector ["out.complete.tsv"] in
let t_choices_plot = t_choices_dir / selector ["out.pdf"] in
let t_choices_condensed_plot = t_choices_dir / selector ["out_condensed.pdf"] in
let rp_plot = t_choices_dir / selector ["out_condensed.recall_precision_ok.pdf"] in
let tree_prefix = ha_PCOC.tree_prefix in
Some {t_choices_max; t_choices_recall09; t_choices_complete ; t_choices_plot; t_choices_condensed_plot; tree_prefix}
Some {t_choices_max; t_choices_recall09; t_choices_complete ; t_choices_plot; t_choices_condensed_plot; rp_plot; tree_prefix}
| _ -> None
......@@ -424,6 +426,7 @@ let repo_of_post_analyses_res ~prefix ~post_analyses_res =
Repo.[
item [prefix ^ ".t_choices.pdf"] w.t_choices_plot ;
item [prefix ^ ".t_choices.condensed.pdf"] w.t_choices_condensed_plot ;
item [prefix ^ ".recall_precision_ok.pdf"] w.rp_plot ;
]
] |> List.concat
);
......
......@@ -388,6 +388,44 @@ plot_out = function(df_out, df_d , df_recall_sup09_per_meth, meths = NULL, suffi
base_aspect_ratio = 1,
limitsize = FALSE
)
plot = ggplot(subset(df_out,couple%in%c("H0/HaPC NeG5", "H0/HaPC NeG5_NeC_div2" , "H0/HaPC NeG5_NeC_x2")), aes(x=sensitivity, y=precision98_02, col = methode))
plot = plot + theme_bw()
plot = plot + labs(x="Sensitivity (= Recall)", y="Precision (98/2)")
plot = plot + theme(legend.position="top")
plot = plot + ylim(c(0,1)) + xlim(c(0,1))
plot = plot + guides(fill=FALSE)
plot = plot + scale_color_manual(values=colors2)
#plot = plot + geom_point(size=1, alpha=alpha)
plot = plot + geom_step(direction="vh", size=0.7, alpha=alpha)
plot = plot + geom_hline( aes(yintercept = 0.9), col="black" , size = 0.5, show.legend = NA,linetype="dashed")
plot = plot + facet_grid(. ~ couple )
plot = plot + theme(axis.text.x = element_text(angle = 45, hjust = 1))
legend_PR = get_legend(plot + theme(legend.position="top",
legend.text = element_text(size=10),
legend.title = element_text(size=0),
plot.margin = unit(x = c(0, 0, 0, 0), units = "cm"),
legend.background = element_rect(fill="white", size=0.5, linetype="solid",
colour ="black")) +
guides( colour = guide_legend(override.aes = list(alpha = 1,size = 3), nrow=2))
)
plot_recall_precision_papier_ok = plot_grid(legend_PR, plot + theme(legend.position="none"), ncol = 1, scale = 1,
labels = c("",""),
align = "h",
rel_heights = c( 0.3, 0.7),
hjust = 0, vjust = 1)
save_plot(paste0(opt$out,suffix,".recall_precision_ok.pdf"),
plot_recall_precision_papier_ok,
ncol = 0.6* 3,
nrow = 1,
base_aspect_ratio = 1,
limitsize = FALSE
)
} else {
plot_recall_precision_papier = NULL
}
......
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