simulateRicker: Simulate one or more time series from the stochastic Ricker...

Description Usage Arguments Details Value Examples

Description

Simulate one or more time series from the stochastic Ricker model

Usage

1
simulateRicker(r, K, sigma2, N0, time_horizon, nsims = 1, keep_N0 = FALSE)

Arguments

r

the low-density growth rate

K

the equilibrium density

sigma2

the variance of the log growth rates

N0

initial population size

time_horizon

number of years to simulate

nsims

number of simulations

keep_N0

(logical) whether to return N0 at the beginning of the time series

Details

This function simulates the Ricker discrete time model with environmental stochasticity, defined as follows

N(t+1) = N(t) exp[r (1 - N(t)/K) + ε(t)]

, where

ε(t) ~ Normal(0, σ)

Value

If nsim=1 (the default), a vector of length time_horizon (or time_horizon+1 if keep_N0=TRUE) with the simulated abundances. Otherwise, a matrix with time_horizon (or time_horizon+1 if keep_N0=TRUE) rows and nsim columns.

Examples

1
2
simulateRicker(2, 100, 0.01, 20, 10)
simulateRicker(2, 100, 0.01, 20, 10, 5, keep_N0=TRUE)

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