simulateSOI {miaSim} | R Documentation |
Generate time-series with The Self-Organised Instability (SOI) model. Implements a K-leap method for accelerating stochastic simulation.
simulateSOI(n.species, I, A, k = 5, com = NULL, tend, norm = FALSE)
n.species |
Integer: number of species |
I |
Integer: community size, number of available sites (individuals) |
A |
interaction matrix |
k |
Integer: the number of transition events that are allowed to take
place during one leap. (default: |
com |
a vector: initial community abundances
If (default: |
tend |
Integer: number of timepoints to be returned in the time series (number of generations) |
norm |
Logical: indicates whether the time series should be returned
with the abundances as proportions ( |
SummarizedExperiment object containing abundance matrix consisting of species abundance as rows and time points as columns
A <- miaSim::powerlawA(10, alpha = 1.2) ExampleSOI <- simulateSOI(n.species = 10, I = 1000, A, k=5, com = NULL, tend = 150, norm = TRUE)