tests/testthat/test-bfactor-bin-vs-multinom.R

In pedro-teles-fonseca/polya: Bayesian Hypothesis Testing for Categorical Data

context("Bayes factors - Bin vs Multinom match")

test_that("Austria log10(BF), unif prior", {
  expect_equal(
    round(
      log10(
        bf_binomial(
          data = austria_bl1,
          success = 1,
          null_par = theta_benford(1)[1],
          hyper_par = c(1, 1)
          )
        ),
      2),
    round(
      log10(
        bf_multinomial(
          data = ifelse(austria_bl1 == 1, "aa", "bb"),
          categories = c("aa", "bb"),
          null_par = c(theta_benford(1)[1], 1 - theta_benford(1)[1]),
          hyper_par = c(1, 1)
          )
        ),
      2)
    )
})

test_that("Austria log10(BF), unif prior", {
  expect_equal(
    round(
      log10(
        bf_binomial(
          data = austria_bl1,
          success = 1,
          null_par = theta_benford(1)[1],
          hyper_par = c(1, 1)
          )
        ),
      2),
    0.07)
})

test_that("Austria log10(BF), unif prior", {
  expect_equal(
    round(
      log10(
        bf_multinomial(
          data = ifelse(austria_bl1 == 1, "bb", "aa"),
          categories = c("bb", "aa"),
          null_par = c(theta_benford(1)[1], 1 - theta_benford(1)[1]))
        ),
      2),
    0.07)
})


test_that("Austria log10(BF), unif prior", {
  expect_equal(
    round(
      log10(
        bf_multinomial(
          data = ifelse(austria_bl1==1, "bb", "aa"),
          categories = c("aa", "bb"),
          null_par = c(1-theta_benford(1)[1], theta_benford(1)[1]))),
      2),
    0.07)
})

test_that("Austria log10(BF), unif prior", {
  expect_equal(
    round(
      log10(
        bf_multinomial(
          data = ifelse(austria_bl1==1, "bb", "aa"),
          categories = c("bb", "aa"),
          null_par = c(theta_benford(1)[1], 1-theta_benford(1)[1]),
          hyper_par = c(1, 1))),
      2),
    0.07)
})

test_that("Austria log10(BF), unif prior", {
  expect_equal(
    round(
      log10(
        bf_multinomial(
          data = ifelse(austria_bl1==1, "aa", "bb"),
          categories = c("aa", "bb"),
          null_par = c(theta_benford(1)[1], 1-theta_benford(1)[1]))),
      2),
    0.07)
})

test_that("Austria log10(BF), unif prior", {
  expect_equal(
    round(
      log10(
        bf_multinomial(
          data = ifelse(austria_bl1==1, "bb", "aa"),
          categories = c("bb", "aa"),
          null_par = c(theta_benford(1)[1], 1-theta_benford(1)[1]))
        ),
      2),
    0.07)
})

test_that("Austria log10(BF), Haldane prior", {
  expect_equal(
    round(
      log10(
        bf_binomial(
          data = austria_bl1,
          success = 1,
          haldane = TRUE,
          null_par = theta_benford(1)[1]
        )
      ),
      2),
    round(
      log10(
        bf_multinomial(
          data = ifelse(austria_bl1 == 1, "aa", "bb"),
          categories = c("aa", "bb"),
          haldane = TRUE,
          null_par = c(theta_benford(1)[1], 1 - theta_benford(1)[1])
        )
      ),
      2)
  )
})