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)
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) |
A table showing frequency of stable complex systems found
sim <- Evo_rand_gen_var(max_sp = 2, iters = 1);#> [1] 2