Sample size of output population

[Zuxi Cui]

Hi there,

I have a quick question about non-WF simulation. I followed 16.12 in the manual and did several

simulations with the setting "defineConstant("K", 200);" in the script of creating mating.txt and

"g = p1.genomes;g.outputVCF("./GBR/sim1.vcf");" in the script of output vcf files.

However, the sample size of simulated outputs varied a little among my parallel simulations from ~180 to ~215. Can someone help explain the reason behind this?

Thanks,

Terry 

[Ben Haller]

Hi Terry.  In nonWF models population size is often an emergent property of the model, not a top-down parameter, and can vary stochastically around the carrying capacity (just as it does in real biological populations), because it is often modeled as a consequence of the balance between stochastic birth and death processes.  See, e.g., sections 1.6 and 16.1.  If you want your nonWF model to behave more like a WF model, with a strictly constant population size, you can do that; all of this is simply up to the script that you write, not in any way hard-coded into SLiM itself.  Section 16.15 has examples and discussion.

Cheers,
-B.

Benjamin C. Haller
Messer Lab
Cornell University


Zuxi Cui wrote on 9/10/22 7:17 PM:

--
SLiM forward genetic simulation: http://messerlab.org/slim/
---
You received this message because you are subscribed to the Google Groups "slim-discuss" group.
To unsubscribe from this group and stop receiving emails from it, send an email to slim-discuss...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/slim-discuss/CAB03vn-ccQg_NbFM7qS4ZT2q3b554bJDgRppnsopD05cKdbyjg%40mail.gmail.com.