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

about RCCR value > 1 #62

Open
cnxiaobo94 opened this issue Jul 15, 2024 · 2 comments
Open

about RCCR value > 1 #62

cnxiaobo94 opened this issue Jul 15, 2024 · 2 comments

Comments

@cnxiaobo94
Copy link

Dear @stschiff,

I have a dataset that consists of three populations as (pop1, (pop2, pop3)). After running msmc2, I have three paired RCCR values, but I met some problems when plotting them. As the picture shows below, some values of time points are out of rates 0-1. I can observe the divergence times of pop1&2 (orange line) and pop1&3 (green line), but the divergence time (blue line) between pop2 and pop3 is hard to know.
I'd appreciate it if you could help with this problem.
rccr_1

Best,
Bo Xiao

@stschiff
Copy link
Owner

Looks to me like pop2 and pop3 are hardly diverged at all. Have you checked with other statistics, like FST, or looked at whether the population sizes start to differ somewhere?

@cnxiaobo94
Copy link
Author

Looks to me like pop2 and pop3 are hardly diverged at all. Have you checked with other statistics, like FST, or looked at whether the population sizes start to differ somewhere?
Thanks for your reply! Sorry for this denied message. I checked the popo]ulations siez of both MSMC and PSMC (as below). The divergence time of pop2 and pop3 is calculated as between 100-200 ka. The weighted Fst estimate of pop2&3 is 0.25, while that of pop1&2 (pop1&3) is 0.56.
image
image

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