forked from phylo42/PEWO
-
Notifications
You must be signed in to change notification settings - Fork 0
/
eval_likelihood.smk
67 lines (59 loc) · 1.88 KB
/
eval_likelihood.smk
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
"""
This is a workflow to evaluate likelihood of extended trees, given the parameters
set in "config.yaml". This snakefile loads all necessary modules and builds
the evaluation workflow itself based on the setup defined in the config file.
"""
__author__ = "Nikolai Romashchenko"
__license__ = "MIT"
configfile: "config.yaml"
config["mode"] = "likelihood"
# Explicitly set config as if there was a single pruning which in fact represents the full (NOT pruned) tree.
# This allow to use the same config file for both 'accuracy' and 'resources' modes of PEWO worflow
# NOTE: this statement MUST be set BEFORE the "includes"
config["pruning_count"] = 1
config["read_length"] = [0]
# Explicitly set config to not repeat binary executions,
# which is an option that should be considered only in 'resource' evaluation mode.
# this allow to use the same config file for both 'accuracy' and 'resources' modes of PEWO worflow
# NOTE: this statement MUST be set BEFORE the "includes"
config["repeats"] = 1
# utils
include:
"rules/utils/workflow.smk"
include:
"rules/utils/etc.smk"
# Tree prunings
include:
"rules/op/operate_prunings.smk"
# Tree optimisation
include:
"rules/op/operate_optimisation.smk"
# phylo-kmer placement, e.g.: rappas
include:
"rules/op/ar.smk"
include:
"rules/placement/rappas.smk"
#alignment (for distance-based and ML approaches)
include:
"rules/alignment/hmmer.smk"
# ML-based placements, e.g.: epa, epang, pplacer
include:
"rules/placement/epa.smk"
include:
"rules/placement/pplacer.smk"
include:
"rules/placement/epang.smk"
# Distance-based placements, e.g.: apples
include:
"rules/placement/apples.smk"
# Results evaluation and plots
include:
"rules/op/operate_likelihood.smk"
include:
"rules/op/operate_plots.smk"
rule all:
"""
top snakemake rule, necessary to launch the workflow
"""
input:
build_likelihood_workflow()