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

Loci length limit? #207

Open
ckastall opened this issue Mar 4, 2020 · 0 comments
Open

Loci length limit? #207

ckastall opened this issue Mar 4, 2020 · 0 comments

Comments

@ckastall
Copy link

ckastall commented Mar 4, 2020

Thank you for this very useful package.

Trying to simulate long sequences I have the following error message:

The sum of the partition lengths in the treefile does not equal
the specified number of sites.

This seems to happen only if I define a model with loci_length > 99998.

Here is the exact commands I am running:

# No problem with this model:
model <- coal_model(c(20, 1), 2, loci_length = 99998) +
        feat_outgroup(2) +
        feat_pop_merge(0.5, 2, 1) +
        feat_mutation(5, model = "GTR", gtr_rates = 1:6)

res <- simulate(model)

# This one won't run:
model <- coal_model(c(20, 1), 2, loci_length = 99999) +
        feat_outgroup(2) +
        feat_pop_merge(0.5, 2, 1) +
        feat_mutation(5, model = "GTR", gtr_rates = 1:6)

res <- simulate(model)

I have tried with both simulators scrm and ms, the problem seems to be the same.

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