Abstract:Simulations are an essential tool in all areas of population genetic research, used in tasks such as the validation of theoretical analysis and the study of complex evolutionary models. Forward-in-time simulations are especially flexible, allowing for various types of natural selection, complex genetic architectures, and non-Wright-Fisher dynamics. However, their intense computational requirements can be prohibitive to simulating large populations and genomes. A popular method to alleviate this burden is to scale down the population size by some scaling factor while scaling up the mutation rate, selection coefficients, and recombination rate by the same factor. However, this rescaling approach may in some cases bias simulation results. To investigate the manner and degree to which rescaling impacts simulation outcomes, we carried out simulations with different demographic histories and distributions of fitness effects using several values of the rescaling factor, , and compared the deviation of key outcomes (fixation times, allele frequencies, linkage disequilibrium, and the fraction of mutations that fix during the simulation) between the scaled and unscaled simulations. Our results indicate that scaling introduces substantial biases to each of these measured outcomes, even at small values of . Moreover, the nature of these effects depends on the evolutionary model and scaling factor being examined. While increasing the scaling factor tends to increase the observed biases, this relationship is not always straightforward, thus it may be difficult to know the impact of scaling on simulation outcomes . However, it appears that for most models, only a small number of replicates was needed to accurately quantify the bias produced by rescaling for a given . In summary, while rescaling forward-in-time simulations may be necessary in many cases, researchers should be aware of the rescaling procedure’s impact on simulation outcomes and consider investigating its magnitude in smaller scale simulations of the desired model(s) before selecting an appropriate value of .

How big does a population need to be before demographers can ignore individual-level randomness in demographic events?

Modeling the Age Pattern of Fertility: An Individual-Level Approach

Demographic inference using genetic data from a single individual: Separating population size variation from population structure

Population Size in Evolutionary Biology Is More Than the Effective Size

Strong dependence of demographic stochasticity on a small effective number of breeders in iteroparous species

The random model for simulation of the growth of small populations

Practical considerations regarding estimating Ne when generations overlap

Population aging caused by rise in sex ratio at birth

Population size rescaling significantly biases outcomes of forward-in-time population genetic simulations

Discreteness Effects in Population Dynamics

Some deterministic structured population models which are limit of stochastic individual based models

Population Aging Caused by a Rise in the Sex Ratio at Birth

Selection Against Demographic Stochasticity in Age-Structured Populations

The Critical Domain Size of Stochastic Population Models

Do estimates of contemporary effective population size tell us what we want to know?

Can Life History Predict the Effect of Demographic Stochasticity on Extinction Risk?

Finite-size effects, demographic noise, and ecosystem dynamics

A Large Population Size Can Be Unhelpful in Evolutionary Algorithms

An inevitably aging world -- Analysis on the evolutionary pattern of age structure in 200 countries

Exact confidence intervals for population growth rate, longevity and generation time

Population growth in discrete time: a renewal equation oriented survey