-
Notifications
You must be signed in to change notification settings - Fork 0
/
inputs.qmd
61 lines (31 loc) · 3.42 KB
/
inputs.qmd
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
# Model inputs: Setting the stage
Having just explored the components of the model, we will now review the input parameters that "set the stage" by defining both how these components are initialized and how they change over the course of the simulation.
We will also review the default configurations of parameters available in roleR and their meanings and assumptions.
## Initialization param
init_type is the model of colonization used to initialize the simulation, either "oceanic_island" or "bridge_island"
The bridge island model has the initial individuals in the local community arriving through a land bridge, while the oceanic has no bridge and is populated by a single dispersal event. Thus, for oceanic island models, all individuals are of a SINGLE species sampled proportional to meta community species abundances, while in bridge island models, individuals are sampled of MANY species proportional to their abundances.
## Community shape params
individuals_local is the number of individuals in the local community
individuals_meta is the number of individuals in the metacommunity
species_meta is the number of species in the metacommunity
## Birth vs. dispersal param
dispersal_prob is probability of dispersal (immigration) occurring from the metacommunity to the local community at each time step
Every time step, either birth or immigration ALWAYS happens, so the probability of birth is just 1 minus the dispersal_prob.
## Speciation and extinction params
speciation_local is the probability of speciation occurring in the local community at each time step.
The new local species can be either from a birth in the local community or an immigration from the metacommunity - the aforementioned dispersal_prob param determines the chance of each of these scenarios.
speciation_meta is the rate of speciation in the metacommunity.
This is used during the simulation of the start phylogeny for metacommunity species. In addition, the sum of speciation_meta and extinction_meta is equal to the average lifetime of phylo branches, and the larger this value the less new individual traits will deviate from the parent's individual trait (if parent is in the local comm) or the metacommunity mean (if parent in the metacommunity).
extinction_meta is the rate of extinction in the metacommunity. Like speciation_meta, it is used in the initial phylogeny simulation and in relation to trait deviation.
## Traits, filtering, and competition params
trait_sigma is the rate of Brownian trait evolution in the metacommunity. It determines how much the trait of a new individual deviates from its parent, i.e. how fast traits change.
env_sigma is the selectivity of the environmental filter, i.e. how strongly the environment selects which trait values are a good match for it. The larger the value, the less chance an individual will survive if it's far from the trait optimum (which is 0).
comp_sigma is the size of the competition kernel: how strongly distance from the traits of other species selects which trait values are likely to survive. The larger the value, the greater the chance an individual will survive if it is close in trait space to others
## Genetics params
mutation_rate
num_basepairs
## Meta params
niter is the number of iteration time steps to run the simulation for.
niterTimestep is the frequency (in numbers of iterations) at which the model state is snapshot and saved in a model's model steps object.
## Pre-set configurations
untbParams