test_that("Plotting functions execute without error", {
  set.seed(123)
  sim <- rbitl(30, xi1 = 1.5, xi2 = 1.2, delta = 0.4)
  bayes_fit <- bayes_bitl(sim[, "x"], sim[, "y"], prior = "flat",
                          n_iter = 200, burn_in = 50, thin = 1)

  pdf(NULL)
  on.exit(dev.off())

  expect_silent(plot_bitl(xi1 = 1.5, xi2 = 1.2, delta = 0.4, type = "pdf", plot_type = "persp"))
  expect_silent(plot_bitl(xi1 = 1.5, xi2 = 1.2, delta = 0.4, type = "survival", plot_type = "contour"))
  expect_silent(plot_ks_bitl(sim[, "x"], sim[, "y"], xi1 = 1.5, xi2 = 1.2))
  expect_silent(plot_ttt_bitl(sim[, "x"], sim[, "y"]))
  expect_silent(plot_mcmc_bitl(bayes_fit, diag_type = "trace"))
})