glv: Simulate time series with the generalized Lotka-Volterra...

In hallucigenia-sparsa/seqtime: Time Series Analysis of Sequencing Data

Simulate time series with the generalized Lotka-Volterra model

Description

Simulate a community time series using the generalized Lotka-Volterra model, defined as \frac{dx}{dt}=x(b+Ax), where x is the vector of species abundances, A is the interaction matrix and b the vector of growth rates.

Usage

glv(N = 4, A, b = runif(N), y = runif(N), tstart = 0, tend = 100,
  tstep = 0.1, perturb = NULL)

Arguments

N	species number
A	interaction matrix
b	growth rates
y	initial abundances
tstart	initial time point
tend	final time point
tstep	time step
perturb	perturbation object describing growth changes

Value

a matrix with species abundances as rows and time points as columns, column names give time points

See Also

ricker for the Ricker model

Examples

tsplot(glv(N=4,generateA(4)),header="gLV")

hallucigenia-sparsa/seqtime documentation built on Jan. 9, 2023, 11:53 p.m.