swh:1:snp:2c68a6c5a8af2f06ac2c0225927f25b54fd1f9d0
Raw File
Tip revision: 428249f43a9c6fd0c425b28deb5fee51a9525d69 authored by Dominique Makowski on 18 September 2022, 01:46:03 UTC
version 0.13.0
Tip revision: 428249f
test-p_direction.R
if (requiet("testthat") && requiet("bayestestR") && requiet("rstanarm") && requiet("brms")) {
  test_that("p_direction", {
    set.seed(333)
    x <- bayestestR::distribution_normal(10000, 1, 1)
    pd <- bayestestR::p_direction(x)
    expect_equal(as.numeric(pd), 0.842, tolerance = 0.1)
    expect_equal(as.numeric(p_direction(x, method = "kernel")), 0.842, tolerance = 0.1)
    expect_equal(nrow(p_direction(data.frame(replicate(4, rnorm(100))))), 4)
    expect_s3_class(pd, "p_direction")
    expect_equal(tail(capture.output(print(pd)), 1), "Probability of Direction: 0.84")
  })


  .runThisTest <- Sys.getenv("RunAllbayestestRTests") == "yes"
  if (.runThisTest) {
    if (require("insight")) {
      m <- insight::download_model("stanreg_merMod_5")
      p <- insight::get_parameters(m, effects = "all")

      test_that("p_direction", {
        expect_equal(
          p_direction(m, effects = "all")$pd,
          p_direction(p)$pd,
          tolerance = 1e-3
        )
      })

      m <- insight::download_model("brms_zi_3")
      p <- insight::get_parameters(m, effects = "all", component = "all")

      test_that("p_direction", {
        expect_equal(
          p_direction(m, effects = "all", component = "all")$pd,
          p_direction(p)$pd,
          tolerance = 1e-3
        )
      })
    }
  }
}
back to top