IPMpack (version 1.6)

stochGrowthRateSampleList: Estimating the stochastic population growth rate (lambda_s).

Description

Estimates the stochastic growth rate (lambda_s) by iteration; operates by sampling a list of IPMs. Note that the function stoch.growth.rate in the package popbio does this more efficiently and with more useful output; but may fail for large bin numbers.

Usage

stochGrowthRateSampleList(listIPMmatrix, nRunIn, tMax)

Arguments

listIPMmatrix

a list of IPMmatrix objects corresponding to possible states of the environment.

nRunIn

numeric, the size of the burnin,

tMax

numeric, the total samples desired.

Value

• a numeric converging on high enough log lambda_s for sufficient tMax; note that if the population size declines to zero, this may return NAs because of logging.

References

Caswell, 2001. Matrix population models: analysis, construction and interpretation. 2nd ed. Sinauer. p452-502.

See Also

makeListIPMs

Examples

dff <- generateData()
IPMlist <- makeListIPMs(dataf = dff, nBigMatrix = 10, minSize = -5, 
maxSize = 15,fecConstants=data.frame(1e6), correction="constant")
stochGrowthRateSampleList(listIPMmatrix = IPMlist,nRunIn = 100, tMax = 5000)