Feature: Instantaneous Size Change

Share:

Description

This feature changes the effective population size of one population. The change is performed at a given time point and applies to the time interval further on into the past from this point. The population size is set to a fraction of N0.

Usage

1
feat_size_change(new_size, population = 1, time = "0")

Arguments

new_size

A parameter giving the new size of the population, as a fraction of N0.

population

The number of the population whichs size changes. Can also be set to "all". Then the size changes applies to all populations.

time

The time at which the population's size is changed.

Value

The feature, which can be added to a model using '+'.

See Also

For continuous size changes over time: feat_growth.

Other features: feat_growth, feat_ignore_singletons, feat_migration, feat_mutation, feat_outgroup, feat_pop_merge, feat_recombination, feat_selection, feat_unphased

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
# A model with one smaller population:
model <- coal_model(c(20, 5), 3) +
  feat_size_change(.1, population = 2) +
  feat_mutation(1.0) +
  feat_migration(0.5, 2, 1) +
  sumstat_sfs()
simulate(model)

# A model of one population that experienced a bottleneck:
model <- coal_model(10, 1) +
  feat_size_change(0.1, time = 0.3) +
  feat_size_change(1.0, time = 0.5) +
  feat_mutation(20) +
  sumstat_sfs()
simulate(model)

Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.