sim.delay: Simulate the Stochastic Delay-Ricker Model

In BruceKendall/PVA: Population viability analysis in R

Description

Simulate the Stochastic delay-Ricker Model

Usage

sim.delay(N0, param, Tmax, nsim, method=c("normal","resample"), resids=NULL)

Arguments

N0	Starting population size. A two-element vector, containing population size at times 0 and -1
param	Vector of parameters for the model. First element is a0, second is a1, third is a2, and the fourth (used only if method="normal") is the variance of the epsilons.
Tmax	Number of years to simulate
nsim	Number of replicate simulations
method	Treat epsilon as normally distributed random numbers ("normal") or resample from the observed residuals ("resample"). The default is "normal"
resids	The residuals from which to resample to simulate epsilon. Used only if method="resample"

Details

This is pretty straightforward.

Value

A matrix with Tmax rows and nsim columns. Each column is the time series from a replicate simulation.

Author(s)

Bruce E. Kendall ([e-mail])

References

Morris, W. F. and Doak D. F. (2002) Quantitative Conservation Biology: Theory and Practice of Population Viability Analysis. Sunderland: Sinauer Associates.

Examples

## Generate 10 replicate simulations 50 years long
sim.delay(c(130,100), param=c(0.5,-0.001,-0.0001,.2), Tmax=50, nsim=10)

BruceKendall/PVA documentation built on Jan. 23, 2021, 2:56 a.m.