tests/testthat/test-stanreg_bridge_sampler_basic.R

In quentingronau/bridgesampling: Bridge Sampling for Marginal Likelihoods and Bayes Factors

context('bridge_sampler.stanreg works.')

test_that("stan_bridge_sampler", {

  testthat::skip_on_cran()
  testthat::skip_on_travis()

  if (require(rstanarm)) {

    fit_1 <- stan_glm(mpg ~ wt + qsec + am, data = mtcars,
                     chains = 2, cores = 2, iter = 5000,
                     diagnostic_file = file.path(tempdir(), "df.csv"))
    bridge_norm <- bridge_sampler(fit_1)

    fit_2 <- update(fit_1, formula = . ~ . + cyl)
    bridge_warp <- bridge_sampler(fit_2, method = "warp3")

    expect_true(bridge_norm$logml > bridge_warp$logml)
  }
})


test_that("stan_bridge_sampler in multicore", {
  testthat::skip_on_cran()
  testthat::skip_on_travis()
  #testthat::skip_on_os("windows")
  if (require(rstanarm)) {
    fit_1 <- stan_glm(mpg ~ wt + qsec + am, data = mtcars,
                     chains = 2, cores = 2, iter = 5000,
                     diagnostic_file = file.path(tempdir(), "df.csv"))
    bridge_norm <- bridge_sampler(fit_1, cores = 2)

    fit_2 <- update(fit_1, formula = . ~ . + cyl)
    bridge_warp <- bridge_sampler(fit_2, method = "warp3", cores = 2)

    expect_true(bridge_norm$logml > bridge_warp$logml)
  }
})