devRateIBMparam: Forecast ectotherm phenology as a function of temperature and...

In devRate: Quantify the Relationship Between Development Rate and Temperature in Ectotherms

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

Description

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

Usage

devRateIBMparam(
  tempTS,
  timeStepTS,
  eq,
  myParam,
  numInd = 10,
  stocha,
  timeLayEggs = 1
)

Arguments

tempTS	The temperature time series (a vector).
timeStepTS	The time step of the temperature time series (a numeric with 1 = one day).
eq	The name of the equation (e.g., lactin2_95).
myParam	The known parameters for the equation (a list of list for each life stage).
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).
timeLayEggs	The delay between emergence of adults and the time where females lay eggs in time steps (a numeric).

Details

Please note that this function is experimental and only works for the briere2_99 equation.

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

# with only one life stage
forecastX <- devRateIBMparam(
  tempTS = rnorm(n = 20, mean = 20, sd = 1),
  timeStepTS = 10,
  eq = briere2_99,
  myParam = list(
    list(
      aa = 0.0002,
      Tmin = 10,
      Tmax = 36.1,
      bb = 2.84)
  ),
  numInd = 10,
  stocha = 0.015,
  timeLayEggs = 1
)
# with two life stages
forecastXX <- devRateIBMparam(
  tempTS = rnorm(n = 20, mean = 20, sd = 1),
  timeStepTS = 10,
  eq = briere2_99,
  myParam = list(
    lifeStage01 = list(
      aa = 0.0002,
      Tmin = 10,
      Tmax = 36.1,
      bb = 2.84),
    lifeStage02 = list(
      aa = 0.0004,
      Tmin = 8,
      Tmax = 35,
      bb = 2.8)
  ),
  numInd = 10,
  stocha = 0.015,
  timeLayEggs = 1
)

devRate documentation built on Aug. 24, 2023, 9:07 a.m.