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

rho/theta--> total recombination events? #247

Closed
pneumowidow opened this issue Apr 15, 2019 · 1 comment
Closed

rho/theta--> total recombination events? #247

pneumowidow opened this issue Apr 15, 2019 · 1 comment

Comments

@pneumowidow
Copy link

Hi,

Since rho/theta corresponds to the ratio of the number of recombination events to point mutations on a branch, I just wanted to confirm that in order to calculate the number of recombination events on a branch, we would have to multiply rho/theta with the average no of snps per site from the final tree? For example if my rho/theta for a node is 0.4 and my final tree gives 200 as the average number of substitutions per site (branch length), then the number of recombination events on that branch would be 0.4*200 = 80 ?

Did I get that right or should I use the total snps found on that node (irrespective of whether it arose from recombination or not).

Thank you in advance.

@pneumowidow
Copy link
Author

I finally figured it out. The recombination events are the recombination blocks and the point mutations refers to the SNPs outside of recombination. The division of these 2 values gives the rho/theta value.

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