devRateIBM {devRate}

Forecast ectotherm phenology as a function of temperature and development rate models

Description

Forecast ectotherm phenology as a function of temperature and development rate models

Usage

devRateIBM(tempTS, timeStepTS, models, numInd = 100, stocha, timeLayEggs = 1)

Arguments

tempTS

The temperature time series (a vector).

timeStepTS

The time step of the temperature time series (a numeric in days).

models

The models for development rate (a list with objects of class nls).

numInd

The number of individuals for the simulation (an integer).

stocha

The standard deviation of a Normal distribution centered on development rate to create stochasticity among individuals (a numeric). Either a single number (same stochasticity for all stages) or a vector of length corresponding to the number of models used (different stochasticity for the phenological stages).

timeLayEggs

The delay between emergence of adults and the time where females lay eggs in time steps (a numeric).

Value

A list with three elements: the table of phenology for each individual, the models used (nls objects), and the time series for temperature.

Examples

data(exTropicalMoth)
forecastTsolanivora <- devRateIBM(
   tempTS = rnorm(n = 100, mean = 15, sd = 1),
   timeStepTS = 1,
   models = exTropicalMoth[[2]],
   numInd = 100,
   stocha = c(0.015, 0.005, 0.01),
   timeLayEggs = 1)