Skip to content

Commit a685ef9

Browse files
committed
Fixup: Phylogenetic: Add Frequencies
Switch from wide_bandwidth to proportion_wide Document frequencies parameter choices in the config file
1 parent 8cf43ac commit a685ef9

File tree

2 files changed

+7
-5
lines changed

2 files changed

+7
-5
lines changed

phylogenetic/defaults/config.yaml

Lines changed: 5 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -79,10 +79,12 @@ traits:
7979
]
8080

8181
tip_frequencies:
82-
min_date: "1980-01-01"
82+
# 2000 since there is an increase in WNV at that time
83+
min_date: "2000-01-01"
8384
max_date: "6M"
84-
narrow_bandwidth: 0.2
85-
wide_bandwidth: 0.6
85+
# Quarterly narrow_bandwidth or every 3 months (3 /12.0 = 0.25)
86+
narrow_bandwidth: 0.25
87+
proportion_wide: 0.0
8688

8789
export:
8890
description: "defaults/description.md"

phylogenetic/rules/export.smk

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -69,7 +69,7 @@ rule tip_frequencies:
6969
min_date = config["tip_frequencies"]["min_date"],
7070
max_date = config["tip_frequencies"]["max_date"],
7171
narrow_bandwidth = config["tip_frequencies"]["narrow_bandwidth"],
72-
wide_bandwidth = config["tip_frequencies"]["wide_bandwidth"]
72+
proportion_wide = config["tip_frequencies"]["proportion_wide"]
7373
shell:
7474
r"""
7575
augur frequencies \
@@ -80,6 +80,6 @@ rule tip_frequencies:
8080
--min-date {params.min_date} \
8181
--max-date {params.max_date} \
8282
--narrow-bandwidth {params.narrow_bandwidth} \
83-
--wide-bandwidth {params.wide_bandwidth} \
83+
--proportion-wide {params.proportion_wide} \
8484
--output {output.tip_freq}
8585
"""

0 commit comments

Comments
 (0)