tests/testthat/test-plot.R

In chemdeg: Analysis of Chemical Degradation Kinetic Data

test_that("plot_ord works with weighted data", {
  t <- c(0, 4, 8, 12, 16, 20)
  set.seed(1)
  conc <- (0.5 * (0.05 * t + (10^(-0.5)) / 0.5))^-2 + rnorm(length(t), 0, 0.05)
  err <- c(0.02, 0.05, 0.04, 0.04, 0.03, 0.02)
  dframe <- data.frame(t, conc, err)
  res <- det_order(dframe)
  expect_no_error(plot_ord(res))
})

test_that("plot_ord works with unweighted data", {
  t <- c(0, 4, 8, 12, 16, 20)
  set.seed(1)
  conc <- (0.5 * (0.05 * t + (10^(-0.5)) / 0.5))^-2 + rnorm(length(t), 0, 0.05)
  dframe <- data.frame(t, conc)
  res <- det_order(dframe)
  expect_no_error(plot_ord(res))
})

test_that("plot_ord does not give error with no reaction order determined", {
  t <- c(0, 4, 8, 12, 16, 20)
  conc <- c(1, 0.51, 0.51, 0.7, 0.12, 0.07)
  err <- c(0.02, 0.05, 0.04, 0.04, 0.03, 0.02)
  dframe <- data.frame(t, conc, err)
  res <- det_order(dframe)
  expect_equal(res[[6]], NaN, ignore_attr = TRUE)
  expect_no_error(plot_ord(res))
})

test_that("plot_ord works when 0-order model is determined", {
  t <- c(0, 4, 8, 12, 16, 20)
  conc <- c(1, 0.91, 0.80, 0.69, 0.59, 0.51)
  err <- c(0.02, 0.05, 0.04, 0.04, 0.03, 0.02)
  dframe <- data.frame(t, conc, err)
  res <- det_order(dframe)
  expect_no_error(plot_ord(res))
})

test_that("plot_ord works only with ord_res objects", {
  t <- c(0, 4, 8, 12, 16, 20)
  conc <- c(1, 0.51, 0.51, 0.7, 0.12, 0.07)
  err <- c(0.02, 0.05, 0.04, 0.04, 0.03, 0.02)
  pippo <- lm(conc ~ t)
  pippo2 <- nls(conc ~ exp(-k * t), start = list(k = 0.5), weights = 1 / err^2)
  dframe <- data.frame(t, conc, err)
  expect_error(plot_ord(t))
  expect_error(plot_ord(dframe))
  expect_error(plot_ord(pippo))
  expect_error(plot_ord(pippo2))
})

test_that("If estimated model parameters are not significant a particular
          message is displayed - case order==0", {
  t <- c(0, 1, 2, 3, 4, 5)

  set.seed(1)
  conc <- 12 * exp(-0.01 * t) + rnorm(length(t), 0, 1)
  pippo <- data.frame(t, conc, err = abs(rnorm(length(t), 0, 1)))
  expect_no_error(res <- det_order(pippo))
  expect_no_error(plot_ord(res))
  expect_message(
    plot_ord(res),
    "NB: parameters of the regression are not significant!"
  )
})