test_that("real continuous data: iris works with jsd()", {
  data(iris)

  x <- iris$Sepal.Length[iris$Species == "setosa"]
  y <- iris$Sepal.Length[iris$Species == "versicolor"]

  res <- jsd(x, y, type = "continuous")

  expect_s3_class(res, "jsd_estimate")
  expect_equal(res$type, "continuous")
  expect_true(is.numeric(res$estimate))
  expect_length(res$estimate, 1)
  expect_true(is.finite(res$estimate))
  expect_true(res$estimate >= 0)
})

test_that("real continuous data: iris works with jsd_ci()", {
  data(iris)

  x <- iris$Sepal.Length[iris$Species == "setosa"]
  y <- iris$Sepal.Length[iris$Species == "versicolor"]

  res <- jsd_ci(x, y, type = "continuous", B = 100, seed = 1)

  expect_s3_class(res, "jsd_ci")
  expect_equal(res$type, "continuous")
  expect_true(is.numeric(res$estimate))
  expect_true(is.finite(res$estimate))
  expect_length(res$conf_int, 2)
  expect_true(all(is.finite(res$conf_int)))
  expect_true(res$conf_int[1] <= res$conf_int[2])
  expect_true(is.numeric(res$boot_se))
  expect_true(is.finite(res$boot_se))
  expect_true(res$boot_se >= 0)
})

test_that("real continuous data: iris plot_dist runs", {
  data(iris)

  x <- iris$Sepal.Length[iris$Species == "setosa"]
  y <- iris$Sepal.Length[iris$Species == "versicolor"]

  expect_no_error(
    plot_dist(
      x, y,
      type = "continuous",
      group_names = c("setosa", "versicolor")
    )
  )
})