tests/testthat/test_devRate_nls.R

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

# ------------------------------------------------------------------------------
# Tests for devRateModel
# ------------------------------------------------------------------------------

myT <- 5:15
myDev <- -0.05 + rnorm(n = length(myT), mean = myT, sd = 1) * 0.01
myDf <- data.frame(myT, myDev)

test_that("Returns NLS object",{
  res <- devRateModel(eq = campbell_74, df = myDf)
  expect_equal(
    object = class(res),
    expected = "nls")
  res2 <- devRateModel(eq = campbell_74, temp = myT, devRate = myDev)
  expect_equal(
    object = class(res2),
    expected = "nls")
})

test_that("Returns NLS object with LM algo",{
  res <- devRateModel(eq = campbell_74, df = myDf, algo = "LM")
  expect_equal(
    object = class(res),
    expected = "nls")
  res2 <- devRateModel(eq = campbell_74, temp = myT, devRate = myDev, algo = "LM")
  expect_equal(
    object = class(res2),
    expected = "nls")
})

test_that("NLS stopCode 0 (convergence)",{
  res <- devRateModel(eq = campbell_74, df = myDf)
  expect_equal(
    object = res$convInfo$stopCode,
    expected = 0)
})

# stinner

rawDevEggs <- matrix(c(10, 0.031, 10, 0.039, 15, 0.047, 15, 0.059, 15.5, 0.066,
  13, 0.072, 16, 0.083, 16, 0.100, 17, 0.100, 20, 0.100, 20, 0.143, 25, 0.171,
  25, 0.200, 30, 0.200, 30, 0.180, 35, 0.001), ncol = 2, byrow = TRUE)

test_that("stinner",{
  res <- devRateModel(
    eq = stinner_74,
    temp = rawDevEggs[,1],
    devRate = rawDevEggs[,2],
    startValues = list(
      list(C = 0.1, k1 = 4, k2 = -0.2),
      list(Topt = 30)
    )
  )
  expect_equal(
    object = class(res[[1]]),
    expected = "nls")
  expect_equal(
    object = class(res[[2]]),
    expected = "nls")
  expect_equal(
    object = res[[1]]$convInfo$stopCode,
    expected = 0)
  expect_equal(
    object = res[[2]]$convInfo$stopCode,
    expected = 0)

})

# ------------------------------------------------------------------------------
# Tests for devRatePrint
# ------------------------------------------------------------------------------

myT <- 5:15
myDev <- -0.05 + rnorm(n = length(myT), mean = myT, sd = 1) * 0.01
myDf <- data.frame(myT, myDev)

test_that("print devRateModel",{
  modL <- devRateModel(eq = campbell_74, df = myDf)
  res <- devRatePrint(myNLS = modL)
  expect_equal(
    object = length(res),
    expected = 6
  )
})

test_that("plot devRateModel",{
  modL <- devRateModel(eq = campbell_74, df = myDf)
  res <- devRatePrint(
    myNLS = modL,
    doPlots = TRUE)
  expect_equal(
    object = length(res),
    expected = 6
  )
})