tests/testthat/test-two-sample-bayes.R

In statsExpressions: Tidy Dataframes and Expressions with Statistical Details

test_that(
  desc = "bayesian (independent samples t-test)",
  code = {
    set.seed(123)
    df <- two_sample_test(
      type = "bayes",
      data = ToothGrowth,
      x = supp,
      y = len,
      paired = FALSE,
      conf.level = 0.99
    )

    expect_snapshot(df[["expression"]])
  }
)

test_that(
  desc = "bayesian (paired t-test)",
  code = {
    dat_tidy <- dplyr::filter(bugs_long, condition %in% c("HDLF", "HDHF"))

    set.seed(123)
    df_results <- two_sample_test(
      type = "bayes",
      data = dat_tidy,
      x = "condition",
      y = desire,
      paired = TRUE,
      bf.prior = 0.8
    )

    expect_equal(df_results$bf10[[1]], 40.36079, tolerance = 0.001)

    set.seed(123)
    df <- two_sample_test(
      type = "bayes",
      data = dat_tidy,
      x = "condition",
      y = desire,
      paired = TRUE,
      bf.prior = 0.8
    )

    expect_snapshot(df[["expression"]])
  }
)

test_that(
  desc = "works with subject id",
  code = {
    df <- filter(data_with_subid, condition %in% c(1, 5))

    set.seed(123)
    expr1 <- two_sample_test(
      type = "bayes",
      data = df,
      x = condition,
      y = score,
      subject.id = id,
      paired = TRUE
    )

    set.seed(123)
    expr2 <- two_sample_test(
      type = "bayes",
      data = arrange(df, id),
      x = condition,
      y = score,
      paired = TRUE
    )

    expect_equal(expr1, expr2, ignore_attr = TRUE)
  }
)