Open h-pawar opened 3 months ago
Hello, @h-pawar! Sorry we didn't notice this earlier. I can't dig in deep here right now, but if this is still an issue:
Would you expect this model to take such a long time to run (~6 days)? Or have I specified the introgression pulse incorrectly, does this need a separate call to be added?
Well, how long does the model without adaptive introgression take to run? If it is much shorter than 6 days, probably what you are seeing is that the events you're conditioning on are very unlikely, so it is restarting a lot. I guess the way I would debug this is to start simple - so, remove some of the conditions and see if that helps? I would also do this debugging using a very strong rescaling factor, so it runs in minutes not days.
Is there a way to check that the mutation has introgressed & correctly undergone selection in the recipient population?
The short answer is: no, not an easy way? We do have pretty rigorous checks that these things function as they should, though. You could come up with ad-hoc ways like turn the selection coefficient way up and check that this increases the ending frequency?
Short description of the problem:
Hello, building on the orangutan “TwoSpecies_2L11” demographic model, I would like to add adaptive introgression (from bornean into sumatran). Mutation arises any time from the species split to 1000 generations ago, passes to sumatran & is beneficial with s=0.1.
Following the tutorial - https://github.com/popsim-consortium/stdpopsim/blob/main/docs/selection_example.py, I
Would you expect this model to take such a long time to run (~6 days)? Or have I specified the introgression pulse incorrectly, does this need a separate call to be added?
Is there a way to check that the mutation has introgressed & correctly undergone selection in the recipient population?
How to reproduce the problem: