devRateQlBio function

Biological likelihood of nls fits

Return a table of 5 metrics of development (CTmin, CTmax, Topt, XTmin, XTmax)

devRateQlBio(nlsDR, propThresh = 0.01, eq, interval = c(0, 50))

Arguments

• nlsDR: A list of nls objects.
• propThresh: The proportion of maximal development rate used as a threshold for estimating XTmin and XTmax for asymptotic equations (default value is 0.01)
• eq: A list of equations used for nls fitting.
• interval: A vector containing the lower and upper boundaries of the interval of temperatures in which metrics are searched.

Returns

An object of class data.frame with development metrics (CTmin, Ctmax, Topt, XTmin, XTmax) in columns and nls objects in rows.

Details

NULL is returned when nlsDR or eq are not a list.

Examples

myDf <- data.frame(temp = seq(from = 0, to = 50, by = 10),
 rT = c(0.001, 0.008, 0.02, 0.03, 0.018, 0.004))
myNLS <- list(
 devRateModel(
   eq = janisch_32,
   df = myDf,
   startValues = list(aa = 0.2, bb = 0.1, Dmin = 10, Topt = 30),
   algo = "LM"),
 devRateModel(
   eq = kontodimas_04,
   df = myDf,
   startValues = list(aa = 1, Tmin = 7, Tmax = 40),
   algo = "LM"),
 devRateModel(
   eq = poly2,
   df = myDf,
   startValues = list(a0 = 1, a1 = 1, a2 = 1),
   algo = "LM"))
devRateQlBio(
  nlsDR = myNLS,
  eq = list(janisch_32, kontodimas_04, poly2),
  propThresh = 0.1)