tests/testthat/test-lsa.R

In mrgsim.sa: Sensitivity Analysis with 'mrgsolve'

library(testthat)
library(dplyr)
library(mrgsim.sa)

context("test-lsa")

mod <- mrgsolve::house()

test_that("lsa", {
  mod <- mrgsolve::ev(mod, mrgsolve::ev(amt = 100))
  out <- lsa(mod, par = "CL,VC", var = "CP")
  expect_is(out, "lsa")
  expect_is(out, "tbl_df")
  expect_equal(names(out), c("time", "dv_name", "dv_value", "p_name", "sens"))
  expect_equal(unique(out$dv_name), "CP")
  expect_equal(unique(out$p_name), c("CL", "VC"))
})

test_that("lsa input and output errors", {
  expect_error(lsa(mod, par = "a,b,c"), regex = "invalid parameter") 
  expect_error(
    lsa(mod, par = "CL", var = "a,b,c"), 
    regex = "invalid output name"
  ) 
  fun <- function(p, ...) {
    out <- mrgsim_df(mod)
    out$time <- NULL
    out
  }
  expect_error(
    lsa(mod, par = "CL", var = "CP", fun = fun), 
    regex = "output from `fun` must contain"
  )
})

test_that("lsa plot", {
  out <- lsa(mod, par = "CL", var = "CP")
  ans <- lsa_plot(out)
  expect_is(ans, "gg")
})