5  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.

5.1 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.

5.2 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

5.3 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.

5.4 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.

5.5 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

5.6 Genetics params

mutation_rate

num_basepairs

5.7 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.

5.8 Pre-set configurations

untbParams