devRateQlBio: Biological likelihood of nls fits

View source: R/devRate_qualityFits.R

devRateQlBioR Documentation

Biological likelihood of nls fits

Description

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

Usage

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.

Details

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

Value

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

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)

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