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

Implement tree sequence #126

Open
gregorgorjanc opened this issue Apr 1, 2023 · 4 comments
Open

Implement tree sequence #126

gregorgorjanc opened this issue Apr 1, 2023 · 4 comments

Comments

@gregorgorjanc
Copy link
Contributor

gregorgorjanc commented Apr 1, 2023

This is a big one and a project on its own, but here are some initial ideas

  • use tskit C library to store tree sequence?
  • import tree sequence, split some sites as QTL and SNP array, the rest as tree sequence sites ...
  • licencing issues? AlphaSimR is MIT, tskitC is GPL
@gregorgorjanc
Copy link
Contributor Author

Peter Ralph directed me to this implementation doc for SLiM https://github.com/MesserLab/SLiM/blob/master/treerec/implementation.md

@gaynorr
Copy link
Owner

gaynorr commented Oct 28, 2023

Need to change recombination tracking in AlphaSimR to keep track of all recombinations and not just observed recombinations. This can be accomplished by not processing recombination output from findBivalent and only processing the output within bivalent for the purposes of copying data.

The full history can then be fowarded along to recombination tracking. Processing will be needed in the current recombination history list, because it depends on this data being processed. Perhaps we need to recombination histories, or some adjustments to the current IBD tracking to process as needed.

@gaynorr
Copy link
Owner

gaynorr commented Oct 28, 2023

Mutations need to be tracked for tree sequences. This can be added to SimParam.

@gaynorr
Copy link
Owner

gaynorr commented Nov 1, 2023

https://www.slendr.net/articles/vignette-05-tree-sequences.html

Look at this article for an example implementation in R.

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

2 participants