Script (R-code) for simulating negative frequency-dependence and ecological drift (Fig. 4)

This file contains the R-code needed to run the simulations in Fig. 4A-F and associated Supplemental figures. The code can be modified to vary community size (J), the strength of negative frequency-dependence (NF) and fitness ratio (FR). Each simulation run follows only 10 populations to speed up computing power, but the code can easily be modified to follow more populations, as well as the number of generations to follow. For further details, see Material and Methods in our paper and Vellend (2016; see the reference list in the paper). The original code that we have modified in this study can be found here: https://github.com/aammd/ecotheory/blob/master/vignettes/Box2.Rmd [MacDonald & Vellend. R code accompanying Vellend (2016) The Theory of Ecological Communities. Princeton University Press]