Skip to content

Commit

Permalink
[cattle-outbreak] add genome clock filter
Browse files Browse the repository at this point in the history
This should help remove outliers unrelated to the cattle-outbreak,
especially as the clock signal is worse than it might otherwise be for
these outliers due to our enforcing of the root strain.

See related discussion on slack <https://bedfordlab.slack.com/archives/CD84ELG0N/p1733773815454369>
  • Loading branch information
jameshadfield committed Dec 9, 2024
1 parent 2fac740 commit 871a0de
Show file tree
Hide file tree
Showing 2 changed files with 5 additions and 1 deletion.
4 changes: 3 additions & 1 deletion Snakefile
Original file line number Diff line number Diff line change
Expand Up @@ -203,7 +203,9 @@ def refine_clock_rates(w):
return f"--clock-rate {info[w.segment][0]} --clock-std-dev {info[w.segment][1]}"

def refine_clock_filter(w):
filter = get_config('refine', 'clock_filter_iqd', w)
filter = get_config('refine', 'genome_clock_filter_iqd', w) \
if w.segment=='genome' \
else get_config('refine', 'clock_filter_iqd', w)
return f"--clock-filter-iqd {filter}" if filter else ""


Expand Down
2 changes: 2 additions & 0 deletions config/h5n1-cattle-outbreak.yaml
Original file line number Diff line number Diff line change
Expand Up @@ -75,6 +75,8 @@ refine:
coalescent: const
date_inference: marginal

genome_clock_filter_iqd:
FALLBACK: 3
clock_filter_iqd:
FALLBACK: false

Expand Down

0 comments on commit 871a0de

Please sign in to comment.