test_that("stats work", { p2 <- study |> tidyplot(treatment, score, color = treatment) |> add_mean_bar(alpha = 0.4) |> add_sem_errorbar() p2 |> add_test_pvalue() |> vdiffr::expect_doppelganger("add stats pvalue no ref.group", fig = _) p2 |> add_test_pvalue(ref.group = 1) |> vdiffr::expect_doppelganger("add stats pvalue ref.group", fig = _) p2 |> add_test_pvalue(comparisons = list(c(1,4),c(2,3))) |> vdiffr::expect_doppelganger("add stats pvalue comparisons", fig = _) p2 |> add_test_pvalue(ref.group = 1, p.adjust.method = "bonferroni") |> vdiffr::expect_doppelganger("add stats pvalue ref.group bonferroni", fig = _) p2 |> add_test_asterisks(ref.group = 1, p.adjust.method = "bonferroni") |> vdiffr::expect_doppelganger("add stats asterisks ref.group bonferroni", fig = _) p2 |> add_test_asterisks(comparisons = list(c(1,4),c(2,3))) |> vdiffr::expect_doppelganger("add stats asterisks comparisons", fig = _) p3 <- study |> tidyplot(x = dose, y = score, color = group) |> add_mean_bar(alpha = 0.4) |> add_sem_errorbar() |> add_data_points_beeswarm() p3 |> add_test_pvalue() |> vdiffr::expect_doppelganger("add grouped stats pvalue", fig = _) p3 |> add_test_pvalue(p.adjust.method = "bonferroni") |> vdiffr::expect_doppelganger("add grouped stats pvalue bonferroni", fig = _) p3 |> add_test_asterisks(p.adjust.method = "bonferroni") |> vdiffr::expect_doppelganger("add grouped stats asterisks bonferroni", fig = _) }) test_that("paired stats work", { x <- c(2.3, 4.5, 6.3, 3.4, 7.8, 6.7) df <- data.frame( x = c(x, x + c(0.8, 0.75)), group = paste0("g", rep(c(1, 2), each = 6)), batch = paste0("b", c(1:6, 1:6)), control = paste0("c", c(1:6, 6:1)) ) df |> tidyplot(group, x, color = group) |> add_boxplot() |> add_data_points_beeswarm() |> add_test_pvalue(paired_by = control) |> add_line(group = control, color = "black") |> vdiffr::expect_doppelganger("add paired stats pvalue", fig = _) df |> tidyplot(group, x, color = group) |> add_boxplot() |> add_data_points_beeswarm() |> add_test_pvalue(paired_by = batch) |> add_line(group = batch, color = "black") |> vdiffr::expect_doppelganger("add paired stats pvalue control", fig = _) df |> tidyplot(group, x, color = group) |> add_boxplot() |> add_data_points_beeswarm() |> add_test_asterisks(paired_by = control) |> add_line(group = control, color = "black") |> vdiffr::expect_doppelganger("add paired stats asterisks", fig = _) df |> tidyplot(group, x, color = group) |> add_boxplot() |> add_data_points_beeswarm() |> add_test_asterisks(paired_by = batch) |> add_line(group = batch, color = "black") |> vdiffr::expect_doppelganger("add paired stats asterisks control", fig = _) })