tests/testthat/test-parset.R

In kylebmetrum/optimhelp: Helper functions for optimizing in R

context("test-parset")

test_that("basics functions", {
  foo <- log_par("foo", 1)
  foo2 <- ident_par("foo2", 2)
  p <- parset(foo,foo2)
  expect_true(is.par(foo))
  expect_true(is.parset(p))
  expect_false(is.parset(list(a=1)))
})

test_that("log_par transforms correctly", {
  foo <- log_par("foo", 1)
  p <- parset(foo)
  i <- initials(p)
  tr <- get_trans(p)
  expect_equal(i["foo"],c(foo=log(1)))
  expect_equal(tr,log(1))
})

test_that("ident_par transforms correctly", {
  foo <- ident_par("foo", 1)
  p <- parset(foo)
  i <- initials(p)
  tr <- get_trans(p)
  expect_equal(i["foo"],c(foo=1))
  expect_equal(tr,1)
})

.logit <- function(x) log(x/(1-x))
test_that("logit_par transforms correctly", {
  foo <- logit_par("foo", 0.75)
  p <- parset(foo)
  i <- initials(p)
  tr <- get_trans(p)
  expect_equal(i["foo"],c(foo=.logit(0.75)))
  expect_equal(tr,.logit(0.75))
})

test_that("fixed parameter", {
  foo <- fixed_par("foo",value=3.3)
  foo2 <- log_par("foo2", value=3.3)
  p <- parset(foo,foo2)
  i <- initials(p)
  expect_equal(length(i),1)
  
  utr <- as.list(untrans(p))
  expect_equal(length(utr),2)
  expect_equal(utr$foo,3.3)
  expect_equal(utr$foo2,3.3)
  
  tr <- as.list(trans(p))
  expect_equal(length(tr),2)
  expect_equal(tr$foo,3.3)
  expect_equal(tr$foo2,(3.3))
})

test_that("duplicate parameter names", {
  foo1 <- log_par("foo", 3.3)
  foo2 <- log_par("foo", 4.4)
  expect_error(parset(foo1,foo2))
})

test_that("graft parameters onto parset", {
   p <- parset(log_par("foo", 2),log_par("foo2",10))
   x <- 0.2
   par <- initials(p) + x
   pp <- as.list(graft(p,par))
   expect_equal(pp$foo,exp(log(2) + x))
   expect_equal(pp$foo2,exp(log(10) + x))
})