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

Improve ancestor building by using more information #911

Open
hyanwong opened this issue Mar 21, 2024 · 0 comments
Open

Improve ancestor building by using more information #911

hyanwong opened this issue Mar 21, 2024 · 0 comments

Comments

@hyanwong
Copy link
Member

It looks like #903 is being caused by bad choices of ancestral vs derived states in the ancestor builder. We do have more information for building, e.g. we know after the first inference (or if we use the PBWT) whether we think that many of the haplotypes contributing to the consensus form a cluster of recent close relatives.

We should think about how to use this information to improve ancestor reconstruction. @petrelharp suggested building a "decision tree" rightwards and leftwards from each site, branching whenever we make a choice of 0 or 1 for an adjacent site, and picking e.g. the longest path. He pointed out that this can be (mostly) kept between ancestors, so should be efficient.

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

No branches or pull requests

1 participant