Outer function for investigating the potential for a complex system to be stabilised by component response rates. For a stretch of values of component number (i.e., system size, S), the function runs some number of iterations, and in each iteration, if a system is not initially found to be stable with a set of component response times, then a genetic algorithm is called in attempt to find a stable solution.

Evo_rand_gen_var(max_sp, iters, int_type = 0, rmx = 0.4, C = 1)

Arguments

max_sp

The maximum system component size (S) to simulate

iters

The number of iterations (complex systems) tried for each S

int_type

The type of interaction, including random (0), competitive (1), mutualist (2), or predator-prey (3)

rmx

The standard devation of interaction (off-diagonal) elements in a complex system

C

The connectedness of a complex system (i.e., the probability that an off-diagonal element does not equal zero)

Value

A table showing frequency of stable complex systems found

Examples

sim <- Evo_rand_gen_var(max_sp = 2, iters = 1);
#> [1] 2