Simulation plots.

<p>Simulation results showing stochastic behavior of the average (over 5 loci) of total-population and subpopulation Shannon entropies for <i>N</i> = 10000, <i>n</i> = 4, <i>μ</i> = 0.005%, <i>m</i> = 0.1% in the simulation. The horizontal line in each panel represents the theoretical equilibrium value. The initial condition was set to be just one allele (all shared) in each subpopulation. (a) The stochastic pattern for the total-population entropy ¹<i>H</i>_<i>T</i> is shown in black curve, and the red curve is ¹<i>H</i>_<i>T</i> = <i>ψ</i>[1/(1−²<i>H</i>_<i>T</i>)]+0.5772, which is the ¹<i>H</i>_<i>T</i> value calculated from a function of heterozygosity under IAM-FIM. (b) The pattern for subpopulation entropy ¹<i>H</i>_<i>S</i> is shown in black curve, and the red curve is obtained via a link from heterozygosity (see Eq. D7 in <a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0125471#pone.0125471.s004" target="_blank">S4 Appendix</a>) under IAM-FIM. In both (a) and (b), the processes converge roughly after 40000 generations, but the two lines become close before equilibrium (around 20000 generations). (c) The stochastic pattern for total- population entropy ¹<i>H</i>_<i>T</i> under SMM-FIM is shown in black curve, and the red curve is log{[1+²<i>H</i>_<i>T</i>−(²<i>H</i>_<i>T</i>)²]/(1−²<i>H</i>_<i>T</i>)}, which is the ¹<i>H</i>_<i>T</i> value calculated from a function of heterozygosity. (d) The pattern for subpopulation entropy ¹<i>H</i>_<i>S</i> is shown in black curve, and the red curve is obtained via a link from heterozygosity (see <a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0125471#pone.0125471.s004" target="_blank">S4 Appendix</a> for the link). The relationship between heterozygosity and Shannon entropy exists in all stages of the stochastic process under SMM-FIM.