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

Incorporate inferred tree sequences to the software #30

Open
3 tasks done
shz9 opened this issue Jul 22, 2020 · 4 comments
Open
3 tasks done

Incorporate inferred tree sequences to the software #30

shz9 opened this issue Jul 22, 2020 · 4 comments
Assignees
Labels
enhancement New feature or request

Comments

@shz9
Copy link
Collaborator

shz9 commented Jul 22, 2020

The ultimate goal of this project is to align inferred tree sequences to the pedigree graph. Therefore, we need to create tools to infer tree sequences from genetic samples.

I talked to @kobrica and she will work on providing this functionality. This task is composed of 3 steps.

  • Given genetic samples (e.g. genotype matrix), infer the corresponding tree sequences using tsinfer (tutorial).
  • Convert the inferred tree sequences to Traversal objects.
  • Given msprime simulation results, obtain the corresponding inferred tree sequences. This will involve 2 steps: First, obtain the genetic samples from the msprime simulation, and then run tsinfer on those samples, as in the tutorial above.
@shz9 shz9 added the enhancement New feature or request label Jul 22, 2020
@andjelatodorovic
Copy link
Collaborator

I have completed all of the above tasks, and added some samples of genotype matrices in csv format in data folder.
This is easier than manually adding row by row to sample data.

The tree sequence to Traversal may need to be double-checked, I'm not quite sure the Traversal object has all of the properties. Please refer to commit 66a2702

@shz9
Copy link
Collaborator Author

shz9 commented Jul 28, 2020

Can you please modify the infer_from_msprime function so that it takes a simulation object (like the one we get from generate_msprime_simulations)? In other words, you can assume you already have the result of the simulation.

@shz9 shz9 reopened this Jul 28, 2020
@andjelatodorovic
Copy link
Collaborator

@shz9 Sure! Should it take both simulation results and traversal object as parameters?

@shz9
Copy link
Collaborator Author

shz9 commented Jul 28, 2020

I think just the simulation results are needed.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
enhancement New feature or request
Projects
None yet
Development

No branches or pull requests

2 participants