library(vdiffr)
library(ggplot2)
library(dplyr)
library(distributional)

df <- data.frame(
  trt = factor(c(1, 1, 2, 2)),
  resp = c(1, 5, 3, 4),
  group = factor(c(1, 2, 1, 2)),
  upper = c(1.1, 5.3, 3.3, 4.2),
  lower = c(0.8, 4.6, 2.4, 3.6)
)

uncertain_df <- df |>
  group_by(trt, group) |>
  mutate(resp = dist_normal(resp, runif(1,0,0.1)),
         upper = dist_normal(upper, runif(1,0,0.2)),
         lower = dist_normal(lower, runif(1,0,0.2))
  )

suppressWarnings({
  set.seed(422)
  test_that("geom_pointrange_sample tests", {
    p1 <- ggplot(uncertain_df, aes(trt, resp, colour = group))
    #+ geom_pointrange_sample(aes(ymin = lower, ymax = upper)) 
    expect_doppelganger("Example 1", p1)
  })
})