Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

2y build targets #14

Merged
merged 12 commits into from
Apr 23, 2024
Merged

2y build targets #14

merged 12 commits into from
Apr 23, 2024

Conversation

trvrb
Copy link
Member

@trvrb trvrb commented Apr 23, 2024

This is mostly @lmoncla's work to implement a new 2y vs all-time wildcard to specify the timespan for individual build targets. This allows a much larger fraction of recent H5N1 sequences to make it into the final subsampled analysis.

I had run a version of this AWS using the h5nx-2-year-test branch. The results of this run are available on staging with URLS like:

Everything here seemed to be in order. However, I couldn't just merge h5nx-2-year-test into master as there were a conflicts due to PR #12 and other work on master. To compensate I manually updated this PR branch 2y-build-targets to bring it up to speed with master in commit 7549b7d.

I tested merging 2y-build-targets into master locally and it went smoothly with just a copy small resolutions during the merge.

@lmoncla: if this looks good to you, I'm happy to help with the merge.

lhmoncla and others added 12 commits April 3, 2024 17:36
…ut the clock rate and clock rate standard deviation for the 2 year build
…ild and an all-time build. This commit introduces this new wildcard structure, as well as new functions for assigning strict clock rates for the 3 year trees
This switches from the older Augur subsampling strategy of --sequences-per-group to the newer --subsampling-max-sequences. By doing this, we can just specific once the target sequence count and don't have to play around with logic for different subtypes / time-spans. This also makes it easier to rapidly test entire set of builds by locally swapping from "3000" to "100" in just one line.
I'd suggest using "3y" or "2y" in the URL wildcard rather than "3-year" or "2-year" as this will be consistent with existing timespan wildcards for seasonal flu (6m, 2y, 3y, 12y) and for ncov (1m, 2m, 6m, all-time).

I'd suggest doing a 2 year target rather than a 3 year target just because there seems to be well sufficient sequence counts in GISAID to pull >3000 genomes for a 2 year span for H5N1. If I look now, I see 5571 H5N1 viruses collected from 2022-04-22 to today.

I moved the wildcard dictionary parameters from being all in one-line to being separated by subtype to improve readability.
This commit pulls in changes on master branch that were not present on 2y branch
@trvrb trvrb requested a review from lmoncla April 23, 2024 01:59
@lmoncla
Copy link
Collaborator

lmoncla commented Apr 23, 2024

Thank you @trvrb for working on this! This looks good to me. Please go ahead and merge!

@trvrb trvrb merged commit 6cc834e into master Apr 23, 2024
6 checks passed
@trvrb trvrb deleted the 2y-build-targets branch April 23, 2024 17:50
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

2 participants