Skip to content
GitLab
Projects
Groups
Snippets
Help
Loading...
Help
Help
Support
Community forum
Keyboard shortcuts
?
Submit feedback
Contribute to GitLab
Sign in / Register
Toggle navigation
Open sidebar
VEBER Philippe
codepi
Commits
bda12e50
Commit
bda12e50
authored
Sep 24, 2018
by
Carine Rey
Browse files
add more hyps in plots
parent
d1645dc4
Changes
3
Expand all
Hide whitespace changes
Inline
Side-by-side
Showing
3 changed files
with
235 additions
and
34 deletions
+235
-34
lib/pipeline.ml
lib/pipeline.ml
+1
-0
lib/post_analyses.ml
lib/post_analyses.ml
+107
-16
lib/scripts/calc_t_per_meth.R
lib/scripts/calc_t_per_meth.R
+127
-18
No files found.
lib/pipeline.ml
View file @
bda12e50
...
...
@@ -94,6 +94,7 @@ let derive_from_tree ~tree_dir ~tree ~profile ~preview ~use_concat ~ns ~no_Ne ~n
Variable
(
1
.,
4
.
);
Variable
(
4
.,
1
.
);
Variable
(
4
.,
8
.
);
Variable
(
8
.,
4
.
);
Variable
(
16
.,
1
.
);
Variable
(
1
.,
8
.
);
Variable
(
8
.,
1
.
);
]
in
let
models
=
h0_hapc_forall
nes_list
in
...
...
lib/post_analyses.ml
View file @
bda12e50
...
...
@@ -73,8 +73,8 @@ let build_cmd_t_choices (opt_name : string) mr_option =
|
None
->
[]
let
make_t_choices
?
(
tree_prefix
=
""
)
~
haPCOC_mr
~
h0_NeG1_mr
~
h0_NeG2_mr
~
h0_NeG3_mr
~
h0_NeG4_mr
~
h0_NeG5_mr
~
h0_NeG6_mr
~
h0_NeG8_mr
~
haPC_NeG1_mr
~
haPC_NeG2_mr
~
haPC_NeG3_mr
~
haPC_NeG4_mr
~
haPC_NeG5_mr
~
haPC_NeG6_mr
~
haPC_NeG8_mr
~
h0_NeG1_mr
~
h0_NeG2_mr
~
h0_NeG3_mr
~
h0_NeG4_mr
~
h0_NeG5_mr
~
h0_NeG6_mr
~
h0_NeG8_mr
~
h0_NeG16_mr
~
haPC_NeG1_mr
~
haPC_NeG2_mr
~
haPC_NeG3_mr
~
haPC_NeG4_mr
~
haPC_NeG5_mr
~
haPC_NeG6_mr
~
haPC_NeG8_mr
~
haPC_NeG16_mr
~
haPC_NeG5_NeC_div2_mr
~
haPC_NeG5_NeC_x2_mr
~
h0_NeG5_NeC_div2_mr
...
...
@@ -84,12 +84,22 @@ let make_t_choices ?(tree_prefix="") ~haPCOC_mr
~
haPC_NeG1_NeC_5_mr
~
h0_NeG1_NeC_5_mr
~
h0_NeG5_NeC_1_mr
~
haPC_NeG1_NeC_8_mr
~
haPC_NeG8_NeC_1_mr
~
h0_NeG8_NeC_1_mr
~
h0_NeG1_NeC_8_mr
~
haPC_NeG1_NeC_16_mr
~
haPC_NeG16_NeC_1_mr
~
h0_NeG16_NeC_1_mr
~
h0_NeG1_NeC_16_mr
~
haPC_NeG4_NeC_1_mr
~
haPC_NeG1_NeC_4_mr
~
haPC_NeG4_NeC_8_mr
~
haPC_NeG8_NeC_4_mr
~
h0_NeG1_NeC_4_mr
~
h0_NeG4_NeC_1_mr
~
h0_NeG4_NeC_8_mr
~
h0_NeG8_NeC_4_mr
~
haPC_NeG4_NeC_6_mr
~
haPC_NeG4_NeC_2_mr
~
haPC_NeG2_NeC_4_mr
...
...
@@ -108,6 +118,7 @@ let make_t_choices ?(tree_prefix="") ~haPCOC_mr
(
"--H0_NeG5"
,
h0_NeG5_mr
);
(
"--H0_NeG6"
,
h0_NeG6_mr
);
(
"--H0_NeG8"
,
h0_NeG8_mr
);
(
"--H0_NeG16"
,
h0_NeG16_mr
);
(
"--HaPC_NeG1"
,
haPC_NeG1_mr
);
(
"--HaPC_NeG2"
,
haPC_NeG2_mr
);
(
"--HaPC_NeG3"
,
haPC_NeG3_mr
);
...
...
@@ -115,16 +126,27 @@ let make_t_choices ?(tree_prefix="") ~haPCOC_mr
(
"--HaPC_NeG5"
,
haPC_NeG5_mr
);
(
"--HaPC_NeG6"
,
haPC_NeG6_mr
);
(
"--HaPC_NeG8"
,
haPC_NeG8_mr
);
(
"--HaPC_NeG16"
,
haPC_NeG16_mr
);
(
"--HaPC_NeG5_NeC_div2_res"
,
haPC_NeG5_NeC_div2_mr
);
(
"--HaPC_NeG5_NeC_x2_res"
,
haPC_NeG5_NeC_x2_mr
);
(
"--H0_NeG5_NeC_div2_res"
,
h0_NeG5_NeC_div2_mr
);
(
"--H0_NeG5_NeC_x2_res"
,
h0_NeG5_NeC_x2_mr
);
(
"--HaPC_NeG4_NeC_1_res"
,
haPC_NeG4_NeC_1_mr
);
(
"--H0_NeG4_NeC_1_res"
,
h0_NeG4_NeC_1_mr
);
(
"--HaPC_NeG8_NeC_4_res"
,
haPC_NeG4_NeC_8_mr
);
(
"--H0_NeG8_NeC_4_res"
,
h0_NeG4_NeC_8_mr
);
(
"--HaPC_NeG4_NeC_8_res"
,
haPC_NeG4_NeC_8_mr
);
(
"--H0_NeG4_NeC_8_res"
,
h0_NeG4_NeC_8_mr
);
(
"--HaPC_NeG1_NeC_4_res"
,
haPC_NeG1_NeC_4_mr
);
(
"--H0_NeG1_NeC_4_res"
,
h0_NeG1_NeC_4_mr
);
(
"--HaPC_NeG16_NeC_1_res"
,
haPC_NeG16_NeC_1_mr
);
(
"--H0_NeG16_NeC_1_res"
,
h0_NeG16_NeC_1_mr
);
(
"--HaPC_NeG1_NeC_16_res"
,
haPC_NeG1_NeC_16_mr
);
(
"--H0_NeG1_NeC_16_res"
,
h0_NeG1_NeC_16_mr
);
(
"--HaPC_NeG8_NeC_1_res"
,
haPC_NeG8_NeC_1_mr
);
(
"--H0_NeG8_NeC_1_res"
,
h0_NeG8_NeC_1_mr
);
(
"--HaPC_NeG1_NeC_8_res"
,
haPC_NeG1_NeC_8_mr
);
(
"--H0_NeG1_NeC_8_res"
,
h0_NeG1_NeC_8_mr
);
(
"--HaPC_NeG5_NeC_1_res"
,
haPC_NeG5_NeC_1_mr
);
(
"--H0_NeG5_NeC_1_res"
,
h0_NeG5_NeC_1_mr
);
(
"--HaPC_NeG1_NeC_5_res"
,
haPC_NeG1_NeC_5_mr
);
...
...
@@ -233,6 +255,7 @@ type res_all_hyp = {
h0_NeG5_res
:
dataset_res
option
;
h0_NeG6_res
:
dataset_res
option
;
h0_NeG8_res
:
dataset_res
option
;
h0_NeG16_res
:
dataset_res
option
;
haPC_NeG1_res
:
dataset_res
option
;
haPC_NeG2_res
:
dataset_res
option
;
haPC_NeG3_res
:
dataset_res
option
;
...
...
@@ -240,18 +263,32 @@ type res_all_hyp = {
haPC_NeG5_res
:
dataset_res
option
;
haPC_NeG6_res
:
dataset_res
option
;
haPC_NeG8_res
:
dataset_res
option
;
haPC_NeG16_res
:
dataset_res
option
;
haPC_NeG5_NeC_div2_res
:
dataset_res
option
;
haPC_NeG5_NeC_x2_res
:
dataset_res
option
;
h0_NeG5_NeC_div2_res
:
dataset_res
option
;
h0_NeG5_NeC_x2_res
:
dataset_res
option
;
haPC_NeG4_NeC_1_res
:
dataset_res
option
;
haPC_NeG4_NeC_8_res
:
dataset_res
option
;
haPC_NeG1_NeC_4_res
:
dataset_res
option
;
h0_NeG4_NeC_1_res
:
dataset_res
option
;
h0_NeG4_NeC_8_res
:
dataset_res
option
;
h0_NeG1_NeC_4_res
:
dataset_res
option
;
h0_NeG4_NeC_8_res
:
dataset_res
option
;
haPC_NeG4_NeC_8_res
:
dataset_res
option
;
h0_NeG8_NeC_4_res
:
dataset_res
option
;
haPC_NeG8_NeC_4_res
:
dataset_res
option
;
haPC_NeG16_NeC_1_res
:
dataset_res
option
;
haPC_NeG1_NeC_16_res
:
dataset_res
option
;
h0_NeG16_NeC_1_res
:
dataset_res
option
;
h0_NeG1_NeC_16_res
:
dataset_res
option
;
haPC_NeG1_NeC_8_res
:
dataset_res
option
;
haPC_NeG8_NeC_1_res
:
dataset_res
option
;
h0_NeG1_NeC_8_res
:
dataset_res
option
;
h0_NeG8_NeC_1_res
:
dataset_res
option
;
haPC_NeG5_NeC_1_res
:
dataset_res
option
;
haPC_NeG1_NeC_5_res
:
dataset_res
option
;
h0_NeG5_NeC_1_res
:
dataset_res
option
;
...
...
@@ -269,24 +306,24 @@ type res_all_hyp = {
}
let
make_t_choices_per_couple
~
tree_prefix
{
ha_PCOC_res
;
h0_NeG1_res
;
h0_NeG2_res
;
h0_NeG3_res
;
h0_NeG4_res
;
h0_NeG5_res
;
h0_NeG6_res
;
h0_NeG8_res
;
h0_NeG1_res
;
h0_NeG2_res
;
h0_NeG3_res
;
h0_NeG4_res
;
h0_NeG5_res
;
h0_NeG6_res
;
h0_NeG8_res
;
h0_NeG16_res
;
haPC_NeG1_res
;
haPC_NeG2_res
;
haPC_NeG3_res
;
haPC_NeG4_res
;
haPC_NeG5_res
;
haPC_NeG6_res
;
haPC_NeG8_res
;
haPC_NeG5_res
;
haPC_NeG6_res
;
haPC_NeG8_res
;
haPC_NeG16_res
;
haPC_NeG5_NeC_div2_res
;
haPC_NeG5_NeC_x2_res
;
h0_NeG5_NeC_div2_res
;
h0_NeG5_NeC_x2_res
;
haPC_NeG4_NeC_1_res
;
haPC_NeG4_NeC_8_res
;
haPC_NeG8_NeC_4_res
;
haPC_NeG1_NeC_4_res
;
h0_NeG4_NeC_1_res
;
h0_NeG4_NeC_8_res
;
h0_NeG8_NeC_4_res
;
h0_NeG1_NeC_4_res
;
haPC_NeG5_NeC_1_res
;
haPC_NeG1_NeC_5_res
;
h0_NeG4_NeC_6_res
;
h0_NeG4_NeC_2_res
;
h0_NeG2_NeC_4_res
;
...
...
@@ -294,7 +331,17 @@ let make_t_choices_per_couple ~tree_prefix {ha_PCOC_res;
haPC_NeG4_NeC_2_res
;
haPC_NeG2_NeC_4_res
;
h0_NeG5_NeC_1_res
;
h0_NeG1_NeC_5_res
;
h0_NeG1_NeC_5_res
;
haPC_NeG5_NeC_1_res
;
haPC_NeG1_NeC_5_res
;
h0_NeG16_NeC_1_res
;
h0_NeG1_NeC_16_res
;
haPC_NeG16_NeC_1_res
;
haPC_NeG1_NeC_16_res
;
h0_NeG8_NeC_1_res
;
h0_NeG1_NeC_8_res
;
haPC_NeG8_NeC_1_res
;
haPC_NeG1_NeC_8_res
;
h0_NeG5_indel_res
;
haPC_NeG5_indel_res
;
}
=
...
...
@@ -308,6 +355,7 @@ let make_t_choices_per_couple ~tree_prefix {ha_PCOC_res;
let
h0_NeG5_mr
=
get_merged_results_opt
h0_NeG5_res
in
let
h0_NeG6_mr
=
get_merged_results_opt
h0_NeG6_res
in
let
h0_NeG8_mr
=
get_merged_results_opt
h0_NeG8_res
in
let
h0_NeG16_mr
=
get_merged_results_opt
h0_NeG16_res
in
let
haPC_NeG1_mr
=
get_merged_results_opt
haPC_NeG1_res
in
let
haPC_NeG2_mr
=
get_merged_results_opt
haPC_NeG2_res
in
let
haPC_NeG3_mr
=
get_merged_results_opt
haPC_NeG3_res
in
...
...
@@ -315,6 +363,7 @@ let make_t_choices_per_couple ~tree_prefix {ha_PCOC_res;
let
haPC_NeG5_mr
=
get_merged_results_opt
haPC_NeG5_res
in
let
haPC_NeG6_mr
=
get_merged_results_opt
haPC_NeG6_res
in
let
haPC_NeG8_mr
=
get_merged_results_opt
haPC_NeG8_res
in
let
haPC_NeG16_mr
=
get_merged_results_opt
haPC_NeG16_res
in
let
haPC_NeG5_NeC_div2_mr
=
get_merged_results_opt
haPC_NeG5_NeC_div2_res
in
let
haPC_NeG5_NeC_x2_mr
=
get_merged_results_opt
haPC_NeG5_NeC_x2_res
in
...
...
@@ -323,15 +372,27 @@ let make_t_choices_per_couple ~tree_prefix {ha_PCOC_res;
let
haPC_NeG4_NeC_1_mr
=
get_merged_results_opt
haPC_NeG4_NeC_1_res
in
let
haPC_NeG4_NeC_8_mr
=
get_merged_results_opt
haPC_NeG4_NeC_8_res
in
let
haPC_NeG8_NeC_4_mr
=
get_merged_results_opt
haPC_NeG8_NeC_4_res
in
let
haPC_NeG1_NeC_4_mr
=
get_merged_results_opt
haPC_NeG1_NeC_4_res
in
let
h0_NeG1_NeC_4_mr
=
get_merged_results_opt
h0_NeG1_NeC_4_res
in
let
h0_NeG4_NeC_1_mr
=
get_merged_results_opt
h0_NeG4_NeC_1_res
in
let
h0_NeG4_NeC_8_mr
=
get_merged_results_opt
h0_NeG4_NeC_8_res
in
let
h0_NeG8_NeC_4_mr
=
get_merged_results_opt
h0_NeG8_NeC_4_res
in
let
haPC_NeG5_NeC_1_mr
=
get_merged_results_opt
haPC_NeG5_NeC_1_res
in
let
haPC_NeG1_NeC_5_mr
=
get_merged_results_opt
haPC_NeG1_NeC_5_res
in
let
h0_NeG1_NeC_5_mr
=
get_merged_results_opt
h0_NeG1_NeC_5_res
in
let
h0_NeG5_NeC_1_mr
=
get_merged_results_opt
h0_NeG5_NeC_1_res
in
let
haPC_NeG8_NeC_1_mr
=
get_merged_results_opt
haPC_NeG8_NeC_1_res
in
let
haPC_NeG1_NeC_8_mr
=
get_merged_results_opt
haPC_NeG1_NeC_8_res
in
let
h0_NeG1_NeC_8_mr
=
get_merged_results_opt
h0_NeG1_NeC_8_res
in
let
h0_NeG8_NeC_1_mr
=
get_merged_results_opt
h0_NeG8_NeC_1_res
in
let
haPC_NeG16_NeC_1_mr
=
get_merged_results_opt
haPC_NeG16_NeC_1_res
in
let
haPC_NeG1_NeC_16_mr
=
get_merged_results_opt
haPC_NeG1_NeC_16_res
in
let
h0_NeG1_NeC_16_mr
=
get_merged_results_opt
h0_NeG1_NeC_16_res
in
let
h0_NeG16_NeC_1_mr
=
get_merged_results_opt
h0_NeG16_NeC_1_res
in
let
haPC_NeG4_NeC_6_mr
=
get_merged_results_opt
haPC_NeG4_NeC_6_res
in
let
haPC_NeG4_NeC_2_mr
=
get_merged_results_opt
haPC_NeG4_NeC_2_res
in
...
...
@@ -351,6 +412,7 @@ let make_t_choices_per_couple ~tree_prefix {ha_PCOC_res;
~
h0_NeG5_mr
~
h0_NeG6_mr
~
h0_NeG8_mr
~
h0_NeG16_mr
~
haPC_NeG1_mr
~
haPC_NeG2_mr
~
haPC_NeG3_mr
...
...
@@ -358,6 +420,7 @@ let make_t_choices_per_couple ~tree_prefix {ha_PCOC_res;
~
haPC_NeG5_mr
~
haPC_NeG6_mr
~
haPC_NeG8_mr
~
haPC_NeG16_mr
~
haPC_NeG5_NeC_div2_mr
~
haPC_NeG5_NeC_x2_mr
~
h0_NeG5_NeC_div2_mr
...
...
@@ -366,12 +429,12 @@ let make_t_choices_per_couple ~tree_prefix {ha_PCOC_res;
~
haPC_NeG5_indel_mr
~
haPC_NeG4_NeC_1_mr
~
haPC_NeG4_NeC_8_mr
~
haPC_NeG8_NeC_4_mr
~
haPC_NeG1_NeC_4_mr
~
h0_NeG1_NeC_4_mr
~
h0_NeG4_NeC_8_mr
~
h0_NeG8_NeC_4_mr
~
h0_NeG4_NeC_1_mr
~
haPC_NeG5_NeC_1_mr
~
haPC_NeG1_NeC_5_mr
~
h0_NeG4_NeC_6_mr
~
h0_NeG4_NeC_2_mr
~
h0_NeG2_NeC_4_mr
...
...
@@ -380,6 +443,16 @@ let make_t_choices_per_couple ~tree_prefix {ha_PCOC_res;
~
haPC_NeG2_NeC_4_mr
~
h0_NeG1_NeC_5_mr
~
h0_NeG5_NeC_1_mr
~
haPC_NeG5_NeC_1_mr
~
haPC_NeG1_NeC_5_mr
~
h0_NeG1_NeC_8_mr
~
h0_NeG8_NeC_1_mr
~
haPC_NeG8_NeC_1_mr
~
haPC_NeG1_NeC_8_mr
~
h0_NeG1_NeC_16_mr
~
h0_NeG16_NeC_1_mr
~
haPC_NeG16_NeC_1_mr
~
haPC_NeG1_NeC_16_mr
()
...
...
@@ -394,6 +467,7 @@ let get_t_choices ~tree_prefix ~(dataset_results_l: dataset_res list) : t_choice
let
h0_NeG5_res
=
List
.
find
dataset_results_l
(
is_hyp
~
hyp
:
"H0_NeG5"
)
in
let
h0_NeG6_res
=
List
.
find
dataset_results_l
(
is_hyp
~
hyp
:
"H0_NeG6"
)
in
let
h0_NeG8_res
=
List
.
find
dataset_results_l
(
is_hyp
~
hyp
:
"H0_NeG8"
)
in
let
h0_NeG16_res
=
List
.
find
dataset_results_l
(
is_hyp
~
hyp
:
"H0_NeG16"
)
in
let
haPC_NeG1_res
=
List
.
find
dataset_results_l
(
is_hyp
~
hyp
:
"HaPC_NeG1"
)
in
let
haPC_NeG2_res
=
List
.
find
dataset_results_l
(
is_hyp
~
hyp
:
"HaPC_NeG2"
)
in
let
haPC_NeG3_res
=
List
.
find
dataset_results_l
(
is_hyp
~
hyp
:
"HaPC_NeG3"
)
in
...
...
@@ -401,6 +475,7 @@ let get_t_choices ~tree_prefix ~(dataset_results_l: dataset_res list) : t_choice
let
haPC_NeG5_res
=
List
.
find
dataset_results_l
(
is_hyp
~
hyp
:
"HaPC_NeG5"
)
in
let
haPC_NeG6_res
=
List
.
find
dataset_results_l
(
is_hyp
~
hyp
:
"HaPC_NeG6"
)
in
let
haPC_NeG8_res
=
List
.
find
dataset_results_l
(
is_hyp
~
hyp
:
"HaPC_NeG8"
)
in
let
haPC_NeG16_res
=
List
.
find
dataset_results_l
(
is_hyp
~
hyp
:
"HaPC_NeG16"
)
in
let
haPC_NeG5_NeC_div2_res
=
List
.
find
dataset_results_l
(
is_hyp
~
hyp
:
"HaPC_NeG5_NeC_div2"
)
in
let
haPC_NeG5_NeC_x2_res
=
List
.
find
dataset_results_l
(
is_hyp
~
hyp
:
"HaPC_NeG5_NeC_x2"
)
in
...
...
@@ -408,16 +483,29 @@ let get_t_choices ~tree_prefix ~(dataset_results_l: dataset_res list) : t_choice
let
h0_NeG5_NeC_x2_res
=
List
.
find
dataset_results_l
(
is_hyp
~
hyp
:
"H0_NeG5_NeC_x2"
)
in
let
haPC_NeG4_NeC_1_res
=
List
.
find
dataset_results_l
(
is_hyp
~
hyp
:
"HaPC_NeG4_NeC_1"
)
in
let
haPC_NeG4_NeC_8_res
=
List
.
find
dataset_results_l
(
is_hyp
~
hyp
:
"HaPC_NeG4_NeC_8"
)
in
let
haPC_NeG1_NeC_4_res
=
List
.
find
dataset_results_l
(
is_hyp
~
hyp
:
"HaPC_NeG1_NeC_4"
)
in
let
h0_NeG4_NeC_1_res
=
List
.
find
dataset_results_l
(
is_hyp
~
hyp
:
"H0_NeG4_NeC_1"
)
in
let
h0_NeG4_NeC_8_res
=
List
.
find
dataset_results_l
(
is_hyp
~
hyp
:
"H0_NeG4_NeC_8"
)
in
let
h0_NeG1_NeC_4_res
=
List
.
find
dataset_results_l
(
is_hyp
~
hyp
:
"H0_NeG1_NeC_4"
)
in
let
h0_NeG4_NeC_8_res
=
List
.
find
dataset_results_l
(
is_hyp
~
hyp
:
"H0_NeG4_NeC_8"
)
in
let
haPC_NeG4_NeC_8_res
=
List
.
find
dataset_results_l
(
is_hyp
~
hyp
:
"HaPC_NeG4_NeC_8"
)
in
let
h0_NeG8_NeC_4_res
=
List
.
find
dataset_results_l
(
is_hyp
~
hyp
:
"H0_NeG8_NeC_4"
)
in
let
haPC_NeG8_NeC_4_res
=
List
.
find
dataset_results_l
(
is_hyp
~
hyp
:
"HaPC_NeG8_NeC_4"
)
in
let
haPC_NeG5_NeC_1_res
=
List
.
find
dataset_results_l
(
is_hyp
~
hyp
:
"HaPC_NeG5_NeC_1"
)
in
let
haPC_NeG1_NeC_5_res
=
List
.
find
dataset_results_l
(
is_hyp
~
hyp
:
"HaPC_NeG1_NeC_5"
)
in
let
h0_NeG5_NeC_1_res
=
List
.
find
dataset_results_l
(
is_hyp
~
hyp
:
"H0_NeG5_NeC_1"
)
in
let
h0_NeG1_NeC_5_res
=
List
.
find
dataset_results_l
(
is_hyp
~
hyp
:
"H0_NeG1_NeC_5"
)
in
let
haPC_NeG8_NeC_1_res
=
List
.
find
dataset_results_l
(
is_hyp
~
hyp
:
"HaPC_NeG8_NeC_1"
)
in
let
haPC_NeG1_NeC_8_res
=
List
.
find
dataset_results_l
(
is_hyp
~
hyp
:
"HaPC_NeG1_NeC_8"
)
in
let
h0_NeG8_NeC_1_res
=
List
.
find
dataset_results_l
(
is_hyp
~
hyp
:
"H0_NeG8_NeC_1"
)
in
let
h0_NeG1_NeC_8_res
=
List
.
find
dataset_results_l
(
is_hyp
~
hyp
:
"H0_NeG1_NeC_8"
)
in
let
haPC_NeG16_NeC_1_res
=
List
.
find
dataset_results_l
(
is_hyp
~
hyp
:
"HaPC_NeG16_NeC_1"
)
in
let
haPC_NeG1_NeC_16_res
=
List
.
find
dataset_results_l
(
is_hyp
~
hyp
:
"HaPC_NeG1_NeC_16"
)
in
let
h0_NeG16_NeC_1_res
=
List
.
find
dataset_results_l
(
is_hyp
~
hyp
:
"H0_NeG16_NeC_1"
)
in
let
h0_NeG1_NeC_16_res
=
List
.
find
dataset_results_l
(
is_hyp
~
hyp
:
"H0_NeG1_NeC_16"
)
in
let
haPC_NeG4_NeC_6_res
=
List
.
find
dataset_results_l
(
is_hyp
~
hyp
:
"HaPC_NeG4_NeC_6"
)
in
let
haPC_NeG4_NeC_2_res
=
List
.
find
dataset_results_l
(
is_hyp
~
hyp
:
"HaPC_NeG4_NeC_2"
)
in
...
...
@@ -433,11 +521,14 @@ let get_t_choices ~tree_prefix ~(dataset_results_l: dataset_res list) : t_choice
match
(
h0_NeG4_res
)
with
|
(
Some
h0_NeG4
)
->
let
t_choices_dir
=
make_t_choices_per_couple
~
tree_prefix
{
ha_PCOC_res
;
h0_NeG1_res
;
h0_NeG2_res
;
h0_NeG3_res
;
h0_NeG4_res
;
h0_NeG5_res
;
h0_NeG6_res
;
h0_NeG8_res
;
haPC_NeG1_res
;
haPC_NeG2_res
;
haPC_NeG3_res
;
haPC_NeG4_res
;
haPC_NeG5_res
;
haPC_NeG6_res
;
haPC_NeG8_res
;
h0_NeG1_res
;
h0_NeG2_res
;
h0_NeG3_res
;
h0_NeG4_res
;
h0_NeG5_res
;
h0_NeG6_res
;
h0_NeG8_res
;
h0_NeG16_res
;
haPC_NeG1_res
;
haPC_NeG2_res
;
haPC_NeG3_res
;
haPC_NeG4_res
;
haPC_NeG5_res
;
haPC_NeG6_res
;
haPC_NeG8_res
;
haPC_NeG16_res
;
haPC_NeG5_NeC_div2_res
;
haPC_NeG5_NeC_x2_res
;
h0_NeG5_NeC_div2_res
;
h0_NeG5_NeC_x2_res
;
haPC_NeG4_NeC_1_res
;
haPC_NeG4_NeC_8_res
;
haPC_NeG1_NeC_4_res
;
h0_NeG1_NeC_4_res
;
h0_NeG4_NeC_1_res
;
h0_NeG4_NeC_8_res
;
haPC_NeG4_NeC_1_res
;
haPC_NeG1_NeC_4_res
;
h0_NeG1_NeC_4_res
;
h0_NeG4_NeC_1_res
;
haPC_NeG4_NeC_8_res
;
h0_NeG4_NeC_8_res
;
haPC_NeG8_NeC_4_res
;
h0_NeG8_NeC_4_res
;
haPC_NeG5_NeC_1_res
;
haPC_NeG1_NeC_5_res
;
h0_NeG1_NeC_5_res
;
h0_NeG5_NeC_1_res
;
haPC_NeG8_NeC_1_res
;
haPC_NeG1_NeC_8_res
;
h0_NeG1_NeC_8_res
;
h0_NeG8_NeC_1_res
;
haPC_NeG16_NeC_1_res
;
haPC_NeG1_NeC_16_res
;
h0_NeG1_NeC_16_res
;
h0_NeG16_NeC_1_res
;
haPC_NeG4_NeC_6_res
;
haPC_NeG4_NeC_2_res
;
haPC_NeG2_NeC_4_res
;
h0_NeG4_NeC_6_res
;
h0_NeG4_NeC_2_res
;
h0_NeG2_NeC_4_res
;
h0_NeG5_indel_res
;
haPC_NeG5_indel_res
}
in
let
t_choices_max
=
t_choices_dir
/
selector
[
"out.max_MCC_per_meth.tsv"
]
in
...
...
lib/scripts/calc_t_per_meth.R
View file @
bda12e50
This diff is collapsed.
Click to expand it.
Write
Preview
Markdown
is supported
0%
Try again
or
attach a new file
.
Attach a file
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Cancel
Please
register
or
sign in
to comment