tests/testthat/test-CEriskav.R

In giabaio/BCEA: Bayesian Cost Effectiveness Analysis

# risk aversion

# library(dplyr)
# library(reshape2)

##TODO:...

test_that("zero risk aversion is equal to default", {
  
  load(test_path("ce_smoking.RData"))
  
  treats <- c("No intervention", "Self-help", "Individual counselling", "Group counselling")
  bcea_smoke <- bcea(eff, cost, ref = 4, interventions = treats, Kmax = 500)
  
  r <- 0
  CEriskav(bcea_smoke) <- r
  tol <- 0.0001
  
  expect_equivalent(bcea_smoke$U[,,1], bcea_smoke$Ur[,,1,1], tolerance = tol)
  expect_equivalent(bcea_smoke$U[,,2], bcea_smoke$Ur[,,2,1], tolerance = tol)
  expect_equivalent(bcea_smoke$U[,,3], bcea_smoke$Ur[,,3,1], tolerance = tol)
  expect_equivalent(bcea_smoke$U[,,4], bcea_smoke$Ur[,,4,1], tolerance = tol)
})