This is the repository for a program simulating the evolution of an annual grass species in a semi-arid patchy environment.
This is an agent-based simulation where sessile individuals dwel in a metapopulation consisting of multiple demes, and each deme contains a facilitated and an unfactilitated habitat patch (faciltiated patches are basically the undergrowth of nurse shrubs). The individuals possess genes that determine some level of tolerance against the harsh conditions encountered outside of the shrubs, which trade off with reproductive output. See manuscript for more details.
You can use the CMake files provided here to build the program, for example by using this setup.
This program is a command-line interface program. Assuming the name of the executable is brachypode, use the following command in the terminal from the working directory to run a simulation with default parameter values:
./brachypodeTo run a simulation with custom parameters, use:
./brachypode parameters.txtwhere parameters.txt is a text file containing a list of parameter names and values.
The file parameters.txt should look like:
type 1
popsize 10
pgood 3 0.8 0.6 0.5 0.3 0.1
maxgrowths 1 2
stress 4 0
zwidths 1 2
capacities 100 10000
steep 2
dispersal 0.1
mutation 0.0001
nchrom 1
nloci 50
effect 0.1
allfreq 0.5
xmax 5
ymax 5
tradeoff 1
selfing 0.95
recombination 1
tend 100
tsave 20
seed 12345
sow 1
loadarch 0
savepars 1
savelog 0
savearch 1
talkative 0
choose 0
and should be present within the working directory where the simulation is launched. Click here for a detailed description of the parameters and their meaning.
For parameters that have to do with the genetic architecture of the evolving traits in this simulation, see here.
The user can choose which variables to save from the simulation. Those have to be encoded in the whattosave.txt file provided that parameter choose is set to 1 within the parameters.txt file, as explained here. Output data are saved into binary files with the .dat extension, which can be read by the custom-made R package brachypoder.
- brachypoder: package to read the simulations in R
- brachypode-results: generated data and analyses of the results
- brachypode-ms: manuscript
- brachypode-nb: Mathematica notebooks used for some derivations
Copyright (c) Raphaël Scherrer, 2024 (open source license will be added upon publication). This code comes with no guarantee whatsoever.