tests/testthat/test-effective_sample.R

test_that("effective_sample", {
  skip_if_not_installed("curl")
  skip_if_offline()
  skip_if_not_installed("httr2")
  skip_if_not_or_load_if_installed("rstanarm")
  skip_if_not_or_load_if_installed("brms")
  skip_if_not_or_load_if_installed("rstan")

  brms_1 <- insight::download_model("brms_1")
  skip_if(is.null(brms_1))

  res <- effective_sample(brms_1)
  expect_equal(
    res,
    data.frame(
      Parameter = c("b_Intercept", "b_wt", "b_cyl"),
      ESS = c(5283, 2120, 2001),
      ESS_tail = c(3255, 2003, 2227),
      stringsAsFactors = FALSE
    ),
    ignore_attr = TRUE,
    tolerance = 1e-1
  )

  brms_null_1 <- insight::download_model("brms_null_1")
  skip_if(is.null(brms_null_1))

  res <- effective_sample(brms_null_1)
  expect_equal(
    res,
    data.frame(
      Parameter = "b_Intercept",
      ESS = 2912,
      ESS_tail = 2388,
      stringsAsFactors = FALSE
    ),
    ignore_attr = TRUE,
    tolerance = 1e-1
  )

  brms_null_2 <- insight::download_model("brms_null_2")
  skip_if(is.null(brms_null_2))

  res <- effective_sample(brms_null_2)
  expect_equal(
    res,
    data.frame(
      Parameter = "b_Intercept",
      ESS = 1098,
      ESS_tail = 954,
      stringsAsFactors = FALSE
    ),
    ignore_attr = TRUE,
    tolerance = 1e-1
  )
})

Try the bayestestR package in your browser

Any scripts or data that you put into this service are public.

bayestestR documentation built on June 8, 2025, 10:18 a.m.