sim.Beverton: Simulate the Stochastic Beverton-Holt Model

In BruceKendall/PVA: Population viability analysis in R

Description

Simulate the Stochastic Beverton-Holt Model

Usage

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

Arguments

N0	Starting population size
param	Vector of parameters for the model. First element is R, second is a, and the third (used only if method="normal") is the variance of the epsilons. If method="DS", then the third parameter is the intercept and the fourth is the slope of the regression of squared residuals on 1/N.
Tmax	Number of years to simulate
nsim	Number of replicate simulations
method	Treat epsilon as normally distributed random numbers ("normal"), resample from the observed residuals ("resample"), or have the variance be inversely related to population size to simulate demographic stochasticity (method="DS"). 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.Beverton(100, param=c(2,0.002,.1), Tmax=50, nsim=10)

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