tests/testthat/test-simulate.R

In investr: Inverse Estimation/Calibration Functions

context("Simulation")

test_that("simulate.nls works as expected", {
  
  # Nonlinear least squares fit
  mod <- nls(weight ~ theta1/(1 + exp(theta2 + theta3 * log(conc))),
             start = list(theta1 = 1000, theta2 = -1, theta3 = 1),
             data = nasturtium)
  
  # Simulate new response values
  sim1 <- simulate(mod, nsim = 10L, seed = 101)
  sim2 <- simulate(mod, nsim = 10L, seed = 101)
  
  # Expectations
  expect_is(sim1, "data.frame")
  expect_equal(dim(sim1), c(42, 10))
  expect_identical(sim1, sim2)
  
})