Genetic drift as a function of population size in the DDH model.
For all panels, the carrying capacity is Ck = 10,000 and the intrinsic growth rate is r = 2. (a) When N increases, E(K) decreases as modeled in Eq. (5). The z value of Eq. (5) (0.1, 1.5 and 3) determines the strength of N regulation, indicated by the slope of E(K) near Ck = 10,000. (b) Depending on the strength of N regulation near Ck, genetic drift can indeed decrease, increase or stay nearly constant as the population size increases. Thus, the conventional view of Ne being positively dependent on N is true only when the regulation of N is weak (the green line). At an intermediate strength (the red line), Ne is nearly independent of N. When the regulation becomes even stronger at z = 3, Ne becomes negatively dependent on N. (c-e) V(K)/E(K) of Eq. (6)) is shown as a function of N. The results of panel (b) are based on a constant V(K)/E(K) shown in panel (c). Interestingly, the results of panel (b) would not be perceptibly changed when V(K)/E(K) varies, as shown in panel (d) and (e).