selsim
trvrb / projects

This is a simulation of the turnover of sequences within a population due to the effects of natural selection and genetic drift. It adheres to most of the standard assumptions of population genetics. There is a population of sequences that differ in their fitnesses. The population size is constrained to remain constant through time. Each generation, the population is reconstituted by sampling with replacement from the pool of sequences. If a sequence has 2× the fitness of another sequence, then it is twice as likely to be represented in the next generation. Mutations occur that change the sequences and their fitnesses.

I used this program to study the effects of selection on genealogy shape. Source code is provided. Please cite this paper if you find the program useful.

Current version: 1.0
Released 8/23/2011

Download:       source