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

Tracer plots for LDhat run #24

Open
Homap opened this issue Jun 14, 2022 · 1 comment
Open

Tracer plots for LDhat run #24

Homap opened this issue Jun 14, 2022 · 1 comment

Comments

@Homap
Copy link

Homap commented Jun 14, 2022

Hello,

I ran interval on 2000 SNPs of a region about 300 Kb. I have 10 diploid individuals. I ran interval with block penalty of 5, 10 million iterations and sampled every 2000 steps. I use a subset of likelihood table for 20 haplotypes and a theta of 0.001. The pi in my samples is about 0.002 and theta is very close to it.

Screenshot 2022-06-14 at 10 56 14

Screenshot 2022-06-14 at 10 56 21

Screenshot 2022-06-14 at 10 56 25

The effective sample sizes (ESS) are very low for likelihood and blocks and there is a large difference between their ESS and the one of map.

Sorry as this is a very general questions but I wonder if you have any suggestions on how I could improve this? Is it the matter of running the chain longer or some other potential issue? Could the likelihood table be the problem? Is there a difference in the likelihood table when using genotypes vs. haplotypes?

Thank you!
Homa

@Homap
Copy link
Author

Homap commented Jun 14, 2022

Very interesting, I just converted the length of region to Kb, keeping everything else constant. The result changed drastically:

Screenshot 2022-06-14 at 11 44 52

Screenshot 2022-06-14 at 11 45 00

Screenshot 2022-06-14 at 11 45 03

Why would this change the output? It's really interesting.

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