devRateIBM function

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

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).

Returns

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)