growth: Population growth

growthR Documentation

Population growth


This function applies negative or positive growth to the population using matrix multiplication. Stochasticity can be added to cell-based transition matrices or globally. Users can also specify a built-in or custom function to modify the transition matrices throughout a simulation. Please see the tutorial vignette titled "Creating custom *steps* functions" for information on how to write custom functions for use in simulations.


  global_stochasticity = 0,
  local_stochasticity = 0,
  transition_function = NULL,
  transition_order = c("fecundity", "survival"),
  two_sex = FALSE



A symmetrical age-based (Leslie) or stage-based (Lefkovitch) population structure matrix.

global_stochasticity, local_stochasticity

Either scalar values or matrices (with the same dimension as transition_matrix) specifying the variability in the transition matrix either for populations in all grid cells (global_stochasticity) or for each grid cell population separately (local_stochasticity). Values supplied here are the standard deviation of a truncated normal distribution where the mean is the value supplied in the transition matrix.


A function to specify or modify life-stage transitions at each timestep. See transition_function.


Order of transitions performed in growth function. This behaviour is only applied when demographic stochasticity is set to "full" (default) and transitions are applied sequentially. By default "fecundity" is performed first (calculating the number of new individuals to be added to the populations), then "survival" is applied. The final population is the sum of these. Users should be cautious of specifying "survival" to be performed first as typically survival of reproductive stages will already be accounted for in the fecundity values of the transition matrix.


Does the transition matrix include life stages for two sexes (i.e. male and female)? Default is FALSE which assumes a single sex matrix (e.g. females only).


# Example of a growth function that changes the populations based on a transition matrix that
# is subject to global stochasticity. 

## Not run: 
stoch_growth <- growth(transition_matrix = egk_mat, global_stochasticity = egk_mat_stoch)

ls <- landscape(population = egk_pop, suitability = NULL, carrying_capacity = NULL)

pd <- population_dynamics(change = stoch_growth)

simulation(landscape = ls, population_dynamics = pd, habitat_dynamics = NULL, timesteps = 20)

## End(Not run)

steps documentation built on Oct. 5, 2022, 1:06 a.m.