msms icon indicating copy to clipboard operation
msms copied to clipboard

-Sc argument not working as expected

Open ornobalam opened this issue 2 years ago • 0 comments

Hi!

I would expect the following pieces of code to be equivalent:

A: -SI 0.01 1 0.001 -Sc 0 1 400 400 0

B: -SI 0.01 1 0.001 -SAA 400 -SaA 400

But A produces a frequency trajectory (across replicates) that starts rising right after t=0.01, while B only rises very recently. Am I understanding the -Sc option incorrectly?

To provide some more context, my goal is to simulate selection on standing variation but in a way that the previously neutral variation extends into the past and I can retrieve their frequencies. I thought I could do this by setting the selection coefficient to zero using -Sc before a certain time.

Thanks! Ornob

ornobalam avatar Nov 14 '22 16:11 ornobalam